Index: /issm/trunk/test/NightlyRun/test1501.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1501.m	(revision 8998)
+++ /issm/trunk/test/NightlyRun/test1501.m	(revision 8999)
@@ -28,6 +28,6 @@
 smb=[ [ones(399,1)*(smb')]' smb  tooth tooth];
 
-md.forcings.surface_mass_balance= smb;
-md.forcings.surface_mass_balance(end+1,:)=[1:2000];
+md.surface_mass_balance= smb;
+md.surface_mass_balance(end+1,:)=[1:2000];
 
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test1502.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1502.m	(revision 8998)
+++ /issm/trunk/test/NightlyRun/test1502.m	(revision 8999)
@@ -29,6 +29,6 @@
 smb=[ [ones(399,1)*(smb')]' smb  tooth tooth];
 
-md.forcings.surface_mass_balance= smb;
-md.forcings.surface_mass_balance(end+1,:)=[1:2000];
+md.surface_mass_balance= smb;
+md.surface_mass_balance(end+1,:)=[1:2000];
 
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test255.m
===================================================================
--- /issm/trunk/test/NightlyRun/test255.m	(revision 8998)
+++ /issm/trunk/test/NightlyRun/test255.m	(revision 8999)
@@ -13,6 +13,6 @@
 smb=[ smb smb*-1 ];
 
-md.forcings.surface_mass_balance= smb;
-md.forcings.surface_mass_balance(end+1,:)=[1.5 3];
+md.surface_mass_balance= smb;
+md.surface_mass_balance(end+1,:)=[1.5 3];
 md.isthermal=0;
 
Index: /issm/trunk/test/NightlyRun/test256.m
===================================================================
--- /issm/trunk/test/NightlyRun/test256.m	(revision 8998)
+++ /issm/trunk/test/NightlyRun/test256.m	(revision 8999)
@@ -13,6 +13,6 @@
 smb=[ smb smb*-1 ];
 
-md.forcings.surface_mass_balance= smb;
-md.forcings.surface_mass_balance(end+1,:)=[1.5 3];
+md.surface_mass_balance= smb;
+md.surface_mass_balance(end+1,:)=[1.5 3];
 md.isthermal=0;
 
Index: /issm/trunk/test/NightlyRun/test257.m
===================================================================
--- /issm/trunk/test/NightlyRun/test257.m	(revision 8998)
+++ /issm/trunk/test/NightlyRun/test257.m	(revision 8999)
@@ -13,6 +13,6 @@
 smb=[ smb smb*2 ];
 
-md.forcings.surface_mass_balance= smb;
-md.forcings.surface_mass_balance(end+1,:)=[1.5 3];
+md.surface_mass_balance= smb;
+md.surface_mass_balance(end+1,:)=[1.5 3];
 md.isthermal=0;
 
Index: /issm/trunk/test/NightlyRun/test258.m
===================================================================
--- /issm/trunk/test/NightlyRun/test258.m	(revision 8998)
+++ /issm/trunk/test/NightlyRun/test258.m	(revision 8999)
@@ -13,6 +13,6 @@
 smb=[ smb smb*2 ];
 
-md.forcings.surface_mass_balance= smb;
-md.forcings.surface_mass_balance(end+1,:)=[1.5 3];
+md.surface_mass_balance= smb;
+md.surface_mass_balance(end+1,:)=[1.5 3];
 md.isthermal=0;
 
Index: /issm/trunk/test/NightlyRun/test259.m
===================================================================
--- /issm/trunk/test/NightlyRun/test259.m	(revision 8998)
+++ /issm/trunk/test/NightlyRun/test259.m	(revision 8999)
@@ -14,6 +14,6 @@
 smb=[ smb smb*-1 ];
 
-md.forcings.surface_mass_balance= smb;
-md.forcings.surface_mass_balance(end+1,:)=[1.5 3];
+md.surface_mass_balance= smb;
+md.surface_mass_balance(end+1,:)=[1.5 3];
 md.isthermal=0;
 
Index: /issm/trunk/test/NightlyRun/test260.m
===================================================================
--- /issm/trunk/test/NightlyRun/test260.m	(revision 8998)
+++ /issm/trunk/test/NightlyRun/test260.m	(revision 8999)
@@ -14,6 +14,6 @@
 smb=[ smb smb*-1 ];
 
-md.forcings.surface_mass_balance= smb;
-md.forcings.surface_mass_balance(end+1,:)=[1.5 3];
+md.surface_mass_balance= smb;
+md.surface_mass_balance(end+1,:)=[1.5 3];
 md.isthermal=0;
 
Index: /issm/trunk/test/NightlyRun/test261.m
===================================================================
--- /issm/trunk/test/NightlyRun/test261.m	(revision 8998)
+++ /issm/trunk/test/NightlyRun/test261.m	(revision 8999)
@@ -14,6 +14,6 @@
 smb=[ smb smb*2 ];
 
-md.forcings.surface_mass_balance= smb;
-md.forcings.surface_mass_balance(end+1,:)=[1.5 3];
+md.surface_mass_balance= smb;
+md.surface_mass_balance(end+1,:)=[1.5 3];
 md.isthermal=0;
 
Index: /issm/trunk/test/NightlyRun/test262.m
===================================================================
--- /issm/trunk/test/NightlyRun/test262.m	(revision 8998)
+++ /issm/trunk/test/NightlyRun/test262.m	(revision 8999)
@@ -14,6 +14,6 @@
 smb=[ smb smb*2 ];
 
-md.forcings.surface_mass_balance= smb;
-md.forcings.surface_mass_balance(end+1,:)=[1.5 3];
+md.surface_mass_balance= smb;
+md.surface_mass_balance(end+1,:)=[1.5 3];
 md.isthermal=0;
 
Index: /issm/trunk/test/Par/79North.par
===================================================================
--- /issm/trunk/test/Par/79North.par	(revision 8998)
+++ /issm/trunk/test/Par/79North.par	(revision 8999)
@@ -29,5 +29,4 @@
 md.basal_melting_rate(find(pos))=10;
 md.surface_mass_balance=15*ones(md.numberofnodes,1);
-md.forcings.surface_mass_balance=[15*ones(md.numberofnodes,1);1];
 
 %Numerical parameters
Index: /issm/trunk/test/Par/RoundSheetEISMINT.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 8998)
+++ /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 8999)
@@ -28,5 +28,4 @@
 rel=450*1000; %m
 md.surface_mass_balance=min(smb_max,sb*(rel-radius));
-md.forcings.surface_mass_balance=[min(smb_max,sb*(rel-radius));1];
 
 disp('      creating velocities');
Index: /issm/trunk/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetShelf.par	(revision 8998)
+++ /issm/trunk/test/Par/RoundSheetShelf.par	(revision 8999)
@@ -47,5 +47,4 @@
 %Surface mass balance and basal melting
 md.surface_mass_balance=-10*ones(md.numberofnodes,1);
-md.forcings.surface_mass_balance=[-10*ones(md.numberofnodes,1);1];
 md.basal_melting_rate=zeros(md.numberofnodes,1);
 pos=find(md.nodeoniceshelf);md.basal_melting_rate(pos)=10;
Index: /issm/trunk/test/Par/RoundSheetStaticEISMINT.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetStaticEISMINT.par	(revision 8998)
+++ /issm/trunk/test/Par/RoundSheetStaticEISMINT.par	(revision 8999)
@@ -34,5 +34,4 @@
 rel=450*1000; %m
 md.surface_mass_balance=min(smb_max,sb*(rel-radius));
-md.forcings.surface_mass_balance=[min(smb_max,sb*(rel-radius));1];
 
 disp('      creating velocities');
Index: /issm/trunk/test/Par/SquareEISMINT.par
===================================================================
--- /issm/trunk/test/Par/SquareEISMINT.par	(revision 8998)
+++ /issm/trunk/test/Par/SquareEISMINT.par	(revision 8999)
@@ -28,5 +28,4 @@
 disp('      creating surface mass balance');
 md.surface_mass_balance=0.2*ones(md.numberofnodes,1); %0m/a
-md.forcings.surface_mass_balance=[0.2*ones(md.numberofnodes,1);1]; %0m/a
 md.basal_melting_rate=0*ones(md.numberofnodes,1); %0m/a
 
Index: /issm/trunk/test/Par/SquareSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/SquareSheetShelf.par	(revision 8998)
+++ /issm/trunk/test/Par/SquareSheetShelf.par	(revision 8999)
@@ -29,5 +29,4 @@
 %Accumulation and melting
 md.surface_mass_balance=10*ones(md.numberofnodes,1);
-md.forcings.surface_mass_balance=[10*ones(md.numberofnodes,1);1];
 md.basal_melting_rate=5*ones(md.numberofnodes,1);
 
Index: /issm/trunk/test/Par/SquareShelfConstrained.par
===================================================================
--- /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 8998)
+++ /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 8999)
@@ -26,5 +26,4 @@
 %Surface mass balance and basal melting
 md.surface_mass_balance=10*ones(md.numberofnodes,1);
-md.forcings.surface_mass_balance=[10*ones(md.numberofnodes,1);1];
 md.basal_melting_rate=5*ones(md.numberofnodes,1);
 
Index: /issm/trunk/test/Par/SquareThermal.par
===================================================================
--- /issm/trunk/test/Par/SquareThermal.par	(revision 8998)
+++ /issm/trunk/test/Par/SquareThermal.par	(revision 8999)
@@ -33,5 +33,4 @@
 disp('      creating surface_mass_balance');
 md.surface_mass_balance=ones(md.numberofnodes,1)/md.yts; %1m/a
-md.forcings.surface_mass_balance=[ones(md.numberofnodes,1)/md.yts;1]; %1m/a
 md.basal_melting_rate=0*ones(md.numberofnodes,1)/md.yts; %1m/a