From: mpa <mpa@opencascade.com>
Date: Tue, 19 Nov 2013 12:29:02 +0000 (+0000)
Subject: - Modifing Geometry and Mesh Python scripts from SALOME 6 and before
X-Git-Tag: V7_3_0a1~36
X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=98eacab61c81c951c8d7662ee279c6604842dc43;p=modules%2Fsmesh.git

- Modifing Geometry and Mesh Python scripts from SALOME 6 and before
---

diff --git a/src/Tools/Verima/Scripts/a.py b/src/Tools/Verima/Scripts/a.py
index ac79c4d54..56ff1da42 100644
--- a/src/Tools/Verima/Scripts/a.py
+++ b/src/Tools/Verima/Scripts/a.py
@@ -34,12 +34,11 @@ l_cas = [ "tout", "syme" ]
 ### GEOM component
 ###
 import GEOM
-import geompy
+from salome.geom import geomBuilder
+geompy = geomBuilder.New(theStudy)
 import math
 import SALOMEDS
 #
-geompy.init_geom(theStudy)
-#
 # 0. Les conventions d'orientation des axes sont les suivantes :
 #    . l'axe Z s'enfonce dans le plan du dessin, son origine étant dans le plan
 #    . l'axe Y est vers le haut
@@ -199,7 +198,7 @@ TRX = 0.0
 TRY = -2.0*LG_ARETE_BLOC + DELTA
 TRZ = - 2.0*DELTA
 geompy.TranslateDXDYDZ(boite_aux_G, TRX, TRY, TRZ)
