- Timestamp:
- 05/17/10 08:06:13 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/Verification/test12_IceSheetNoIceFrontH2d/Square.par
r2331 r3757 16 16 17 17 md.drag_type=2; %0 none 1 plastic 2 viscous 18 md.drag =200*ones(md.numberofgrids,1); %q=1.18 md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1. 19 19 %Take care of iceshelves: no basal drag 20 20 pos=find(md.elementoniceshelf); 21 md.drag (md.elements(pos,:))=0;22 md. p=ones(md.numberofelements,1);23 md. q=ones(md.numberofelements,1);21 md.drag_coefficient(md.elements(pos,:))=0; 22 md.drag_p=ones(md.numberofelements,1); 23 md.drag_q=ones(md.numberofelements,1); 24 24 md.viscosity_overshoot=0.3; 25 25 26 26 md.observed_temperature=(273-20)*ones(md.numberofgrids,1); 27 27 28 md. B=paterson(md.observed_temperature);29 md. n=3*ones(md.numberofelements,1);28 md.rheology_B=paterson(md.observed_temperature); 29 md.rheology_n=3*ones(md.numberofelements,1); 30 30 31 31 %Deal with boundary conditions:
Note:
See TracChangeset
for help on using the changeset viewer.