Index: /issm/trunk-jpl/test/NightlyRun/test237.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test237.py	(revision 25374)
+++ /issm/trunk-jpl/test/NightlyRun/test237.py	(revision 25375)
@@ -27,5 +27,6 @@
 # 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,4 +64,8 @@
     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.)
Index: /issm/trunk-jpl/test/NightlyRun/test239.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test239.py	(revision 25374)
+++ /issm/trunk-jpl/test/NightlyRun/test239.py	(revision 25375)
@@ -25,5 +25,6 @@
 # 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,4 +46,8 @@
 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')