Index: /issm/trunk-jpl/src/m/classes/SMBcomponents.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBcomponents.m	(revision 26208)
+++ /issm/trunk-jpl/src/m/classes/SMBcomponents.m	(revision 26209)
@@ -12,4 +12,5 @@
 		averaging = 0;
 		requested_outputs     = {};
+		isclimatology;
 	end
 	methods
Index: /issm/trunk-jpl/src/m/classes/SMBforcing.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBforcing.m	(revision 26208)
+++ /issm/trunk-jpl/src/m/classes/SMBforcing.m	(revision 26209)
@@ -10,4 +10,5 @@
 		requested_outputs = {};
 		averaging         = 0;
+		isclimatology;
 	end
 	methods
Index: /issm/trunk-jpl/src/m/classes/SMBgemb.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBgemb.m	(revision 26208)
+++ /issm/trunk-jpl/src/m/classes/SMBgemb.m	(revision 26209)
@@ -23,4 +23,5 @@
 		isturbulentflux;
 		isconstrainsurfaceT;
+		isclimatology;
 
 		%inputs:
Index: /issm/trunk-jpl/src/m/classes/SMBmeltcomponents.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBmeltcomponents.m	(revision 26208)
+++ /issm/trunk-jpl/src/m/classes/SMBmeltcomponents.m	(revision 26209)
@@ -13,4 +13,5 @@
 		averaging = 0;
 		requested_outputs      = {};
+		isclimatology;
 	end
 	methods
Index: /issm/trunk-jpl/src/m/classes/model.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.m	(revision 26208)
+++ /issm/trunk-jpl/src/m/classes/model.m	(revision 26209)
@@ -168,9 +168,11 @@
 			%2019 Mar 28, updated 2021 April 23
 			if isa(md.smb,'SMBcomponents') | isa(md.smb,'SMBmeltcomponents') | isa(md.smb,'SMBforcing') | isa(md.smb,'SMBgemb') 
-				if isa(md.smb.isclimatology,'double')
-					if prod(size(md.smb.isclimatology)) ~= 1
-						md.smb.isclimatology = 0;
+				if any(strcmp(fieldnames(md.smb),'isclimatology'))
+					if isa(md.smb.isclimatology,'double')
+						if prod(size(md.smb.isclimatology)) ~= 1
+							md.smb.isclimatology = 0;
+						end
+						md.timestepping.cycle_forcing=md.smb.isclimatology;
 					end
-					md.timestepping.cycle_forcing=md.smb.isclimatology;
 				end
 			end
