Changeset 7359


Ignore:
Timestamp:
02/05/11 20:56:49 (14 years ago)
Author:
Eric.Larour
Message:

Rewrote grounding line test

Location:
issm/trunk/test
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/test/NightlyRun/test446.m

    r7075 r7359  
    22shelfextent=2e5;
    33resolution=.5e5;
    4 
    54%mesh {{{1
    65md=roundmesh(model,radius,resolution);
     
    2322md=setelementstype(md,'macayeal','all');
    2423md=SetParallel(md,3);
    25 md=solve(md,'analysis_type',Transient2DSolutionEnum);
     24%plug holes into the ice sheet, to test for grounding line migration. {{{1
     25di=md.rho_ice/md.rho_water;
     26rad=sqrt((md.x).*md.x+(md.y).*md.y);
     27pos=find(rad<200000);
     28md.thickness(pos)=100;
     29md.bed(pos)=-di*md.thickness(pos)-20;
     30md.surface(pos)=md.bed(pos)+md.thickness(pos);
     31
     32pos=find(md.x<.2*1e6 & md.x>-.2*1e6 & md.y>0);
     33md.thickness(pos)=100;
     34md.bed(pos)=-di*md.thickness(pos)-20;
     35md.surface(pos)=md.bed(pos)+md.thickness(pos);
     36
     37pos=find(md.x<.1*1e6 & md.x>-.1*1e6 & md.y<-.5*1e6 & md.y>-.6*1e6);
     38md.thickness(pos)=100;
     39md.bed(pos)=-di*md.thickness(pos)-20;
     40md.surface(pos)=md.bed(pos)+md.thickness(pos);
     41%}}}
     42
     43%test different grounding line dynamics.
     44%grounding line:
     45md.gl_migration=AgressiveMigrationEnum;
     46md=solve(md,'analysis_type',GroundingLineMigration2DSolutionEnum);
     47element_on_iceshelf_agressive=PatchToVec(md.results.GroundingLineMigration2DSolution.ElementOnIceShelf);
     48
     49md.gl_migration=SoftMigrationEnum;
     50md=solve(md,'analysis_type',GroundingLineMigration2DSolutionEnum);
     51element_on_iceshelf_soft=PatchToVec(md.results.GroundingLineMigration2DSolution.ElementOnIceShelf);
    2652
    2753%Fields and tolerances to track changes
    28 field_names     ={'ElementOnIceShelf','Vx50','Vy50','Vel50'};
    29 field_tolerances={1e-13,1e-13,1e-13,1e-13};
    30 field_values={...
    31         PatchToVec(md.results.Transient2DSolution(50).ElementOnIceShelf),...
    32         PatchToVec(md.results.Transient2DSolution(50).Vx),...
    33         PatchToVec(md.results.Transient2DSolution(50).Vy),...
    34         PatchToVec(md.results.Transient2DSolution(50).Vel),...
    35         };
     54field_names     ={'ElementOnIceShelfAgressive','ElementOnIceShelfSoft'};
     55field_tolerances={1e-13,1e-13};
     56field_values={element_on_iceshelf_agressive,element_on_iceshelf_soft};
  • issm/trunk/test/Par/RoundSheetShelf.par

    r7329 r7359  
    11%Start defining model parameters here
    22
     3di=md.rho_ice/md.rho_water;
    34rad=1e6;
    45shelfextent=2e5;
     
    1516md.bed(pos)=md.bed(pos)-300*(radius(pos)-(rad-shelfextent))/(rad-shelfextent);
    1617md.surface=md.bed+md.thickness;
     18
     19pos=find(radius<200000);
     20md.thickness(pos)=100;
     21md.bed(pos)=-di*md.thickness(pos)-20;
     22md.surface(pos)=md.bed(pos)+md.thickness(pos);
     23
     24
     25pos=find(md.x<.2*1e6 & md.x>-.2*1e6 & md.y>0);
     26md.thickness(pos)=100;
     27md.bed(pos)=-di*md.thickness(pos)-20;
     28md.surface(pos)=md.bed(pos)+md.thickness(pos);
     29
     30pos=find(md.x<.1*1e6 & md.x>-.1*1e6 & md.y<-.5*1e6 & md.y>-.6*1e6);
     31md.thickness(pos)=100;
     32md.bed(pos)=-di*md.thickness(pos)-20;
     33md.surface(pos)=md.bed(pos)+md.thickness(pos);
    1734
    1835%Initial velocity
Note: See TracChangeset for help on using the changeset viewer.