Index: /issm/trunk-jpl/src/m/mech/newforcing.m
===================================================================
--- /issm/trunk-jpl/src/m/mech/newforcing.m	(revision 20186)
+++ /issm/trunk-jpl/src/m/mech/newforcing.m	(revision 20186)
@@ -0,0 +1,35 @@
+function forcing=newforcing(t0,t1,deltaT,f0,f1,nodes)
+%FUNCTION NEWFORCING Build forcing that extends temporally from t0 to t1, and in magnitude from f0 to f1. Equal time 
+%                    and magnitude spacing. 
+%
+%       Usage: forcing=newforcing(t0,t1,deltaT,f0,f1,nodes);  
+%       Where: 
+%          t0:t1: time interval. 
+%          deltaT: time step
+%          f0:f1: magnitude interval.
+%          nodes: number of vertices where we have a temporal forcing
+%
+%       Example: 
+%           md.smb.mass_balance=newforcing(md.timestepping.start_time,md.timestepping.final_time,...
+%                                          md.timestepping.time_step,-1,+2,md.mesh.numberofvertices);
+%
+
+	%some cheks:
+	if nargin~=6, error('newforcing error message: should have 5 input arguments'); end
+
+	%Number of time steps: 
+	nsteps=(t1-t0)/deltaT+1;
+
+	%delta forcing:
+	deltaf=(f1-f0)/(nsteps-1);
+
+	%creates times:
+	times=t0:deltaT:t1;
+
+	%create forcing:
+	forcing=(f0:deltaf:f1);
+
+	%replicate for all nodes
+	forcing=repmat(forcing,nodes+1,1);
+	forcing(end,:)=times;
+