source: issm/oecreview/Archive/16554-17801/ISSM-17197-17198.diff

../trunk-jpl/test/NightlyRun/test333.m

@@ -25 +25 @@
 
 md=solve(md,HydrologySolutionEnum());
 
-%store=md.constants.g*md.hydrology.sediment_porosity* ...
-%                       md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility)
-
-%sed=ones(1,size(md.results.HydrologySolution,2));
-%epl=ones(1,size(md.results.HydrologySolution,2));
-%res=ones(1,size(md.results.HydrologySolution,2));
-%input=ones(1,size(md.results.HydrologySolution,2));
-%for i= 1:size(md.results.HydrologySolution,2)
-%       sed(i)=mean(md.results.HydrologySolution(i).SedimentHead);
-%       res(i)=mean(md.results.HydrologySolution(i).SedimentHeadResidual);
-%       epl(i)=mean(md.results.HydrologySolution(i).EplHead);
-%       input(i)=2.0*(i*0.2);
+%store=md.constants.g*md.hydrology.sediment_porosity*md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility);
+%sedstore=20.0*store;
+%for i=1:10
+%diff=(mean(md.results.HydrologySolution(i).EplHead)*store+ ...
+%                       mean(md.results.HydrologySolution(i).SedimentHead)*sedstore-0.4*i)
 %end
 
 %Fields and tolerances to track changes