Fix bug resulting in successive additions of vturb_add (#33)

ITA-Solar · Oct 5, 2021 · 65105a9 · 65105a9
1 parent 43ccf74
commit 65105a9
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/rh15d/hdf5atmos.c b/rh15d/hdf5atmos.c
@@ -238,6 +238,10 @@ void readAtmos_hdf5(int xi, int yi, Atmosphere *atmos, Geometry *geometry,
     setTcut(atmos, geometry, input.p15d_tmax);
   } else {
     mpi.zcut = 0;
+    /*Free vturb in order to avoid successive addition of vturb_add. 
+      This is done in realloc_ndep (called by setTcut) if 15D_DEFINE_ZCUT = TRUE*/
+    free(atmos->vturb);
+    atmos->vturb  = (double *) calloc(atmos->Nspace , sizeof(double));
   }
 
   /* Memory dataspace, redefine for Nspace */