Changeset 26906
- Timestamp:
- 02/25/22 00:30:49 (3 years ago)
- Location:
- issm/trunk-jpl/test/NightlyRun
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test247.py
r26402 r26906 1 #Test Name: SquareShelfTranIspddIsdeltaO18pdNoInterpSSA2d 1 #Test Name: SquareShelfTranIspddIsdeltaO18pdNoInterpSSA2d 2 2 import numpy as np 3 3 from model import * … … 57 57 tim1 = np.linspace(1,12,12)/12 58 58 59 md.smb.temperatures_reconstructed[md.mesh.numberofvertices,0:12] = tim1 59 md.smb.temperatures_reconstructed[md.mesh.numberofvertices,0:12] = tim1 60 60 md.smb.temperatures_reconstructed[md.mesh.numberofvertices,12:24] = tim1+3 61 61 md.smb.temperatures_reconstructed[md.mesh.numberofvertices,24:36] = tim1+5 62 62 63 63 md.smb.precipitations_reconstructed[md.mesh.numberofvertices,0:12] = tim1 64 md.smb.precipitations_reconstructed[md.mesh.numberofvertices,12:24] = tim1+3 64 md.smb.precipitations_reconstructed[md.mesh.numberofvertices,12:24] = tim1+3 65 65 md.smb.precipitations_reconstructed[md.mesh.numberofvertices,24:36] = tim1+5 66 66 … … 68 68 # spc 69 69 md.thermal.spctemperature=np.mean(md.smb.temperatures_presentday[0:md.mesh.numberofvertices,0:12],axis=1)-10. 70 md.initialization.temperature=md.thermal.spctemperature 70 md.initialization.temperature=md.thermal.spctemperature 71 71 72 72 md.smb.s0p = np.maximum(md.geometry.surface.reshape(-1, 1),np.zeros((md.mesh.numberofvertices,1))) … … 89 89 md.thermal.isdrainicecolumn = 0 90 90 91 md =solve(md,'thermal')91 md = solve(md, 'thermal') 92 92 93 93 md.initialization.temperature = md.results.ThermalSolution.Temperature … … 121 121 122 122 #Fields and tolerances to track changes 123 field_names = ['Vx1', 'Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','Enthalpy1','SmbMassBalance1',124 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','Enthalpy2','SmbMassBalance2',125 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','Enthalpy3','SmbMassBalance3',126 'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','Temperature4','Enthalpy4','SmbMassBalance4']127 field_tolerances=[1e-09, 1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,128 1e-09, 1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,129 1e-09, 1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,130 1e-09, 1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13]123 field_names = ['Vx1', 'Vy1', 'Vz1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'Temperature1', 'Enthalpy1', 'SmbMassBalance1', 124 'Vx2', 'Vy2', 'Vz2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'Temperature2', 'Enthalpy2', 'SmbMassBalance2', 125 'Vx3', 'Vy3', 'Vz3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'Temperature3', 'Enthalpy3', 'SmbMassBalance3', 126 'Vx4', 'Vy4', 'Vz4', 'Vel4', 'Pressure4', 'Bed4', 'Surface4', 'Thickness4', 'Temperature4', 'Enthalpy4', 'SmbMassBalance4'] 127 field_tolerances=[1e-09, 1e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-13, 128 1e-09, 1e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-13, 129 1e-09, 1e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-13, 130 1e-09, 1e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-13] 131 131 field_values = [md.results.TransientSolution[0].Vx, 132 md.results.TransientSolution[0].Vy,133 md.results.TransientSolution[0].Vz,134 md.results.TransientSolution[0].Vel,135 md.results.TransientSolution[0].Pressure,136 md.results.TransientSolution[0].Base,137 md.results.TransientSolution[0].Surface,138 md.results.TransientSolution[0].Thickness,139 md.results.TransientSolution[0].Temperature,140 md.results.TransientSolution[0].Enthalpy,141 md.results.TransientSolution[0].SmbMassBalance,142 md.results.TransientSolution[1].Vx,143 md.results.TransientSolution[1].Vy,144 md.results.TransientSolution[1].Vz,145 md.results.TransientSolution[1].Vel,146 md.results.TransientSolution[1].Pressure,147 md.results.TransientSolution[1].Base,148 md.results.TransientSolution[1].Surface,149 md.results.TransientSolution[1].Thickness,150 md.results.TransientSolution[1].Temperature,151 md.results.TransientSolution[1].Enthalpy,152 md.results.TransientSolution[1].SmbMassBalance,153 md.results.TransientSolution[2].Vx,154 md.results.TransientSolution[2].Vy,155 md.results.TransientSolution[2].Vz,156 md.results.TransientSolution[2].Vel,157 md.results.TransientSolution[2].Pressure,158 md.results.TransientSolution[2].Base,159 md.results.TransientSolution[2].Surface,160 md.results.TransientSolution[2].Thickness,161 md.results.TransientSolution[2].Temperature,162 md.results.TransientSolution[2].Enthalpy,163 md.results.TransientSolution[2].SmbMassBalance,164 md.results.TransientSolution[3].Vx,165 md.results.TransientSolution[3].Vy,166 md.results.TransientSolution[3].Vz,167 md.results.TransientSolution[3].Vel,168 md.results.TransientSolution[3].Pressure,169 md.results.TransientSolution[3].Base,170 md.results.TransientSolution[3].Surface,171 md.results.TransientSolution[3].Thickness,172 md.results.TransientSolution[3].Temperature,173 md.results.TransientSolution[3].Enthalpy,174 md.results.TransientSolution[3].SmbMassBalance]132 md.results.TransientSolution[0].Vy, 133 md.results.TransientSolution[0].Vz, 134 md.results.TransientSolution[0].Vel, 135 md.results.TransientSolution[0].Pressure, 136 md.results.TransientSolution[0].Base, 137 md.results.TransientSolution[0].Surface, 138 md.results.TransientSolution[0].Thickness, 139 md.results.TransientSolution[0].Temperature, 140 md.results.TransientSolution[0].Enthalpy, 141 md.results.TransientSolution[0].SmbMassBalance, 142 md.results.TransientSolution[1].Vx, 143 md.results.TransientSolution[1].Vy, 144 md.results.TransientSolution[1].Vz, 145 md.results.TransientSolution[1].Vel, 146 md.results.TransientSolution[1].Pressure, 147 md.results.TransientSolution[1].Base, 148 md.results.TransientSolution[1].Surface, 149 md.results.TransientSolution[1].Thickness, 150 md.results.TransientSolution[1].Temperature, 151 md.results.TransientSolution[1].Enthalpy, 152 md.results.TransientSolution[1].SmbMassBalance, 153 md.results.TransientSolution[2].Vx, 154 md.results.TransientSolution[2].Vy, 155 md.results.TransientSolution[2].Vz, 156 md.results.TransientSolution[2].Vel, 157 md.results.TransientSolution[2].Pressure, 158 md.results.TransientSolution[2].Base, 159 md.results.TransientSolution[2].Surface, 160 md.results.TransientSolution[2].Thickness, 161 md.results.TransientSolution[2].Temperature, 162 md.results.TransientSolution[2].Enthalpy, 163 md.results.TransientSolution[2].SmbMassBalance, 164 md.results.TransientSolution[3].Vx, 165 md.results.TransientSolution[3].Vy, 166 md.results.TransientSolution[3].Vz, 167 md.results.TransientSolution[3].Vel, 168 md.results.TransientSolution[3].Pressure, 169 md.results.TransientSolution[3].Base, 170 md.results.TransientSolution[3].Surface, 171 md.results.TransientSolution[3].Thickness, 172 md.results.TransientSolution[3].Temperature, 173 md.results.TransientSolution[3].Enthalpy, 174 md.results.TransientSolution[3].SmbMassBalance] -
issm/trunk-jpl/test/NightlyRun/test317.py
r23793 r26906 20 20 21 21 #Fields and tolerances to track changes 22 field_names = ['Vx1', 'Vy1', 'Vz1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'Temperature1', 'BasalforcingsGroundediceMeltingRate1', 'GroundedArea1', 'FloatingArea1', 'TotalFloatingBmb1', 'TotalGroundedBmb1' 22 field_names = ['Vx1', 'Vy1', 'Vz1', 'Vel1', 'Pressure1', 'Bed1', 'Surface1', 'Thickness1', 'Temperature1', 'BasalforcingsGroundediceMeltingRate1', 'GroundedArea1', 'FloatingArea1', 'TotalFloatingBmb1', 'TotalGroundedBmb1', 23 23 'Vx2', 'Vy2', 'Vz2', 'Vel2', 'Pressure2', 'Bed2', 'Surface2', 'Thickness2', 'Temperature2', 'BasalforcingsGroundediceMeltingRate2', 'GroundedArea2', 'FloatingArea2', 'TotalFloatingBmb2', 'TotalGroundedBmb2', 24 24 'Vx3', 'Vy3', 'Vz3', 'Vel3', 'Pressure3', 'Bed3', 'Surface3', 'Thickness3', 'Temperature3', 'BasalforcingsGroundediceMeltingRate3', 'GroundedArea3', 'FloatingArea3', 'TotalFloatingBmb2', 'TotalGroundedBmb2'] -
issm/trunk-jpl/test/NightlyRun/test427.py
r24862 r26906 2 2 from model import * 3 3 from socket import gethostname 4 5 from setmask import *6 from parameterize import *7 from setflowequation import *8 from solve import *4 from triangle import triangle 5 from setmask import setmask 6 from parameterize import parameterize 7 from setflowequation import setflowequation 8 from solve import solve 9 9 10 10 md = triangle(model(), '../Exp/Square.exp', 150000.) -
issm/trunk-jpl/test/NightlyRun/test436.py
r24214 r26906 3 3 from model import * 4 4 from socket import gethostname 5 from triangle import *6 from setmask import *7 from parameterize import *8 from setflowequation import *9 from solve import *5 from triangle import triangle 6 from setmask import setmask 7 from parameterize import parameterize 8 from setflowequation import setflowequation 9 from solve import solve 10 10 11 11 md = triangle(model(), '../Exp/Square.exp', 150000.) … … 26 26 for i in ['LliboutryDuval', 'CuffeyTemperate']: 27 27 print(' ') 28 print(' == ==== Testing rheology law: ' + i + ' = ')28 print(' ====== Testing rheology law: ' + i + ' = ') 29 29 30 30 md.materials.rheology_law = i -
issm/trunk-jpl/test/NightlyRun/test443.py
r24862 r26906 3 3 from model import * 4 4 from socket import gethostname 5 6 from triangle import * 7 from setmask import * 8 from parameterize import * 9 from setflowequation import * 10 from solve import * 5 from triangle import triangle 6 from setmask import setmask 7 from parameterize import parameterize 8 from setflowequation import setflowequation 9 from solve import solve 11 10 from frictioncoulomb import frictioncoulomb 12 11 from ContourToMesh import ContourToMesh -
issm/trunk-jpl/test/NightlyRun/test612.py
r25023 r26906 2 2 import numpy as np 3 3 from model import * 4 from setmask import * 4 from setmask import setmask 5 from setflowequation import setflowequation 5 6 from triangle import triangle 6 7 from parameterize import parameterize -
issm/trunk-jpl/test/NightlyRun/test701.py
r24862 r26906 66 66 md.flowequation.fe_FS = i 67 67 md = solve(md, 'Stressbalance') 68 field_names = field_names + [ ['Vx' + i], ['Vy' + i], ['Vel' + i], ['Pressure' + i]]68 field_names = field_names + ['Vx' + i, 'Vy' + i, 'Vel' + i, 'Pressure' + i] 69 69 field_tolerances = field_tolerances + [9e-5, 9e-5, 9e-5, 1e-10] 70 70 field_values = field_values + [md.results.StressbalanceSolution.Vx,
Note:
See TracChangeset
for help on using the changeset viewer.