Changeset 22864


Ignore:
Timestamp:
06/22/18 01:02:47 (7 years ago)
Author:
bdef
Message:

BUG: random fixes in nightlies

Location:
issm/trunk-jpl/test/NightlyRun
Files:
11 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/test/NightlyRun/test236.py

    r21408 r22864  
    4040
    4141# creating initialization and spc temperatures initialization and spc
    42 md.thermal.spctemperature=np.mean(md.smb.temperatures_lgm[0:md.mesh.numberofvertices,:],axis=1).reshape(-1,1)    #-10*ones(md.mesh.numberofvertices,1)
    43 md.thermal.spctemperature=np.tile(md.thermal.spctemperature,(1,md.timestepping.final_time/md.timestepping.time_step))
     42md.thermal.spctemperature=np.mean(md.smb.temperatures_lgm[0:md.mesh.numberofvertices,:],axis=1)    #-10*ones(md.mesh.numberofvertices,1)
     43md.thermal.spctemperature=np.tile(md.thermal.spctemperature,(int(md.timestepping.final_time/md.timestepping.time_step),1)).T
    4444itemp=np.arange(0,md.timestepping.final_time,md.timestepping.time_step)
    45 md.thermal.spctemperature=np.vstack((md.thermal.spctemperature,itemp.reshape(1,-1)))
     45md.thermal.spctemperature=np.vstack((md.thermal.spctemperature,itemp))
    4646
    4747md.initialization.temperature=md.smb.temperatures_lgm[0:md.mesh.numberofvertices,0]    #*ones(md.mesh.numberofvertices,1)
    … …  
    5252md.smb.precipitations_lgm=np.zeros((md.mesh.numberofvertices+1,12))
    5353for imonth in xrange(0,12):
    54     md.smb.precipitations_presentday[0:md.mesh.numberofvertices,imonth]=-0.4*10**(-6)*md.mesh.y+0.5
    55     md.smb.precipitations_presentday[md.mesh.numberofvertices,imonth]=((float(imonth)+1.)/12.)
    56     md.smb.precipitations_lgm[0:md.mesh.numberofvertices,imonth]=-0.4*10**(-6)*md.mesh.y+0.5
    57     md.smb.precipitations_lgm[md.mesh.numberofvertices,imonth]=((float(imonth)+1.)/12.)
     54        md.smb.precipitations_presentday[0:md.mesh.numberofvertices,imonth]=-0.4*10**(-6)*md.mesh.y+0.5
     55        md.smb.precipitations_presentday[md.mesh.numberofvertices,imonth]=((float(imonth)+1.)/12.)
     56        md.smb.precipitations_lgm[0:md.mesh.numberofvertices,imonth]=-0.4*10**(-6)*md.mesh.y+0.5
     57        md.smb.precipitations_lgm[md.mesh.numberofvertices,imonth]=((float(imonth)+1.)/12.)
    5858
    5959# Interpolation factors
    60 md.smb.Tdiff=0.5*np.ones((2,md.timestepping.final_time))
    61 md.smb.sealev=0.5*np.ones((2,md.timestepping.final_time))
     60md.smb.Tdiff=0.5*np.ones((2,int(md.timestepping.final_time)))
     61md.smb.sealev=0.5*np.ones((2,int(md.timestepping.final_time)))
    6262# Year of each data point
    63 md.smb.Tdiff[1,1:md.timestepping.final_time]=np.arange(1.,md.timestepping.final_time)
    64 md.smb.sealev[1,1:md.timestepping.final_time]=np.arange(1.,md.timestepping.final_time)
     63md.smb.Tdiff[1,1:int(md.timestepping.final_time)]=np.arange(1.,int(md.timestepping.final_time))
     64md.smb.sealev[1,1:int(md.timestepping.final_time)]=np.arange(1.,int(md.timestepping.final_time))
    6565
    6666# time steps and resolution
    … …  
    6969
    7070
    71 # 
     71#
    7272md.transient.requested_outputs=['default','SmbMonthlytemperatures']
    7373md=setflowequation(md,'SSA','all')
    … …  
    7676
    7777#Fields and tolerances to track changes
    78 field_names     =['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMonthlytemperatures1','SmbMassBalance1',\
    79                             'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMonthlytemperatures2','SmbMassBalance2',\
    80                                       'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMonthlytemperatures3','SmbMassBalance3']
    81 field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
    82                             1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
    83                                  1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
    84 field_values=[\
    85         md.results.TransientSolution[0].Vx,\
    86         md.results.TransientSolution[0].Vy,\
    87         md.results.TransientSolution[0].Vel,\
    88         md.results.TransientSolution[0].Pressure,\
    89         md.results.TransientSolution[0].Base,\
    90         md.results.TransientSolution[0].Surface,\
    91         md.results.TransientSolution[0].Thickness,\
    92         md.results.TransientSolution[0].SmbMonthlytemperatures,\
    93         md.results.TransientSolution[0].SmbMassBalance,\
    94         md.results.TransientSolution[1].Vx,\
    95         md.results.TransientSolution[1].Vy,\
    96         md.results.TransientSolution[1].Vel,\
    97         md.results.TransientSolution[1].Pressure,\
    98         md.results.TransientSolution[1].Base,\
    99         md.results.TransientSolution[1].Surface,\
    100         md.results.TransientSolution[1].Thickness,\
    101         md.results.TransientSolution[1].SmbMonthlytemperatures,\
    102         md.results.TransientSolution[1].SmbMassBalance,\
    103         md.results.TransientSolution[2].Vx,\
    104         md.results.TransientSolution[2].Vy,\
    105         md.results.TransientSolution[2].Vel,\
    106         md.results.TransientSolution[2].Pressure,\
    107         md.results.TransientSolution[2].Base,\
    108         md.results.TransientSolution[2].Surface,\
    109         md.results.TransientSolution[2].Thickness,\
    110         md.results.TransientSolution[2].SmbMonthlytemperatures,\
    111         md.results.TransientSolution[2].SmbMassBalance,\
    112         ]
     78field_names     =['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMonthlytemperatures1','SmbMassBalance1',
     79                                                                        'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMonthlytemperatures2','SmbMassBalance2',
     80                                                                        'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMonthlytemperatures3','SmbMassBalance3']
     81field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,
     82                                                                        1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,
     83                                                                        1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
     84field_values=[md.results.TransientSolution[0].Vx,
     85                                                        md.results.TransientSolution[0].Vy,
     86                                                        md.results.TransientSolution[0].Vel,
     87                                                        md.results.TransientSolution[0].Pressure,
     88                                                        md.results.TransientSolution[0].Base,
     89                                                        md.results.TransientSolution[0].Surface,
     90                                                        md.results.TransientSolution[0].Thickness,
     91                                                        md.results.TransientSolution[0].SmbMonthlytemperatures,
     92                                                        md.results.TransientSolution[0].SmbMassBalance,
     93                                                        md.results.TransientSolution[1].Vx,
     94                                                        md.results.TransientSolution[1].Vy,
     95                                                        md.results.TransientSolution[1].Vel,
     96                                                        md.results.TransientSolution[1].Pressure,
     97                                                        md.results.TransientSolution[1].Base,
     98                                                        md.results.TransientSolution[1].Surface,
     99                                                        md.results.TransientSolution[1].Thickness,
     100                                                        md.results.TransientSolution[1].SmbMonthlytemperatures,
     101                                                        md.results.TransientSolution[1].SmbMassBalance,
     102                                                        md.results.TransientSolution[2].Vx,
     103                                                        md.results.TransientSolution[2].Vy,
     104                                                        md.results.TransientSolution[2].Vel,
     105                                                        md.results.TransientSolution[2].Pressure,
     106                                                        md.results.TransientSolution[2].Base,
     107                                                        md.results.TransientSolution[2].Surface,
     108                                                        md.results.TransientSolution[2].Thickness,
     109                                                        md.results.TransientSolution[2].SmbMonthlytemperatures,
     110                                                        md.results.TransientSolution[2].SmbMassBalance]

  • issm/trunk-jpl/test/NightlyRun/test237.py

    r21408 r22864  
    3838
    3939# creating initialization and spc temperatures initialization and spc
    40 md.thermal.spctemperature=np.mean(md.smb.temperatures_lgm[0:md.mesh.numberofvertices,:],axis=1).reshape(-1,1)    #-10*ones(md.mesh.numberofvertices,1)
    41 md.thermal.spctemperature=np.tile(md.thermal.spctemperature,(1,md.timestepping.final_time/md.timestepping.time_step))
     40md.thermal.spctemperature=np.mean(md.smb.temperatures_lgm[0:md.mesh.numberofvertices,:],axis=1)    #-10*ones(md.mesh.numberofvertices,1)
     41md.thermal.spctemperature=np.tile(md.thermal.spctemperature,(int(md.timestepping.final_time/md.timestepping.time_step),1)).T
    4242itemp=np.arange(0,md.timestepping.final_time,md.timestepping.time_step)
    43 md.thermal.spctemperature=np.vstack((md.thermal.spctemperature,itemp.reshape(1,-1)))
     43md.thermal.spctemperature=np.vstack((md.thermal.spctemperature,itemp))
    4444
    45 md.initialization.temperature=md.smb.temperatures_lgm[0:md.mesh.numberofvertices,0].reshape(-1,1)   #*ones(md.mesh.numberofvertices,1)
     45md.initialization.temperature=md.smb.temperatures_lgm[0:md.mesh.numberofvertices,0]   #*ones(md.mesh.numberofvertices,1)
    4646md.smb.initialize(md)
    4747

  • issm/trunk-jpl/test/NightlyRun/test292.py

    r22341 r22864  
    2424# Fields and tolerances to track changes
    2525
    26 field_names     =['Vx1','Vy1','Vel1','Pressure1','Base1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Base2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Base3','Surface3','Thickness3']
    27 field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
    28 field_values=[\
    29         md.results.TransientSolution[0].Vx,\
    30         md.results.TransientSolution[0].Vy,\
    31         md.results.TransientSolution[0].Vel,\
    32         md.results.TransientSolution[0].Pressure,\
    33         md.results.TransientSolution[0].Base,\
    34         md.results.TransientSolution[0].Surface,\
    35         md.results.TransientSolution[0].Thickness,\
    36         md.results.TransientSolution[1].Vx,\
    37         md.results.TransientSolution[1].Vy,\
    38         md.results.TransientSolution[1].Vel,\
    39         md.results.TransientSolution[1].Pressure,\
    40         md.results.TransientSolution[1].Base,\
    41         md.results.TransientSolution[1].Surface,\
    42         md.results.TransientSolution[1].Thickness,\
    43         md.results.TransientSolution[2].Vx,\
    44         md.results.TransientSolution[2].Vy,\
    45         md.results.TransientSolution[2].Vel,\
    46         md.results.TransientSolution[2].Pressure,\
    47         md.results.TransientSolution[2].Base,\
    48         md.results.TransientSolution[2].Surface,\
    49         md.results.TransientSolution[2].Thickness,\
    50         ]
     26field_names     =['Vx1','Vy1','Vel1','Pressure1','Base1','Surface1','Thickness1',
     27                                                                        'Vx2','Vy2','Vel2','Pressure2','Base2','Surface2','Thickness2',
     28                                                                        'Vx3','Vy3','Vel3','Pressure3','Base3','Surface3','Thickness3']
     29field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,
     30                                                                        1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,
     31                                                                        1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
     32field_values=[md.results.TransientSolution[0].Vx,
     33                                                        md.results.TransientSolution[0].Vy,
     34                                                        md.results.TransientSolution[0].Vel,
     35                                                        md.results.TransientSolution[0].Pressure,
     36                                                        md.results.TransientSolution[0].Base,
     37                                                        md.results.TransientSolution[0].Surface,
     38                                                        md.results.TransientSolution[0].Thickness,
     39                                                        md.results.TransientSolution[1].Vx,
     40                                                        md.results.TransientSolution[1].Vy,
     41                                                        md.results.TransientSolution[1].Vel,
     42                                                        md.results.TransientSolution[1].Pressure,
     43                                                        md.results.TransientSolution[1].Base,
     44                                                        md.results.TransientSolution[1].Surface,
     45                                                        md.results.TransientSolution[1].Thickness,
     46                                                        md.results.TransientSolution[2].Vx,
     47                                                        md.results.TransientSolution[2].Vy,
     48                                                        md.results.TransientSolution[2].Vel,
     49                                                        md.results.TransientSolution[2].Pressure,
     50                                                        md.results.TransientSolution[2].Base,
     51                                                        md.results.TransientSolution[2].Surface,
     52                                                        md.results.TransientSolution[2].Thickness]

  • issm/trunk-jpl/test/NightlyRun/test501.py

    r21408 r22864  
    1919field_names     =['Vx','Vy','Vel','Pressure']
    2020field_tolerances=[1e-12,2e-12,2e-12,1e-13]
    21 field_values=[\
    22         md.results.StressbalanceSolution.Vx,\
    23         md.results.StressbalanceSolution.Vy,\
    24         md.results.StressbalanceSolution.Vel,\
    25         md.results.StressbalanceSolution.Pressure,\
    26         ]
     21field_values=[md.results.StressbalanceSolution.Vx,
     22                                                        md.results.StressbalanceSolution.Vy,
     23                                                        md.results.StressbalanceSolution.Vel,
     24                                                        md.results.StressbalanceSolution.Pressure]

  • issm/trunk-jpl/test/NightlyRun/test514.py

    r22216 r22864  
    4646
    4747#refine existing mesh 4
    48 md2=bamg(copy.deepcopy(md),'field',md.inversion.vy_obs,'hmin',1000.,'hmax',20000.,'gradation',3.,'Hessiantype',0,'err',np.array([[1.]]))
     48md2=bamg(copy.deepcopy(md),'field',md.inversion.vy_obs,'hmin',1000.,'hmax',20000.,'gradation',3.,'Hessiantype',0,'err',np.array([1.]))
    4949x6=md2.mesh.x
    5050y6=md2.mesh.y
    5151
    5252#refine existing mesh 5
    53 md2=bamg(copy.deepcopy(md),'field',np.hstack((md.inversion.vy_obs,md.geometry.thickness)),'hmin',1000.,'hmax',20000.,'gradation',3.,'Hessiantype',1,'err',np.array([[10.,100.]]))
     53md2=bamg(copy.deepcopy(md),'field',np.vstack((md.inversion.vy_obs,md.geometry.thickness)).T,'hmin',1000.,'hmax',20000.,'gradation',3.,'Hessiantype',1,'err',np.array([[10.,100.]]))
    5454x7=md2.mesh.x
    5555y7=md2.mesh.y

  • issm/trunk-jpl/test/NightlyRun/test515.py

    r21408 r22864  
    2424
    2525# Fields and tolerances to track changes
    26 field_names     =['Temperature1','BasalforcingsGroundediceMeltingRate1', \
    27                                       'Temperature2','BasalforcingsGroundediceMeltingRate2']
     26field_names     =['Temperature1','BasalforcingsGroundediceMeltingRate1',
     27                                                                        'Temperature2','BasalforcingsGroundediceMeltingRate2']
    2828field_tolerances=[1e-13,1e-8,1e-13,5e-8]
    29 field_values=[\
    30         md.results.TransientSolution[0].Temperature,\
    31         md.results.TransientSolution[0].BasalforcingsGroundediceMeltingRate,\
    32         md.results.TransientSolution[1].Temperature,\
    33         md.results.TransientSolution[1].BasalforcingsGroundediceMeltingRate,\
    34         ]
     29field_values=[md.results.TransientSolution[0].Temperature,
     30                                                        md.results.TransientSolution[0].BasalforcingsGroundediceMeltingRate,
     31                                                        md.results.TransientSolution[1].Temperature,
     32                                                        md.results.TransientSolution[1].BasalforcingsGroundediceMeltingRate]

  • issm/trunk-jpl/test/NightlyRun/test601.py

    r21408 r22864  
    1919field_names     =['Thickness']
    2020field_tolerances=[1e-13]
    21 field_values=[\
    22         md.results.MasstransportSolution.Thickness,\
    23         ]
     21field_values=[md.results.MasstransportSolution.Thickness]

  • issm/trunk-jpl/test/NightlyRun/test611.py

    r21408 r22864  
    3535field_names     =['Gradient','Misfits','BalancethicknessThickeningRate','Thickness']
    3636field_tolerances=[1e-12,1e-12,1e-12,1e-12,1e-12,1e-12]
    37 field_values=[\
    38         md.results.BalancethicknessSolution.Gradient1,\
    39         md.results.BalancethicknessSolution.J,\
    40         md.results.BalancethicknessSolution.BalancethicknessThickeningRate,\
    41         md.results.BalancethicknessSolution.Thickness
    42 ]
     37field_values=[md.results.BalancethicknessSolution.Gradient1,
     38                                                        md.results.BalancethicknessSolution.J,
     39                                                        md.results.BalancethicknessSolution.BalancethicknessThickeningRate,
     40                                                        md.results.BalancethicknessSolution.Thickness]

  • issm/trunk-jpl/test/NightlyRun/test701.py

    r22268 r22864  
    3636
    3737#Boundary conditions
    38 md.stressbalance.referential = float('NaN')*np.ones((md.mesh.numberofvertices,6))
     38md.stressbalance.referential = np.nan*np.ones((md.mesh.numberofvertices,6))
    3939md.stressbalance.loadingforce = 0. * np.ones((md.mesh.numberofvertices,3))
    40 md.stressbalance.spcvx = float('NaN') * np.ones((md.mesh.numberofvertices,))
    41 md.stressbalance.spcvy = float('NaN') * np.ones((md.mesh.numberofvertices,))
    42 md.stressbalance.spcvz = float('NaN') * np.ones((md.mesh.numberofvertices,))
     40md.stressbalance.spcvx = np.nan * np.ones((md.mesh.numberofvertices,))
     41md.stressbalance.spcvy = np.nan * np.ones((md.mesh.numberofvertices,))
     42md.stressbalance.spcvz = np.nan * np.ones((md.mesh.numberofvertices,))
    4343md.stressbalance.spcvx[np.where(md.mesh.vertexflags(4))] = 0.
    4444md.stressbalance.spcvy[np.where(md.mesh.vertexflags(4))] = 0.
    … …  
    4646#Misc
    4747md = setflowequation(md,'FS','all')
    48 md.stressbalance.abstol = float('NaN')
     48md.stressbalance.abstol = np.nan
    4949#md.stressbalance.reltol = 10**-16
    5050md.stressbalance.FSreconditioning = 1.
    … …  
    6767        field_names = field_names + [['Vx'+i],['Vy'+i],['Vel'+i],['Pressure'+i]]
    6868        field_tolerances = field_tolerances + [9e-5,9e-5,9e-5,1e-10]
    69         field_values = field_values + [
    70                 md.results.StressbalanceSolution.Vx,
    71                 md.results.StressbalanceSolution.Vy,
    72                 md.results.StressbalanceSolution.Vel,
    73                 md.results.StressbalanceSolution.Pressure
    74                 ]
     69        field_values = field_values + [md.results.StressbalanceSolution.Vx,
     70                                                                                                                                 md.results.StressbalanceSolution.Vy,
     71                                                                                                                                 md.results.StressbalanceSolution.Vel,
     72                                                                                                                                 md.results.StressbalanceSolution.Pressure]

  • issm/trunk-jpl/test/NightlyRun/test703.py

    r22857 r22864  
    3232#mesh domain
    3333md = bamgflowband(model(),x,b+h,b,'hmax',150.)
    34 print md.mesh.numberofvertices
    3534
    3635#geometry
    … …  
    4039
    4140#mask
    42 md.mask.ice_levelset = -np.ones((md.mesh.numberofvertices,))
     41md.mask.ice_levelset = -np.ones((md.mesh.numberofvertices))
    4342md.mask.ice_levelset[np.where(md.mesh.vertexflags(2))] = 0
    44 md.mask.groundedice_levelset = np.zeros((md.mesh.numberofvertices,)) - 0.5
     43md.mask.groundedice_levelset = -0.5*np.ones((md.mesh.numberofvertices))
     44md.mask.groundedice_levelset[np.where(md.mesh.x<0)]=0.5
    4545
    4646#materials
    47 md.initialization.temperature = (273. - 20.) * np.ones((md.mesh.numberofvertices,))
     47md.initialization.temperature = (273. - 20.) * np.ones((md.mesh.numberofvertices))
    4848md.materials.rheology_B = paterson(md.initialization.temperature)
    49 md.materials.rheology_n = 3 * np.ones((md.mesh.numberofelements,))
     49md.materials.rheology_n = 3 * np.ones((md.mesh.numberofelements))
    5050
    5151#friction
    52 md.friction.coefficient = np.zeros((md.mesh.numberofvertices,))
     52md.friction.coefficient = np.zeros((md.mesh.numberofvertices))
    5353md.friction.coefficient[np.where(md.mesh.vertexflags(1))] = 20
    54 md.friction.p = np.ones((md.mesh.numberofelements,))
    55 md.friction.q = np.ones((md.mesh.numberofelements,))
     54md.friction.p = np.ones((md.mesh.numberofelements))
     55md.friction.q = np.ones((md.mesh.numberofelements))
    5656
    5757#boundary conditions
    58 md.stressbalance.spcvx = np.nan * np.ones((md.mesh.numberofvertices,))
    59 md.stressbalance.spcvy = np.nan * np.ones((md.mesh.numberofvertices,))
    60 md.stressbalance.spcvz = np.nan * np.ones((md.mesh.numberofvertices,))
     58md.stressbalance.spcvx = np.nan * np.ones((md.mesh.numberofvertices))
     59md.stressbalance.spcvy = np.nan * np.ones((md.mesh.numberofvertices))
     60md.stressbalance.spcvz = np.nan * np.ones((md.mesh.numberofvertices))
    6161md.stressbalance.referential = np.nan * np.ones((md.mesh.numberofvertices,6))
    6262md.stressbalance.loadingforce = 0 * np.ones((md.mesh.numberofvertices,3))
    6363md.stressbalance.spcvx[np.where(md.mesh.vertexflags(4))] = 800.
    6464md.stressbalance.spcvy[np.where(md.mesh.vertexflags(4))] = 0.
    65 
    66 #print md.stressbalance.referential
    67 #print md.stressbalance.loadingforce
    68 #print md.stressbalance.spcvx
    69 #print md.stressbalance.spcvy
    70 #print md.stressbalance.spcvz
    71 #print np.shape(md.stressbalance.referential)
    72 #print np.shape(md.stressbalance.loadingforce)
    73 #print np.shape(md.stressbalance.spcvx)
    74 #print np.shape(md.stressbalance.spcvy)
    75 #print np.shape(md.stressbalance.spcvz)
    7665
    7766#Misc
    … …  
    10392md.cluster = generic('np',3)
    10493md.stressbalance.shelf_dampening = 1
    105 md1 = solve(md,'Transient')
     94md1=copy.deepcopy(md)
     95md1 = solve(md1,'Transient')
    10696
    10797md.stressbalance.shelf_dampening = 0
    … …  
    10999
    110100#Fields and tolerances to track changes
    111 field_names      = [
    112         'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1',
    113         'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2',
    114         'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3',
    115         'Vx1_damp','Vy1_damp','Vel1_damp','Pressure1_damp','Bed1_damp','Surface1_damp','Thickness1_damp',
    116         'Vx2_damp','Vy2_damp','Vel2_damp','Pressure2_damp','Bed2_damp','Surface2_damp','Thickness2_damp',
    117         'Vx3_damp','Vy3_damp','Vel3_damp','Pressure3_damp','Bed3_damp','Surface3_damp','Thickness3_damp']
    118 field_tolerances = [
    119         2e-08,2e-08,2e-08,1e-08,1e-10,1e-10,1e-10,
    120         2e-08,2e-08,2e-08,1e-08,1e-10,1e-10,1e-10,
    121         2e-08,2e-08,2e-08,1e-08,1e-10,1e-10,1e-10,
    122         5e-08,5e-08,5e-08,1e-08,1e-10,1e-10,1e-10,
    123         5e-08,5e-08,5e-08,1e-08,1e-10,1e-10,1e-10,
    124         5e-08,5e-08,5e-08,1e-08,1e-10,1e-10,1e-10]
    125 field_values = [
    126         md.results.TransientSolution[0].Vx,
    127         md.results.TransientSolution[0].Vy,
    128         md.results.TransientSolution[0].Vel,
    129         md.results.TransientSolution[0].Pressure,
    130         md.results.TransientSolution[0].Base,
    131         md.results.TransientSolution[0].Surface,
    132         md.results.TransientSolution[0].Thickness,
    133         md.results.TransientSolution[1].Vx,
    134         md.results.TransientSolution[1].Vy,
    135         md.results.TransientSolution[1].Vel,
    136         md.results.TransientSolution[1].Pressure,
    137         md.results.TransientSolution[1].Base,
    138         md.results.TransientSolution[1].Surface,
    139         md.results.TransientSolution[1].Thickness,
    140         md.results.TransientSolution[2].Vx,
    141         md.results.TransientSolution[2].Vy,
    142         md.results.TransientSolution[2].Vel,
    143         md.results.TransientSolution[2].Pressure,
    144         md.results.TransientSolution[2].Base,
    145         md.results.TransientSolution[2].Surface,
    146         md.results.TransientSolution[2].Thickness,
    147         md1.results.TransientSolution[0].Vx,
    148         md1.results.TransientSolution[0].Vy,
    149         md1.results.TransientSolution[0].Vel,
    150         md1.results.TransientSolution[0].Pressure,
    151         md1.results.TransientSolution[0].Base,
    152         md1.results.TransientSolution[0].Surface,
    153         md1.results.TransientSolution[0].Thickness,
    154         md1.results.TransientSolution[1].Vx,
    155         md1.results.TransientSolution[1].Vy,
    156         md1.results.TransientSolution[1].Vel,
    157         md1.results.TransientSolution[1].Pressure,
    158         md1.results.TransientSolution[1].Base,
    159         md1.results.TransientSolution[1].Surface,
    160         md1.results.TransientSolution[1].Thickness,
    161         md1.results.TransientSolution[2].Vx,
    162         md1.results.TransientSolution[2].Vy,
    163         md1.results.TransientSolution[2].Vel,
    164         md1.results.TransientSolution[2].Pressure,
    165         md1.results.TransientSolution[2].Base,
    166         md1.results.TransientSolution[2].Surface,
    167         md1.results.TransientSolution[2].Thickness,
    168         ]
     101field_names      = ['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1',
     102                                                                                'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2',
     103                                                                                'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3',
     104                                                                                'Vx1_damp','Vy1_damp','Vel1_damp','Pressure1_damp','Bed1_damp','Surface1_damp','Thickness1_damp',
     105                                                                                'Vx2_damp','Vy2_damp','Vel2_damp','Pressure2_damp','Bed2_damp','Surface2_damp','Thickness2_damp',
     106                                                                                'Vx3_damp','Vy3_damp','Vel3_damp','Pressure3_damp','Bed3_damp','Surface3_damp','Thickness3_damp']
     107field_tolerances = [2e-08,2e-08,2e-08,1e-08,1e-10,1e-10,1e-10,
     108                                                                                2e-08,2e-08,2e-08,1e-08,1e-10,1e-10,1e-10,
     109                                                                                2e-08,2e-08,2e-08,1e-08,1e-10,1e-10,1e-10,
     110                                                                                5e-08,5e-08,5e-08,1e-08,1e-10,1e-10,1e-10,
     111                                                                                5e-08,5e-08,5e-08,1e-08,1e-10,1e-10,1e-10,
     112                                                                                5e-08,5e-08,5e-08,1e-08,1e-10,1e-10,1e-10]
     113field_values = [md.results.TransientSolution[0].Vx,
     114                                                                md.results.TransientSolution[0].Vy,
     115                                                                md.results.TransientSolution[0].Vel,
     116                                                                md.results.TransientSolution[0].Pressure,
     117                                                                md.results.TransientSolution[0].Base,
     118                                                                md.results.TransientSolution[0].Surface,
     119                                                                md.results.TransientSolution[0].Thickness,
     120                                                                md.results.TransientSolution[1].Vx,
     121                                                                md.results.TransientSolution[1].Vy,
     122                                                                md.results.TransientSolution[1].Vel,
     123                                                                md.results.TransientSolution[1].Pressure,
     124                                                                md.results.TransientSolution[1].Base,
     125                                                                md.results.TransientSolution[1].Surface,
     126                                                                md.results.TransientSolution[1].Thickness,
     127                                                                md.results.TransientSolution[2].Vx,
     128                                                                md.results.TransientSolution[2].Vy,
     129                                                                md.results.TransientSolution[2].Vel,
     130                                                                md.results.TransientSolution[2].Pressure,
     131                                                                md.results.TransientSolution[2].Base,
     132                                                                md.results.TransientSolution[2].Surface,
     133                                                                md.results.TransientSolution[2].Thickness,
     134                                                                md1.results.TransientSolution[0].Vx,
     135                                                                md1.results.TransientSolution[0].Vy,
     136                                                                md1.results.TransientSolution[0].Vel,
     137                                                                md1.results.TransientSolution[0].Pressure,
     138                                                                md1.results.TransientSolution[0].Base,
     139                                                                md1.results.TransientSolution[0].Surface,
     140                                                                md1.results.TransientSolution[0].Thickness,
     141                                                                md1.results.TransientSolution[1].Vx,
     142                                                                md1.results.TransientSolution[1].Vy,
     143                                                                md1.results.TransientSolution[1].Vel,
     144                                                                md1.results.TransientSolution[1].Pressure,
     145                                                                md1.results.TransientSolution[1].Base,
     146                                                                md1.results.TransientSolution[1].Surface,
     147                                                                md1.results.TransientSolution[1].Thickness,
     148                                                                md1.results.TransientSolution[2].Vx,
     149                                                                md1.results.TransientSolution[2].Vy,
     150                                                                md1.results.TransientSolution[2].Vel,
     151                                                                md1.results.TransientSolution[2].Pressure,
     152                                                                md1.results.TransientSolution[2].Base,
     153                                                                md1.results.TransientSolution[2].Surface,
     154                                                                md1.results.TransientSolution[2].Thickness]

  • issm/trunk-jpl/test/NightlyRun/test808.py

    r22267 r22864  
    4444
    4545#Fields and tolerances to track changes
    46 field_names = [
    47         'Vx1','Vy1','Vel1','Pressure1','Thickness1','Surface1','MaskIceLevelset1'
    48         'Vx2','Vy2','Vel2','Pressure2','Thickness2','Surface2','MaskIceLevelset2'
    49         'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3']
    50 field_tolerances = [
    51         1e-8,1e-8,1e-8,1e-9,1e-9,1e-9,3e-9,
    52         1e-8,1e-8,1e-8,1e-9,1e-9,1e-9,3e-9,
    53         1e-8,1e-8,1e-8,1e-9,1e-9,1e-9,3e-9]
    54 field_values = [
    55         md.results.TransientSolution[0].Vx,
    56         md.results.TransientSolution[0].Vy,
    57         md.results.TransientSolution[0].Vel,
    58         md.results.TransientSolution[0].Pressure,
    59         md.results.TransientSolution[0].Thickness,
    60         md.results.TransientSolution[0].Surface,
    61         md.results.TransientSolution[0].MaskIceLevelset,
    62         md.results.TransientSolution[1].Vx,
    63         md.results.TransientSolution[1].Vy,
    64         md.results.TransientSolution[1].Vel,
    65         md.results.TransientSolution[1].Pressure,
    66         md.results.TransientSolution[1].Thickness,
    67         md.results.TransientSolution[1].Surface,
    68         md.results.TransientSolution[1].MaskIceLevelset,
    69         md.results.TransientSolution[2].Vx,
    70         md.results.TransientSolution[2].Vy,
    71         md.results.TransientSolution[2].Vel,
    72         md.results.TransientSolution[2].Pressure,
    73         md.results.TransientSolution[2].Thickness,
    74         md.results.TransientSolution[2].Surface,
    75         md.results.TransientSolution[2].MaskIceLevelset
    76         ]
     46field_names = [ 'Vx1','Vy1','Vel1','Pressure1','Thickness1','Surface1','MaskIceLevelset1'
     47                                                                'Vx2','Vy2','Vel2','Pressure2','Thickness2','Surface2','MaskIceLevelset2'
     48                                                                'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3']
     49field_tolerances = [1e-8,1e-8,1e-8,1e-9,1e-9,1e-9,3e-9,
     50                                                                                1e-8,1e-8,1e-8,1e-9,1e-9,1e-9,3e-9,
     51                                                                                1e-8,1e-8,1e-8,1e-9,1e-9,1e-9,3e-9]
     52field_values = [md.results.TransientSolution[0].Vx,
     53                                                                md.results.TransientSolution[0].Vy,
     54                                                                md.results.TransientSolution[0].Vel,
     55                                                                md.results.TransientSolution[0].Pressure,
     56                                                                md.results.TransientSolution[0].Thickness,
     57                                                                md.results.TransientSolution[0].Surface,
     58                                                                md.results.TransientSolution[0].MaskIceLevelset,
     59                                                                md.results.TransientSolution[1].Vx,
     60                                                                md.results.TransientSolution[1].Vy,
     61                                                                md.results.TransientSolution[1].Vel,
     62                                                                md.results.TransientSolution[1].Pressure,
     63                                                                md.results.TransientSolution[1].Thickness,
     64                                                                md.results.TransientSolution[1].Surface,
     65                                                                md.results.TransientSolution[1].MaskIceLevelset,
     66                                                                md.results.TransientSolution[2].Vx,
     67                                                                md.results.TransientSolution[2].Vy,
     68                                                                md.results.TransientSolution[2].Vel,
     69                                                                md.results.TransientSolution[2].Pressure,
     70                                                                md.results.TransientSolution[2].Thickness,
     71                                                                md.results.TransientSolution[2].Surface,
     72                                                                md.results.TransientSolution[2].MaskIceLevelset]
Note: See TracChangeset for help on using the changeset viewer.