Salome HOME
Copyright update 2022
[modules/homard.git] / doc / files / tutorial_3.py
index dd0554888cd3158b415a3770f72296ec52c751e6..a94287ae4c2d5a41b0180c053a9751b560abdb36 100755 (executable)
@@ -1,7 +1,7 @@
-#!/usr/bin/env python
+#!/usr/bin/env python3
 # -*- coding: utf-8 -*-
 
-# Copyright (C) 2011-2015  CEA/DEN, EDF R&D
+# Copyright (C) 2011-2022  CEA/DEN, EDF R&D
 #
 # This library is free software; you can redistribute it and/or
 # modify it under the terms of the GNU Lesser General Public
 
 """
 Exemple de couplage HOMARD-Salome
-Copyright EDF-R&D 1996, 2010, 2014
+Copyright EDF 1996, 2010, 2019
 """
-__revision__ = "V1.9"
+__revision__ = "V3.03"
 #
 import os
 import sys
 #
 # ==================================
-pathHomard = os.getenv('HOMARD_ROOT_DIR')
+PATH_HOMARD = os.getenv("HOMARD_ROOT_DIR")
 # Repertoire des donnees du tutorial
-data_dir = os.path.join(pathHomard, "share", "doc", "salome", "gui", "HOMARD", "fr", "_downloads")
-data_dir = os.path.normpath(data_dir)
-sys.path.append(data_dir)
-from tutorial_util import gzip_gunzip
+DATA_TUTORIAL = os.path.join(PATH_HOMARD, "share", "doc", "salome", "gui", "HOMARD", "fr", "_downloads")
+DATA_TUTORIAL = os.path.normpath(DATA_TUTORIAL)
+sys.path.append(DATA_TUTORIAL)
 from tutorial_util import creation_dircase
 # ==================================
-dircase = creation_dircase(3)
-gzip_gunzip(data_dir, 3, -1)
+DIRCASE = creation_dircase(3)
 # ==================================
 #
 import salome
@@ -47,85 +45,84 @@ salome.salome_init()
 import HOMARD
 #
 homard = salome.lcc.FindOrLoadComponent("FactoryServer", "HOMARD")
-study_main = salome.myStudyManager.NewStudy("HOMARD")
-homard.SetCurrentStudy(salome.myStudy)
+homard.UpdateStudy()
 #
-# Hypothese "Hypo_0vers1"
+#============================= Début des commandes =============================
+#
+# Hypothese "hypo_0vers1"
 # =======================
-Hypo_0vers1 = homard.CreateHypothesis('Hypo_0vers1')
+hypo_0vers1 = homard.CreateHypothesis('hypo_0vers1')
 # Characterization of the field
-Hypo_0vers1.SetField('SOLU_0__QIRE_ELEM_SIGM__________')
-Hypo_0vers1.SetUseComp(0)
-Hypo_0vers1.AddComp('ERREST          ')
-Hypo_0vers1.SetRefinThr(3, 1.0)
-Hypo_0vers1.SetTypeFieldInterp(2)
-Hypo_0vers1.AddFieldInterp('SOLU_0__DEPL____________________')
-Hypo_0vers1.AddFieldInterp('SOLU_0__ERRE_ELEM_SIGM__________')
-#
-# Hypothese "Hypo_1vers2"
+hypo_0vers1.SetField('SOLU_0__QIRE_ELEM_SIGM__________')
+hypo_0vers1.SetUseComp(0)
+hypo_0vers1.AddComp('ERREST          ')
+hypo_0vers1.SetRefinThr(3, 1.0)
+hypo_0vers1.SetTypeFieldInterp(2)
+hypo_0vers1.AddFieldInterp('SOLU_0__DEPL____________________')
+hypo_0vers1.AddFieldInterp('SOLU_0__ERRE_ELEM_SIGM__________')
+#
+# Hypothese "hypo_1vers2"
 # =======================
-Hypo_1vers2 = homard.CreateHypothesis('Hypo_1vers2')
+hypo_1vers2 = homard.CreateHypothesis('hypo_1vers2')
 # Characterization of the field
-Hypo_1vers2.SetField('SOLU_1__QIRE_ELEM_SIGM__________')
-Hypo_1vers2.SetUseComp(0)
-Hypo_1vers2.AddComp('ERREST          ')
-Hypo_1vers2.SetRefinThr(3, 1.5)
-Hypo_1vers2.SetUnRefThr(3, 6.)
-Hypo_1vers2.SetTypeFieldInterp(2)
-Hypo_1vers2.AddFieldInterp('SOLU_1__DEPL____________________')
-Hypo_1vers2.AddFieldInterp('SOLU_1__QIRE_ELEM_SIGM__________')
-#
-# Hypothese "Hypo_1vers2_bis"
+hypo_1vers2.SetField('SOLU_1__QIRE_ELEM_SIGM__________')
+hypo_1vers2.SetUseComp(0)
+hypo_1vers2.AddComp('ERREST          ')
+hypo_1vers2.SetRefinThr(3, 1.5)
+hypo_1vers2.SetUnRefThr(3, 6.)
+hypo_1vers2.SetTypeFieldInterp(2)
+hypo_1vers2.AddFieldInterp('SOLU_1__DEPL____________________')
+hypo_1vers2.AddFieldInterp('SOLU_1__QIRE_ELEM_SIGM__________')
+#
+# Hypothese "hypo_1vers2_bis"
 # ===========================
