Index: /issm/trunk-jpl/src/m/classes/flowequation.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/flowequation.py	(revision 15654)
+++ /issm/trunk-jpl/src/m/classes/flowequation.py	(revision 15655)
@@ -52,4 +52,8 @@
 		#}}}
 	def setdefaultparameters(self): # {{{
+
+		#MINI element for FS by default
+		self.fe_FS=3
+
 		return self
 	#}}}
@@ -64,5 +68,5 @@
 			md = checkfield(md,'flowequation.fe_SSA','numel',[1],'values',[0,1])
 			md = checkfield(md,'flowequation.fe_HO','numel',[1],'values',numpy.arange(0,3+1))
-			md = checkfield(md,'flowequation.fe_FS','numel',[1],'values',[0])
+			md = checkfield(md,'flowequation.fe_FS','numel',[1],'values',numpy.arange(0,3+1))
 			md = checkfield(md,'flowequation.borderSSA','size',[md.mesh.numberofvertices],'values',[0,1])
 			md = checkfield(md,'flowequation.borderHO','size',[md.mesh.numberofvertices],'values',[0,1])
Index: /issm/trunk-jpl/src/m/solve/WriteData.m
===================================================================
--- /issm/trunk-jpl/src/m/solve/WriteData.m	(revision 15654)
+++ /issm/trunk-jpl/src/m/solve/WriteData.m	(revision 15655)
@@ -39,4 +39,8 @@
 		data  = scale.*data;
 	end
+end
+if(size(data,1)==forcinglength),
+	yts=365.0*24.0*3600.0;
+	data(end,:) = data(end,:)*yts;
 end
 
Index: /issm/trunk-jpl/src/m/solve/WriteData.py
===================================================================
--- /issm/trunk-jpl/src/m/solve/WriteData.py	(revision 15654)
+++ /issm/trunk-jpl/src/m/solve/WriteData.py	(revision 15655)
@@ -53,4 +53,8 @@
 		else:
 			data  = scale*data
+	if numpy.size(data) > 1 :
+		if numpy.size(data,0)==forcinglength:
+			yts=365.0*24.0*3600.0
+			data[-1,:] = yts*data[-1,:]
 
 	#Step 1: write the enum to identify this record uniquely