Changeset 24624

Timestamp:

03/02/20 14:55:12 (5 years ago)

Author:

Mathieu Morlighem

Message:

CHG: use reconstructed field now

Location:

issm/trunk-jpl/test

Files:

• Archives/Archive247.arch (modified) ( previous)
• NightlyRun/test247.m (modified) (5 diffs)

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

@@ -1 @@
-%Test Name: SquareShelfTranIspddIsdeltaO18pdNoInterpSSA2d(
+%Test Name: SquareShelfTranIspddIsdeltaO18pdNoInterpSSA2d
 md=triangle(model(),'../Exp/Square.exp',60000.);
 md=setmask(md,'all','');
@@ -14 +14 @@
 load '../Data/delta18o.data'
 md.smb.delta18o=delta18o;
+md.smb.istemperaturescaled = 0;
+md.smb.isprecipscaled = 0;
 
 % creating Present day  temperatures
@@ -23 +25 @@
 md.smb = initialize(md.smb,md);
 
+% creating precipitation
+for imonth=0:11
+        md.smb.precipitations_presentday(1:md.mesh.numberofvertices,imonth+1)=-0.4*10^(-6)*md.mesh.y+0.5;
+end
+
+% 3 total years of input
+md.smb.temperatures_reconstructed = nan(md.mesh.numberofvertices+1,12*3);
+md.smb.precipitations_reconstructed = nan(md.mesh.numberofvertices+1,12*3);
+
+md.smb.temperatures_reconstructed(1:end-1,1:12) = md.smb.temperatures_presentday(1:end,:);
+md.smb.temperatures_reconstructed(1:end-1,13:24) = md.smb.temperatures_presentday(1:end,:)+1.2;
+md.smb.temperatures_reconstructed(1:end-1,25:36) = md.smb.temperatures_presentday(1:end,:)-0.8;
+
+md.smb.precipitations_reconstructed(1:end-1,1:12) = md.smb.precipitations_presentday(1:end,:);
+md.smb.precipitations_reconstructed(1:end-1,13:24) = md.smb.precipitations_presentday(1:end,:)+0.1;
+md.smb.precipitations_reconstructed(1:end-1,25:36) = md.smb.precipitations_presentday(1:end,:)-0.1;
+
+tim1 = [1/12:1/12:1];
+
+md.smb.temperatures_reconstructed(end,1:12) = tim1; md.smb.temperatures_reconstructed(end,13:24) = tim1+3; md.smb.temperatures_reconstructed(end,25:36) = tim1+5;
+md.smb.precipitations_reconstructed(end,1:12) = tim1; md.smb.precipitations_reconstructed(end,13:24) = tim1+3; md.smb.precipitations_reconstructed(end,25:36) = tim1+5;
+
 % creating initialization and spc temperatures initialization and
 % spc
@@ -28 +52 @@
 md.initialization.temperature=md.thermal.spctemperature; %md.smb.temperatures_presentday(1:md.mesh.numberofvertices,1);
 
-% creating precipitation
-for imonth=0:11
-        md.smb.precipitations_presentday(1:md.mesh.numberofvertices,imonth+1)=-0.4*10^(-6)*md.mesh.y+0.5;
-end
 
 md.smb.s0p = (max(md.geometry.surface,0));
@@ -73 +93 @@
 md.friction.pressure_adjusted_temperature = temperature - Tm;
 md.friction.gamma= 5;
-
 
 % time steps and resolution