le maillage généré avec MedCoupling ne contenait pas les groupes

diff --git a/doc/salome/tutorial/interpolationZ.rst b/doc/salome/tutorial/interpolationZ.rst
index 0538697f8c3c6d574c920dc4939f03f8f78a60ab..ce2477663282c5e96ec6da2f01a27404d75021fd 100644 (file)
--- a/doc/salome/tutorial/interpolationZ.rst
+++ b/doc/salome/tutorial/interpolationZ.rst
@@ -92,8 +92,7 @@ avec des suffixes différents, rangés dans le répertoire du fichier d'origine
 * garonne_1.med  : fichier d'origine (coordonnée z = 0)
 * garonne_1.xyz  : fichier xyz (ASCII) des altitudes aux noeuds
 * garonne_1Z.med : coordonnée Z à sa valeur calculée 
-* garonne_1F.med : identique au précédent
-* garonne_1L.med : coordonnée Z à sa valeur calculée, et champ "BOTTOM" avec la valeur Z aux noeuds
+* garonne_1F.med : coordonnée Z à sa valeur calculée, et champ "BOTTOM" avec la valeur Z aux noeuds
 
 Il faut recopier le script et l'adapter en fonction des noms utilisés dans le cas de calcul 
 et pour le maillage.

diff --git a/src/HYDROTools/interpolZ.py b/src/HYDROTools/interpolZ.py
index efe490a266b6685de2c808cff673992c3d7bd74c..d803355bc49c057f5a8325e2bea57168e0e54213 100644 (file)
--- a/src/HYDROTools/interpolZ.py
+++ b/src/HYDROTools/interpolZ.py
@@ -109,7 +109,7 @@ def createZfield1(fichierMaillage):
 
 # -----------------------------------------------------------------------------
 
-from MEDLoader import MEDLoader, MEDCouplingFieldDouble, ON_NODES, DataArrayDouble
+from MEDLoader import MEDLoader, MEDCouplingFieldDouble, ON_NODES, DataArrayDouble, MEDFileMesh
 
 def createZfield2(fichierMaillage):
   """
@@ -117,24 +117,24 @@ def createZfield2(fichierMaillage):
   Add a field on nodes, named "BOTTOM", of type double, containing z coordinates of nodes.
   createZfield2 is used after interpolZ. createZfield1 is base on MEDLoader interface.
   There is an alternate method based on Med file, equivalent (createZfield1).
-  The file <fichierMaillage>F.med produced by interpolz must exist, and is modified.
   fichierMaillage : 2D (x,y) med file produced by SMESH and used by interpolZ.
-  return <fichierMaillage>L.med : med file containing the field "BOTTOM"
+  return <fichierMaillage>F.med : med file containing the field "BOTTOM"
   """
   
   noms = string.split(fichierMaillage,'.')
   basename = string.join(noms[:-1], '.')
   fichierZMaillage = basename + 'Z.med'
-  fichierLMaillage = basename + 'L.med'
-  print fichierLMaillage
+  fichierFMaillage = basename + 'F.med'
+  print fichierFMaillage
 
   mymesh = MEDLoader.ReadUMeshFromFile(fichierZMaillage,0)
   fieldOnNodes=MEDCouplingFieldDouble.New(ON_NODES)
   fieldOnNodes.setName("BOTTOM")
   fieldOnNodes.setMesh(mymesh)
   fieldOnNodes.setArray(mymesh.getCoords()[:,2])
-
-  MEDLoader.WriteField(fichierLMaillage,fieldOnNodes,True)
+  mm=MEDFileMesh.New(fichierZMaillage)
+  mm.write(fichierFMaillage,2)
+  MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fichierFMaillage,fieldOnNodes)
 
 # -----------------------------------------------------------------------------