Index: /issm/trunk-jpl/src/c/classes/Elements/Element.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 22472)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 22473)
@@ -274,4 +274,5 @@
 		virtual int        NumberofNodesVelocity(void)=0;
 		virtual void       PicoUpdateBoxid(int* pmax_boxid_basin)=0;
+		virtual void       PicoUpdateFirstBox(void){_error_("not implemented");};
 		virtual void       PotentialUngrounding(Vector<IssmDouble>* potential_sheet_ungrounding)=0;
 		virtual int        PressureInterpolation()=0;
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 22472)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 22473)
@@ -2788,5 +2788,5 @@
 	rhow			= this->GetMaterialParameter(MaterialsRhoSeawaterEnum);
 	earth_grav  = this->GetMaterialParameter(ConstantsGEnum);
-	rho_star		= 1033;             //kg/m^3
+	rho_star		= 1033.;             //kg/m^3
 	nu				= rhoi/rhow;
 	latentheat	= this->GetMaterialParameter(MaterialsLatentheatEnum);
@@ -2801,6 +2801,6 @@
 	this->parameters->FindParam(&gamma_T,BasalforcingsPicoGammaTEnum);
 	this->parameters->FindParam(&overturning_coeff,BasalforcingsPicoOverturningCoeffEnum);
-	this->parameters->FindParam(&t_farocean,BasalforcingsPicoFarOceanTemperatureEnum);
-	this->parameters->FindParam(&s_farocean,BasalforcingsPicoFarOceanSalinityEnum);
+	this->parameters->FindParam(&t_farocean,BasalforcingsPicoFarOceansalinityEnum);
+	this->parameters->FindParam(&s_farocean,BasalforcingsPicoFarOceansalinityEnum);
 
 	this->inputs->GetInputValue(&basinid,BasalforcingsPicoBasinIdEnum);
@@ -2808,4 +2808,7 @@
 	this->inputs->GetInputValue(&maxbox,BasalforcingsPicoMaxboxcountEnum);
 	this->inputs->GetInputValue(&thickness,ThicknessEnum);
+
+	_error_("to be continued");
+	/*
 
 	toc_farocean = t_farocean[basinid];
@@ -2827,4 +2830,5 @@
 
 	this->AddInput(BasalforcingsGroundediceMeltingRateEnum,basalmeltingrates,this->GetElementType());
+	*/
 
 }
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 22472)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 22473)
@@ -111,4 +111,5 @@
 		int         NumberofNodesVelocity(void);
 		void        PicoUpdateBoxid(int* pmax_boxid_basin);
+		void        PicoUpdateFirstBox(void);
 		void        PotentialUngrounding(Vector<IssmDouble>* potential_sheet_ungrounding);
 		int         PressureInterpolation();
Index: /issm/trunk-jpl/src/c/modules/FloatingiceMeltingRatex/FloatingiceMeltingRatex.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/FloatingiceMeltingRatex/FloatingiceMeltingRatex.cpp	(revision 22472)
+++ /issm/trunk-jpl/src/c/modules/FloatingiceMeltingRatex/FloatingiceMeltingRatex.cpp	(revision 22473)
@@ -107,5 +107,5 @@
 
 	//2 Get area of the boxes
-	boxareas=new Vector<IssmDouble>(num_basins*maxbox);
+	boxareas=xNew<IssmDouble>(num_basins*maxbox);
 	for(i=0;i<num_basins*maxbox;i++){boxareas[i]=0.;}
 	for(i=0;i<femmodel->elements->Size();i++){
@@ -114,5 +114,6 @@
 		element->inputs->GetInputValue(&basinid,BasalforcingsPicoBasinIdEnum);
 		element->inputs->GetInputValue(&boxid,BasalforcingsPicoBoxIdEnum);
-		boxareas[basinid*maxbox+boxid]+=element->GetArea();
+		_error_("fix thi");
+		//boxareas[basinid*maxbox+boxid]+=element->GetArea();
 	}
 
@@ -138,4 +139,3 @@
 	xDelete<IssmDouble>(t_farocean);
    xDelete<IssmDouble>(s_farocean);
-	xDelete<IssmDouble>(element);
 }/*}}}*/
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp	(revision 22472)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp	(revision 22473)
@@ -213,7 +213,7 @@
 				parameters->AddObject(iomodel->CopyConstantObject("md.basalforcings.gamma_T",BasalforcingsPicoGammaTEnum));
 				iomodel->FetchData(&transparam,&N,&M,"md.basalforcings.farocean_temperature");
-				parameters->AddObject(new TransientArrayParam(BasalforcingsPicoFarOceanTemperatureEnum,&transparam[0],&transparam[M*(N-1)],interp,M,N));
+				parameters->AddObject(new TransientArrayParam(BasalforcingsPicoFarOceantemperatureEnum,&transparam[0],&transparam[M*(N-1)],interp,M,N));
 				iomodel->FetchData(&transparam,&N,&M,"md.basalforcings.farocean_salinity");
-				parameters->AddObject(new TransientArrayParam(BasalforcingsPicoFarOceanSalinityEnum,&transparam[0],&transparam[M*(N-1)],interp,M,N));
+				parameters->AddObject(new TransientArrayParam(BasalforcingsPicoFarOceansalinityEnum,&transparam[0],&transparam[M*(N-1)],interp,M,N));
 			break;
 		default: