Index: /issm/trunk-jpl/test/NightlyRun/test247.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test247.m	(revision 24623)
+++ /issm/trunk-jpl/test/NightlyRun/test247.m	(revision 24624)
@@ -1,3 +1,3 @@
-%Test Name: SquareShelfTranIspddIsdeltaO18pdNoInterpSSA2d(
+%Test Name: SquareShelfTranIspddIsdeltaO18pdNoInterpSSA2d
 md=triangle(model(),'../Exp/Square.exp',60000.);
 md=setmask(md,'all','');
@@ -14,4 +14,6 @@
 load '../Data/delta18o.data'
 md.smb.delta18o=delta18o;
+md.smb.istemperaturescaled = 0;
+md.smb.isprecipscaled = 0;
 
 % creating Present day  temperatures
@@ -23,4 +25,26 @@
 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,8 +52,4 @@
 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,5 +93,4 @@
 md.friction.pressure_adjusted_temperature = temperature - Tm;
 md.friction.gamma= 5;
-
 
 % time steps and resolution