Changeset 24760 for issm/trunk-jpl/test/NightlyRun/test243.py

Timestamp:

04/29/20 14:07:15 (5 years ago)

Author:

schlegel

Message:

CHG: Gemb maintain initial layer depth profile throughout the simulation

File:

• issm/trunk-jpl/test/NightlyRun/test243.py (modified) (2 diffs)

issm/trunk-jpl/test/NightlyRun/test243.py

@@ -21 +21 @@
 md.smb = SMBgemb()
 md.smb.setdefaultparameters(md.mesh, md.geometry)
-md.smb.dsnowIdx = 0
+md.smb.dsnowIdx = 1
 
 #load hourly surface forcing date from 1979 to 2009:
@@ -60 +60 @@
 md = solve(md, 'Transient')
 
-nlayers=np.size(md.results.TransientSolution[-1].SmbT,1)
+nlayers=np.size(md.results.TransientSolution[0].SmbT,1)
+for i in range(1, np.size(md.results.TransientSolution,0)):
+    nlayers=np.minimum(np.size(md.results.TransientSolution[i].SmbT,1),nlayers)
 
 #Fields and tolerances to track changes
 field_names = ['Layers','SmbDz', 'SmbT', 'SmbD', 'SmbRe', 'SmbGdn', 'SmbGsp', 'SmbA', 'SmbEC', 'SmbMassBalance', 'SmbMAdd', 'SmbDzAdd', 'SmbFAC', 'SmbMeanSHF', 'SmbMeanLHF', 'SmbMeanULW', 'SmbNetLW', 'SmbNetSW']
-field_tolerances = [1e-12,1e-12,1e-12,1e-11,1e-11,1e-11,1e-11,1e-12,1e-11,1e-12,1e-12,1e-12,1e-11,2e-11,2e-11,1e-11,9e-10,2e-11] 
+field_tolerances = [1e-12,2e-12,1e-12,1e-11,1e-11,2e-11,1e-11,1e-12,1e-11,1e-12,1e-12,1e-12,1e-11,2e-11,2e-11,1e-11,9e-10,2e-11] 
 #shape is different in python solution (fixed using reshape) which can cause test failure:
 field_values = [nlayers,md.results.TransientSolution[-1].SmbDz[0, 0:nlayers].reshape(1, -1),