source: issm/oecreview/Archive/15392-16133/ISSM-15993-15994.diff@ 16134

Index: ../trunk-jpl/test/Par/RoundSheetShelf.par
===================================================================
--- ../trunk-jpl/test/Par/RoundSheetShelf.par   (revision 15993)
+++ ../trunk-jpl/test/Par/RoundSheetShelf.par   (revision 15994)
@@ -12,7 +12,7 @@
 md.geometry.thickness=hmax+(hmin-hmax)*(radius-ymin)/(ymax-ymin);
 md.geometry.bed=-md.materials.rho_ice/md.materials.rho_water*md.geometry.thickness;
 
-pos=find(md.mask.vertexongroundedice);
+pos=find(md.mask.groundedice_levelset>0.);
 md.geometry.bed(pos)=md.geometry.bed(pos)-300.*(radius(pos)-(rad-shelfextent))/(rad-shelfextent);
 md.geometry.surface=md.geometry.bed+md.geometry.thickness;
 
@@ -63,12 +63,19 @@
 %Surface mass balance and basal melting
 md.surfaceforcings.mass_balance=-10.*ones(md.mesh.numberofvertices,1);
 md.basalforcings.melting_rate=zeros(md.mesh.numberofvertices,1);
-pos=find(md.mask.vertexonfloatingice);md.basalforcings.melting_rate(pos)=10.;
+pos=find(md.mask.groundedice_levelset>0.);md.basalforcings.melting_rate(pos)=10.;
 md.basalforcings.geothermalflux=ones(md.mesh.numberofvertices,1);
 
 %Friction
+radius=1.e6;
+shelfextent=2.e5;
 md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
-md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
+xelem=md.mesh.x(md.mesh.elements)*[1;1;1]/3.;
+yelem=md.mesh.y(md.mesh.elements)*[1;1;1]/3.;
+rad=sqrt(xelem.^2+yelem.^2);
+flags=zeros(md.mesh.numberofelements,1);
+pos=find(rad>=(radius-shelfextent));
+md.friction.coefficient(md.mesh.elements(pos,:))=0.;
 md.friction.p=ones(md.mesh.numberofelements,1);
 md.friction.q=ones(md.mesh.numberofelements,1);
 
@@ -88,7 +95,7 @@
 %bathymetry and grounding line migration:
 md.groundingline.migration='AgressiveMigration';
 md.geometry.bathymetry=md.geometry.bed;
-pos=find(md.mask.vertexonfloatingice); md.geometry.bathymetry(pos)=md.geometry.bed(pos)-900.;
+pos=find(md.mask.groundedice_levelset<0.); md.geometry.bathymetry(pos)=md.geometry.bed(pos)-900.;
 
 %Deal with boundary conditions:
 md.stressbalance.spcvx=NaN*ones(md.mesh.numberofvertices,1);