Changeset 10282

Timestamp:

10/25/11 08:25:05 (13 years ago)

Author:

seroussi

Message:

groundingline test

File:

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

issm/trunk/test/NightlyRun/test446.m

@@ -5 +5 @@
 md=roundmesh(model,radius,resolution);
 %fix center node to 0,0
-rad=sqrt((md.mesh.x).*md.mesh.x+(md.mesh.y).*md.mesh.y);
+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
@@ -12 +12 @@
 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).*xelem+(yelem).*yelem);
+rad=sqrt(xelem.^2+yelem.^2);
 flags=zeros(md.mesh.numberofelements,1);
 pos=find(rad>=(radius-shelfextent));
 flags(pos)=1;
-flags2=zeros(md.mesh.numberofelements,1);
-md=setmask(md,flags,flags2); 
+md=setmask(md,flags,''); 
 %}}}
 md=parameterize(md,'../Par/RoundSheetShelf.par');
@@ -24 +23 @@
 %plug holes into the ice sheet, to test for grounding line migration. {{{1
 di=md.materials.rho_ice/md.materials.rho_water;
-rad=sqrt((md.mesh.x).*md.mesh.x+(md.mesh.y).*md.mesh.y);
+rad=sqrt(md.mesh.x.^2+md.mesh.y.^2);
 pos=find(rad<200000);
 md.geometry.thickness(pos)=100;
@@ -42 +41 @@
 
 %test different grounding line dynamics.
-%grounding line: 
+md.transient.isthermal=0;
+md.transient.isprognostic=0;
+md.transient.isdiagnostic=0;
+md.transient.isgroundingline=1;
 md.groundingline.migration='AgressiveMigration';
-md=solve(md,GroundinglineMigration2dSolutionEnum);
-element_on_iceshelf_agressive=PatchToVec(md.results.GroundinglineMigration2dSolution.MaskElementonfloatingice);
+md=solve(md,TransientSolutionEnum);
+element_on_iceshelf_agressive=PatchToVec(md.results.TransientSolution.MaskElementonfloatingice);
 
 md.groundingline.migration='SoftMigration';
-md=solve(md,GroundinglineMigration2dSolutionEnum);
-element_on_iceshelf_soft=PatchToVec(md.results.GroundinglineMigration2dSolution.MaskElementonfloatingice);
+md=solve(md,TransientSolutionEnum);
+element_on_iceshelf_soft=PatchToVec(md.results.TransientSolution.MaskElementonfloatingice);
 
 %Fields and tolerances to track changes