Index: /issm/trunk-jpl/src/m/classes/miscellaneous.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/miscellaneous.m	(revision 27628)
+++ /issm/trunk-jpl/src/m/classes/miscellaneous.m	(revision 27629)
@@ -25,4 +25,8 @@
 		end % }}}
 		function self = setdefaultparameters(self) % {{{
+
+			%Add some information about the model for future reference
+			self.notes = ['Model created on ' date ' by ' char(java.lang.System.getProperty('user.name')) ', ISSM v' num2str(issmversion()) ' (path: ' pwd() ')'];
+			self.name  = '';
 
 		end % }}}
Index: /issm/trunk-jpl/src/m/classes/spatiallinearbasalforcings.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/spatiallinearbasalforcings.m	(revision 27628)
+++ /issm/trunk-jpl/src/m/classes/spatiallinearbasalforcings.m	(revision 27629)
@@ -21,13 +21,15 @@
 				case 1
 					lb=varargin{1};
-					if strcmpi(class(lb),'linearbasalforcings'),
+					if isa(lb,'linearbasalforcings');
 						nvertices=length(lb.groundedice_melting_rate);
 						self.groundedice_melting_rate=lb.groundedice_melting_rate;
 						self.geothermalflux=lb.geothermalflux;
-						self.deepwater_elevation=lb.deepwater_elevation*ones(nvertices,1);
-						self.deepwater_melting_rate=lb.deepwater_melting_rate*ones(nvertices,1);
-						self.upperwater_melting_rate=lb.upperwater_melting_rate*ones(nvertices,1);
-						self.upperwater_elevation=lb.upperwater_elevation*ones(nvertices,1);
-						self.perturbation_melting_rate=lb.perturbation_melting_rate*ones(nvertices,1);
+						self.deepwater_elevation       = lb.deepwater_elevation*ones(nvertices,1);
+						self.deepwater_melting_rate    = lb.deepwater_melting_rate*ones(nvertices,1);
+						self.upperwater_melting_rate   = lb.upperwater_melting_rate*ones(nvertices,1);
+						self.upperwater_elevation      = lb.upperwater_elevation*ones(nvertices,1);
+						if ~isnan(lb.perturbation_melting_rate)
+							self.perturbation_melting_rate = lb.perturbation_melting_rate*ones(nvertices,1);
+						end
 					else 
 						self=structtoobj(spatiallinearbasalforcings(),varargin{1});