Changeset 10290

Timestamp:

10/25/11 10:09:29 (13 years ago)

Author:

seroussi

Message:

minor

File:

• issm/trunk/test/NightlyRun/test446.m (modified) (3 diffs)

issm/trunk/test/NightlyRun/test446.m

@@ -1 +1 @@
 radius=1e6;
 shelfextent=2e5;
-resolution=.5e5;
-%mesh {{{1
-md=roundmesh(model,radius,resolution);
 %fix center node to 0,0
 rad=sqrt(md.mesh.x.^2+md.mesh.y.^2);
 pos=find(rad==min(rad));
 md.mesh.x(pos)=0; md.mesh.y(pos)=0; %the closest node to the center is changed to be exactly at the center
-%}}}
-%put our grounding line 'shelfextent' meters from the icefront {{{1
 xelem=md.mesh.x(md.mesh.elements)*[1;1;1]/3;
 yelem=md.mesh.y(md.mesh.elements)*[1;1;1]/3;
@@ -17 +12 @@
 flags(pos)=1;
 md=setmask(md,flags,''); 
-%}}}
 md=parameterize(md,'../Par/RoundSheetShelf.par');
 md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
-%plug holes into the ice sheet, to test for grounding line migration. {{{1
+%plug holes into the ice sheet, to test for grounding line migration. 
 di=md.materials.rho_ice/md.materials.rho_water;
-rad=sqrt(md.mesh.x.^2+md.mesh.y.^2);
 pos=find(rad<200000);
 md.geometry.thickness(pos)=100;
@@ -38 +31 @@
 md.geometry.bed(pos)=-di*md.geometry.thickness(pos)-20;
 md.geometry.surface(pos)=md.geometry.bed(pos)+md.geometry.thickness(pos);
-%}}}
 
 %test different grounding line dynamics.