Changeset 19270

Timestamp:

04/08/15 07:45:44 (10 years ago)

Author:

lemorzad

Message:

updating PDD nightly run 237 and 240

Location:

issm/trunk-jpl/test/NightlyRun

Files:

• test237.m (modified) (2 diffs)
• test240.m (modified) (2 diffs)

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

@@ -11 +11 @@
 md.surfaceforcings.ismungsm=1;
 
-if md.surfaceforcings.isdelta18o==0 & md.surfaceforcings.ismungsm==0
-    md.surfaceforcings.precipitation=zeros(md.mesh.numberofvertices,1);
-    md.surfaceforcings.monthlytemperatures=273*ones(md.mesh.numberofvertices,1);
-end
-    
-
-% Add temperature, precipitation and delta18o needed to measure the surface mass balance
-% % creating delta18o
-% load '../Data/delta18o.data'
-% md.surfaceforcings.delta18o=delta18o;
-% % creating delta18oSurface
-% md.surfaceforcings.delta18o_surface(1,1:(length(delta18o))) = 0;
-% md.surfaceforcings.delta18o_surface(2,:) = delta18o(2,:);
+% time steps and resolution
+md.timestepping.time_step=20;
+md.settings.output_frequency=1;
+md.timestepping.final_time=60;
 
 % creating Present day and lgm temperatures
@@ -53 +44 @@
 end
 
-md.surfaceforcings.Pfac(1,1:md.timestepping.final_time)=0.5;
-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.Pfac(2,1:md.timestepping.final_time)=1:1:md.timestepping.final_time;
-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
-md.timestepping.time_step=20;
-md.settings.output_frequency=1;
-md.timestepping.final_time=60;
+for iint=1:(md.timestepping.final_time/md.timestepping.time_step)+1
+    md.surfaceforcings.Pfac(1,iint)=0.15*iint;
+    md.surfaceforcings.Tdiff(1,iint)=0.15*iint;
+    md.surfaceforcings.sealev(1,iint)=0.15*iint;
+    % Year of each data point
+    md.surfaceforcings.Pfac(2,iint)=(iint-1)*20;
+    md.surfaceforcings.Tdiff(2,iint)=(iint-1)*20;
+    md.surfaceforcings.sealev(2,iint)=(iint-1)*20;
+end
 
 %

issm/trunk-jpl/test/NightlyRun/test240.m

@@ -12 +12 @@
 %md.surfaceforcings.monthlytemperatures(1:md.mesh.numberofvertices,1:12)=273;
 
+% time steps and resolution
+md.timestepping.time_step=20;
+md.settings.output_frequency=1;
+md.timestepping.final_time=60;
+md.timestepping.interp_forcings=0;
+
 % Add temperature, precipitation and delta18o needed to measure the surface mass balance
 %  creating delta18o
 load '../Data/delta18o.data'
-md.surfaceforcings.delta18o=delta18o;
+md.surfaceforcings.delta18o=delta18o(:,1:(md.timestepping.final_time/md.timestepping.time_step)+1);
 
 % creating Present day  temperatures
@@ -39 +45 @@
 end
 
-% time steps and resolution
-md.timestepping.time_step=0.5;
-md.settings.output_frequency=1;
-md.timestepping.final_time=2;
-md.timestepping.interp_forcings=0;
 
 %