Changeset 28193

Timestamp:

03/29/24 11:39:58 (12 months ago)

Author:

Mathieu Morlighem

Message:

BUG: avoid created NaNs if basin is not present (would divide by 0, which makes CoDiPack trip)

File:

• issm/trunk-jpl/src/c/modules/FloatingiceMeltingRatex/FloatingiceMeltingRatex.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/modules/FloatingiceMeltingRatex/FloatingiceMeltingRatex.cpp

@@ -193 +193 @@
                 ISSM_MPI_Allreduce(areas_summed,areas_summed_cpu,num_basins,ISSM_MPI_DOUBLE,ISSM_MPI_SUM,IssmComm::GetComm());
 
+                /*Make sure Area is not zero to avoid dividing by 0 if a basin is not present in the model*/
+                for(int k=0;k<num_basins;k++) if(areas_summed_cpu[k]==0.) areas_summed_cpu[k] = 1.;
+
                 /*Compute weighted means and save*/
                 for(int k=0;k<num_basins;k++){tf_weighted_avg_cpu[k] = tf_weighted_avg_cpu[k]/areas_summed_cpu[k];}