Index: /issm/trunk-jpl/src/c/classes/Elements/Element.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 24013)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 24014)
@@ -2178,5 +2178,5 @@
 	int         basinid,num_basins,M,N;
 	IssmDouble  tf,gamma0,base,delta_t_basin,mean_tf_basin,absval;
-	IssmDouble  basalmeltrate[NUM_VERTICES];
+	IssmDouble* basalmeltrate = xNew<IssmDouble>(NUM_VERTICES);
 	bool        islocal;
 	IssmDouble* delta_t = NULL;
@@ -2682,8 +2682,8 @@
 	IssmDouble g1               = area_boxi*gamma_T;
 
-	IssmDouble basalmeltrates_shelf[NUM_VERTICES];
-	IssmDouble potential_pressure_melting_point[NUM_VERTICES];
-	IssmDouble Tocs[NUM_VERTICES];
-	IssmDouble Socs[NUM_VERTICES];
+	IssmDouble* basalmeltrates_shelf 				= xNew<IssmDouble>(NUM_VERTICES);
+	IssmDouble* potential_pressure_melting_point	= xNew<IssmDouble>(NUM_VERTICES);
+	IssmDouble* Tocs 								= xNew<IssmDouble>(NUM_VERTICES);
+	IssmDouble* Socs								= xNew<IssmDouble>(NUM_VERTICES);
 
 	/* First box calculations */
@@ -2694,6 +2694,6 @@
 		this->parameters->FindParam(&toc_farocean, basinid, time, BasalforcingsPicoFarOceantemperatureEnum);
 		this->parameters->FindParam(&soc_farocean, basinid, time, BasalforcingsPicoFarOceansalinityEnum);
-		IssmDouble s1 = soc_farocean/(nu*lambda);
-		IssmDouble overturnings[NUM_VERTICES];
+		IssmDouble 	s1 				= soc_farocean/(nu*lambda);
+		IssmDouble* overturnings 	= xNew<IssmDouble>(NUM_VERTICES);
 
 		/* Start looping on the number of verticies and calculate ocean vars */
@@ -2793,5 +2793,5 @@
 
 	/*Define arrays*/
-	IssmDouble basalmeltrates_shelf[NUM_VERTICES];  //Basal melt-rate
+	IssmDouble* basalmeltrates_shelf = xNew<IssmDouble>(NUM_VERTICES); //Basal melt-rate
 
 	/*Polynomial coefficients*/