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
RevLine 
[13670]1md=triangle(model(),'../Exp/Square.exp',150000.);
[11591]2md=setmask(md,'../Exp/SquareShelf.exp','');
3md=parameterize(md,'../Par/SquareSheetShelf.par');
[13671]4md=extrude(md,3,2.);
[15565]5md=setflowequation(md,'HO','all');
[11591]6md.cluster=generic('name',oshostname(),'np',3);
[13681]7md.timestepping.time_step=0.;
[11591]8md.thermal.isenthalpy=1;
9md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
[16039]10md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1);
[13663]11md=solve(md,SteadystateSolutionEnum());
[11591]12
13%Fields and tolerances to track changes
[11604]14field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','Water fraction','Enthalpy'};
[15607]15field_tolerances={2e-09,1e-09,1e-09,1e-09,1e-13,1e-10,1e-10,1e-10};
[11591]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),...
[11604]23 (md.results.SteadystateSolution.Waterfraction),...
24 (md.results.SteadystateSolution.Enthalpy),...
[11591]25 };
Note: See TracBrowser for help on using the repository browser.