Changeset 18968 for issm/trunk-jpl/test/NightlyRun/test236.py

Timestamp:

12/30/14 08:26:56 (10 years ago)

Author:

lemorzad

Message:

updating pdd module and NR 236 and 237. Updated python scripts may not work as I am not familiar with python.

File:

• issm/trunk-jpl/test/NightlyRun/test236.py (modified) (3 diffs)

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

@@ -16 +16 @@
 md.surfaceforcings = SMBpdd();
 md.surfaceforcings.isdelta18o=1
+md.surfaceforcings.ismungsm=0
 
 # Add temperature, precipitation and delta18o needed to measure the surface mass balance
@@ -50 +51 @@
     md.surfaceforcings.precipitations_presentday[0:md.mesh.numberofvertices,imonth]=-0.4*10**(-6)*md.mesh.y+0.5
     md.surfaceforcings.precipitations_presentday[md.mesh.numberofvertices,imonth]=((float(imonth)+1.)/12.)
+    md.surfaceforcings.precipitations_lgm[0:md.mesh.numberofvertices,imonth]=-0.4*10**(-6)*md.mesh.y+0.5
+    md.surfaceforcings.precipitations_lgm[md.mesh.numberofvertices,imonth]=((float(imonth)+1.)/12.)
+
+# Interpolation factors
+md.surfaceforcings.Tdiff[1,1:md.timestepping.final_time]=0.5;
+md.surfaceforcings.sealev[1,1:md.timestepping.final_time]=0.5;
+# Year of each data point
+md.surfaceforcings.Tdiff[2,1:md.timestepping.final_time]=1:1:md.timestepping.final_time;
+md.surfaceforcings.sealev[2,1:md.timestepping.final_time]=1:1:md.timestepping.final_time;
 
 # time steps and resolution
@@ -55 +65 @@
 md.timestepping.final_time=60.
 
-md.surfaceforcings.Pfac=[1,1]
-md.surfaceforcings.Tdiff=[1,1]
-md.surfaceforcings.sealev=[1,1]
 
 #