Index: /issm/trunk-jpl/src/m/classes/stochasticforcing.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/stochasticforcing.m	(revision 26616)
+++ /issm/trunk-jpl/src/m/classes/stochasticforcing.m	(revision 26617)
@@ -89,5 +89,5 @@
             md = checkfield(md,'fieldname','stochasticforcing.defaultdimension','numel',1,'NaN',1,'Inf',1,'>',0);
             md = checkfield(md,'fieldname','stochasticforcing.default_id','Inf',1,'>=',0,'<=',self.defaultdimension,'size',[md.mesh.numberofelements,1]);
-         end
+			end
 		end % }}}
 		function disp(self) % {{{
@@ -142,4 +142,8 @@
 					end
 				end
+				%Set dummy default_id vector if defaults not used
+				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 26616)
+++ /issm/trunk-jpl/src/m/classes/stochasticforcing.py	(revision 26617)
@@ -132,4 +132,7 @@
                     for col in inds: # scale columns corresponding to scaled field
                         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);
             WriteData(fid, prefix, 'data', num_fields, 'name', 'md.stochasticforcing.num_fields', 'format', 'Integer')
             WriteData(fid, prefix, 'object', self, 'fieldname', 'fields', 'format', 'StringArray')