Changeset 7359
- Timestamp:
- 02/05/11 20:56:49 (14 years ago)
- Location:
- issm/trunk/test
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/NightlyRun/test446.m
r7075 r7359 2 2 shelfextent=2e5; 3 3 resolution=.5e5; 4 5 4 %mesh {{{1 6 5 md=roundmesh(model,radius,resolution); … … 23 22 md=setelementstype(md,'macayeal','all'); 24 23 md=SetParallel(md,3); 25 md=solve(md,'analysis_type',Transient2DSolutionEnum); 24 %plug holes into the ice sheet, to test for grounding line migration. {{{1 25 di=md.rho_ice/md.rho_water; 26 rad=sqrt((md.x).*md.x+(md.y).*md.y); 27 pos=find(rad<200000); 28 md.thickness(pos)=100; 29 md.bed(pos)=-di*md.thickness(pos)-20; 30 md.surface(pos)=md.bed(pos)+md.thickness(pos); 31 32 pos=find(md.x<.2*1e6 & md.x>-.2*1e6 & md.y>0); 33 md.thickness(pos)=100; 34 md.bed(pos)=-di*md.thickness(pos)-20; 35 md.surface(pos)=md.bed(pos)+md.thickness(pos); 36 37 pos=find(md.x<.1*1e6 & md.x>-.1*1e6 & md.y<-.5*1e6 & md.y>-.6*1e6); 38 md.thickness(pos)=100; 39 md.bed(pos)=-di*md.thickness(pos)-20; 40 md.surface(pos)=md.bed(pos)+md.thickness(pos); 41 %}}} 42 43 %test different grounding line dynamics. 44 %grounding line: 45 md.gl_migration=AgressiveMigrationEnum; 46 md=solve(md,'analysis_type',GroundingLineMigration2DSolutionEnum); 47 element_on_iceshelf_agressive=PatchToVec(md.results.GroundingLineMigration2DSolution.ElementOnIceShelf); 48 49 md.gl_migration=SoftMigrationEnum; 50 md=solve(md,'analysis_type',GroundingLineMigration2DSolutionEnum); 51 element_on_iceshelf_soft=PatchToVec(md.results.GroundingLineMigration2DSolution.ElementOnIceShelf); 26 52 27 53 %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 }; 54 field_names ={'ElementOnIceShelfAgressive','ElementOnIceShelfSoft'}; 55 field_tolerances={1e-13,1e-13}; 56 field_values={element_on_iceshelf_agressive,element_on_iceshelf_soft}; -
issm/trunk/test/Par/RoundSheetShelf.par
r7329 r7359 1 1 %Start defining model parameters here 2 2 3 di=md.rho_ice/md.rho_water; 3 4 rad=1e6; 4 5 shelfextent=2e5; … … 15 16 md.bed(pos)=md.bed(pos)-300*(radius(pos)-(rad-shelfextent))/(rad-shelfextent); 16 17 md.surface=md.bed+md.thickness; 18 19 pos=find(radius<200000); 20 md.thickness(pos)=100; 21 md.bed(pos)=-di*md.thickness(pos)-20; 22 md.surface(pos)=md.bed(pos)+md.thickness(pos); 23 24 25 pos=find(md.x<.2*1e6 & md.x>-.2*1e6 & md.y>0); 26 md.thickness(pos)=100; 27 md.bed(pos)=-di*md.thickness(pos)-20; 28 md.surface(pos)=md.bed(pos)+md.thickness(pos); 29 30 pos=find(md.x<.1*1e6 & md.x>-.1*1e6 & md.y<-.5*1e6 & md.y>-.6*1e6); 31 md.thickness(pos)=100; 32 md.bed(pos)=-di*md.thickness(pos)-20; 33 md.surface(pos)=md.bed(pos)+md.thickness(pos); 17 34 18 35 %Initial velocity
Note:
See TracChangeset
for help on using the changeset viewer.