Changeset 14805
- Timestamp:
- 04/30/13 08:56:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test332.m
r14772 r14805 1 md=triangle(model(),'../Exp/Square.exp',10000 .);1 md=triangle(model(),'../Exp/Square.exp',100000.); 2 2 md=setmask(md,'',''); 3 3 md=parameterize(md,'../Par/SquareSheetConstrained.par'); … … 6 6 md.hydrology=(hydrologydc); 7 7 md.hydrology.isefficientlayer=0; 8 md.hydrology.sedimentlimit_flag=1; 9 md.hydrology.sedimentlimit=8000.0; 8 10 md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1); 9 11 md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1); … … 18 20 19 21 %Fields and tolerances to track changes 20 %you cqn also compqre with an analitic solution 22 %you can also compare with an analitic solution, but it is exact 23 %only if no limits are applied 21 24 %analitic=(md.mesh.y.^2-2*md.mesh.y*1.0e6)*(-2.0/(2*md.constants.yts*md.hydrology.sediment_transmitivity)) 22 25 field_names ={'SedimentWaterHead'};
Note:
See TracChangeset
for help on using the changeset viewer.