add some conditionals to prevent segfault

tduretz · Feb 13, 2024 · 435ba89 · 435ba89
1 parent cc60f31
commit 435ba89
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/MDLIB/Setup.c b/MDLIB/Setup.c
@@ -195,9 +195,9 @@ void SetBCs(SetBCs_ff setBCs, MdoodzInput *instance, grid *mesh, markers *topo_c
   double VxWestSum = 0.0;
   double VxEastSum = 0.0;
 
-  if (instance->topo_height) {
+  if (instance->topo_height && topo_chain->Nb_part && topo_chain->z[0] && topo_chain->z[topo_chain->Nb_part - 1]) {
     instance->topo_height->west = topo_chain->z[0];
-    instance->topo_height->east = topo_chain->z[topo_chain->Nb_part-1];
+    instance->topo_height->east = topo_chain->z[topo_chain->Nb_part - 1];
   }
 
   for (int l = 0; l < mesh->Nz + 1; l++) {