Index: /issm/trunk/test/Par/79North.par
===================================================================
--- /issm/trunk/test/Par/79North.par	(revision 8367)
+++ /issm/trunk/test/Par/79North.par	(revision 8368)
@@ -29,5 +29,5 @@
 md.melting_rate(find(pos))=10;
 md.accumulation_rate=15*ones(md.numberofnodes,1);
-md.forcing.accumulation_rate=[15*ones(md.numberofnodes,1);1];
+md.forcings.accumulation_rate=[15*ones(md.numberofnodes,1);1];
 
 %Numerical parameters
Index: /issm/trunk/test/Par/RoundSheetEISMINT.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 8367)
+++ /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 8368)
@@ -28,5 +28,5 @@
 rel=450*1000; %m
 md.accumulation_rate=min(acc_max,sb*(rel-radius));
-md.forcing.accumulation_rate=[min(acc_max,sb*(rel-radius));1];
+md.forcings.accumulation_rate=[min(acc_max,sb*(rel-radius));1];
 
 disp('      creating velocities');
Index: /issm/trunk/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetShelf.par	(revision 8367)
+++ /issm/trunk/test/Par/RoundSheetShelf.par	(revision 8368)
@@ -47,5 +47,5 @@
 %Accumulation and melting
 md.accumulation_rate=-10*ones(md.numberofnodes,1);
-md.forcing.accumulation_rate=[-10*ones(md.numberofnodes,1);1];
+md.forcings.accumulation_rate=[-10*ones(md.numberofnodes,1);1];
 md.melting_rate=zeros(md.numberofnodes,1);
 pos=find(md.nodeoniceshelf);md.melting_rate(pos)=10;
Index: /issm/trunk/test/Par/RoundSheetStaticEISMINT.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetStaticEISMINT.par	(revision 8367)
+++ /issm/trunk/test/Par/RoundSheetStaticEISMINT.par	(revision 8368)
@@ -34,5 +34,5 @@
 rel=450*1000; %m
 md.accumulation_rate=min(acc_max,sb*(rel-radius));
-md.forcing.accumulation_rate=[min(acc_max,sb*(rel-radius));1];
+md.forcings.accumulation_rate=[min(acc_max,sb*(rel-radius));1];
 
 disp('      creating velocities');
Index: /issm/trunk/test/Par/SquareEISMINT.par
===================================================================
--- /issm/trunk/test/Par/SquareEISMINT.par	(revision 8367)
+++ /issm/trunk/test/Par/SquareEISMINT.par	(revision 8368)
@@ -28,5 +28,5 @@
 disp('      creating accumulation rates');
 md.accumulation_rate=0.2*ones(md.numberofnodes,1); %0m/a
-md.forcing.accumulation_rate=[0.2*ones(md.numberofnodes,1);1]; %0m/a
+md.forcings.accumulation_rate=[0.2*ones(md.numberofnodes,1);1]; %0m/a
 md.melting_rate=0*ones(md.numberofnodes,1); %0m/a
 
Index: /issm/trunk/test/Par/SquareSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/SquareSheetShelf.par	(revision 8367)
+++ /issm/trunk/test/Par/SquareSheetShelf.par	(revision 8368)
@@ -29,5 +29,5 @@
 %Accumulation and melting
 md.accumulation_rate=10*ones(md.numberofnodes,1);
-md.forcing.accumulation_rate=[10*ones(md.numberofnodes,1);1];
+md.forcings.accumulation_rate=[10*ones(md.numberofnodes,1);1];
 md.melting_rate=5*ones(md.numberofnodes,1);
 
Index: /issm/trunk/test/Par/SquareShelfConstrained.par
===================================================================
--- /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 8367)
+++ /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 8368)
@@ -26,5 +26,5 @@
 %Accumulation and melting
 md.accumulation_rate=10*ones(md.numberofnodes,1);
-md.forcing.accumulation_rate=[10*ones(md.numberofnodes,1);1];
+md.forcings.accumulation_rate=[10*ones(md.numberofnodes,1);1];
 md.melting_rate=5*ones(md.numberofnodes,1);
 
Index: /issm/trunk/test/Par/SquareThermal.par
===================================================================
--- /issm/trunk/test/Par/SquareThermal.par	(revision 8367)
+++ /issm/trunk/test/Par/SquareThermal.par	(revision 8368)
@@ -35,5 +35,5 @@
 disp('      creating accumulation rates');
 md.accumulation_rate=ones(md.numberofnodes,1)/md.yts; %1m/a
-md.forcing.accumulation_rate=[ones(md.numberofnodes,1)/md.yts;1]; %1m/a
+md.forcings.accumulation_rate=[ones(md.numberofnodes,1)/md.yts;1]; %1m/a
 md.melting_rate=0*ones(md.numberofnodes,1)/md.yts; %1m/a