Index: /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 27009)
+++ /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 27010)
@@ -146,5 +146,7 @@
 		case FrontalForcingsDefaultEnum:
 			iomodel->FetchDataToInput(inputs,elements,"md.frontalforcings.meltingrate",CalvingMeltingrateEnum);
-			iomodel->FetchDataToInput(inputs,elements,"md.frontalforcings.ablationrate",CalvingAblationrateEnum);
+			if (calvinglaw == CalvingParameterizationEnum) {
+				iomodel->FetchDataToInput(inputs,elements,"md.frontalforcings.ablationrate",CalvingAblationrateEnum);
+			}
 			break;
 		case FrontalForcingsRignotEnum:
Index: /issm/trunk-jpl/src/m/classes/frontalforcings.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/frontalforcings.m	(revision 27009)
+++ /issm/trunk-jpl/src/m/classes/frontalforcings.m	(revision 27010)
@@ -42,5 +42,7 @@
 
 			md = checkfield(md,'fieldname','frontalforcings.meltingrate','NaN',1,'Inf',1,'timeseries',1,'>=',0);
-			md = checkfield(md,'fieldname','frontalforcings.ablationrate','NaN',1,'Inf',1,'timeseries',1);
+			if ~isnan(md.frontalforcings.ablationrate)
+				md = checkfield(md,'fieldname','frontalforcings.ablationrate','Inf',1,'timeseries',1);
+			end
 
 		end % }}}
@@ -54,5 +56,7 @@
 			WriteData(fid,prefix,'name','md.frontalforcings.parameterization','data',1,'format','Integer');
 			WriteData(fid,prefix,'object',self,'fieldname','meltingrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1./yts);
-			WriteData(fid,prefix,'object',self,'fieldname','ablationrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1./yts);
+			if ~isnan(md.frontalforcings.ablationrate)
+				WriteData(fid,prefix,'object',self,'fieldname','ablationrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1./yts);
+			end
 		end % }}}
 	end