Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 26910)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 26911)
@@ -5500,7 +5500,8 @@
 	/*Get inputs*/
 	IssmDouble flux = 0.;
-	IssmDouble calvingratex,calvingratey,vx,vy,vel,meltingrate,meltingratex,meltingratey,thickness,Jdet;
+	IssmDouble calvingratex,calvingratey,vx,vy,vel,meltingrate,meltingratex,meltingratey,thickness,Jdet,groundedice;
 	IssmDouble rho_ice=FindParam(MaterialsRhoIceEnum);
-	Input* thickness_input=this->GetInput(ThicknessEnum); _assert_(thickness_input);
+	Input* thickness_input = this->GetInput(ThicknessEnum); _assert_(thickness_input);
+	Input *gr_input        = this->GetInput(MaskOceanLevelsetEnum); _assert_(gr_input);
 	Input* calvingratex_input=NULL;
 	Input* calvingratey_input=NULL;
@@ -5526,8 +5527,10 @@
 		calvingratex_input->GetInputValue(&calvingratex,gauss);
 		calvingratey_input->GetInputValue(&calvingratey,gauss);
+		gr_input->GetInputValue(&groundedice,gauss);
 		vx_input->GetInputValue(&vx,gauss);
 		vy_input->GetInputValue(&vy,gauss);
 		vel=vx*vx+vy*vy;
 		meltingrate_input->GetInputValue(&meltingrate,gauss);
+		if(groundedice<0) meltingrate = 0.;
 		meltingratex=meltingrate*vx/(sqrt(vel)+1.e-14);
 		meltingratey=meltingrate*vy/(sqrt(vel)+1.e-14);