-l_aux = geompy.GetShapesOnBox (boite_aux_G, B_G1, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+l_aux = geompy.GetShapesOnBox (boite_aux_G, B_G1, geompy.ShapeType["SOLID"], GEOM.ST_IN )
 #print "l_aux =", l_aux
 B_G2 = geompy.MakeCut(l_aux[0], OUTIL_G_1)
 # . Repérage du volume interne inférieur
@@ -207,7 +206,7 @@ TRX = 0.0
 TRY = -RAYON_H - DELTA
 TRZ = 0.0
 geompy.TranslateDXDYDZ(boite_aux_G, TRX, TRY, TRZ)
-l_aux = geompy.GetShapesOnBox (boite_aux_G, B_G1, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+l_aux = geompy.GetShapesOnBox (boite_aux_G, B_G1, geompy.ShapeType["SOLID"], GEOM.ST_IN )
 B_G3 = geompy.MakeCut(l_aux[0], OUTIL_G_2)
 #
 # 3.3.3. Partie droite
@@ -224,7 +223,7 @@ TRX = -DX
 TRY = -2.0*LG_ARETE_BLOC + DELTA
 TRZ = - 2.0*DELTA
 geompy.TranslateDXDYDZ(boite_aux_D, TRX, TRY, TRZ)
-l_aux = geompy.GetShapesOnBox (boite_aux_D, B_D1, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+l_aux = geompy.GetShapesOnBox (boite_aux_D, B_D1, geompy.ShapeType["SOLID"], GEOM.ST_IN )
 if ( verbose_max ) :
   print "3.3.3. supérieur l_aux =", l_aux
 B_D2 = geompy.MakeCut(l_aux[0], OUTIL_D_1)
@@ -233,7 +232,7 @@ TRX = 0.0
 TRY = -RAYON_H - DELTA
 TRZ = 0.0
 geompy.TranslateDXDYDZ(boite_aux_D, TRX, TRY, TRZ)
-l_aux = geompy.GetShapesOnBox (boite_aux_D, B_D1, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+l_aux = geompy.GetShapesOnBox (boite_aux_D, B_D1, geompy.ShapeType["SOLID"], GEOM.ST_IN )
 if ( verbose_max ) :
   print "3.3.3. inférieur l_aux =", l_aux
 B_D3 = geompy.MakeCut(l_aux[0], OUTIL_D_2)
@@ -393,7 +392,7 @@ for cas in l_cas :
   TRY = -0.5*LG_ARETE_MASSIF - DELTA
   TRZ = - DELTA
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_solid_id = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+  l_solid_id = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["SOLID"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.1. l_solid_id =", l_solid_id
   ROCHE_G = geompy.CreateGroup(MASSIF_G, geompy.ShapeType["SOLID"])
@@ -411,7 +410,7 @@ for cas in l_cas :
   TRY = -0.5*LG_ARETE_MASSIF - DELTA
   TRZ = - DELTA
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.1.1. Gauche l_aux =", l_aux
   GAUCHE_G = geompy.CreateGroup ( MASSIF_G, geompy.ShapeType["FACE"] )
@@ -426,7 +425,7 @@ for cas in l_cas :
     TRY = 0.0
     TRZ = 0.0
     geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-    l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+    l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
     if ( verbose_max ) :
       print "6.2.1.2. Droite l_aux =", l_aux
     DROITE_G = geompy.CreateGroup ( MASSIF_G, geompy.ShapeType["FACE"] )
@@ -438,7 +437,7 @@ for cas in l_cas :
     TRY = 0.0
     TRZ = 0.0
     geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-    l_face_id_syme  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+    l_face_id_syme  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
     if ( verbose_max ) :
       print "6.2.1.2. Symetrie - tout l_face_id_syme =", l_face_id_syme
     l_aux = []
@@ -463,7 +462,7 @@ for cas in l_cas :
       if ( verbose_max ) :
         print "6.2.1.2. Cavite %02d" % (iaux+1)
       # Les id des faces bordant la cavite courante
-      l_aux_2 = geompy.GetShapesOnBoxIDs (boite_2, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+      l_aux_2 = geompy.GetShapesOnBoxIDs (boite_2, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
       if ( verbose_max ) :
         print ".. l_aux_2 =", l_aux_2
       for face_id in l_aux_2 :
@@ -485,7 +484,7 @@ for cas in l_cas :
   TRY =  0.5*LG_ARETE_MASSIF - DELTA
   TRZ = - DELTA
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.2. Haut l_aux =", l_aux
   HAUT_G = geompy.CreateGroup ( MASSIF_G, geompy.ShapeType["FACE"] )
@@ -495,7 +494,7 @@ for cas in l_cas :
   TRY = -LG_ARETE_MASSIF
   TRZ = 0.0
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.2. Bas l_aux =", l_aux
   BAS_G = geompy.CreateGroup ( MASSIF_G, geompy.ShapeType["FACE"] )
@@ -512,7 +511,7 @@ for cas in l_cas :
   TRY = -0.5*LG_ARETE_BLOC - DELTA
   TRZ = - DELTA
   geompy.TranslateDXDYDZ(boite_aux_1, TRX, TRY, TRZ)
-  l_face_avant_tunnel  = geompy.GetShapesOnBoxIDs (boite_aux_1, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_face_avant_tunnel  = geompy.GetShapesOnBoxIDs (boite_aux_1, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.3.1. Devant l_face_avant_tunnel =", l_face_avant_tunnel
   # Création du groupe de la face avant
@@ -530,7 +529,7 @@ for cas in l_cas :
   TRY = -0.5*LG_ARETE_MASSIF - DELTA
   TRZ = - DELTA
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_aux_2  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_aux_2  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.3.2. Devant l_aux_2  =", l_aux_2
 # Pour le groupe, on retire les faces correspondant a l'empreinte du tunnel
@@ -549,7 +548,7 @@ for cas in l_cas :
   TRY = 0.0
   TRZ = LG_ARETE_MASSIF
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.3.3. Derriere l_aux =", l_aux
   DERRIERE_G = geompy.CreateGroup ( MASSIF_G, geompy.ShapeType["FACE"] )
@@ -567,7 +566,7 @@ for cas in l_cas :
   DY = -0.5*LG_ARETE_BLOC
   point = geompy.MakeVertex(DX, DY, DZ)
   #        Les faces posees sur ce plan
-  liste_face_tunnel_base = geompy.GetShapesOnPlaneWithLocationIDs (MASSIF_G, geompy.ShapeType["FACE"], normale, point, geompy.GEOM.ST_ON )
+  liste_face_tunnel_base = geompy.GetShapesOnPlaneWithLocationIDs (MASSIF_G, geompy.ShapeType["FACE"], normale, point, GEOM.ST_ON )
   if ( verbose_max ) :
     print "6.3.1. liste_face_tunnel_base =", liste_face_tunnel_base
 #        Création du groupe associe
@@ -595,7 +594,7 @@ for cas in l_cas :
     geompy.TranslateDXDYDZ(boite, TRX, TRY, TRZ)
     l_cyl_supp.append((boite, cle))
   # Reperage des faces
-    l_aux_1 = geompy.GetShapesOnShapeIDs (boite, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+    l_aux_1 = geompy.GetShapesOnShapeIDs (boite, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
     if ( verbose_max ) :
       print ".. l_aux_1 =", l_aux_1
     l_aux = []
@@ -648,7 +647,7 @@ for cas in l_cas :
       geompy.UnionIDs ( GR_CAV_G, l_cav_toutes_p_id )
       l_groupe_cav_e_g.append((GR_CAV_G, "PAROI_%02d" % (iaux+1)))
     # Les id des solides dans la cavite courante
-    l_aux_1 = geompy.GetShapesOnBoxIDs (boite_cav_v, MASSIF_G, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+    l_aux_1 = geompy.GetShapesOnBoxIDs (boite_cav_v, MASSIF_G, geompy.ShapeType["SOLID"], GEOM.ST_IN )
     if ( verbose_max ) :
       print ".. l_aux_1 =", l_aux_1
     # Création du groupe solide de la cavite courante
@@ -672,7 +671,7 @@ for cas in l_cas :
     geompy.UnionIDs ( GR_CAV_G, l_solid_id_1 )
     l_groupe_roche_g.append((GR_CAV_G, "ROCHE_%02d" % (iaux+1)))
     # Les id des faces du fond de la cavite courante
-    l_aux_1 = geompy.GetShapesOnBoxIDs (boite_cav_f, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+    l_aux_1 = geompy.GetShapesOnBoxIDs (boite_cav_f, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
     if ( verbose_max ) :
       print ".. l_aux_1 =", l_aux_1
     # Création du groupe du fond de la cavite
@@ -680,7 +679,7 @@ for cas in l_cas :
     geompy.UnionIDs ( GR_CAV_G, l_aux_1 )
     l_groupe_cav_f_g.append((GR_CAV_G, "FOND_%02d" % (iaux+1)))
     # Les id des faces dans la cavite courante
-    l_aux_2 = geompy.GetShapesOnBoxIDs (boite_cav_v, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+    l_aux_2 = geompy.GetShapesOnBoxIDs (boite_cav_v, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
     if ( verbose_max ) :
       print ".. l_aux_2 =", l_aux_2
     # Création du groupe du pourtour de la cavite : penser a retirer les fonds et l'eventuel plan de symetrie !
@@ -803,11 +802,11 @@ for cas in l_cas :
 ### SMESH component
 ###=======================================================================
 
-import smesh, SMESH, SALOMEDS
-
-smesh.SetCurrentStudy(theStudy)
-import BLSURFPlugin
-import GHS3DPlugin
+import SMESH, SALOMEDS
+from salome.smesh import smeshBuilder
+smesh = smeshBuilder.New(theStudy)
+from salome.BLSURFPlugin import BLSURFPluginBuilder
+from salome.GHS3DPlugin import GHS3DPluginBuilder
 #
 for cas in l_cas :
 #
@@ -823,19 +822,19 @@ for cas in l_cas :
 # 2. Parametres du maillage volumique
 #
   if cas == "tout" :
-    GHS3D_3D = MASSIF_M.Tetrahedron(algo=smesh.GHS3D)
+    GHS3D_3D = MASSIF_M.Tetrahedron(algo=smeshBuilder.GHS3D)
     GHS3D_Parameters = smesh.CreateHypothesis('GHS3D_Parameters', 'GHS3DEngine')
     # Niveau d'optimisation : 3 ==> standard +
     GHS3D_Parameters.SetOptimizationLevel( 3 )
   else :
     isdone = MASSIF_M.AddHypothesis(GHS3D_Parameters)
-    GHS3D_3D_1 = MASSIF_M.Tetrahedron(algo=smesh.GHS3D)
+    GHS3D_3D_1 = MASSIF_M.Tetrahedron(algo=smeshBuilder.GHS3D)
 #
 # 3. Parametres du maillage surfacique
 #
   if cas == "tout" :
     #
-    BLSURF = MASSIF_M.Triangle(algo=smesh.BLSURF)
+    BLSURF = MASSIF_M.Triangle(algo=smeshBuilder.BLSURF)
     BLSURF_Parameters = BLSURF.Parameters()
     # Geometrical mesh - if set to "Custom", allows user input in Angle Mesh S, Angle Mesh C and Gradation fields.
     # These fields control computation of the element size, so called geometrical size, conform to the
@@ -852,7 +851,7 @@ for cas in l_cas :
     BLSURF_Parameters.SetPhySize( 25. )
   else :
     isdone = MASSIF_M.AddHypothesis(BLSURF_Parameters)
-    BLSURF_1 = MASSIF_M.Triangle(algo=smesh.BLSURF)
+    BLSURF_1 = MASSIF_M.Triangle(algo=smeshBuilder.BLSURF)
 #
 # 4. Calcul
 #
diff --git a/src/Tools/Verima/Scripts/bb.py b/src/Tools/Verima/Scripts/bb.py
index 581fed896..43bd23180 100644
--- a/src/Tools/Verima/Scripts/bb.py
+++ b/src/Tools/Verima/Scripts/bb.py
@@ -1,11 +1,10 @@
 # Get Information About Mesh by GetMeshIn
 import salome
-import smesh
-import SMESH, SALOMEDS
-
 salome.salome_init()
-theStudy = salome.myStudy
-smesh.SetCurrentStudy(theStudy)
+
+import SMESH, SALOMEDS
+from salome.smesh import smeshBuilder
+smesh = smeshBuilder.New(salome.myStudy)
 
 aMesh,aResult = smesh.CreateMeshesFromMED('/tmp/MASSIF.new.med')
 print aResult, aMesh
diff --git a/src/Tools/Verima/Scripts/excavation_6.6.py b/src/Tools/Verima/Scripts/excavation_6.6.py
index 583b1a210..85fa6aa22 100644
--- a/src/Tools/Verima/Scripts/excavation_6.6.py
+++ b/src/Tools/Verima/Scripts/excavation_6.6.py
@@ -36,9 +36,6 @@ l_cas = [ "tout", "syme" ]
 import GEOM
 import math
 import SALOMEDS
-#
-#import geompy
-#geompy.init_geom(theStudy)
 from salome.geom import geomBuilder
 geompy = geomBuilder.New(theStudy) 
 #
@@ -201,7 +198,7 @@ TRX = 0.0
 TRY = -2.0*LG_ARETE_BLOC + DELTA
 TRZ = - 2.0*DELTA
 geompy.TranslateDXDYDZ(boite_aux_G, TRX, TRY, TRZ)
-l_aux = geompy.GetShapesOnBox (boite_aux_G, B_G1, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+l_aux = geompy.GetShapesOnBox (boite_aux_G, B_G1, geompy.ShapeType["SOLID"], GEOM.ST_IN )
 #print "l_aux =", l_aux
 B_G2 = geompy.MakeCut(l_aux[0], OUTIL_G_1)
 # . Repérage du volume interne inférieur
@@ -209,7 +206,7 @@ TRX = 0.0
 TRY = -RAYON_H - DELTA
 TRZ = 0.0
 geompy.TranslateDXDYDZ(boite_aux_G, TRX, TRY, TRZ)
-l_aux = geompy.GetShapesOnBox (boite_aux_G, B_G1, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+l_aux = geompy.GetShapesOnBox (boite_aux_G, B_G1, geompy.ShapeType["SOLID"], GEOM.ST_IN )
 B_G3 = geompy.MakeCut(l_aux[0], OUTIL_G_2)
 #
 # 3.3.3. Partie droite
@@ -226,7 +223,7 @@ TRX = -DX
 TRY = -2.0*LG_ARETE_BLOC + DELTA
 TRZ = - 2.0*DELTA
 geompy.TranslateDXDYDZ(boite_aux_D, TRX, TRY, TRZ)
-l_aux = geompy.GetShapesOnBox (boite_aux_D, B_D1, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+l_aux = geompy.GetShapesOnBox (boite_aux_D, B_D1, geompy.ShapeType["SOLID"], GEOM.ST_IN )
 if ( verbose_max ) :
   print "3.3.3. supérieur l_aux =", l_aux
 B_D2 = geompy.MakeCut(l_aux[0], OUTIL_D_1)
@@ -235,7 +232,7 @@ TRX = 0.0
 TRY = -RAYON_H - DELTA
 TRZ = 0.0
 geompy.TranslateDXDYDZ(boite_aux_D, TRX, TRY, TRZ)
-l_aux = geompy.GetShapesOnBox (boite_aux_D, B_D1, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+l_aux = geompy.GetShapesOnBox (boite_aux_D, B_D1, geompy.ShapeType["SOLID"], GEOM.ST_IN )
 if ( verbose_max ) :
   print "3.3.3. inférieur l_aux =", l_aux
 B_D3 = geompy.MakeCut(l_aux[0], OUTIL_D_2)
@@ -395,7 +392,7 @@ for cas in l_cas :
   TRY = -0.5*LG_ARETE_MASSIF - DELTA
   TRZ = - DELTA
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_solid_id = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+  l_solid_id = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["SOLID"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.1. l_solid_id =", l_solid_id
   ROCHE_G = geompy.CreateGroup(MASSIF_G, geompy.ShapeType["SOLID"])
@@ -413,7 +410,7 @@ for cas in l_cas :
   TRY = -0.5*LG_ARETE_MASSIF - DELTA
   TRZ = - DELTA
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.1.1. Gauche l_aux =", l_aux
   GAUCHE_G = geompy.CreateGroup ( MASSIF_G, geompy.ShapeType["FACE"] )
@@ -428,7 +425,7 @@ for cas in l_cas :
     TRY = 0.0
     TRZ = 0.0
     geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-    l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+    l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
     if ( verbose_max ) :
       print "6.2.1.2. Droite l_aux =", l_aux
     DROITE_G = geompy.CreateGroup ( MASSIF_G, geompy.ShapeType["FACE"] )
@@ -440,7 +437,7 @@ for cas in l_cas :
     TRY = 0.0
     TRZ = 0.0
     geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-    l_face_id_syme  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+    l_face_id_syme  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
     if ( verbose_max ) :
       print "6.2.1.2. Symetrie - tout l_face_id_syme =", l_face_id_syme
     l_aux = []
@@ -465,7 +462,7 @@ for cas in l_cas :
       if ( verbose_max ) :
         print "6.2.1.2. Cavite %02d" % (iaux+1)
       # Les id des faces bordant la cavite courante
-      l_aux_2 = geompy.GetShapesOnBoxIDs (boite_2, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+      l_aux_2 = geompy.GetShapesOnBoxIDs (boite_2, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
       if ( verbose_max ) :
         print ".. l_aux_2 =", l_aux_2
       for face_id in l_aux_2 :
@@ -487,7 +484,7 @@ for cas in l_cas :
   TRY =  0.5*LG_ARETE_MASSIF - DELTA
   TRZ = - DELTA
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.2. Haut l_aux =", l_aux
   HAUT_G = geompy.CreateGroup ( MASSIF_G, geompy.ShapeType["FACE"] )
@@ -497,7 +494,7 @@ for cas in l_cas :
   TRY = -LG_ARETE_MASSIF
   TRZ = 0.0
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.2. Bas l_aux =", l_aux
   BAS_G = geompy.CreateGroup ( MASSIF_G, geompy.ShapeType["FACE"] )
@@ -514,7 +511,7 @@ for cas in l_cas :
   TRY = -0.5*LG_ARETE_BLOC - DELTA
   TRZ = - DELTA
   geompy.TranslateDXDYDZ(boite_aux_1, TRX, TRY, TRZ)
-  l_face_avant_tunnel  = geompy.GetShapesOnBoxIDs (boite_aux_1, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_face_avant_tunnel  = geompy.GetShapesOnBoxIDs (boite_aux_1, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.3.1. Devant l_face_avant_tunnel =", l_face_avant_tunnel
   # Création du groupe de la face avant
@@ -532,7 +529,7 @@ for cas in l_cas :
   TRY = -0.5*LG_ARETE_MASSIF - DELTA
   TRZ = - DELTA
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_aux_2  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_aux_2  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.3.2. Devant l_aux_2  =", l_aux_2
 # Pour le groupe, on retire les faces correspondant a l'empreinte du tunnel
@@ -551,7 +548,7 @@ for cas in l_cas :
   TRY = 0.0
   TRZ = LG_ARETE_MASSIF
   geompy.TranslateDXDYDZ(boite_aux, TRX, TRY, TRZ)
-  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+  l_aux  = geompy.GetShapesOnBoxIDs (boite_aux, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
   if ( verbose_max ) :
     print "6.2.3.3. Derriere l_aux =", l_aux
   DERRIERE_G = geompy.CreateGroup ( MASSIF_G, geompy.ShapeType["FACE"] )
@@ -569,7 +566,7 @@ for cas in l_cas :
   DY = -0.5*LG_ARETE_BLOC
   point = geompy.MakeVertex(DX, DY, DZ)
   #        Les faces posees sur ce plan
-  liste_face_tunnel_base = geompy.GetShapesOnPlaneWithLocationIDs (MASSIF_G, geompy.ShapeType["FACE"], normale, point, geompy.GEOM.ST_ON )
+  liste_face_tunnel_base = geompy.GetShapesOnPlaneWithLocationIDs (MASSIF_G, geompy.ShapeType["FACE"], normale, point, GEOM.ST_ON )
   if ( verbose_max ) :
     print "6.3.1. liste_face_tunnel_base =", liste_face_tunnel_base
 #        Création du groupe associe
@@ -597,7 +594,7 @@ for cas in l_cas :
     geompy.TranslateDXDYDZ(boite, TRX, TRY, TRZ)
     l_cyl_supp.append((boite, cle))
   # Reperage des faces
-    l_aux_1 = geompy.GetShapesOnShapeIDs (boite, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+    l_aux_1 = geompy.GetShapesOnShapeIDs (boite, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
     if ( verbose_max ) :
       print ".. l_aux_1 =", l_aux_1
     l_aux = []
@@ -650,7 +647,7 @@ for cas in l_cas :
       geompy.UnionIDs ( GR_CAV_G, l_cav_toutes_p_id )
       l_groupe_cav_e_g.append((GR_CAV_G, "PAROI_%02d" % (iaux+1)))
     # Les id des solides dans la cavite courante
-    l_aux_1 = geompy.GetShapesOnBoxIDs (boite_cav_v, MASSIF_G, geompy.ShapeType["SOLID"], geompy.GEOM.ST_IN )
+    l_aux_1 = geompy.GetShapesOnBoxIDs (boite_cav_v, MASSIF_G, geompy.ShapeType["SOLID"], GEOM.ST_IN )
     if ( verbose_max ) :
       print ".. l_aux_1 =", l_aux_1
     # Création du groupe solide de la cavite courante
@@ -674,7 +671,7 @@ for cas in l_cas :
     geompy.UnionIDs ( GR_CAV_G, l_solid_id_1 )
     l_groupe_roche_g.append((GR_CAV_G, "ROCHE_%02d" % (iaux+1)))
     # Les id des faces du fond de la cavite courante
-    l_aux_1 = geompy.GetShapesOnBoxIDs (boite_cav_f, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+    l_aux_1 = geompy.GetShapesOnBoxIDs (boite_cav_f, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
     if ( verbose_max ) :
       print ".. l_aux_1 =", l_aux_1
     # Création du groupe du fond de la cavite
@@ -682,7 +679,7 @@ for cas in l_cas :
     geompy.UnionIDs ( GR_CAV_G, l_aux_1 )
     l_groupe_cav_f_g.append((GR_CAV_G, "FOND_%02d" % (iaux+1)))
     # Les id des faces dans la cavite courante
-    l_aux_2 = geompy.GetShapesOnBoxIDs (boite_cav_v, MASSIF_G, geompy.ShapeType["FACE"], geompy.GEOM.ST_IN )
+    l_aux_2 = geompy.GetShapesOnBoxIDs (boite_cav_v, MASSIF_G, geompy.ShapeType["FACE"], GEOM.ST_IN )
     if ( verbose_max ) :
       print ".. l_aux_2 =", l_aux_2
     # Création du groupe du pourtour de la cavite : penser a retirer les fonds et l'eventuel plan de symetrie !
@@ -805,11 +802,12 @@ for cas in l_cas :
 ### SMESH component
 ###=======================================================================
 
-import smesh, SMESH, SALOMEDS
+import SMESH, SALOMEDS
+from salome.smesh import smeshBuilder
+smesh = smeshBuilder.New(theStudy)
 
-smesh.SetCurrentStudy(theStudy)
-import BLSURFPlugin
-import GHS3DPlugin
+from salome.BLSURFPlugin import BLSURFPluginBuilder
+from salome.GHS3DPlugin import GHS3DPluginBuilder
 #
 for cas in l_cas :
 #
@@ -825,19 +823,19 @@ for cas in l_cas :
 # 2. Parametres du maillage volumique
 #
   if cas == "tout" :
-    GHS3D_3D = MASSIF_M.Tetrahedron(algo=smesh.GHS3D)
+    GHS3D_3D = MASSIF_M.Tetrahedron(algo=smeshBuilder.GHS3D)
     GHS3D_Parameters = smesh.CreateHypothesis('GHS3D_Parameters', 'GHS3DEngine')
     # Niveau d'optimisation : 3 ==> standard +
     GHS3D_Parameters.SetOptimizationLevel( 3 )
   else :
     isdone = MASSIF_M.AddHypothesis(GHS3D_Parameters)
-    GHS3D_3D_1 = MASSIF_M.Tetrahedron(algo=smesh.GHS3D)
+    GHS3D_3D_1 = MASSIF_M.Tetrahedron(algo=smeshBuilder.GHS3D)
 #
 # 3. Parametres du maillage surfacique
 #
   if cas == "tout" :
     #
-    BLSURF = MASSIF_M.Triangle(algo=smesh.BLSURF)
+    BLSURF = MASSIF_M.Triangle(algo=smeshBuilder.BLSURF)
     BLSURF_Parameters = BLSURF.Parameters()
     # Geometrical mesh - if set to "Custom", allows user input in Angle Mesh S, Angle Mesh C and Gradation fields.
     # These fields control computation of the element size, so called geometrical size, conform to the
@@ -854,7 +852,7 @@ for cas in l_cas :
     BLSURF_Parameters.SetPhySize( 25. )
   else :
     isdone = MASSIF_M.AddHypothesis(BLSURF_Parameters)
-    BLSURF_1 = MASSIF_M.Triangle(algo=smesh.BLSURF)
+    BLSURF_1 = MASSIF_M.Triangle(algo=smeshBuilder.BLSURF)
 #
 # 4. Calcul
 #
diff --git a/src/Tools/Verima/Scripts/toto.py b/src/Tools/Verima/Scripts/toto.py
index 94137d754..52358fd3a 100644
--- a/src/Tools/Verima/Scripts/toto.py
+++ b/src/Tools/Verima/Scripts/toto.py
@@ -1,5 +1,11 @@
-import geompy
-import smesh
+import GEOM
+from salome.geom import geomBuilder
+geompy = geomBuilder.New(salome.myStudy)
+
+import SMESH
+from salome.smesh import smeshBuilder
+smesh = smeshBuilder.New(salome.myStudy)
+
 import sys
 
 # create a box