Changeset 9641 for issm/trunk/test/Par/RoundSheetShelf.par

Timestamp:

09/07/11 08:46:46 (14 years ago)

Author:

Mathieu Morlighem

Message:

Added mask object

File:

• issm/trunk/test/Par/RoundSheetShelf.par (modified) (5 diffs)

issm/trunk/test/Par/RoundSheetShelf.par

@@ -13 +13 @@
 md.bed=-md.materials.rho_ice/md.materials.rho_water*md.thickness;
 
-pos=find(md.nodeonicesheet);
+pos=find(md.mask.vertexongroundedice);
 md.bed(pos)=md.bed(pos)-300*(radius(pos)-(rad-shelfextent))/(rad-shelfextent);
 md.surface=md.bed+md.thickness;
@@ -47 +47 @@
 md.surfaceforcings.mass_balance=-10*ones(md.numberofnodes,1);
 md.basalforcings.melting_rate=zeros(md.numberofnodes,1);
-pos=find(md.nodeoniceshelf);md.basalforcings.melting_rate(pos)=10;
+pos=find(md.mask.vertexonfloatingice);md.basalforcings.melting_rate(pos)=10;
 
 %Friction
-pos=find(md.elementoniceshelf);
+pos=find(md.mask.elementonfloatingice);
 md.friction.coefficient=20*ones(md.numberofnodes,1);
 md.friction.coefficient(md.elements(pos,:))=0;
@@ -59 +59 @@
 md.viscosity_overshoot=0.0;
 md.artificial_diffusivity=1;
-md.thermal.stabilisation=1;
+md.thermal.stabilization=1;
 md.verbose=verbose(0);
 md.waitonlock=30;
@@ -71 +71 @@
 md.groundingline.migration='AgressiveMigration';
 md.bathymetry=md.bed;
-pos=find(md.nodeoniceshelf); md.bathymetry(pos)=md.bed(pos)-900;
+pos=find(md.mask.vertexonfloatingice); md.bathymetry(pos)=md.bed(pos)-900;
 
 %Deal with boundary conditions:
@@ -83 +83 @@
 
 pressureload=md.segments;
-pressureload=[pressureload 1*md.elementoniceshelf(pressureload(:,end)) + 0*md.elementonicesheet(pressureload(:,end))];
+pressureload=[pressureload 1*md.mask.elementonfloatingice(pressureload(:,end)) + 0*md.mask.elementongroundedice(pressureload(:,end))];
 md.pressureload=pressureload;
 md.spcthickness=NaN*ones(md.numberofnodes,1);