preliminaires.rst
sinusX_Format.rst
streams.rst
+ completeCase.rst
+ copyMesh.rst
)
SET(SPHINXOPTS )
SALOME_CONFIGURE_FILE(conf.py.in conf.py)
ADD_CUSTOM_TARGET(BUILD_HTML ALL
- COMMAND ${CMAKE_COMMAND} -E make_directory html
+ COMMAND ${CMAKE_COMMAND} -E make_directory html
COMMAND ${CMAKE_COMMAND} -E make_directory doctrees
COMMAND ${SPHINX_EXECUTABLE} -c ${CMAKE_BINARY_DIR}/doc/salome/tutorial -b html ${ALLSPHINXOPTS} html
DEPENDS ${RSTFILES}
WORKING_DIRECTORY ${CMAKE_CURRENT_BINARY_DIR}
)
-
+
ADD_CUSTOM_TARGET(BUILD_PDF ALL
- COMMAND ${CMAKE_COMMAND} -E make_directory latex
+ COMMAND ${CMAKE_COMMAND} -E make_directory latex
COMMAND ${SPHINX_EXECUTABLE} -c ${CMAKE_BINARY_DIR}/doc/salome/tutorial -b latex ${ALLSPHINXOPTS} latex
COMMAND make -C latex all-pdf
DEPENDS ${RSTFILES}
WORKING_DIRECTORY ${CMAKE_CURRENT_BINARY_DIR}
)
-INSTALL(DIRECTORY ${CMAKE_CURRENT_BINARY_DIR}/html/
+INSTALL(DIRECTORY ${CMAKE_CURRENT_BINARY_DIR}/html/
DESTINATION ${SALOME_INSTALL_DOC}/gui/HYDRO_tutorial
PATTERN ".buildinfo" EXCLUDE
)
|createCalculationCase2|
-Le dialogue propose de sélectionner des groupes à conserver :
+Le dialogue propose de sélectionner des groupes à conserver :
il s'agit des contours des différents objets pris en compte dans le cas.
Ces groupes peuvent se révéler utiles par la suite pour définir des conditions aux limites, ou contrôler le maillage.
Le fait de les conserver permet de les retrouver sous leur nom dans les différentes étapes du calcul.
Le dialogue suivant permet d'affecter des coefficients de frottement à différentes zones du domaine,
à l'aide d'une table de coefficients de Strickler, et d'un découpage en zones `land cover map`.
Il faut avoir préalablement importé et/ou édité cette table de Strickler et ce `land cover map`.
-Pour cet exemple, la `land cover map` est optionnelle. Pour l'introduire, on peut se référer au
+Pour cet exemple, la `land cover map` est optionnelle. Pour l'introduire, on peut se référer au
paragraphe *Land Cover Map* plus bas.
Il faut valider cette étape en appuyant sur le bouton *next*.
Quand la résolution des conflits est terminée, il ne doit plus rester de zone rouge.
Il est possible de régénérer les couleurs des zones (bouton *regenerate colors*) si nécessaire.
-
+
|createCalculationCase7|
Regroupement des zones en régions
Le cas est publié dans l'arbre d'étude.Nous sélectionnons chaque région pour lui donner un nom significatif (menu contextuel *Rename*).
Pour identifier la région, il faut s'aider en examinant les zones qu'elle contient.
-La vue graphique affiche les zones. En développant l'arbre d'étude au niveau de chaque région,
+La vue graphique affiche les zones. En développant l'arbre d'étude au niveau de chaque région,
il est possible de sélectionner les zones pour les mettre en surbrillance.
En principe, il n'est pas nécessaire de renommer les zones elles mêmes, mais seulement les régions.
Dans le module *HYDROSOLVER*, on utilise le dialogue *change Coordinates* accessible depuis le menu *Hydro*
ou l'icone |changes_coords|.
-Si on laisse décochée la case *Change coordinates system*, Le dialogue permet une simple translation du maillage,
+Si on laisse décochée la case *Change coordinates system*, Le dialogue permet une simple translation du maillage,
**sans perte des champs que le fichier peut contenir**. Une fois sélectionné le fichier du maillage d'origine,
le dialogue propose un nom de fichier pour le maillage translaté en ajoutant un suffixe. Il est possible de prendre le nom de son choix.
Il faut indiquer les coordonnées du repère local d'origine, et celles d'arrivée.
--- /dev/null
+..
+ Copyright (C) 2015-2016 EDF
+
+ This file is part of SALOME HYDRO module.
+
+ SALOME HYDRO module is free software: you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ SALOME HYDRO module is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with SALOME HYDRO module. If not, see <http://www.gnu.org/licenses/>.
+
+#########################################
+Compléter un cas de calcul
+#########################################
+
+.. |CompleteCase| image:: /_static/CompleteCase.png
+ :align: middle
+
+.. |CompleteCase1| image:: /_static/CompleteCase1.png
+ :align: middle
+
+Une fois que le maillage a été réalisé et la topographie/bathymétrie interpolée, il peut être nécessaire de rajouter des polylignes dans la géométrie (lignes de crête de digue, par exemple). Ces polylignes sont des lignes de contraintes que l’on souhaite intégrer au maillage par la suite.
+
+Il faut, dans un premier temps, rajouter sa polyline dans la rubrique POLYLINES.
+Puis sélectionner le cas de calcul à modifier et faire clic droit / Complete calculation case.
+
+ |CompleteCase|
+
+Rajouter alors la polyline dans Included objects et cliquer sur Next jusqu’à la dernière page.
+
+ |CompleteCase1|
+
+Cette option *Complete calculation case* permet de ne pas devoir récréer le cas de calcul de zéro et de ne pas devoir réassembler les zones qui auraient été préalablement déjà regroupées. Seules les zones coupées par la polyline rajoutée devront être regroupées selon le souhait de l'utilisateur.
+Cliquer sur Finish et Exporter vers GEOM. On pourra alors appliquer un sous-maillage sur cette polyline.
+
+.. only:: html
+
+ :ref:`ref_casParticuliers`
+
--- /dev/null
+..
+ Copyright (C) 2015-2016 EDF
+
+ This file is part of SALOME HYDRO module.
+
+ SALOME HYDRO module is free software: you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ SALOME HYDRO module is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with SALOME HYDRO module. If not, see <http://www.gnu.org/licenses/>.
+
+#####################################################
+Copier un maillage existant sur une autre géométrie
+#####################################################
+
+.. |CopyMesh| image:: /_static/CopyMesh.png
+ :align: middle
+
+.. |CopyMesh1| image:: /_static/CopyMesh1.png
+ :align: middle
+
+.. |CopyMesh2| image:: /_static/CopyMesh2.png
+ :align: middle
+
+.. |CopyMesh3| image:: /_static/CopyMesh3.png
+ :align: middle
+
+.. |CopyMesh4| image:: /_static/CopyMesh4.png
+ :align: middle
+
+.. |CopyMesh5| image:: /_static/CopyMesh5.png
+ :align: middle
+
+.. |CopyMesh6| image:: /_static/CopyMesh6.png
+ :align: middle
+
+.. |CopyMesh7| image:: /_static/CopyMesh7.png
+ :align: middle
+
+.. |CopyMesh8| image:: /_static/CopyMesh8.png
+ :align: middle
+
+On dispose d’une géométrie A et d’un maillage A s’appuyant sur la géométrie A. On souhaite réaliser un maillage B de la géométrie B avec les mêmes hypothèses que le maillage A, sachant que la géométrie B diffère très peu de la A.
+
+Pour avoir la géométrie B, on a utilisé l’option *Complete calculation case* du module *HYDRO*. Une fois le cas de calcul complété, exporté dans GEOM, penser à créer le groupe d’edge correspondant à la polyline dans le module GEOM. La polyline apparaîtra alors dans le maillage dans SMESH et l’on pourra réaliser un sous-maillage s’appuyant dessus.
+
+ |CopyMesh|
+
+Dans le module SMESH, aller dans *Mesh / Copy Mesh*.
+
+ |CopyMesh1|
+
+Renseigner alors l’interface suivante :
+
+ |CopyMesh2|
+
+Dans Mesh : sélectionner le maillage à copier (maillage A).
+Dans New Mesh Name, donner un nom au futur maillage (maillage B).
+Dans New Geometry, sélectionner la nouvelle géométrie sur laquelle on souhaite réaliser un nouveau maillage (géométrie B).
+Cliquer sur *Apply and Close*.
+
+Voici le résultat que l’on obtient. La partie comprenant la nouvelle polyline, dans notre cas, repose sur une partie de la géométrie qui a été modifiée : elle n’est donc pas maillée immédiatement.
+
+ |CopyMesh3|
+
+Si on souhaite appliquer un sous-maillage à la polyline que l’on a rajoutée, cliquer sur *Create sub-mesh* et renseigner les différents champs (exemple possible de choix sur la figure ci-dessous).
+
+ |CopyMesh4|
+
+On obtient alors le maillage suivant :
+
+ |CopyMesh5|
+
+On remarque que la face qui vient d’être maillée n’est pas orientée dans le même sens que les autres faces (couleur foncée au lieu de claire).
+
+Pour réorienter la face, aller dans *Modification / Orientation*. Dans Group, sélectionner la face à réorienter en cliquant directement dans l’Object Browser sur le groupe de face concerné. Cliquer sur Add, la face devient jaune.
+Cliquer sur *Apply and Close*.
+
+ |CopyMesh6|
+
+ |CopyMesh7|
+
+Voici le maillage final :
+
+ |CopyMesh8|
+
+.. only:: html
+
+ :ref:`ref_casParticuliers`
+
.. |paravis08| image:: /_static/paravis08.png
:align: middle
-.. |paravis08| image:: /_static/paravis08.png
- :align: middle
-
.. |paravis09| image:: /_static/paravis09.png
:align: middle
.. |paravis27| image:: /_static/paravis27.png
:align: middle
+.. |paravis28| image:: /_static/paravis28.png
+ :align: middle
+
+.. |paravis29| image:: /_static/paravis29.png
+ :align: middle
+
+.. |paravis30| image:: /_static/paravis30.png
+ :align: middle
+
+.. |paravis31| image:: /_static/paravis31.png
+ :align: middle
+
+.. |paravis32| image:: /_static/paravis32.png
+ :align: middle
+
+.. |paravis33| image:: /_static/paravis33.png
+ :align: middle
+
+.. |paravis34| image:: /_static/paravis34.png
+ :align: middle
+
+.. |paravis35| image:: /_static/paravis35.png
+ :align: middle
+
+.. |paravis36| image:: /_static/paravis36.png
+ :align: middle
+
Paravis correspond à l'adaptation du produit open source Paraview à SALOME.
-Paravis peut être utilisé de de façons. Soit comme module depuis l'interface graphique de SALOME,
+Paravis peut être utilisé de plusieur façons. Soit comme module depuis l'interface graphique de SALOME,
|paravisDansSalome|
ftp://ftp.pleiade.edf.fr/projets/salome/Releases/
-Une fois SALOME installé, 2 dossiers sont créés :
+Une fois SALOME installé, 2 dossiers sont créés :
* appli_xxx
* Salome_xxx
-Lancer Paraview par le salome shell de SALOME.
+Se placer dans le dossier appli_xxx et Lancer Paraview par le salome shell de SALOME.
La commande est : ::
<chemin appli_xxx>/salome shell paraview
-
+
Exemple ::
~/salome-hydro/SALOME_DEV/14-10-2018/appli_DEV/salome shell paraview
Aller dans *Filters / Hydraulics*. C’est là que sont répertoriés les filtres dédiés à l’hydraulique tels que :
-* Profil spatial
-* Profil temporel
-* Calcul de débit à travers une section
-* Calcul de sédiments érodés / déposés dans une zone
-* Evolution temporelle d’une variable dans une colonne d’eau (en un x et y fixé,
+* Profil spatial `Filtre Spatial Profile`_
+* :ref:Profil temporel `Filtre Temporal on Point`_
+* :ref:Calcul de débit à travers une section `Filtre Rate Of Flow Through Section`_
+* :ref:Calcul de sédiments érodés / déposés dans une zone `Filtre Sediment Deposit`_
+* :ref:Evolution temporelle d’une variable dans une colonne d’eau `Filtre Depth Vs Time On Point`_ (en un x et y fixé,
évolution d’une variable comme la température, par exemple, en fonction du temps
en abscisse et de la profondeur en ordonnées)
-
+
+
|hydraulicsFilters|
-
+
Filtre Temporal on Point
------------------------
Pour l’instant, se placer en 3D. Il existe un bug pour déplacer la croix si on reste en 2D.
- |paravis02|
-
Déplacer la croix blanche à l’endroit souhaité. La croix blanche est située initialement au centre de la boîte encadrant le modèle.
+ |paravis02|
+
+Jouer sur le Radius (encadré en rouge) pour diminuer le rayon de la sphère en fonction de la taille des mailles.
+
|paravis03|
Cliquer sur *Apply*.
|paravis05|
-Déplacer les 2 extrémités de cette ligne à l’endroit où l’on souhaite commencer la polyligne.
+Déplacer les 2 extrémités de cette ligne à l’endroit où l’on souhaite commencer la polyligne :
+pour cela, placer la souris où l’on souhaite déplacer l’extrémité sur le fichier résultat et taper 1 ou Ctrl+1. Faire de même pour la deuxième extrémité en faisant 2 ou Ctrl+2.
|paravis06|
Adapter l’échelle de couleur si besoin : pour cela se placer sur le fichier résultat et cliquer sur :
|paravis07a|
-
+
|paravis07|
Se replacer sur le filtre SpatialProfile1.
|paravis08|
-Utiliser les commandes de l’encadré rouge :
+Utiliser les commandes de l’encadré rouge :
on conseille d’utiliser la dernière commande *Alt + Left click* pour ajouter des points à la fin de la polyligne.
|paravis09|
Cliquer sur *Apply*.
-Sélectionner la variable souhaitée, par exemple *Curv Abscissa*.
+Sélectionner la variable souhaitée, par exemple *Free Surface*.
+Dans **X Axis Parameters**, décocher *Use Index For Axis* et sélectionner *Curv Abscissa* dans *X Array Name*.
|paravis11|
* **Astuce :** La polyligne a disparu à ce moment sur la fenêtre de gauche RenderView.
- Pour la faire réapparaître, se placer sur la fenêtre RenderView montrant le fichier résultat.
+ Pour la faire réapparaître, se placer sur la fenêtre RenderView montrant le fichier résultat et cliquer sur *Show Spline*.
|paravis12|
-* **Remarque importante :**
+Penser à **sélectionner l’option Subdivide Input Polyline**. Elle permet d’avoir une interpolation plus ou moins précise entre les points que l’on a cliqués à l’écran. C’est à l’utilisateur de bouger le curseur de façon à ce que la courbe ne change plus. 50 suffit la plupart du temps mais il faut vérifier qu’il n’y a plus de changement pour des valeurs supérieures (dans certains cas c’est 400).
+
+ |paravis28|
+
+* **Remarque :**
+
+Un exemple de comparaison pour un profil en travers d’un cours d’eau si on ne sélectionne pas Subdivide Input Polyline (à gauche) et si on le sélectionne et le met égal à 50 (à droite) :
- Au 16-11-2018, il existe un bug sur l’interpolation le long de cette polyligne.
- Le tracé de la polyligne fonctionne bien mais les points rajoutés ne sont pas pris en compte dans l’interpolation.
+ |paravis29|
-Filtre Spatial Profile On Points
+ |paravis30|
+
+ |paravis31|
+
+Poly Line Source
--------------------------------
-Charger un fichier résultat.
+Les trois filtres qui suivent peuvent être appliqués sur des couches shape, sinusx mais aussi sur des polylines qui seraient dessinées à main levée sur le fichier résultat.
+
+Pour cela, aller dans Sources et sélectionner Poly Line Source.
+
+ |paravis32|
+
+Voici ce qu’on obtient une fois la source Poly Line Source sélectionnée :
+
+ |paravis33|
+
+ |paravis34|
+
+Par défaut la polyline est placée en 0, 0. Ce qui, dans la plupart des cas, sera très loin de notre fichier résultat. Pour ramener la polyline facilement sur le fichier résultat, utiliser les 2 commandes suivantes :
+
+* Déplacer les 2 extrémités de cette ligne à l’endroit où l’on souhaite commencer la polyligne : pour cela, placer la souris où l’on souhaite déplacer l’extrémité de la ligne sur le fichier résultat et taper 1 ou Ctrl+1.
+* Faire de même pour la deuxième extrémité en faisant 2 ou Ctrl+2.
+
+Puis, pour dessiner la polyline :
+
+* Sélectionner un des deux points extrémité de la ligne
+* Alt + Left Clic permet de rajouter des points à la suite de l’extrémité sélectionnée préalablement.
+* Ctrl + Left click permet d’insérer un nouveau point sur la ligne au niveau de la position de la souris.
+* Shift + Left Click permet d’enlever un point au niveau de la position de la souris.
+
+ |paravis35|
Cliquer sur *Apply*.
-Charger une polyligne déjà créée. Cette polyligne peut avoir 2 formats : shape ou sinusx.
-(créée dans Qgis par exemple au format shape).
+On peut ensuite lui appliquer les filtres Spatial Profile With Source ou Sediment Deposit ou Rate Of Flow Through Section. Pour Sediment Deposit, penser à fermer la polyline en cochant Closed (entouré en rouge sur la Figure ci-dessus).
+
+
+Filtre Spatial Profile With Source
+----------------------------------
+
+Charger un fichier résultat.
+
+Cliquer sur *Apply*.
+
+Charger une polyligne déjà créée. Cette polyligne peut avoir 2 formats : shape ou sinusx (créée dans Qgis, par exemple, au format shape).
Cliquer sur *Apply*.
Penser aux possibles translations qui ont pu être effectuées lors de la construction du modèle.
Si la polyligne est construite dans Qgis sur fond de carte, elle sera dans le référentiel global.
-Il va donc falloir la translater dans Paraview.
+Il va donc falloir la translater dans Paraview.
Pour cela une fois la polyligne chargée, sélectionner le filtre *Transform* et renseigner
la translation dans l’encadré rouge.
|paravis14|
-Lancer le filtre Spatial *Profil On Points*.
+Lancer le filtre *Spatial Profil With Source*.
-Une fenêtre s’ouvre :
+Une fenêtre s’ouvre :
-* Renseigner dans un premier temps l’input : il s’agit du fichier résultat.
-* Renseigner dans un deuxième temps la source : il s’agit de la polyligne chargée (ou de sa translation s’il y a lieu).
+* Renseigner dans un premier temps l’input : il s’agit du fichier résultat.
+* Renseigner dans un deuxième temps la source : il s’agit de la polyligne (chargée (ou de sa translation s’il y a lieu) ou dessinée).
|paravis15| |paravis16|
|paravis17|
-* **Remarque :** Un point important est de choisir le nombre de subdivisions de la polyligne (chiffre variant entre 1 et 100).
- La valeur conseillée est 50. La courbe traçant la variable sera alors plus précise, plus détaillée.
-
- Pour cela, sélectionner la polyligne dans *Pipeline Browser*.
-
- Changer le nombre de subdivisions : ici on passe de 1 à 50.
+* **Remarque :** Comme pour le filtre Spatial Profile, il est important de cocher Subdivide Input Polyline (qui apparaît une fois le filtre activé et sélectionné) afin d’avoir une interpolation précise. La valeur conseillée est 50 (mais ne pas hésiter à jouer sur la valeur pour vérifier que la courbe ne bouge plus). La courbe traçant la variable sera alors plus précise, plus détaillée.
+
+Pour cela, sélectionner la polyligne dans *Pipeline Browser* et cocher *Subdivide Input Polyline*.
|paravis18|
L’interpolation est plus précise.
+On peut également charger des couches shape comprenant plusieurs polylines, ce qui permettra de comparer des profils en travers par exemple :
+
+ |paravis36|
+
+* **Astuce :** Si vous avez une polyline au format SinusX et qu’elle est chargée sans message d’erreur dans PARAVIEW mais que vous ne la voyez pas, pensez à utiliser la commande suivante : ::
+
+ sed -i 's/\r//g' nom_du_fichier
+
+Elle permet de Transformer un fichier au format SinusX créé sous Windows en format SinusX lisible sous Linux.
+
+
Filtre Rate Of Flow Through Section
-----------------------------------
Cliquer sur *Apply*.
Charger un fichier polyligne décrivant une section (à travers laquelle on souhaite connaître le débit transitant).
-Cette polyligne peut avoir 2 formats : shape ou sinusx.
+Cette polyligne peut avoir 2 formats : shape ou sinusx.
+
+Ou dessiner directement sur le résultat une polyline en utilisant Sources / Poly Line Source.
Cliquer sur *Apply*.
Lancer le filtre *Rate Of Flow Through Section*.
-Une fenêtre s’ouvre :
+Une fenêtre s’ouvre :
-* Renseigner dans un premier temps l’input : il s’agit du fichier résultat.
-* Renseigner dans un deuxième temps la source : il s’agit de la polyligne chargée décrivant
+* Renseigner dans un premier temps l’input : il s’agit du fichier résultat.
+* Renseigner dans un deuxième temps la source : il s’agit de la polyligne (chargée ou dessinée) décrivant
la section à travers laquelle on souhaite connaître le débit transitant.
|paravis20| |paravis21|
|paravis22|
-* **Remarque :** Dans ce cas aussi, on peut changer le nombre de subdivisions de la polyligne chargée.
- Le résultat change en fonction du nombre de subdivisions. Le probleme est est en cours de traitement par l'équipe de Paraview.
-
- Dans notre cas, en fonction du nombre de subdivisions (1, 50 ou 100), le débit change :
-
-* 1 => 1487 m3/s – 1428 m3/s
-* 50 => 1603 m3/s – 1531 m3/s
-* 100 => 1592 m3/s – 1525 m3/s
-
Filtre Sediment Deposit
-----------------------
Cliquer sur *Apply*
Charger un fichier contour (polygone au niveau duquel on souhaite connaître les évolutions du fond (érosion / dépôt)).
-Ce polygone est au format shape.
+Ce polygone est au format shape.
+
+Ou dessiner directement sur le résultat une polyline en utilisant Sources / Poly Line Source. Penser à la fermer.
Cliquer sur *Apply*.
Lancer le filtre *Sediment Deposit*.
-Une fenêtre s’ouvre :
+Une fenêtre s’ouvre :
-* Renseigner dans un premier temps l’input : il s’agit du fichier résultat.
-* Renseigner dans un deuxième temps la source : il s’agit du polygone chargé décrivant le contour au niveau
+* Renseigner dans un premier temps l’input : il s’agit du fichier résultat.
+* Renseigner dans un deuxième temps la source : il s’agit du polygone (chargé ou dessiné) décrivant le contour au niveau
duquel on souhaite connaître les évolutions du fond.
|paravis24| |paravis25|
Cliquer sur *Apply*.
|paravis26|
-
+
On a alors 3 courbes tracées que l’on peut sélectionner comme on le souhaite :
|paravis27|
-----------------------------
Description à venir.
-
-
-.. only:: html
-
- :ref:`ref_exempleInondation`
.. |casCalculExtensionSimple| image:: /_static/casCalculExtensionSimple.png
:align: middle
-.. |casCalculExtensionSimple_0| image:: /_static/casCalculExtensionSimple_0.png
- :align: middle
-
-.. |casCalculExtensionSimple_1| image:: /_static/casCalculExtensionSimple_1.png
- :align: middle
-
.. |geometrieExtensionSimple| image:: /_static/geometrieExtensionSimple.png
:align: middle
tout en gardant si possible le maillage existant.
Nous supposons ici que nous ne disposons que d'un maillage au format MED, avec éventuellement des groupes à conserver.
Ce maillage doit être remis dans le système de coordonnées de la nouvelle étude, ici Lambert 93.
-Nous pouvons garder un décalage de l'origine.
Le script *<appli_xxx>/bin/salome/test/HYDRO/g022_extensionSimpleComplete.py* permet de générer ce maillage, pour les besoins du tutoriel.
-Ce maillage est alors disponible ici : *<appli_xxx>/bin/salome/test/HYDRO/tmp_test/garonne_2.med*.
-Il est en coordonnées Lambert 93, avec une origine à (430000, 6350000).
|maillageExistant|
-Une fois ce maillage généré, nous poursuivons dans une nouvelle étude SALOME, pour repartir d'un module HYDRO vierge.
-
Afin de pouvoir facilement dessiner les contours de l'extension dans Qgis ou avec le module HYDRO,
nous générons des shapefiles correspondants aux frontières du domaine existant,
-et des groupes d'edges présents dans le maillage.
+et des groupes présents dans le maillage.
Pour cela, nous utilisons le dialogue *Mesh edges to shapes* du module HYDRO (menu *HYDRO / Python plugins / Mesh edges to shapes*).
|dialogMeshEdgesToShapes|
Les coordonnées de l'origine du maillage doivent être correctement renseignées pour un positionnement correct des shapes générées.
-Comme indiqué précédemment, le maillage est en coordonnées Lambert 93, avec une origine à (430000, 6350000).
+Ici, le maillage est en coordonnées Lambert 93, avec une origine à (430000, 6350000).
Le répertoire de sortie contiendra les shapes et le maillage, enrichi d'un groupe correspondant aux mailles de bord (*FreeBorders*).
-Dans Qgis, nous pouvons importer les shapefiles *garonne_2_brd_FreeBorders.shp* et *garonne_2_brd_FreeBorders_pts.shp*
+Dans Qgis, nous importons les shapefiles *garonne_2_brd_FreeBorders.shp* et *garonne_2_brd_FreeBorders_pts.shp*
correspondants aux arêtes et noeuds des limites du maillage.
|qgisDomaineInitial|
|contourExtensionSimple|
-Il faut tout de suite faire le changement de repère local dans la nouvelle étude HYDRO : origine à (430000, 6350000).
-Ce repère local doit être **strictement identique** au repère local du maillage existant, pour permettre la création d'une extension.
-
Il est souvent pratique d'importer la shapefile du contour du domaine d'origine en tant que spline (éditer la polyligne pour la transformer en spline).
Nous ne tiendrons pas compte des éventuelles îles (régions insubmersibles) du domaine initial.
-* **Remarque:** La transformation du contour en spline facilite la sélection de groupe d'edges dans GEOM, si on a besoin d'en faire
- (une seule ligne à selectionner plutot que beaucoup de segments), mais, l'algorithme de maillage dans lequel on associe
- la ligne de contour au maillage d'origine échoue parfois avec certaines splines. Il faut alors reprendre le problème
- en changeant la nature de la polyligne.
-
-Nous créons deux objets naturels, le domaine étendu, et le domaine initial, **marqué en tant que région insubmersible**.
-Il ne faut pas tenir compte de l'ilôt dans le maillage existant.
-
- |casCalculExtensionSimple|
-
-Nous créons un cas de calcul, en mode manuel, en ne prenant que les deux objets naturels ci-dessus.
-
- |casCalculExtensionSimple_0|
-
+Nous créons un cas de calcul avec deux objets naturels, le domaine étendu, et le domaine initial, **marqué en tant que région insubmersible**.
Lors de la création du cas, il faut inclure les groupes de bord générés automatiquement, pour ne pas avoir à les recréer lors de l'étape de géométrie.
- |casCalculExtensionSimple_1|
+ |casCalculExtensionSimple|
Après export dans GEOM, la géométrie ressemble à ceci :
|maillageExtensionSimpleResult|
-Le maillage obtenu est conforme, il contient les groupes du maillage d'origine. Selon le sens de dessin des contours,
-les nouvelles mailles peuvent être à l'envers, comme dans l'image ci-dessus. Il reste alors à inverser les nouvelles mailles,
+Le maillage obtenu est conforme, il contient les groupes du maillage d'origine. Il reste à inverser les nouvelles mailles,
détecter et corriger les triangles de bord surcontraints, et, bien sûr, introduire d'éventuels raffinements dans l'extension
pour avoir un maillage plus réaliste.
|limiteExtensionAmont|
-**Attention**, il vaut mieux que la section correspondant au maillage d'origine soit transformée en spline pour faciliter sa manipulation par la suite
-*(voir la remarque du paragraphe précédent sur l'intérêt et les inconvénients de cette transformation)*.
+**Attention**, il vaut mieux que la section correspondant au maillage d'origine soit transformée en spline pour faciliter sa manipulation par la suite.
Nous pouvons alors créer un cas de calcul très simple avec un seul objet naturel construit sur le contour amont que l'on vient de produire.
Après export dans GEOM, on recrée un groupe correspondant à la limite du maillage d'origine, en vert ci-dessous.
Après import de la shapefile ajustée au premier maillage, nous réutilisons le dialogue *Adjust shapefile to mesh* avec le contour du deuxième maillage
(et la shapefile ajustée), en décochant toujours les options *"split mesh edges shapefile"* et *"Split shapefile to adjust"*.
-Nous pouvons importer la shapefile ajustée au deuxième maillage. **Dans cet exemple, nous ne transformons aucune des polylignes importées en spline**,
-une des splines posant problème à l'algorithme de maillage dans ce cas précis (il est toujours possible de refaire un test avec spline).
+Nous pouvons importer la shapefile ajustée au deuxième maillage.
Nous créons alors trois objets naturels avec la shapefile ajustée du raccord des domaines, et avec les deux domaines d'origine,
**marqués en tant que régions insubmersibles**.
+Notons que, dans cet exemple, aucune polyligne importée n'a été transformée en spline.
|troisRegions|
Ici,
- * *meshToCut* : le chemin complet du maillage à découper, translaté dans le même repère que le polygone de découpe, exemple :"/home/user/meshOrig.med"
+ * *meshToCut* : le chemin complet du maillage à découper, translaté dans le même repère que le polygone de découpe,
* *cutShape* : le chemin complet de la shapefile du polygone de découpe,
* *meshCut* : le chemin complet du maillage découpé, retranslaté.
* *offsetX, offsetY* : les coordonnées du repère local du maillage découpé.
objetsArtificiels.rst
manipsPolys.rst
landCoverMap.rst
+ completeCase.rst
+ copyMesh.rst
.. _ref_outilsReprise:
.. |genereInterpolz| image:: /_static/genereInterpolz.png
:align: middle
-.. |BottomImage| image:: /_static/Bottom.png
+.. |Bottom| image:: /_static/Bottom.png
:align: middle
.. |Capture_meshZ| image:: /_static/Capture_meshZ.png
Il faut sélectionner le cas de calcul dans la rubrique *HYDRO / CALCULATION CASE* de l'arbre d'étude.
Son nom apparaît dans la première ligne *Calculation cas* du dialogue.
-Le chemin complet du script à créer doit être renseigné dans la rubrique *Output path* (le nom du script doit se terminer par .py).
+Le chemin complet du script à créer doit être renseigné dans la rubrique *Output path*.
Il faut désigner le fichier du maillage de départ construit à l'étape précédente dans rubrique *MED file*.
Le champ s'affiche dans la vue 3D. Le menu contextuel de la vue 3D propose la commande *Representation / Surface with Edges*
- |BottomImage|
+ |Bottom|
Visualisation dans le module SMESH
----------------------------------
.. |BottomFriction| image:: /_static/BottomFriction.png
:align: middle
+.. |GenerateKs| image:: /_static/GenerateKs.png
+ :align: middle
+
+.. |InterfaceScriptKs| image:: /_static/InterfaceScriptKs.png
+ :align: middle
+
Import d'une Land Cover Map
===========================
Telemac exploite un champ des coefficients de Strickler aux noeuds du maillage. Ce champ a pour nom *BOTTOM FRICTION*.
Ce champ est ajouté au fichier MED du maillage, comme le champ d'altitude aux noeuds.
-Il faut adapter manuellement le script ci-dessous :
+Aller dans HYDROSOLVER / Hydro et cliquer sur Generate assignStrickler.py.
-Il faut recopier le script ci-dessous et l'adapter en fonction du cas de calcul
-et des noms de fichiers en entrée et en sortie.
-Il est possible d'utiliser le même fichier MED en entrée et en sortie.
+ |GenerateKs|
-.. literalinclude:: interpolStrickler.py
- :lines: 1-
-
+ |InterfaceScriptKs|
+
+Il faut sélectionner le cas de calcul dans la rubrique HYDRO / CALCULATION CASE de l’arbre d’étude. Son nom apparaît dans la première ligne Calculation case du dialogue.
+
+Le chemin complet du script à créer doit être renseigné dans la rubrique Output path.
+
+Il faut désigner le fichier du maillage de départ construit à l’étape précédente dans rubrique MED file.
+
Nous pouvons voir le champ résultat avec le module MED, comme pour le champ d'altitude.
|BottomFriction|
.. |initialFieldDialog| image:: /_static/initialFieldDialog.png
:align: middle
+.. |CheckBoundaryConditions| image:: /_static/CheckBoundaryConditions.png
+ :align: middle
+
+.. |CheckBoundaryConditions1| image:: /_static/CheckBoundaryConditions1.png
+ :align: middle
+
+.. |CheckBoundaryConditions2| image:: /_static/CheckBoundaryConditions2.png
+ :align: middle
+
+.. |CheckBoundaryConditions3| image:: /_static/CheckBoundaryConditions3.png
+ :align: middle
+
Une fois le maillage généré avec l'altimétrie, il reste à définir la nature des zones de conditions limites,
les valeurs des conditions limites de débit et de hauteur d'eau au cours du temps, et l'ensemble des
paramètres physiques et numériques du calcul.
qui liste l'ensemble des noeuds de conditions limites avec les types associés. Ce fichier est généré
automatiquement.
+Vérification que le contour est bien fermé et affichage des conditions limites
+------------------------------------------------------------------------------
+
+Il est possible de vérifier que le contour du modèle est bien fermé.
+
+Dans *HYDROSOLVER / Hydro*, cliquer sur *Check boundary conditions*.
+
+ |CheckBoundaryConditions|
+
+ |CheckBoundaryConditions1|
+
+Il faut désigner le fichier de maillage MED et le fichier des conditions limites.
+
+Puis on a la possibilité d’afficher les frontières solides et liquides du modèle :
+
+ |CheckBoundaryConditions2|
+
+Ou alors d’afficher les frontières de type Télémac (débit imposé, cote imposée…) :
+
+ |CheckBoundaryConditions3|
+
Valeurs des conditions limites de débit et hauteur d'eau
========================================================
Ici,
- * *meshFileIn* : le chemin complet du maillage d'origine avec son champ 'BOTTOM' à modifier, exemple :"/home/user/meshOrig.med"
+ * *meshFileIn* : le chemin complet du maillage d'origine avec son champ 'BOTTOM' à modifier,
* *meshFileOut* : le chemin complet du maillage avec le nouveau champ 'BOTTOM',
* *newBathy* : le chemin complet du nouveau fichier de bathymétrie,
* *'litMineur'* : le nom du groupe de mailles dont on modifie la bathymétrie,
* *Xoffset, Yoffset* : les coordonnées du repère local, optionnelles, si non nulles.
-
.. only:: html
:ref:`ref_outilsReprise`
