Changeset 17198 for issm/trunk-jpl/test/NightlyRun/test333.m
- Timestamp:
- 01/29/14 15:20:11 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test333.m
r17060 r17198 26 26 md=solve(md,HydrologySolutionEnum()); 27 27 28 %store=md.constants.g*md.hydrology.sediment_porosity* ... 29 % md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility) 30 31 %sed=ones(1,size(md.results.HydrologySolution,2)); 32 %epl=ones(1,size(md.results.HydrologySolution,2)); 33 %res=ones(1,size(md.results.HydrologySolution,2)); 34 %input=ones(1,size(md.results.HydrologySolution,2)); 35 %for i= 1:size(md.results.HydrologySolution,2) 36 % sed(i)=mean(md.results.HydrologySolution(i).SedimentHead); 37 % res(i)=mean(md.results.HydrologySolution(i).SedimentHeadResidual); 38 % epl(i)=mean(md.results.HydrologySolution(i).EplHead); 39 % input(i)=2.0*(i*0.2); 28 %store=md.constants.g*md.hydrology.sediment_porosity*md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility); 29 %sedstore=20.0*store; 30 %for i=1:10 31 %diff=(mean(md.results.HydrologySolution(i).EplHead)*store+ ... 32 % mean(md.results.HydrologySolution(i).SedimentHead)*sedstore-0.4*i) 40 33 %end 41 34
Note:
See TracChangeset
for help on using the changeset viewer.