Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 23803)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 23804)
@@ -2526,5 +2526,5 @@
 		TF_input->GetInputValue(&TF,gauss);
 
-		if(basinid[iv]==0 | basin_icefront_area[reCast<int>(basinid[iv])-1]==0.) meltrates[iv]=0.;
+		if(basinid[iv]==0 || basin_icefront_area[reCast<int>(basinid[iv])-1]==0.) meltrates[iv]=0.;
 		else{
 			/* change the unit of qsg (m^3/d -> m/d) with ice front area */
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 23803)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 23804)
@@ -3448,5 +3448,5 @@
 		TF_input->GetInputValue(&TF,gauss);
 
-		if(basinid[iv]==0 | basin_icefront_area[reCast<int>(basinid[iv])-1]==0.) meltrates[iv]=0.;
+		if(basinid[iv]==0 || basin_icefront_area[reCast<int>(basinid[iv])-1]==0.) meltrates[iv]=0.;
 		else{
 			/* change the unit of qsg (m^3/d -> m/d) with ice front area */