Changeset 25375

Timestamp:

08/11/20 17:43:14 (5 years ago)

Author:

schlegel

Message:

CHG: update python tests for pdd melt

Location:

issm/trunk-jpl/test/NightlyRun

Files:

• test237.py (modified) (2 diffs)
• test239.py (modified) (2 diffs)

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

@@ -27 +27 @@
 # creating Present day and lgm temperatures
 # Same temperature over the all region:
-tmonth = np.ones(12) * (238.15 + 20.)
+curve=np.sin(np.linspace(0, np.pi, 12))
+tmonth = np.ones(12) * (238.15 + 20.) + 20.*curve
 md.smb.temperatures_presentday = np.zeros((md.mesh.numberofvertices, 12))
 md.smb.temperatures_lgm = np.zeros((md.mesh.numberofvertices, 12))
@@ -63 +64 @@
     md.smb.Tdiff[1, iint] = (float(iint)) * 20
     md.smb.sealev[1, iint] = (float(iint)) * 20
+
+md.smb.issetpddfac=1.
+md.smb.pddfac_snow=2.
+md.smb.pddfac_ice=2.
 
 md.extrude(3, 1.)

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

@@ -25 +25 @@
 # creating Present day temperatures
 # Same temperature over the all region:
-tmonth = np.ones(12) * (238.15 + 20.)
+curve=np.sin(np.linspace(0, np.pi, 12))
+tmonth = np.ones(12) * (238.15 + 20.) + 20.*curve
 md.smb.temperatures_presentday = np.zeros((md.mesh.numberofvertices, 12))
 for imonth in range(0, 12):
@@ -45 +46 @@
 md.settings.output_frequency = 1
 md.timestepping.final_time = 2
+
+md.smb.issetpddfac = 1.
+md.smb.pddfac_snow = 3.
+md.smb.pddfac_ice = 3.
 
 md = setflowequation(md, 'SSA', 'all')