Index: /issm/trunk-jpl/src/c/analyses/FreeSurfaceBaseAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/FreeSurfaceBaseAnalysis.cpp	(revision 27278)
+++ /issm/trunk-jpl/src/c/analyses/FreeSurfaceBaseAnalysis.cpp	(revision 27279)
@@ -89,7 +89,9 @@
 	/*Get what we need for ocean-induced basal melting*/
 	bool isstochastic;
-	int basalforcing_model;
-	iomodel->FindConstant(&basalforcing_model,"md.basalforcings.model");
-	iomodel->FindConstant(&isstochastic,"md.stochasticforcing.isstochasticforcing");
+   int basalforcing_model;
+   int melt_parameterization;
+   iomodel->FindConstant(&basalforcing_model,"md.basalforcings.model");
+   iomodel->FindConstant(&isstochastic,"md.stochasticforcing.isstochasticforcing");
+   iomodel->FindConstant(&melt_parameterization,"md.frontalforcings.parameterization");
 	switch(basalforcing_model){
 		case FloatingMeltRateEnum:
@@ -132,5 +134,5 @@
             iomodel->FetchDataToInput(inputs,elements,"md.basalforcings.ocean_temp",BasalforcingsOceanTempEnum);
          }
-         else{
+         else if(melt_parameterization!=FrontalForcingsRignotarmaEnum){
             iomodel->FetchDataToInput(inputs,elements,"md.basalforcings.ocean_thermalforcing",FrontalForcingsAndBasalforcingsThermalForcingEnum);
          }
Index: /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp	(revision 27278)
+++ /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp	(revision 27279)
@@ -163,6 +163,8 @@
 	bool isstochastic;
 	int basalforcing_model;
+	int melt_parameterization;
 	iomodel->FindConstant(&basalforcing_model,"md.basalforcings.model");
 	iomodel->FindConstant(&isstochastic,"md.stochasticforcing.isstochasticforcing");
+	iomodel->FindConstant(&melt_parameterization,"md.frontalforcings.parameterization");
 	switch(basalforcing_model){
 		case FloatingMeltRateEnum:
@@ -222,5 +224,5 @@
 				iomodel->FetchDataToInput(inputs,elements,"md.basalforcings.ocean_temp",BasalforcingsOceanTempEnum);
 			}
-			else{
+			else if(melt_parameterization!=FrontalForcingsRignotarmaEnum){
 				iomodel->FetchDataToInput(inputs,elements,"md.basalforcings.ocean_thermalforcing",FrontalForcingsAndBasalforcingsThermalForcingEnum);
 			}
Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 27278)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 27279)
@@ -817,7 +817,9 @@
 		/*Add basal forcings to compute melt rate*/
 		bool isstochastic;
-		int basalforcing_model;
-		iomodel->FindConstant(&basalforcing_model,"md.basalforcings.model");
-		iomodel->FindConstant(&isstochastic,"md.stochasticforcing.isstochasticforcing");
+	   int basalforcing_model;
+	   int melt_parameterization;
+	   iomodel->FindConstant(&basalforcing_model,"md.basalforcings.model");
+	   iomodel->FindConstant(&isstochastic,"md.stochasticforcing.isstochasticforcing");
+	   iomodel->FindConstant(&melt_parameterization,"md.frontalforcings.parameterization");
 		switch(basalforcing_model){
 			case FloatingMeltRateEnum:
@@ -860,5 +862,5 @@
       	      iomodel->FetchDataToInput(inputs,elements,"md.basalforcings.ocean_temp",BasalforcingsOceanTempEnum);
       	   }
-      	   else{
+      	   else if(melt_parameterization!=FrontalForcingsRignotarmaEnum){
       	      iomodel->FetchDataToInput(inputs,elements,"md.basalforcings.ocean_thermalforcing",FrontalForcingsAndBasalforcingsThermalForcingEnum);
       	   }
Index: /issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp	(revision 27278)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp	(revision 27279)
@@ -120,7 +120,9 @@
 	/*Add basal forcings to compute melt rate*/
 	bool isstochastic;
-	int basalforcing_model;
-	iomodel->FindConstant(&basalforcing_model,"md.basalforcings.model");
-	iomodel->FindConstant(&isstochastic,"md.stochasticforcing.isstochasticforcing");
+   int basalforcing_model;
+   int melt_parameterization;
+   iomodel->FindConstant(&basalforcing_model,"md.basalforcings.model");
+   iomodel->FindConstant(&isstochastic,"md.stochasticforcing.isstochasticforcing");
+   iomodel->FindConstant(&melt_parameterization,"md.frontalforcings.parameterization");
 	switch(basalforcing_model){
 		case FloatingMeltRateEnum:
@@ -162,5 +164,5 @@
             iomodel->FetchDataToInput(inputs,elements,"md.basalforcings.ocean_temp",BasalforcingsOceanTempEnum);
          }
-         else{
+         else if(melt_parameterization!=FrontalForcingsRignotarmaEnum){
             iomodel->FetchDataToInput(inputs,elements,"md.basalforcings.ocean_thermalforcing",FrontalForcingsAndBasalforcingsThermalForcingEnum);
          }