Changeset 18739
- Timestamp:
- 11/05/14 15:52:15 (10 years ago)
- Location:
- issm/trunk-jpl/src/m
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m
r18738 r18739 54 54 if strcmp(elementtype(md.mesh),'Penta'), 55 55 numbernodesfront=4; 56 %I added the next part. If not, there would be no segments !57 md.mesh.vertexconnectivity=NodeConnectivity(md.mesh.elements2d,md.mesh.numberofvertices2d);58 md.mesh.elementconnectivity=ElementConnectivity(md.mesh.elements2d,md.mesh.vertexconnectivity);59 segments2=contourenvelope(md);60 md.mesh.vertexonboundary=zeros(md.mesh.numberofvertices/md.mesh.numberoflayers,1); md.mesh.vertexonboundary(segments2(:,1:2))=1;61 md.mesh.vertexonboundary=repmat(md.mesh.vertexonboundary,md.mesh.numberoflayers,1);62 %Then do it for 3d as usual63 md.mesh.vertexconnectivity=NodeConnectivity(md.mesh.elements,md.mesh.numberofvertices);64 md.mesh.elementconnectivity=ElementConnectivity(md.mesh.elements,md.mesh.vertexconnectivity);65 56 elseif strcmp(elementtype(md.mesh),'Tria'), 66 57 numbernodesfront=2; 67 segments2=md.mesh.segments;68 58 else 69 59 error('mesh type not supported yet'); 70 60 end 71 segmentsfront=md.mask.ice_levelset( segments2(:,1:numbernodesfront))==0;61 segmentsfront=md.mask.ice_levelset(md.mesh.segments(:,1:numbernodesfront))==0; 72 62 segments=find(sum(segmentsfront,2)~=numbernodesfront); 73 63 %Find all nodes for these segments and spc them 74 pos= segments2(segments,1:numbernodesfront);64 pos=md.mesh.segments(segments,1:numbernodesfront); 75 65 md.stressbalance.spcvx(pos(:))=0; 76 66 md.stressbalance.spcvy(pos(:))=0; -
issm/trunk-jpl/src/m/classes/clusters/hpc.m
r18738 r18739 50 50 function md = checkconsistency(cluster,md,solution,analyses) % {{{ 51 51 52 available_queues={'pub64','free64','free48','free*,pub64','free*' ,'free*,free64'};53 queue_requirements_time=[Inf Inf Inf Inf Inf Inf];54 queue_requirements_np=[64 64 48 48 48 48];52 available_queues={'pub64','free64','free48','free*,pub64','free*'}; 53 queue_requirements_time=[Inf Inf Inf Inf Inf]; 54 queue_requirements_np=[64 64 48 48 48]; 55 55 56 56 QueueRequirements(available_queues,queue_requirements_time,queue_requirements_np,cluster.queue,cluster.np,1) -
issm/trunk-jpl/src/m/os/oshostname.m
r18738 r18739 5 5 % 6 6 % 7 hostname='sylvestre';8 return;9 7 10 8 if ispc(), -
issm/trunk-jpl/src/m/solve/waitonlock.m
r18738 r18739 61 61 time=etime(clock,time0)/60; 62 62 else 63 pause( 60);63 pause(5); 64 64 time=etime(clock,time0); 65 65 fprintf('\rchecking for job completion (time: %i min %i sec) ',floor(time/60),floor(rem(time,60)));
Note:
See TracChangeset
for help on using the changeset viewer.