-Hypo_1vers2_bis = homard.CreateHypothesis('Hypo_1vers2_bis')
+hypo_1vers2_bis = homard.CreateHypothesis('hypo_1vers2_bis')
 # Characterization of the field
-Hypo_1vers2_bis.SetField('SOLU_1__DEPL____________________')
-Hypo_1vers2_bis.SetUseComp(1)
-Hypo_1vers2_bis.AddComp('DX')
-Hypo_1vers2_bis.AddComp('DY')
-Hypo_1vers2_bis.AddComp('DZ')
-Hypo_1vers2_bis.SetRefinThr(1, 0.0001)
-Hypo_1vers2_bis.SetUnRefThr(1, 0.000001)
-Hypo_1vers2_bis.SetTypeFieldInterp(0)
+hypo_1vers2_bis.SetField('SOLU_1__DEPL____________________')
+hypo_1vers2_bis.SetUseComp(1)
+hypo_1vers2_bis.AddComp('DX')
+hypo_1vers2_bis.AddComp('DY')
+hypo_1vers2_bis.AddComp('DZ')
+hypo_1vers2_bis.SetRefinThr(1, 0.0001)
+hypo_1vers2_bis.SetUnRefThr(1, 0.000001)
+hypo_1vers2_bis.SetTypeFieldInterp(0)
 #
 # Cas
 # ===
-Case_3 = homard.CreateCase('Case_3', 'G_0', data_dir+'/tutorial_3.00.med')
-Case_3.SetDirName(dircase)
+le_cas = homard.CreateCase('Case_3', 'G_0', os.path.join(DATA_TUTORIAL, "tutorial_3.00.med"))
+le_cas.SetDirName(DIRCASE)
 #
-# Iteration "Iter_3_1"
+# Iteration "iter_3_1"
 # ====================
-Iter_3_1 = Case_3.NextIteration('Iter_3_1')
-Iter_3_1.SetMeshName('H_1')
-Iter_3_1.SetMeshFile(dircase+'/maill.01.med')
-Iter_3_1.SetFieldFile(data_dir+'/tutorial_3.00.med')
-Iter_3_1.SetTimeStepRank( 1, 1)
-Iter_3_1.AssociateHypo('Hypo_0vers1')
-error = Iter_3_1.Compute(1, 2)
-#
-# Iteration "Iter_3_2"
+iter_3_1 = le_cas.NextIteration('iter_3_1')
+iter_3_1.SetMeshName('H_1')
+iter_3_1.SetMeshFile(os.path.join(DIRCASE, "maill.01.med"))
+iter_3_1.SetFieldFile(os.path.join(DATA_TUTORIAL, "tutorial_3.00.med"))
+iter_3_1.SetTimeStepRank( 1, 1)
+iter_3_1.AssociateHypo('hypo_0vers1')
+erreur = iter_3_1.Compute(1, 2)
+#
+# Iteration "iter_3_2"
 # ====================
-Iter_3_2 = Iter_3_1.NextIteration('Iter_3_2')
-Iter_3_2.SetMeshName('H_2')
-Iter_3_2.SetMeshFile(dircase+'/maill.02.med')
-Iter_3_2.SetFieldFile(data_dir+'/tutorial_3.01.med')
-Iter_3_2.SetTimeStepRank(1, 1)
-Iter_3_2.AssociateHypo('Hypo_1vers2')
-error = Iter_3_2.Compute(1, 2)
-#
-# Iteration "Iter_3_2_bis"
+iter_3_2 = iter_3_1.NextIteration('iter_3_2')
+iter_3_2.SetMeshName('H_2')
+iter_3_2.SetMeshFile(os.path.join(DIRCASE, "maill.02.med"))
+iter_3_2.SetFieldFile(os.path.join(DATA_TUTORIAL, "tutorial_3.01.med"))
+iter_3_2.SetTimeStepRank(1, 1)
+iter_3_2.AssociateHypo('hypo_1vers2')
+erreur = iter_3_2.Compute(1, 2)
+#
+# Iteration "iter_3_2_bis"
 # ========================
-Iter_3_2_bis = Iter_3_1.NextIteration('Iter_3_2_bis')
-Iter_3_2_bis.SetMeshName('H_2_bis')
-Iter_3_2_bis.SetMeshFile(dircase+'/maill.02.bis.med')
-Iter_3_2_bis.SetFieldFile(data_dir+'/tutorial_3.01.med')
-Iter_3_2_bis.SetTimeStepRank(1, 1)
-Iter_3_2_bis.AssociateHypo('Hypo_1vers2_bis')
-error = Iter_3_2_bis.Compute(1, 2)
-
-# ==================================
-gzip_gunzip(data_dir, 3, 1)
-# ==================================
-
+iter_3_2_bis = iter_3_1.NextIteration('iter_3_2_bis')
+iter_3_2_bis.SetMeshName('H_2_bis')
+iter_3_2_bis.SetMeshFile(os.path.join(DIRCASE, "maill.02.bis.med"))
+iter_3_2_bis.SetFieldFile(os.path.join(DATA_TUTORIAL, "tutorial_3.01.med"))
+iter_3_2_bis.SetTimeStepRank(1, 1)
+iter_3_2_bis.AssociateHypo('hypo_1vers2_bis')
+erreur = iter_3_2_bis.Compute(1, 2)
+#
+#============================== Fin des commandes ==============================
+#
 if salome.sg.hasDesktop():
-  salome.sg.updateObjBrowser(1)
+  salome.sg.updateObjBrowser()