Index: /issm/trunk-jpl/src/m/classes/stochasticforcing.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/stochasticforcing.m	(revision 26617)
+++ /issm/trunk-jpl/src/m/classes/stochasticforcing.m	(revision 26618)
@@ -143,7 +143,7 @@
 				end
 				%Set dummy default_id vector if defaults not used
-				if isnan(self.default_id)
-					self.default_id = zeros(md.mesh.numberofelements,1);
-				end
+				%if isnan(self.default_id)
+				%	self.default_id = zeros(md.mesh.numberofelements,1);
+				%end
 				WriteData(fid,prefix,'data',num_fields,'name','md.stochasticforcing.num_fields','format','Integer');
 				WriteData(fid,prefix,'object',self,'fieldname','fields','format','StringArray');
Index: /issm/trunk-jpl/src/m/classes/stochasticforcing.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/stochasticforcing.py	(revision 26617)
+++ /issm/trunk-jpl/src/m/classes/stochasticforcing.py	(revision 26618)
@@ -133,6 +133,6 @@
                         tempcovariance[:, col] = 1 / yts * tempcovariance[:, col]
             #Set dummy default_id vector if defaults not used
-            if np.isnan(self.default_id)
-               self.default_id = np.zeros(md.mesh.numberofelements);
+            #if np.isnan(self.default_id)
+              # self.default_id = np.zeros(md.mesh.numberofelements);
             WriteData(fid, prefix, 'data', num_fields, 'name', 'md.stochasticforcing.num_fields', 'format', 'Integer')
             WriteData(fid, prefix, 'object', self, 'fieldname', 'fields', 'format', 'StringArray')