source: issm/trunk-jpl/test/NightlyRun/test142.m@ 11281

Last change on this file since 11281 was 11281, checked in by seroussi, 13 years ago

added NR for enthalpy solution

File size: 1.1 KB
Line 
1md=triangle(model,'../Exp/Square.exp',180000);
2md=setmask(md,'all','');
3md=parameterize(md,'../Par/SquareShelfConstrained.par');
4md=extrude(md,3,1);
5md=setflowequation(md,'macayeal','all');
6md.cluster=generic('name',oshostname(),'np',3);
7md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
8%md=solve(md,ThermalSolutionEnum);
9md=solve(md,EnthalpySolutionEnum);
10
11%Fields and tolerances to track changes
12field_names ={'Enthalpy1','Waterfraction1','Temperature1',...
13 'Enthalpy2','Waterfraction2','Temperature2',...
14 'Enthalpy3','Waterfraction3','Temperature3'};
15field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
16field_values={...
17 (md.results.EnthalpySolution(1).Enthalpy),...
18 (md.results.EnthalpySolution(1).Waterfraction),...
19 (md.results.EnthalpySolution(1).Temperature),...
20 (md.results.EnthalpySolution(2).Enthalpy),...
21 (md.results.EnthalpySolution(2).Waterfraction),...
22 (md.results.EnthalpySolution(2).Temperature),...
23 (md.results.EnthalpySolution(3).Enthalpy),...
24 (md.results.EnthalpySolution(3).Waterfraction),...
25 (md.results.EnthalpySolution(3).Temperature),...
26 };
Note: See TracBrowser for help on using the repository browser.