Changeset 26744 for issm/trunk/examples/shakti/runme.m
- Timestamp:
- 12/22/21 10:39:44 (3 years ago)
- Location:
- issm/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
- Property svn:mergeinfo changed
/issm/trunk-jpl merged: 25837-25866,25868-25993,25995-26330,26332-26733,26736-26739,26741
- Property svn:mergeinfo changed
-
issm/trunk/examples/shakti/runme.m
r23189 r26744 1 1 steps=[1:3]; 2 2 3 if any(steps==1) 3 if any(steps==1) 4 4 disp(' Step 1: Mesh'); 5 5 … … 8 8 9 9 save MoulinMesh md 10 end 11 if any(steps==2) 10 end 11 12 if any(steps==2) 12 13 disp(' Step 2: Parameterization'); 13 14 md=loadmodel('MoulinMesh'); … … 21 22 % Change hydrology class to Sommers' SHaKTI model 22 23 md.hydrology=hydrologyshakti(); 23 24 24 25 % Define initial water head such that water pressure is 50% of ice overburden pressure 25 26 md.hydrology.head = 0.5*md.materials.rho_ice/md.materials.rho_freshwater*md.geometry.thickness + md.geometry.base; … … 48 49 49 50 save MoulinParam md; 50 end 51 if any(steps==3); 51 end 52 53 if any(steps==3) 52 54 disp(' Step 3: Solve!'); 53 55 md=loadmodel('MoulinParam'); … … 80 82 81 83 save MoulinTransient md 82 end 84 end
Note:
See TracChangeset
for help on using the changeset viewer.