Changeset 17441
- Timestamp:
- 03/14/14 15:30:20 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test336.m
r17439 r17441 1 1 md=triangle(model(),'../Exp/Square.exp',100000.); 2 2 md=setmask(md,'',''); 3 md=parameterize(md,'../Par/SquareSheetConstrained.p y');3 md=parameterize(md,'../Par/SquareSheetConstrained.par'); 4 4 md=setflowequation(md,'SSA','all'); 5 5 md.cluster=generic('name',oshostname(),'np',3); … … 14 14 15 15 % init levelset function 16 xmin, xmax=min(md.mesh.x), max(md.mesh.x) 17 ymin, ymax=min(md.mesh.y), max(md.mesh.y) 18 xmed, ymed=(xmax+xmin)/2.,(ymax+ymin)/2. 19 20 md.mask.ice_levelset=md.mesh.y - ymed 16 ymed = (max(md.mesh.y)+min(md.mesh.y))/2; 17 md.mask.ice_levelset=md.mesh.y - ymed; 21 18 22 19 md=solve(md,TransientSolutionEnum()) 23 20 24 21 %Fields and tolerances to track changes 25 field_names ={'Vx','Vy','Vel','Pressure','Thickness','Surface','MaskIceLevelset'} 26 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\ 27 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\ 22 field_names ={'Vx1','Vy1','Vel1','Pressure1','Thickness1','Surface1','MaskIceLevelset1',... 23 'Vx2','Vy2','Vel2','Pressure2','Thickness2','Surface2','MaskIceLevelset2',... 24 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3'} 25 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... 26 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... 28 27 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13} 29 field_values={ \30 md.results.TransientSolution [0].Vx,\31 md.results.TransientSolution [0].Vy,\32 md.results.TransientSolution [0].Vel,\33 md.results.TransientSolution [0].Pressure,\34 md.results.TransientSolution [0].Thickness,\35 md.results.TransientSolution [0].Surface,\36 md.results.TransientSolution [0].MaskIceLevelset,\37 md.results.TransientSolution [1].Vx,\38 md.results.TransientSolution [1].Vy,\39 md.results.TransientSolution [1].Vel,\40 md.results.TransientSolution [1].Pressure,\41 md.results.TransientSolution [1].Thickness,\42 md.results.TransientSolution [1].Surface,\43 md.results.TransientSolution [1].MaskIceLevelset,\44 md.results.TransientSolution [2].Vx,\45 md.results.TransientSolution [2].Vy,\46 md.results.TransientSolution [2].Vel,\47 md.results.TransientSolution [2].Pressure,\48 md.results.TransientSolution [2].Thickness,\49 md.results.TransientSolution [2].Surface,\50 md.results.TransientSolution [2].MaskIceLevelset,\28 field_values={... 29 md.results.TransientSolution(1).Vx,... 30 md.results.TransientSolution(1).Vy,... 31 md.results.TransientSolution(1).Vel,... 32 md.results.TransientSolution(1).Pressure,... 33 md.results.TransientSolution(1).Thickness,... 34 md.results.TransientSolution(1).Surface,... 35 md.results.TransientSolution(1).MaskIceLevelset,... 36 md.results.TransientSolution(2).Vx,... 37 md.results.TransientSolution(2).Vy,... 38 md.results.TransientSolution(2).Vel,... 39 md.results.TransientSolution(2).Pressure,... 40 md.results.TransientSolution(2).Thickness,... 41 md.results.TransientSolution(2).Surface,... 42 md.results.TransientSolution(2).MaskIceLevelset,... 43 md.results.TransientSolution(3).Vx,... 44 md.results.TransientSolution(3).Vy,... 45 md.results.TransientSolution(3).Vel,... 46 md.results.TransientSolution(3).Pressure,... 47 md.results.TransientSolution(3).Thickness,... 48 md.results.TransientSolution(3).Surface,... 49 md.results.TransientSolution(3).MaskIceLevelset... 51 50 }
Note:
See TracChangeset
for help on using the changeset viewer.