Changeset 15610

Timestamp:

07/25/13 09:15:25 (12 years ago)

Author:

seroussi

Message:

FIX: test 423 no front was created, do it in the par file

Location:

issm/trunk-jpl/test/Par

Files:

• RoundSheetShelf.par (modified) (1 diff)
• RoundSheetShelf.py (modified) (1 diff)

issm/trunk-jpl/test/Par/RoundSheetShelf.par

@@ -104 +104 @@
 pressureload=[pressureload 1*md.mask.elementonfloatingice(pressureload(:,end)) + 0*md.mask.elementongroundedice(pressureload(:,end))];
 md.diagnostic.icefront=pressureload;
+md.mask.icelevelset(find(md.mesh.vertexonboundary))=0;
 md.balancethickness.spcthickness=NaN*ones(md.mesh.numberofvertices,1);
 md.prognostic.spcthickness=NaN*ones(md.mesh.numberofvertices,1);

issm/trunk-jpl/test/Par/RoundSheetShelf.py

@@ -113 +113 @@
 pressureload=numpy.hstack((pressureload,(1*md.mask.elementonfloatingice[pressureload[:,-1].astype(int)-1] + 0*md.mask.elementongroundedice[pressureload[:,-1].astype(int)-1]).reshape(-1,1)))
 md.diagnostic.icefront=pressureload
+pos=numpy.nonzero(md.mesh.vertexonboundary)
+md.mask.icelevelset[pos]=0
 md.balancethickness.spcthickness=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
 md.prognostic.spcthickness=float('nan')*numpy.ones((md.mesh.numberofvertices,1))