source: issm/oecreview/Archive/19101-20495/ISSM-19337-19338.diff@ 20498

Last change on this file since 20498 was 20498, checked in by Mathieu Morlighem, 9 years ago

CHG: done with Archive/19101-20495
File size: 15.8 KB

  • ../trunk-jpl/test/NightlyRun/test125.m

     
     1%Test Name: SquareShelfConstrainedRestartTranSSA2d
     2md=triangle(model(),'../Exp/Square.exp',150000.);
     3md=setmask(md,'all','');
     4md=parameterize(md,'../Par/SquareShelfConstrained.par');
     5md=setflowequation(md,'SSA','all');
     6md.cluster=generic('name',oshostname(),'np',1);
     7md.transient.requested_outputs={'IceVolume','TotalSmb'};
     8
     9md.verbose=verbose('solution',true);
     10md.settings.recording_frequency=5;
     11
     12% time steps and resolution
     13md.timestepping.final_time=8;
     14
     15md=solve(md,TransientSolutionEnum());
     16md2=solve(md,TransientSolutionEnum(),'restart',1);
     17
     18%Fields and tolerances to track changes
     19field_names     ={'Vx1','Vy1','Vel1','TotalSmb1','Bed1','Surface1','Thickness1','Volume1','Vx2','Vy2','Vel2','TotalSmb2','Bed2','Surface2','Thickness2','Volume2','Vx3','Vy3','Vel3','TotalSmb3','Bed3','Surface3','Thickness3','Volume3'};
     20field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
     21                                                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
     22                                                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     23field_values={...
     24        (md.results.TransientSolution(6).Vx)-(md2.results.TransientSolution(1).Vx),...
     25        (md.results.TransientSolution(6).Vy)-(md2.results.TransientSolution(1).Vy),...
     26        (md.results.TransientSolution(6).Vel)-(md2.results.TransientSolution(1).Vel),...
     27        (md.results.TransientSolution(6).TotalSmb)-(md2.results.TransientSolution(1).TotalSmb),...
     28        (md.results.TransientSolution(6).Base)-(md2.results.TransientSolution(1).Base),...
     29        (md.results.TransientSolution(6).Surface)-(md2.results.TransientSolution(1).Surface),...
     30        (md.results.TransientSolution(6).Thickness)-(md2.results.TransientSolution(1).Thickness),...
     31        (md.results.TransientSolution(6).IceVolume)-(md2.results.TransientSolution(1).IceVolume),...
     32        (md.results.TransientSolution(7).Vx)-(md2.results.TransientSolution(2).Vx),...
     33        (md.results.TransientSolution(7).Vy)-(md2.results.TransientSolution(2).Vy),...
     34        (md.results.TransientSolution(7).Vel)-(md2.results.TransientSolution(2).Vel),...
     35        (md.results.TransientSolution(7).TotalSmb)-(md2.results.TransientSolution(2).TotalSmb),...
     36        (md.results.TransientSolution(7).Base)-(md2.results.TransientSolution(2).Base),...
     37        (md.results.TransientSolution(7).Surface)-(md2.results.TransientSolution(2).Surface),...
     38        (md.results.TransientSolution(7).Thickness)-(md2.results.TransientSolution(2).Thickness),...
     39        (md.results.TransientSolution(7).IceVolume)-(md2.results.TransientSolution(2).IceVolume),...
     40        (md.results.TransientSolution(8).Vx)-(md2.results.TransientSolution(3).Vx),...
     41        (md.results.TransientSolution(8).Vy)-(md2.results.TransientSolution(3).Vy),...
     42        (md.results.TransientSolution(8).Vel)-(md2.results.TransientSolution(3).Vel),...
     43        (md.results.TransientSolution(8).TotalSmb)-(md2.results.TransientSolution(3).TotalSmb),...
     44        (md.results.TransientSolution(8).Base)-(md2.results.TransientSolution(3).Base),...
     45        (md.results.TransientSolution(8).Surface)-(md2.results.TransientSolution(3).Surface),...
     46        (md.results.TransientSolution(8).Thickness)-(md2.results.TransientSolution(3).Thickness),...
     47        (md.results.TransientSolution(8).IceVolume)-(md2.results.TransientSolution(3).IceVolume),...
     48        };

  • ../trunk-jpl/test/NightlyRun/test126.m

     
     1%Test Name: SquareShelfConstrainedRestartTranHO3d
     2md=triangle(model(),'../Exp/Square.exp',150000.);
     3md=setmask(md,'all','');
     4md=parameterize(md,'../Par/SquareShelfConstrained.par');
     5md.cluster=generic('name',oshostname(),'np',1);
     6md=extrude(md,3,1.);
     7md=setflowequation(md,'HO','all');
     8md.transient.requested_outputs={'IceVolume','TotalSmb'};
     9
     10md.verbose=verbose('solution',true);
     11md.settings.recording_frequency=5;
     12
     13% time steps and resolution
     14md.timestepping.final_time=8;
     15
     16md=solve(md,TransientSolutionEnum());
     17md2=solve(md,TransientSolutionEnum(),'restart',1);
     18
     19%Fields and tolerances to track changes
     20field_names     ={'Vx1','Vy1','Vel1','TotalSmb1','Bed1','Surface1','Thickness1','Volume1','Temperature1','Pressure1','Vx2','Vy2','Vel2','TotalSmb2','Bed2','Surface2','Thickness2','Volume2','Temperature2','Pressure2','Vx3','Vy3','Vel3','TotalSmb3','Bed3','Surface3','Thickness3','Volume3','Temperature3','Pressure3'};
     21field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
     22                                                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
     23                                                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
     24                                                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     25field_values={...
     26        (md.results.TransientSolution(6).Vx)-(md2.results.TransientSolution(1).Vx),...
     27        (md.results.TransientSolution(6).Vy)-(md2.results.TransientSolution(1).Vy),...
     28        (md.results.TransientSolution(6).Vel)-(md2.results.TransientSolution(1).Vel),...
     29        (md.results.TransientSolution(6).TotalSmb)-(md2.results.TransientSolution(1).TotalSmb),...
     30        (md.results.TransientSolution(6).Base)-(md2.results.TransientSolution(1).Base),...
     31        (md.results.TransientSolution(6).Surface)-(md2.results.TransientSolution(1).Surface),...
     32        (md.results.TransientSolution(6).Thickness)-(md2.results.TransientSolution(1).Thickness),...
     33        (md.results.TransientSolution(6).IceVolume)-(md2.results.TransientSolution(1).IceVolume),...
     34        (md.results.TransientSolution(6).Temperature)-(md2.results.TransientSolution(1).Temperature),...
     35        (md.results.TransientSolution(6).Pressure)-(md2.results.TransientSolution(1).Pressure),...
     36        (md.results.TransientSolution(7).Vx)-(md2.results.TransientSolution(2).Vx),...
     37        (md.results.TransientSolution(7).Vy)-(md2.results.TransientSolution(2).Vy),...
     38        (md.results.TransientSolution(7).Vel)-(md2.results.TransientSolution(2).Vel),...
     39        (md.results.TransientSolution(7).TotalSmb)-(md2.results.TransientSolution(2).TotalSmb),...
     40        (md.results.TransientSolution(7).Base)-(md2.results.TransientSolution(2).Base),...
     41        (md.results.TransientSolution(7).Surface)-(md2.results.TransientSolution(2).Surface),...
     42        (md.results.TransientSolution(7).Thickness)-(md2.results.TransientSolution(2).Thickness),...
     43        (md.results.TransientSolution(7).IceVolume)-(md2.results.TransientSolution(2).IceVolume),...
     44        (md.results.TransientSolution(7).Temperature)-(md2.results.TransientSolution(2).Temperature),...
     45        (md.results.TransientSolution(7).Pressure)-(md2.results.TransientSolution(2).Pressure),...
     46        (md.results.TransientSolution(8).Vx)-(md2.results.TransientSolution(3).Vx),...
     47        (md.results.TransientSolution(8).Vy)-(md2.results.TransientSolution(3).Vy),...
     48        (md.results.TransientSolution(8).Vel)-(md2.results.TransientSolution(3).Vel),...
     49        (md.results.TransientSolution(8).TotalSmb)-(md2.results.TransientSolution(3).TotalSmb),...
     50        (md.results.TransientSolution(8).Base)-(md2.results.TransientSolution(3).Base),...
     51        (md.results.TransientSolution(8).Surface)-(md2.results.TransientSolution(3).Surface),...
     52        (md.results.TransientSolution(8).Thickness)-(md2.results.TransientSolution(3).Thickness),...
     53        (md.results.TransientSolution(8).IceVolume)-(md2.results.TransientSolution(3).IceVolume),...
     54        (md.results.TransientSolution(8).Temperature)-(md2.results.TransientSolution(3).Temperature),...
     55        (md.results.TransientSolution(8).Pressure)-(md2.results.TransientSolution(3).Pressure),...
     56        };

  • ../trunk-jpl/test/NightlyRun/test125.py

     
     1#Test Name: SquareShelfConstrainedRestartTranSSA2d
     2from model import *
     3from triangle import *
     4from setmask import *
     5from parameterize import *
     6from setflowequation import *
     7from EnumDefinitions import *
     8from solve import *
     9from MatlabFuncs import *
     10from massfluxatgate import massfluxatgate
     11from generic import generic
     12import copy
     13
     14md=triangle(model(),'../Exp/Square.exp',150000.)
     15md=setmask(md,'all','')
     16md=parameterize(md,'../Par/SquareShelfConstrained.py')
     17md=setflowequation(md,'SSA','all')
     18md.cluster=generic('name',oshostname(),'np',1)
     19md.transient.requested_outputs=['IceVolume','TotalSmb']
     20
     21md.verbose=verbose('solution',1);
     22md.settings.recording_frequency=5
     23
     24# time steps and resolution
     25md.timestepping.final_time=8
     26
     27md=solve(md,TransientSolutionEnum())
     28md2=copy.deepcopy(md)
     29md=solve(md,TransientSolutionEnum(),'restart',1)
     30
     31#Fields and tolerances to track changes
     32field_names     =['Vx1','Vy1','Vel1','TotalSmb1','Bed1','Surface1','Thickness1','Volume1','Vx2','Vy2','Vel2','TotalSmb2','Bed2','Surface2','Thickness2','Volume2','Vx3','Vy3','Vel3','TotalSmb3','Bed3','Surface3','Thickness3','Volume3'];
     33field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
     34                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
     35                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
     36                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
     37field_values=[\
     38                md2.results.TransientSolution[5].Vx-md.results.TransientSolution[0].Vx,\
     39                md2.results.TransientSolution[5].Vy-md.results.TransientSolution[0].Vy,\
     40                md2.results.TransientSolution[5].Vel-md.results.TransientSolution[0].Vel,\
     41                md2.results.TransientSolution[5].TotalSmb-md.results.TransientSolution[0].TotalSmb,\
     42                md2.results.TransientSolution[5].Base-md.results.TransientSolution[0].Base,\
     43                md2.results.TransientSolution[5].Surface-md.results.TransientSolution[0].Surface,\
     44                md2.results.TransientSolution[5].Thickness-md.results.TransientSolution[0].Thickness,\
     45                md2.results.TransientSolution[5].IceVolume-md.results.TransientSolution[0].IceVolume,\
     46                md2.results.TransientSolution[6].Vx-md.results.TransientSolution[1].Vx,\
     47                md2.results.TransientSolution[6].Vy-md.results.TransientSolution[1].Vy,\
     48                md2.results.TransientSolution[6].Vel-md.results.TransientSolution[1].Vel,\
     49                md2.results.TransientSolution[6].TotalSmb-md.results.TransientSolution[1].TotalSmb,\
     50                md2.results.TransientSolution[6].Base-md.results.TransientSolution[1].Base,\
     51                md2.results.TransientSolution[6].Surface-md.results.TransientSolution[1].Surface,\
     52                md2.results.TransientSolution[6].Thickness-md.results.TransientSolution[1].Thickness,\
     53                md2.results.TransientSolution[6].IceVolume-md.results.TransientSolution[1].IceVolume,\
     54                md2.results.TransientSolution[7].Vx-md.results.TransientSolution[2].Vx,\
     55                md2.results.TransientSolution[7].Vy-md.results.TransientSolution[2].Vy,\
     56                md2.results.TransientSolution[7].Vel-md.results.TransientSolution[2].Vel,\
     57                md2.results.TransientSolution[7].TotalSmb-md.results.TransientSolution[2].TotalSmb,\
     58                md2.results.TransientSolution[7].Base-md.results.TransientSolution[2].Base,\
     59                md2.results.TransientSolution[7].Surface-md.results.TransientSolution[2].Surface,\
     60                md2.results.TransientSolution[7].Thickness-md.results.TransientSolution[2].Thickness,\
     61                md2.results.TransientSolution[7].IceVolume-md.results.TransientSolution[2].IceVolume\
     62                ]
     63

  • ../trunk-jpl/test/NightlyRun/test126.py

     
     1#Test Name: SquareShelfConstrainedRestartTranHO3d
     2from model import *
     3from triangle import *
     4from setmask import *
     5from parameterize import *
     6from setflowequation import *
     7from EnumDefinitions import *
     8from solve import *
     9from MatlabFuncs import *
     10from massfluxatgate import massfluxatgate
     11from generic import generic
     12import copy
     13
     14md=triangle(model(),'../Exp/Square.exp',150000.)
     15md=setmask(md,'all','')
     16md=parameterize(md,'../Par/SquareShelfConstrained.py')
     17md.cluster=generic('name',oshostname(),'np',1)
     18md.transient.requested_outputs=['IceVolume','TotalSmb']
     19md.extrude(3,1.)
     20md=setflowequation(md,'HO','all')
     21
     22md.verbose=verbose('solution',1);
     23md.settings.recording_frequency=5
     24
     25# time steps and resolution
     26md.timestepping.final_time=8
     27
     28md=solve(md,TransientSolutionEnum())
     29md2=copy.deepcopy(md)
     30md=solve(md,TransientSolutionEnum(),'restart',1)
     31
     32#Fields and tolerances to track changes
     33field_names     =['Vx1','Vy1','Vel1','TotalSmb1','Bed1','Surface1','Thickness1','Volume1','Temperature1','Pressure1','Vx2','Vy2','Vel2','TotalSmb2','Bed2','Surface2','Thickness2','Volume2','Temperature2','Pressure2','Vx3','Vy3','Vel3','TotalSmb3','Bed3','Surface3','Thickness3','Volume3','Temperature3','Pressure3'];
     34field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
     35                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
     36                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
     37                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
     38                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
     39field_values=[\
     40                md2.results.TransientSolution[5].Vx-md.results.TransientSolution[0].Vx,\
     41                md2.results.TransientSolution[5].Vy-md.results.TransientSolution[0].Vy,\
     42                md2.results.TransientSolution[5].Vel-md.results.TransientSolution[0].Vel,\
     43                md2.results.TransientSolution[5].TotalSmb-md.results.TransientSolution[0].TotalSmb,\
     44                md2.results.TransientSolution[5].Base-md.results.TransientSolution[0].Base,\
     45                md2.results.TransientSolution[5].Surface-md.results.TransientSolution[0].Surface,\
     46                md2.results.TransientSolution[5].Thickness-md.results.TransientSolution[0].Thickness,\
     47                md2.results.TransientSolution[5].IceVolume-md.results.TransientSolution[0].IceVolume,\
     48                md2.results.TransientSolution[5].Temperature-md.results.TransientSolution[0].Temperature,\
     49                md2.results.TransientSolution[5].Pressure-md.results.TransientSolution[0].Pressure,\
     50                md2.results.TransientSolution[6].Vx-md.results.TransientSolution[1].Vx,\
     51                md2.results.TransientSolution[6].Vy-md.results.TransientSolution[1].Vy,\
     52                md2.results.TransientSolution[6].Vel-md.results.TransientSolution[1].Vel,\
     53                md2.results.TransientSolution[6].TotalSmb-md.results.TransientSolution[1].TotalSmb,\
     54                md2.results.TransientSolution[6].Base-md.results.TransientSolution[1].Base,\
     55                md2.results.TransientSolution[6].Surface-md.results.TransientSolution[1].Surface,\
     56                md2.results.TransientSolution[6].Thickness-md.results.TransientSolution[1].Thickness,\
     57                md2.results.TransientSolution[6].IceVolume-md.results.TransientSolution[1].IceVolume,\
     58                md2.results.TransientSolution[6].Temperature-md.results.TransientSolution[1].Temperature,\
     59                md2.results.TransientSolution[6].Pressure-md.results.TransientSolution[1].Pressure,\
     60                md2.results.TransientSolution[7].Vx-md.results.TransientSolution[2].Vx,\
     61                md2.results.TransientSolution[7].Vy-md.results.TransientSolution[2].Vy,\
     62                md2.results.TransientSolution[7].Vel-md.results.TransientSolution[2].Vel,\
     63                md2.results.TransientSolution[7].TotalSmb-md.results.TransientSolution[2].TotalSmb,\
     64                md2.results.TransientSolution[7].Base-md.results.TransientSolution[2].Base,\
     65                md2.results.TransientSolution[7].Surface-md.results.TransientSolution[2].Surface,\
     66                md2.results.TransientSolution[7].Thickness-md.results.TransientSolution[2].Thickness,\
     67                md2.results.TransientSolution[7].IceVolume-md.results.TransientSolution[2].IceVolume,\
     68                md2.results.TransientSolution[7].Temperature-md.results.TransientSolution[2].Temperature,\
     69                md2.results.TransientSolution[7].Pressure-md.results.TransientSolution[2].Pressure\
     70                ]
     71
Note: See TracBrowser for help on using the repository browser.