Changeset 17468 for issm/trunk-jpl/test/NightlyRun/test804.m
- Timestamp:
- 03/18/14 16:29:26 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test804.m
r17464 r17468 1 md=triangle(model(),'../Exp/Square.exp',50000) 2 md=setmask(md,'','') 3 md=parameterize(md,'../Par/ValleyGlacierShelf.py') 4 md=setflowequation(md,'SIA','all') 5 md.cluster=generic('name',oshostname(),'np',3) 1 md=triangle(model(),'../Exp/Square.exp',50000); 2 md=setmask(md,'',''); 3 md=parameterize(md,'../Par/ValleyGlacierShelf.py'); 4 md=setflowequation(md,'SIA','all'); 5 md.cluster=generic('name',oshostname(),'np',3); 6 6 7 7 %Transient 8 md.transient.isstressbalance=1 9 md.transient.islevelset=1 10 md.transient.ismasstransport=1 11 md.transient.isthermal=0 12 md.transient.isgroundingline=1 13 md.transient.isgia=0 8 md.transient.isstressbalance=1; 9 md.transient.islevelset=1; 10 md.transient.ismasstransport=1; 11 md.transient.isthermal=0; 12 md.transient.isgroundingline=1; 13 md.transient.isgia=0; 14 14 15 md.masstransport.calvingrate=1000.*ones(md.mesh.numberofvertices,1) 15 md.masstransport.calvingrate=1000.*ones(md.mesh.numberofvertices,1); 16 16 17 md=solve(md,TransientSolutionEnum()) 17 md=solve(md,TransientSolutionEnum()); 18 18 19 19 %Fields and tolerances to track changes 20 20 field_names ={'Vx1','Vy1','Vel1','Pressure1','Thickness1','Surface1','MaskIceLevelset1',... 21 21 'Vx2','Vy2','Vel2','Pressure2','Thickness2','Surface2','MaskIceLevelset2',... 22 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3'} 22 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3'}; 23 23 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... 24 24 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... 25 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13} 25 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; 26 26 field_values={... 27 27 md.results.TransientSolution(1).Vx,... … … 46 46 md.results.TransientSolution(3).Surface,... 47 47 md.results.TransientSolution(3).MaskIceLevelset,... 48 } 48 };
Note:
See TracChangeset
for help on using the changeset viewer.