casting Cs_w values explicitly to float

Author: erikvansebille

parcels/compilation/codegenerator.py

@@ -830,7 +830,9 @@ def visit_FieldEvalNode(self, node):
         args = self._check_FieldSamplingArguments(node.args.ccode)
         if "croco" in node.field.obj.gridindexingtype and node.field.obj.name != "H" and node.field.obj.name != "Zeta":
             # Get Cs_w values directly from fieldset (since they are 1D in vertical only)
-            Cs_w = [self.fieldset.Cs_w.data[0][zi][0][0] for zi in range(self.fieldset.Cs_w.data.shape[1])]
+            Cs_w = [
+                self.fieldset.Cs_w.data[0][zi][0][0].astype(float) for zi in range(self.fieldset.Cs_w.data.shape[1])
+            ]
             statements_croco = [
                 c.Statement(f"float cs_w[] = {*Cs_w, }".replace("(", "{").replace(")", "}")),
                 c.Statement(

parcels/include/parcels.h

@@ -1253,8 +1253,8 @@ static inline double croco_from_z_to_sigma(CField *U, CField *H, CField *Zeta,
   float *sigma_levels = grid->depth;
   int zdim = grid->zdim;
   float zvec[zdim];
-  status = temporal_interpolation(x, y, 0, time, H, xi, yi, zi, ti, &local_h, LINEAR, CROCO);
-  status = temporal_interpolation(x, y, 0, time, Zeta, xi, yi, zi, ti, &local_zeta, LINEAR, CROCO);
+  status = temporal_interpolation(x, y, 0, time, H, xi, yi, zi, ti, &local_h, LINEAR, CROCO); CHECKSTATUS(status);
+  status = temporal_interpolation(x, y, 0, time, Zeta, xi, yi, zi, ti, &local_zeta, LINEAR, CROCO); CHECKSTATUS(status);
   for (zii = 0; zii < zdim; zii++)  {
     z0 = hc*sigma_levels[zii] + (local_h - hc) *cs_w[zii];
     zvec[zii] = z0 + local_zeta * (1 + z0 / local_h);