![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
source:
issm/oecreview/Archive/16554-17801/ISSM-17467-17468.diff
| Last change on this file was 17802, checked in by , 11 years ago | |
|---|---|
| File size: 11.3 KB | |
-
../trunk-jpl/test/NightlyRun/test801.m
1 md=triangle(model(),'../Exp/Square.exp',50000) 2 md=setmask(md,'','') 3 md=parameterize(md,'../Par/ValleyGlacierShelf.py') 4 md=setflowequation(md,'SSA','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,'SSA','all'); 5 md.cluster=generic('name',oshostname(),'np',3); 6 6 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 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; 14 14 15 md=solve(md,TransientSolutionEnum()) 15 md=solve(md,TransientSolutionEnum()); 16 16 17 17 %Fields and tolerances to track changes 18 18 field_names ={'Vx1','Vy1','Vel1','Pressure1','Thickness1','Surface1','MaskIceLevelset1',... 19 19 'Vx2','Vy2','Vel2','Pressure2','Thickness2','Surface2','MaskIceLevelset2',... 20 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3'} 20 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3'}; 21 21 field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,... 22 22 1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,... 23 1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12} 23 1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12}; 24 24 field_values={... 25 25 md.results.TransientSolution(1).Vx,... 26 26 md.results.TransientSolution(1).Vy,... … … 43 43 md.results.TransientSolution(3).Thickness,... 44 44 md.results.TransientSolution(3).Surface,... 45 45 md.results.TransientSolution(3).MaskIceLevelset,... 46 } 46 }; -
../trunk-jpl/test/NightlyRun/test802.m
1 md=triangle(model(),'../Exp/Square.exp',50000) 2 md=setmask(md,'','') 3 md=parameterize(md,'../Par/ValleyGlacierShelf.py') 4 md=extrude(md,3,2.) 5 md=setflowequation(md,'SSA','all') 6 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=extrude(md,3,2.); 5 md=setflowequation(md,'SSA','all'); 6 md.cluster=generic('name',oshostname(),'np',3); 7 7 8 8 %Thermal model 9 pos_surf=find(md.mesh.vertexonsurface) 10 md.thermal.spctemperature(pos_surf)=md.initialization.temperature(pos_surf) 9 pos_surf=find(md.mesh.vertexonsurface); 10 md.thermal.spctemperature(pos_surf)=md.initialization.temperature(pos_surf); 11 11 12 12 %Transient 13 md.transient.isstressbalance=1 14 md.transient.islevelset=1 15 md.transient.ismasstransport=1 16 md.transient.isthermal=1 17 md.transient.isgroundingline=1 18 md.transient.isgia=0 13 md.transient.isstressbalance=1; 14 md.transient.islevelset=1; 15 md.transient.ismasstransport=1; 16 md.transient.isthermal=1; 17 md.transient.isgroundingline=1; 18 md.transient.isgia=0; 19 19 20 md=solve(md,TransientSolutionEnum()) 20 md=solve(md,TransientSolutionEnum()); 21 21 22 22 %Fields and tolerances to track changes 23 23 field_names ={'Vx1','Vy1','Vel1','Pressure1','Thickness1','Surface1','MaskIceLevelset1','Temperature1',... 24 24 'Vx2','Vy2','Vel2','Pressure2','Thickness2','Surface2','MaskIceLevelset2','Temperature2',... 25 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3','Temperature3'} 25 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3','Temperature3'}; 26 26 field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,... 27 27 1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... 28 1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11} 28 1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11}; 29 29 field_values={... 30 30 md.results.TransientSolution(1).Vx,... 31 31 md.results.TransientSolution(1).Vy,... … … 51 51 md.results.TransientSolution(3).Surface,... 52 52 md.results.TransientSolution(3).MaskIceLevelset,... 53 53 md.results.TransientSolution(3).Temperature,... 54 } 54 }; -
../trunk-jpl/test/NightlyRun/test803.m
1 md=triangle(model(),'../Exp/Square.exp',50000) 2 md=setmask(md,'','') 3 md=parameterize(md,'../Par/ValleyGlacierShelf.py') 4 md=extrude(md,3,2.) 5 md=setflowequation(md,'HO','all') 6 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=extrude(md,3,2.); 5 md=setflowequation(md,'HO','all'); 6 md.cluster=generic('name',oshostname(),'np',3); 7 7 8 8 %Thermal model 9 pos_surf=find(md.mesh.vertexonsurface) 10 md.thermal.spctemperature(pos_surf)=md.initialization.temperature(pos_surf) 11 md.thermal.isenthalpy=1 9 pos_surf=find(md.mesh.vertexonsurface); 10 md.thermal.spctemperature(pos_surf)=md.initialization.temperature(pos_surf); 11 md.thermal.isenthalpy=1; 12 12 13 13 %Transient 14 md.transient.isstressbalance=1 15 md.transient.islevelset=1 16 md.transient.ismasstransport=1 17 md.transient.isthermal=1 18 md.transient.isgroundingline=1 19 md.transient.isgia=0 14 md.transient.isstressbalance=1; 15 md.transient.islevelset=1; 16 md.transient.ismasstransport=1; 17 md.transient.isthermal=1; 18 md.transient.isgroundingline=1; 19 md.transient.isgia=0; 20 20 21 md=solve(md,TransientSolutionEnum()) 21 md=solve(md,TransientSolutionEnum()); 22 22 23 23 %Fields and tolerances to track changes 24 24 field_names ={'Vx1','Vy1','Vel1','Pressure1','Thickness1','Surface1','MaskIceLevelset1','Enthalpy1','Watercolumn1',... 25 25 'Vx2','Vy2','Vel2','Pressure2','Thickness2','Surface2','MaskIceLevelset2','Enthalpy2','Watercolumn2',... 26 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3','Enthalpy3','Watercolumn3'} 26 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3','Enthalpy3','Watercolumn3'}; 27 27 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,... 28 28 1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... 29 1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10} 29 1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10}; 30 30 field_values={... 31 31 md.results.TransientSolution(1).Vx,... 32 32 md.results.TransientSolution(1).Vy,... … … 55 55 md.results.TransientSolution(3).MaskIceLevelset,... 56 56 md.results.TransientSolution(3).Enthalpy,... 57 57 md.results.TransientSolution(3).Watercolumn,... 58 } 58 }; -
../trunk-jpl/test/NightlyRun/test804.m
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,... 28 28 md.results.TransientSolution(1).Vy,... … … 45 45 md.results.TransientSolution(3).Thickness,... 46 46 md.results.TransientSolution(3).Surface,... 47 47 md.results.TransientSolution(3).MaskIceLevelset,... 48 } 48 }; -
../trunk-jpl/test/NightlyRun/test805.m
1 md=triangle(model(),'../Exp/Square.exp',50000) 2 md=setmask(md,'','') 3 md=parameterize(md,'../Par/ValleyGlacierShelf.py') 4 md=extrude(md,3,2.) 5 md=setflowequation(md,'HO','all') 6 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=extrude(md,3,2.); 5 md=setflowequation(md,'HO','all'); 6 md.cluster=generic('name',oshostname(),'np',3); 7 7 8 8 %Thermal model 9 pos_surf=find(md.mesh.vertexonsurface) 10 md.thermal.spctemperature(pos_surf)=md.initialization.temperature(pos_surf) 11 md.thermal.isenthalpy=1 9 pos_surf=find(md.mesh.vertexonsurface); 10 md.thermal.spctemperature(pos_surf)=md.initialization.temperature(pos_surf); 11 md.thermal.isenthalpy=1; 12 12 13 13 %Transient 14 md.transient.isstressbalance=1 15 md.transient.islevelset=1 16 md.transient.ismasstransport=1 17 md.transient.isthermal=1 18 md.transient.isgroundingline=1 19 md.transient.isgia=0 14 md.transient.isstressbalance=1; 15 md.transient.islevelset=1; 16 md.transient.ismasstransport=1; 17 md.transient.isthermal=1; 18 md.transient.isgroundingline=1; 19 md.transient.isgia=0; 20 20 21 md.masstransport.calvingrate=1000.*ones(md.mesh.numberofvertices,1) 21 md.masstransport.calvingrate=1000.*ones(md.mesh.numberofvertices,1); 22 22 23 md=solve(md,TransientSolutionEnum()) 23 md=solve(md,TransientSolutionEnum()); 24 24 25 25 %Fields and tolerances to track changes 26 26 field_names ={'Vx1','Vy1','Vel1','Pressure1','Thickness1','Surface1','MaskIceLevelset1','Enthalpy1','Watercolumn1',... 27 27 'Vx2','Vy2','Vel2','Pressure2','Thickness2','Surface2','MaskIceLevelset2','Enthalpy2','Watercolumn2',... 28 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3','Enthalpy3','Watercolumn3'} 28 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3','Enthalpy3','Watercolumn3'}; 29 29 field_tolerances={1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,... 30 30 1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... 31 1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10} 31 1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10}; 32 32 field_values={... 33 33 md.results.TransientSolution(1).Vx,... 34 34 md.results.TransientSolution(1).Vy,... … … 57 57 md.results.TransientSolution(3).MaskIceLevelset,... 58 58 md.results.TransientSolution(3).Enthalpy,... 59 59 md.results.TransientSolution(3).Watercolumn,... 60 } 60 };
Note:
See TracBrowser
for help on using the repository browser.
