the parallel extrusion was only reading idx files. It now reads also nc

and msh

dgpy/scripts/mesh.py

@@ -224,6 +224,9 @@ def read_msh(filename, vals={}) :
 from os import path
 def read_idx(filename) :
     vals = {}
+    if filename[-3:] == 'msh':
+        read_msh(filename, vals)
+        return vals
     reMerge = re.compile("\s*Merge\s* \"(.*)\"\s*;")
     dirname = path.dirname(filename)
     with open(filename, "r") as f:

dgpy/scripts/slimPre.py

@@ -841,6 +841,11 @@ def get_distance(region, physical_tags):
     return np.array(dist)
 
 def extrude(mesh_file_name, bath_file_name, nb_layers=None, z_layers=None, layers_function=None, mesh_file_name_out='', factor_show=0, periodicity=None):
+
+    if bath_file_name[0][-3:] == ".nc":
+        netcdf_to_msh(mesh_file_name, bath_file_name[0], bath_file_name[1], bath_file_name[0][:-3])
+        bath_file_name = bath_file_name[0][:-3]+'.msh'
+
     if dgpy.Msg.GetCommRank() == 0:
 
       if (nb_layers and z_layers) or (nb_layers and layers_function) or (z_layers and layers_function):