source: issm/trunk-jpl/test/NightlyRun/test432.m@ 16039

Last change on this file since 16039 was 16039, checked in by seroussi, 12 years ago

CHG: Need watercolumn initial values to run enthalpy

File size: 1.0 KB
Line 
1md=triangle(model(),'../Exp/Square.exp',150000.);
2md=setmask(md,'../Exp/SquareShelf.exp','');
3md=parameterize(md,'../Par/SquareSheetShelf.par');
4md=extrude(md,3,2.);
5md=setflowequation(md,'HO','all');
6md.cluster=generic('name',oshostname(),'np',3);
7md.timestepping.time_step=0.;
8md.thermal.isenthalpy=1;
9md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
10md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1);
11md=solve(md,SteadystateSolutionEnum());
12
13%Fields and tolerances to track changes
14field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','Water fraction','Enthalpy'};
15field_tolerances={2e-09,1e-09,1e-09,1e-09,1e-13,1e-10,1e-10,1e-10};
16field_values={...
17 (md.results.SteadystateSolution.Vx),...
18 (md.results.SteadystateSolution.Vy),...
19 (md.results.SteadystateSolution.Vz),...
20 (md.results.SteadystateSolution.Vel),...
21 (md.results.SteadystateSolution.Pressure),...
22 (md.results.SteadystateSolution.Temperature),...
23 (md.results.SteadystateSolution.Waterfraction),...
24 (md.results.SteadystateSolution.Enthalpy),...
25 };
Note: See TracBrowser for help on using the repository browser.