Changeset 24014

Timestamp:

06/12/19 15:10:30 (6 years ago)

Author:

jdquinn

Message:

BUG: Changed all array of IssmDouble initializations from IssmDouble <name>[NUM_VERTICES] to IssmDouble* <name> = xNew<IssmDouble>(NUM_VERTICES)

File:

• issm/trunk-jpl/src/c/classes/Elements/Element.cpp (modified) (4 diffs)

issm/trunk-jpl/src/c/classes/Elements/Element.cpp

@@ -2178 +2178 @@
         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 +2682 @@
         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 +2694 @@
                 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 +2793 @@
 
         /*Define arrays*/
-        IssmDouble basalmeltrates_shelf[NUM_VERTICES];  //Basal melt-rate
+        IssmDouble* basalmeltrates_shelf = xNew<IssmDouble>(NUM_VERTICES); //Basal melt-rate
 
         /*Polynomial coefficients*/