Changeset 7359

Timestamp:

02/05/11 20:56:49 (14 years ago)

Author:

Eric.Larour

Message:

Rewrote grounding line test

Location:

issm/trunk/test

Files:

• Archives/Archive446.mat (modified) ( previous)
• NightlyRun/test446.m (modified) (2 diffs)
• Par/RoundSheetShelf.par (modified) (2 diffs)

issm/trunk/test/NightlyRun/test446.m

@@ -2 +2 @@
 shelfextent=2e5;
 resolution=.5e5;
-
 %mesh {{{1
 md=roundmesh(model,radius,resolution);
@@ -23 +22 @@
 md=setelementstype(md,'macayeal','all');
 md=SetParallel(md,3);
-md=solve(md,'analysis_type',Transient2DSolutionEnum);
+%plug holes into the ice sheet, to test for grounding line migration. {{{1
+di=md.rho_ice/md.rho_water;
+rad=sqrt((md.x).*md.x+(md.y).*md.y);
+pos=find(rad<200000);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
+
+pos=find(md.x<.2*1e6 & md.x>-.2*1e6 & md.y>0);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
+
+pos=find(md.x<.1*1e6 & md.x>-.1*1e6 & md.y<-.5*1e6 & md.y>-.6*1e6);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
+%}}}
+
+%test different grounding line dynamics.
+%grounding line: 
+md.gl_migration=AgressiveMigrationEnum;
+md=solve(md,'analysis_type',GroundingLineMigration2DSolutionEnum);
+element_on_iceshelf_agressive=PatchToVec(md.results.GroundingLineMigration2DSolution.ElementOnIceShelf);
+
+md.gl_migration=SoftMigrationEnum;
+md=solve(md,'analysis_type',GroundingLineMigration2DSolutionEnum);
+element_on_iceshelf_soft=PatchToVec(md.results.GroundingLineMigration2DSolution.ElementOnIceShelf);
 
 %Fields and tolerances to track changes
-field_names     ={'ElementOnIceShelf','Vx50','Vy50','Vel50'};
-field_tolerances={1e-13,1e-13,1e-13,1e-13};
-field_values={...
-        PatchToVec(md.results.Transient2DSolution(50).ElementOnIceShelf),...
-        PatchToVec(md.results.Transient2DSolution(50).Vx),...
-        PatchToVec(md.results.Transient2DSolution(50).Vy),...
-        PatchToVec(md.results.Transient2DSolution(50).Vel),...
-        };
+field_names     ={'ElementOnIceShelfAgressive','ElementOnIceShelfSoft'};
+field_tolerances={1e-13,1e-13};
+field_values={element_on_iceshelf_agressive,element_on_iceshelf_soft};

issm/trunk/test/Par/RoundSheetShelf.par

@@ -1 +1 @@
 %Start defining model parameters here
 
+di=md.rho_ice/md.rho_water;
 rad=1e6;
 shelfextent=2e5;
@@ -15 +16 @@
 md.bed(pos)=md.bed(pos)-300*(radius(pos)-(rad-shelfextent))/(rad-shelfextent);
 md.surface=md.bed+md.thickness;
+
+pos=find(radius<200000);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
+
+
+pos=find(md.x<.2*1e6 & md.x>-.2*1e6 & md.y>0);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
+
+pos=find(md.x<.1*1e6 & md.x>-.1*1e6 & md.y<-.5*1e6 & md.y>-.6*1e6);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
 
 %Initial velocity