Changeset 14748

Timestamp:

04/24/13 16:24:52 (12 years ago)

Author:

bdef

Message:

CHG: added the one layer flag to hydrology test

Location:

issm/trunk-jpl/test/NightlyRun

Files:

• test332.m (modified) (1 diff)
• tmp332.m (modified) (1 diff)

issm/trunk-jpl/test/NightlyRun/test332.m

@@ -5 +5 @@
 md.cluster=generic('name',oshostname(),'np',2);
 md.hydrology=(hydrologydc);
+md.hydrology.isefficientlayer=0;
 md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1);
 md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1);

issm/trunk-jpl/test/NightlyRun/tmp332.m

@@ -3 +3 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=setflowequation(md,'macayeal','all');
+md.cluster=generic('name',oshostname(),'np',2);
 md.hydrology=(hydrologydc);
-md.cluster=generic('name',oshostname(),'np',3);
 md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1);
 md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1);
-pos=find(md.mesh.y==50000.);
+pos=find(md.mesh.x==0);
 md.hydrology.spcsediment_head(pos)=0.0;
+md.basalforcings.melting_rate = 2.0*ones(md.mesh.numberofvertices,1);
+md.timestepping.time_step=0;
+md.timestepping.final_time=1.0;
 md=solve(md,HydrologySolutionEnum());
+
+analitic=(md.mesh.x.^2-2*md.mesh.x*10^6)*(-2.0/(365*3600*24*2.0*md.hydrology.sediment_transmitivity));
+plotmodel(md,'data',md.results.HydrologySolution.SedimentHead-analitic)
+%plotmodel(md,'data',analitic)