ci: Fix pipeline release job.

Stricklers and friction.

.editorconfig

@@ -0,0 +1,12 @@
+# Pamhyr .editorconfig
+
+root = true
+
+[*.py]
+charset = utf-8
+indent_style = space
+indent_size = 4
+tab_width = 8
+end_of_line = lf
+trim_trailing_whitespace = true
+insert_final_newline = true

.gitlab-ci.yml

@@ -15,96 +15,55 @@
 # along with this program.  If not, see <https://www.gnu.org/licenses/>.
 
 stages:
+  - downloads
   - build
   - test
   - package
   - release
 
+#############
+# DOWNLOADS #
+#############
+
+dl-mage-linux:
+  stage: downloads
+  tags:
+    - linux
+  rules:
+    - if: $CI_COMMIT_BRANCH == 'ci-test' || $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
+  script:
+    - curl -L -o mage.tgz https://gitlab.irstea.fr/jean-baptiste.faure/mage/-/releases/Test6/downloads/packages/mage_linux.tgz
+    - mkdir -p mage-linux
+    - cd mage-linux
+    - tar xvf ../mage.tgz
+  artifacts:
+    paths:
+      - mage-linux/mage
+      - mage-linux/mage_extraire
+      - mage-linux/mailleurPF
+
+dl-mage-windows:
+  stage: downloads
+  tags:
+    - linux
+  rules:
+    - if: $CI_COMMIT_BRANCH == 'ci-test' || $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
+  script:
+    - curl -L -o mage.tgz https://gitlab.irstea.fr/jean-baptiste.faure/mage/-/releases/Test6/downloads/packages/mage_windows.tgz
+    - mkdir -p mage-windows
+    - cd mage-windows
+    - tar xvf ../mage.tgz
+  artifacts:
+    paths:
+      - mage-windows/mage.exe
+      - mage-windows/mage_extraire.exe
+      - mage-windows/mailleurPF.exe
+
 #########
 # BUILD #
 #########
 
-build-mage-linux:
+set-version:
-  stage: build
-  tags:
-    - linux
-  rules:
-    - if: $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
-  script:
-    - curl -L -o mage.tgz https://gitlab.irstea.fr/jean-baptiste.faure/mage/-/releases/Test5/downloads/packages/mage_linux.tgz
-    - mkdir -p mage
-    - cd mage
-    - tar xvf ../mage.tgz
-  artifacts:
-    paths:
-      - mage/mage
-      - mage/mage_extraire
-      - mage/mailleurPF
-
-# build-mage-linux:
-#   stage: build
-#   tags:
-#     - linux
-#   rules:
-#     - if: $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
-#   variables:
-#     GIT_SUBMODULE_STRATEGY: recursive
-#     GIT_SUBMODULE_DEPTH: 20
-#   script:
-#     - cd mage/src/
-#     - mkdir -p build
-#     - cd build
-#     - cmake ..
-#     - make
-#   artifacts:
-#     paths:
-#       - mage/src/build/mage
-#       - mage/src/build/mage_extraire
-#       - mage/src/build/mailleurPF
-
-build-mage-windows:
-  stage: build
-  tags:
-    - linux
-  rules:
-    - if: $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
-  script:
-    - curl -L -o mage.tgz https://gitlab.irstea.fr/jean-baptiste.faure/mage/-/releases/Test5/downloads/packages/mage_windows.tgz
-    - mkdir -p mage
-    - cd mage
-    - tar xvf ../mage.tgz
-  artifacts:
-    paths:
-      - mage/mage.exe
-      - mage/mage_extraire.exe
-      - mage/mailleurPF.exe
-
-# build-mage-windows:
-#   stage: build
-#   tags:
-#     - wine
-#   rules:
-#     - if: $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
-#   variables:
-#     GIT_SUBMODULE_STRATEGY: recursive
-#     GIT_SUBMODULE_DEPTH: 20
-#   script:
-#     - cd mage/src/
-#     - sed s/add_compile_definitions(win)/add_compile_definitions(windows)/ CMakeLists.txt > __tmp__
-#     - echo set(CMAKE_Fortran_COMPILER "C:/Program\ Files/gcc/bin/gfortran.exe") > CMakeLists.txt
-#     - echo set(CMAKE_MAKE_PROGRAM "C:/Program Files\ \(x86\)/GnuWin32/bin/make.exe") >> CMakeLists.txt
-#     - type __tmp__ >> CMakeLists.txt
-#     - mkdir build
-#     - cd build
-#     - cmake -G "MinGW Makefiles" ..
-#     - make
-#   artifacts:
-#     paths:
-#       - mage/src/build/mage.exe
-#       - mage/src/build/mage_extraire.exe
-#       - mage/src/build/mailleurPF.exe
-
-build:
   stage: build
   tags:
     - linux
@@ -126,60 +85,214 @@ build-lang:
     paths:
       - src/lang/*.qm
 
+build-users-doc:
+  stage: build
+  tags:
+    - linux
+  script:
+    - cd ./doc/users/
+    - ./build.sh
+  artifacts:
+    paths:
+      - doc/users/documentation.pdf
+
+build-developers-doc:
+  stage: build
+  tags:
+    - linux
+  script:
+    - cd ./doc/dev/
+    - ./build.sh
+  artifacts:
+    paths:
+      - doc/dev/documentation.pdf
+
+build-linux:
+  stage: build
+  tags:
+    - linux
+  rules:
+    - if: $CI_COMMIT_BRANCH == 'ci-test' || $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
+  needs:
+    - job: dl-mage-linux
+      artifacts: true
+    - job: set-version
+      artifacts: true
+    - job: build-lang
+      artifacts: true
+    - job: build-developers-doc
+      artifacts: true
+    - job: build-users-doc
+      artifacts: true
+  script:
+    - mkdir -p linux
+    - cd linux
+    # Setup virtual env
+    - python3 -m venv linux-venv
+    - . linux-venv/bin/activate
+    - pip3 install -r ../requirements.txt
+    # Run Pyinstaller
+    - pyinstaller -y --paths linux-venv/lib/python3.8/site-packages ../src/pamhyr.py
+    # Create directory
+    - mkdir -p pamhyr
+    - mkdir -p pamhyr/pamhyr
+    # Copy Pamhyr info
+    - cp -v ../VERSION pamhyr/pamhyr/
+    - cp -v ../AUTHORS pamhyr/pamhyr/
+    - cp -v ../LICENSE pamhyr/pamhyr/
+    # Copy mage
+    - mkdir -p pamhyr/mage
+    - cp -v ../mage-linux/* pamhyr/mage/
+    # Copy Pamhyr
+    - cp -r dist/pamhyr pamhyr/
+    # Copy Pamhyr resources
+    - mkdir -p pamhyr/pamhyr/View
+    - mkdir -p pamhyr/pamhyr/View/ui
+    - mkdir -p pamhyr/pamhyr/View/ui/Widgets
+    - cp -r ../src/View/ui/*.ui pamhyr/pamhyr/View/ui/
+    - cp -r ../src/View/ui/ressources/ pamhyr/pamhyr/View/ui
+    - cp -r ../src/View/ui/Widgets/*.ui pamhyr/pamhyr/View/ui/Widgets
+    # Copy lang
+    - mkdir -p pamhyr/pamhyr/lang
+    - cp -r ../src/lang/*.qm pamhyr/pamhyr/lang
+    # Copy tests cases
+    - mkdir -p pamhyr/tests_cases/
+    - mkdir -p pamhyr/tests_cases/Saar
+    - cp ../tests_cases/Saar/Saar.pamhyr pamhyr/tests_cases/Saar/
+    # Copy documentations
+    - mkdir -p pamhyr/doc/
+    - cp ../doc/users/documentation.pdf pamhyr/doc/Pamhyr2-users.pdf
+    - cp ../doc/dev/documentation.pdf pamhyr/doc/Pamhyr2-dev.pdf
+    # Create running script
+    - echo "#!/bin/sh" > pamhyr/Pamhyr2
+    - echo "./pamhyr/pamhyr $@" >> pamhyr/Pamhyr2
+    - chmod u+x pamhyr/Pamhyr2
+  artifacts:
+    paths:
+      - linux/pamhyr
+
+build-windows:
+  stage: build
+  tags:
+    - wine
+  rules:
+    - if: $CI_COMMIT_BRANCH == 'ci-test' || $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
+  needs:
+    - job: dl-mage-windows
+      artifacts: true
+    - job: set-version
+      artifacts: true
+    - job: build-lang
+      artifacts: true
+    - job: build-developers-doc
+      artifacts: true
+    - job: build-users-doc
+      artifacts: true
+  script:
+    - mkdir windows
+    - cd windows
+    # Setup env
+    - python -m pip install -r ..\requirements.txt
+    # Run Pyinstaller
+    - pyinstaller -y ..\src\pamhyr.py
+    # Create directory
+    - mkdir pamhyr
+    - dir
+    - copy ..\packages\make-windows-dir.bat .\make-windows-dir.bat
+    - .\make-windows-dir.bat
+  artifacts:
+    paths:
+      - windows/pamhyr
+
 #########
 # TESTS #
 #########
 
-test:
+# test-windows:
-  stage: test
+#   stage: test
-  script:
+#   tags:
-    - echo "TODO PAMHYR tests"
+#     - wine
+#   needs:
+#     - job: set-version
+#       artifacts: true
+#     - job: build-windows
+#       artifacts: true
+#   script:
+#     - cd windows\pamhyr
+#     - pamhyr\pamhyr.exe hello
+
+# test-linux:
+#   stage: test
+#   tags:
+#     - linux
+#   needs:
+#     - job: set-version
+#       artifacts: true
+#     - job: build-linux
+#       artifacts: true
+#   script:
+#     - cd linux/pamhyr
+#     - ./Pamhyr2 hello
 
 ############
 # PACKAGES #
 ############
 
-linux-package:
+linux-package-tar:
   stage: package
   tags:
     - release
     - linux
   needs:
-    - job: build-lang
+    # - job: test-linux
-      artifacts: true
+    - job: build-linux
-    - job: build
-      artifacts: true
-    - job: build-mage-linux
       artifacts: true
   rules:
-    - if: $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
+    - if: $CI_COMMIT_BRANCH == 'ci-test' || $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
   artifacts:
     paths:
-      - packages/pamhyr-gnulinux-amd64.tar.xz
+      - linux/pamhyr-gnulinux.tar.xz
   script:
-    - cd packages
+    - cd linux
-    - ./linux.sh
+    - tar --xz -cf pamhyr-gnulinux.tar.xz ./pamhyr --checkpoint=.100
 
-windows-package:
+windows-package-zip:
+  stage: package
+  tags:
+    - release
+    - linux
+  needs:
+    # - job: test-windows
+    - job: build-windows
+      artifacts: true
+  rules:
+    - if: $CI_COMMIT_BRANCH == 'ci-test' || $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
+  artifacts:
+    paths:
+      - windows/pamhyr-windows.zip
+  script:
+    - cd windows
+    - zip -r pamhyr-windows.zip pamhyr
+
+windows-package-exe:
   stage: package
   tags:
     - release
     - wine
   needs:
-    - job: build-lang
+    # - job: test-windows
-      artifacts: true
+    - job: build-windows
-    - job: build
-      artifacts: true
-    - job: build-mage-windows
       artifacts: true
   rules:
-    - if: $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
+    - if: $CI_COMMIT_BRANCH == 'ci-test' || $CI_COMMIT_BRANCH == 'master' || $CI_COMMIT_TAG
   artifacts:
     paths:
-      - packages/pamhyr-win-amd64.exe
+      - windows/pamhyr-windows.exe
   script:
-    - cd packages
+    - cd windows
-    - ./windows.bat
+    - copy ..\packages\make-windows-exe.bat make-windows-exe.bat
+    - copy ..\packages\pamhyr.nsi pamhyr.nsi
+    - .\make-windows-exe.bat
 
 ###########
 # RELEASE #
@@ -191,16 +304,25 @@ tag-release:
     - release
     - linux
   needs:
-    - job: linux-package
+    - job: linux-package-tar
       artifacts: true
-    - job: windows-package
+    - job: windows-package-zip
+      artifacts: true
+    - job: windows-package-exe
+      artifacts: true
+    - job: build-users-doc
+      artifacts: true
+    - job: build-developers-doc
       artifacts: true
   rules:
     - if: $CI_COMMIT_TAG
   artifacts:
     paths:
-      - packages/pamhyr-gnulinux-amd64.tar.xz
+      - linux/pamhyr-gnulinux.tar.xz
-      - packages/pamhyr-win-amd64.exe
+      - windows/pamhyr-windows.zip
+      - windows/pamhyr-windows.exe
+      - doc/dev/documentation.pdf
+      - doc/users/documentation.pdf
   script:
     - cd packages
   release:
@@ -210,11 +332,23 @@ tag-release:
     ref: '$CI_COMMIT_SHA'
     assets:
       links:
-        - name: 'GNU/Linux amd64 (tar.xz)'
+        - name: 'GNU/Linux (tar.xz)'
-          url: '${CI_PROJECT_URL}/-/jobs/${CI_JOB_ID}/artifacts/raw/packages/pamhyr-gnulinux-amd64.tar.xz'
+          url: '${CI_PROJECT_URL}/-/jobs/${CI_JOB_ID}/artifacts/raw/linux/pamhyr-gnulinux.tar.xz'
-          filepath: '/packages/pamhyr-gnulinux-amd64.tar.xz'
+          filepath: '/linux/pamhyr-gnulinux.tar.xz'
           link_type: 'package'
-        - name: 'Windows amd64 (exe)'
+        - name: 'Windows (exe)'
-          url: '${CI_PROJECT_URL}/-/jobs/${CI_JOB_ID}/artifacts/raw/packages/pamhyr-win-amd64.exe'
+          url: '${CI_PROJECT_URL}/-/jobs/${CI_JOB_ID}/artifacts/raw/windows/pamhyr-windows.exe'
-          filepath: '/packages/pamhyr-win-amd64.exe'
+          filepath: '/windows/pamhyr-windows.exe'
           link_type: 'package'
+        - name: 'Windows (zip)'
+          url: '${CI_PROJECT_URL}/-/jobs/${CI_JOB_ID}/artifacts/raw/windows/pamhyr-windows.zip'
+          filepath: '/windows/pamhyr-windows.zip'
+          link_type: 'package'
+        - name: 'User documentations (pdf)'
+          url: '${CI_PROJECT_URL}/-/jobs/${CI_JOB_ID}/artifacts/raw/doc/users/documentation.pdf'
+          filepath: '/doc/users/documentation.pdf'
+          link_type: 'other'
+        - name: 'Dev documentations (pdf)'
+          url: '${CI_PROJECT_URL}/-/jobs/${CI_JOB_ID}/artifacts/raw/doc/dev/documentation.pdf'
+          filepath: '/doc/dev/documentation.pdf'
+          link_type: 'other'

README.md

@@ -1 +1,50 @@
-# PamHyr (Programme d’Aide à la Modélisation en Hydraulique des Rivières)
+# Pamhyr2
+
+Pamhyr is a free and open source graphical user interface for 1D hydro-sedimentary
+modelling of rivers.
+
+## Features
+
++ Edit a complex river hydraulics network
++ Define the river geometry
++ Define boundary conditions
++ Define lateral contribution
++ Edit river friction
++ Generate initial conditions
++ Edit river sedimentary layers
++ Save the study in a unique file
++ Run Mage 8 on a study and visualize results
+
+## Install
+
+### GNU/Linux
+
+Downloads [last
+realease](https://gitlab.irstea.fr/theophile.terraz/pamhyr/-/releases/permalink/latest)
+GNU/Linux binary archive.
+
+```shell
+tar xf pamhyr-gnulinux-amd64.tar.xz
+cd pamhyr
+./pamhyr
+```
+
+### Windows
+
+Downloads [last
+realease](https://gitlab.irstea.fr/theophile.terraz/pamhyr/-/releases/permalink/latest)
+Windows installer.
+
+### From source code
+
+```shell
+# Get source code
+git clone https://gitlab.irstea.fr/theophile.terraz/pamhyr/
+cd pamhyr
+
+# Install dependencies
+pip3 install -r requirements.txt
+
+# Run pamhyr
+./src/pamhyr.py
+```

doc/dev/build.sh

@@ -0,0 +1,10 @@
+#! /bin/sh
+
+# Build the PDF
+
+emacs -Q -q -nw --batch                                 \
+      --eval "(require 'org)"                           \
+      --eval "(require 'ox-latex)"                      \
+      --eval "(load-file \"../tools/setup.el\")"        \
+      --file ./documentation.org                        \
+      -f org-latex-export-to-pdf

doc/dev/documentation.bib

@@ -0,0 +1 @@
+../tools/ref.bib

doc/dev/documentation.org

@@ -0,0 +1,551 @@
+# documentation.org -- Pamhyr developers documentation
+# Copyright (C) 2023  INRAE
+#
+# This program 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.
+#
+# This program 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 this program.  If not, see <https://www.gnu.org/licenses/>.
+
+# -*- coding: utf-8 -*-
+
+#+STARTUP: indent
+
+#+INCLUDE: ../tools/macro.org
+#+INCLUDE: ../tools/latex.org
+
+#+TITLE: Developers documentation
+#+SUBTITLE: Version: {{{version}}}
+#+AUTHOR: {{{INRAE}}}
+
+#+OPTIONS: toc:t
+#+LANGUAGE: UKenglish
+
+#+NAME: attr_wrap
+#+HEADER: :var width="\\textwidth"
+#+HEADER: :var caption=""
+#+HEADER: :var smallcaption=""
+#+HEADER: :var name=""
+#+HEADER: :var data=""
+#+HEADER: :var float="nil"
+#+BEGIN_SRC sh :results output :exports none
+  echo "#+CAPTION[$smallcaption]: $caption"
+  echo "#+NAME: $name"
+  echo "#+ATTR_LATEX: :width $width :float $float"
+  echo "$data"
+#+END_SRC
+
+#+BEGIN_abstract
+This document is for the use of developers. It describes the project
+architecture, the tools available to assist development and
+debugging. It also describes the procedures for creating packages and
+the configurations required to set up the gitlab runners. Finally,
+this document explains how documentation is written and modified, and
+how to contribute to the project by modifying, improving or adding
+documentation, translations or code.
+#+END_abstract
+
+* Introduction
+
+Pamhyr2 is free and open source software (FOSS) graphical user
+interface (GUI) for 1D hydro-sedimentary modelling of rivers developed
+in Python (with version 3.8). It use PyQt at version 5 and matplotlib
+in version 3.4.1 or later for the user insterface (see
+{{{file(/requirements.txt)}}} for details). The architecture of
+project code follow the Qt Model/View architecture [fn:qt-arch] (see
+details in section [[Architecture]]). Pamhyr2 packages can be build
+manually (see section [[Building packages]]), but there are automatically
+build with the gitlab-ci (see the section [[Setup the CI
+environment]]). Documentation files are written with org-mode[fn:org],
+let see section [[Documentation files]]. Finally, to see the contribution
+rules, see the section [[How to contribute?]].
+
+[fn:qt-arch] Qt Model/View documentation:
+https://doc.qt.io/qt-5/model-view-programming.html (last access
+2023-09-15)
+[fn:org] The org-mode website: https://orgmode.org/ (last access
+2023-09-15)
+
+* TODO Architecture
+
+Pamhyr2 architecture is based on Qt Model/View programming[fn:qt-mv] architecture.
+
+#+name: graph-pipeline
+#+header: :results drawer
+#+header: :exports results
+#+header: :post attr_wrap(width="12cm", data=*this*, name="graph-pipeline", caption="Pamhyr2 solver execution pipeline architecture scheme", float="t")
+#+begin_src dot :file "images/graph-pipeline.png" :cache no
+  digraph {
+        node[colorscheme=set19,shape=box,style="filled",fillcolor=9];
+
+        subgraph cluster0 {
+            label="Pamhyr2"
+            config[label="Config",fillcolor=5];
+            model[label="Model",fillcolor=2];
+            obj[label="Solver",fillcolor=6];
+            results[label="Results",fillcolor=2];
+            view[label="View",fillcolor=1];
+
+            results -> model[style="dotted"];
+            results -> obj[style="dotted"];
+        }
+
+        config -> obj[label=""];
+        obj -> model;
+
+        subgraph cluster1{
+            label="System";
+            in[label="Solver input files",shape=note];
+            out[label="Solver output files",shape=note];
+            bin[label="Solver binary",shape=note];
+        }
+
+        obj -> in[label="Write"];
+        bin -> in[label="Read"];
+        obj -> bin[label="Execute"];
+        bin -> out[label="Write"];
+        obj -> results[label="Create"];
+        results -> out[label="Read"];
+        view -> model;
+        view -> results;
+  }
+#+end_src
+
+#+name: graph-architecture
+#+header: :results drawer
+#+header: :exports results
+#+header: :post attr_wrap(width="7cm", data=*this*, name="graph-architecture", caption="Pamhyr2 Model/View architecture scheme (inspired by Qt Model/View architecture [[https://doc.qt.io/qt-5/model-view-programming.html]])", float="t")
+#+begin_src dot :file "images/graph-architecture.png" :cache no
+  digraph {
+        node[colorscheme=set19,shape=box,style="filled",fillcolor=white];
+
+        save[label="Pamhyr save",fillcolor="9",shape=note];
+        model[label="Model",fillcolor="2"];
+        view[label="View",fillcolor="1"];
+        delegate[label="Delegate",fillcolor="3"];
+        undocommand[label="UndoCommand",fillcolor="4"];
+        user[label="User",shape=ellipse];
+
+        model -> save[label="Save"];
+        save -> model[label="Load"];
+
+        model -> view[label="Rendering"];
+        view -> delegate[label="Rendering"];
+        delegate -> undocommand[label="Create"];
+        undocommand -> model[label="Modify"]
+
+        view -> user[label="Vizualize"];
+        user -> delegate[label="Modify"];
+  }
+#+end_src
+
+#+name: graph-model
+#+header: :results drawer
+#+header: :exports results
+#+header: :post attr_wrap(width="16cm", data=*this*, name="graph-model", caption="Pamhyr2 Model architecture scheme", float="t")
+#+begin_src dot :file "images/graph-model.png" :cache no
+  digraph {
+        node[colorscheme=set19,shape=box,style="filled",fillcolor="2"];
+
+        study[label="Study"];
+        river[label="River"];
+
+        subgraph cluster00 {
+            label="Network"
+            rnode[label="RiverNode"];
+            redge[label="RiverReach"];
+        }
+
+        frictionlist[label="FrictionList"];
+
+        subgraph cluster01 {
+            label="Stricklers";
+            stricklers[label="Stricklers"];
+            stricklerslist[label="StricklersList"];
+        }
+        subgraph cluster02 {
+            label="BoundaryCondition";
+            boundaryconditionlist[label="BoundaryConditionList"];
+            boundarycondition[label="BoundaryCondition"];
+        }
+        subgraph cluster03 {
+            label="LateralContribution";
+            lateralcontributionlist[label="LateralContributionList"];
+            lateralcontribution[label="LateralContribution"];
+        }
+        subgraph cluster04 {
+            label="InitialConditions";
+            initialconditionsdict[label="InitialConditionsDict"];
+            initialconditions[label="InitialConditions"];
+        }
+
+        solverparameterslist[label="SolverParametersList"];
+
+        subgraph cluster05 {
+            label="Sediment";
+            sedimentlayerlist[label="SedimentLayerList"];
+            sedimentlayer[label="SedimentLayer"];
+            layer[label="Layer"];
+        }
+
+        subgraph cluster06 {
+            label="Greometry"
+            georeach[label="Reach"];
+            geocrosssection[label="Cross-section"];
+            geopoint[label="Point"];
+        }
+
+        study -> river;
+        river -> rnode;
+        river -> redge;
+        redge -> rnode;
+        river -> boundaryconditionlist -> boundarycondition -> rnode;
+        river -> lateralcontributionlist -> lateralcontribution -> redge;
+        river -> initialconditionsdict -> initialconditions;
+        initialconditions -> redge;
+        river -> stricklerslist -> stricklers;
+        river -> solverparameterslist;
+        river -> sedimentlayerlist -> sedimentlayer -> layer;
+        redge -> frictionlist -> stricklers;
+        redge -> georeach -> geocrosssection -> geopoint;
+        geocrosssection -> sedimentlayer;
+        geopoint -> sedimentlayer;
+
+  }
+#+end_src
+
+[fn:qt-mv] The Qt Model/View documentation web page:
+https://doc.qt.io/qt-5/model-view-programming.html
+
+** TODO Model
+** TODO View
+*** TODO UI file
+*** TODO Window class
+** TODO UndoCommand
+** TODO Table
+** TODO Plot
+** TODO Solver
+** TODO Unit tests
+** TODO The debug mode
+* TODO Build the project
+
+The project uses gitlab-ci runners to build packages, but it is possible
+to build packages manually.
+
+** Building packages
+
+If you need an hand made package, you can script available in
+{{{file(packages)}}} directory.
+
+*** GNU/Linux
+
+On GNU/Linux building GNU/Linux packages is easy, you just need python
+in version 3.8 must be installed with venv and pyinstaller packages
+(see Listing [[linux-env-deb]] for Debian and derived system). Finally,
+run the {{{file(linux.sh)}}} script (see Listing [[linux-pkg]]).
+
+#+NAME: linux-env-deb
+#+CAPTION: Install environment on GNU/Linux
+#+begin_src shell
+sudo apt install python3.8
+python3 -m pip install venv
+python3 -m pip install pyinstaller
+#+end_src
+
+#+NAME: linux-pkg
+#+CAPTION: Build GNU/Linux package
+#+begin_src shell
+cd packages
+./linux.sh
+#+end_src
+
+*** Windows
+
+To make the Windows packages you have two choice: If you use Windows
+you can use the script {{{file(packages/windows.bat)}}}, other else
+you can use the script {{{file(packages/wine.sh)}}}. Each script need
+a specific software environment.
+
+On windows, you needs python on version 3.8, pyinstaller and
+NSIS[fn:nsis] installed. On GNU/Linux you need wget, wine and
+winetricks installed.
+
+[fn:nsis] The NSIS web site: https://sourceforge.net/projects/nsis/
+
+** TODO Setup the CI environment
+* Documentation files
+
+This document and the user documentation are org files. This text file
+format is formatted so that it can be exported in different formats:
+PDF (with latex), ODT, HTML, etc. It was originally designed for the
+GNUEmacs[fn:emacs] text editor, but can be edited with any text editor. Here we
+take a look at the different features used in these documents.
+
+[fn:org] The org-mode website: https://orgmode.org/ (last access
+2023-09-15)
+[fn:emacs] The GNUEmacs project website: https://gnu.org/s/emacs/
+(last access 2023-09-15)
+
+** Org-mode
+*** Document structure
+
+Org uses the =*= character to define a new document section. To add a
+sub-section, you can add an additional =*= to the current section[fn::
+See document structure documentation:
+https://orgmode.org/org.html#Headlines (last access 2023-09-15)].
+
+#+BEGIN_EXAMPLE
+ * Top level headline
+ ** Second level
+ *** Third level
+     some text
+ *** Third level
+     more text
+ * Another top level headline
+#+END_EXAMPLE
+
+*** Format
+
+Org-mode is a markup file, using markup in the text to modify the
+appearance of a portion of text[fn:: See markup documentation:
+https://orgmode.org/org.html#Emphasis-and-Monospace (last access
+2023-09-15)].
+
+| Markup             | Results          |
+|--------------------+------------------|
+| =*Bolt*=           | *Bolt*           |
+| =/Italic/=         | /Italic/         |
+| =_underline_=      | _underline_      |
+| ==verbatim==       | =verbatim=       |
+| =~code~=           | ~code~           |
+| =+strike-through+= | +strike-through+ |
+
+*** Source code blocks
+
+You can add some code blocks[fn:: See org-mode documentation for
+source code: https://orgmode.org/org.html#Working-with-Source-Code
+(last access 2023-09-15)] in the document.
+
+Here is an example for python source code:
+#+BEGIN_EXAMPLE
+#+CAPTION: Get os type name in Python code
+#+begin_src python
+import os
+
+print(f"Document build on system: {os.name}")
+#+end_src
+#+END_EXAMPLE
+
+If you use GNUEmacs, it is also possible to run the code inside a
+block and export (or not) the reuslts in the document.
+
+#+OPTIONS: float:nil
+#+CAPTION: Get os type name in Python code
+#+begin_src python :python python3 :results output :exports both :noweb yes
+import os
+
+print(f"Document build on system: {os.name}")
+#+end_src
+
+#+RESULTS:
+: Document build on system: posix
+
+*** LaTeX
+
+If we export the file to PDF, org-mode use {{{latex}}}. So we can add
+some piece of {{{latex}}} into the document[fn:: See {{{latex}}} part
+in documentation: https://orgmode.org/org.html#Embedded-LaTeX (last
+access 2023-09-15)]. For exemple, we can add math formula like
+=$E=mc^2$= ($E=mc^2$) or =\[E=mc^2\]=:
+
+\[E=mc^2\]
+
+But we can also add every type of {{{latex}}}:
+
+#+BEGIN_EXAMPLE
+# Add latex in line
+#+LATEX: <my line of latex>
+
+# Add multiple line of LaTeX
+#+BEGIN_EXPORT latex
+<my latex here>
+#+END_EXPORT
+#+END_EXAMPLE
+
+It is also possible to add specific {{{latex}}} file header with
+=#+LATEX_HEADER=. In this document we use the file
+{{{file(doc/tools/latex.org)}}} for all {{{latex}}} headers.
+
+*** Macro
+
+In this document, we use a few macros[fn:: See marcos documentation
+https://orgmode.org/org.html#Macro-Replacement (last access
+2023-09-15)] to simplify writing. They allow you to define sequences
+of text to be replaced, so that the macro name is replaced by its
+value. They are defined in the {{{file(doc/tools/macro.org)}}}
+file. Once defined, they can be used in the document as follows:
+={{{<macro-name>}}}=. You can also have macros with arguments, in this
+case: ={{{<macro-name>(arg1,...)}}}=. Les macros peuvent aussi
+utiliser du code emacs-lisp.
+
+#+BEGIN_EXAMPLE
+# Exemple of macro définition
+
+#+MACRO: toto               tata
+#+MACRO: add                \(($1 + $2)\)
+#+MACRO: emacs-version      (eval (nth 2 (split-string (emacs-version))))
+#+END_EXAMPLE
+
+#+MACRO: toto               tata
+#+MACRO: add                \(($1 + $2)\)
+#+MACRO: emacs-version      (eval (nth 2 (split-string (emacs-version))))
+
+Macro apply:
+- Marco ={{{toto}}}=: {{{toto}}}
+- Marco ={{{add(x,y)}}}=: {{{add(x,y)}}}
+- Marco ={{{emacs-version}}}=: {{{emacs-version}}}
+
+*** Footnotes
+
+Footnote in org-mode is define with marker =[fn:...]=[fn:: Create
+footnotes in org-mode documentation
+https://orgmode.org/org.html#Creating-Footnotes (last access
+2023-09-15)]:
+#+BEGIN_EXAMPLE
+The Org website[fn:1] now looks a lot better than it used to.
+...
+[fn:1] The link is: https://orgmode.org
+#+END_EXAMPLE
+or:
+#+BEGIN_EXAMPLE
+The Org website[fn:: The link is: https://orgmode.org] now looks
+a lot better than it used to.
+...
+#+END_EXAMPLE
+
+*** References
+
+The references use the {{{latex}}} bibtex tools. The bib file is in
+{{{file(/doc/tools/ref.bib)}}} and use for developers and user
+documentation. In document, use ={{{cite(<name>)}}}= to cite a paper.
+
+** Export
+
+To export the files, a {{{file(build.sh)}}} script is available in the org
+files directories. On GNU/Linux system you can build the documentation
+PDF file with the command =./build.sh=. Texlive package must be
+installed, you can install only needed packages or all texlive
+packages, for example on Debian (and some derived system) use command
+Listing [[texlive-install]].
+
+#+NAME: texlive-install
+#+CAPTION: Installation command for texlive full on Debian system
+#+begin_src shell
+  sudo apt install texlive-full
+#+end_src
+
+Some org-mode configuration used in documentations files are define in
+=/doc/tools/=:
+- {{{file(PamhyrDoc.cls)}}}: The {{{latex}}} document class
+- {{{file(macro.org)}}}: Available macro
+- {{{file(latex.org)}}}: {{{latex}}} configutation for documentations
+  files
+- {{{file(setup.el)}}}: GNUEmacs configuration to build documentations
+- {{{file(ref.bib)}}}: Bibtex files for documentations files
+
+* How to contribute?
+
+Pamhyr2 is free software: you can redistribute it and/or modify it
+under the terms of the GNU General Public License[fn:license], either
+version 3 of the License, or any later version.
+
+[fn:license] The GPLv3 web page:
+https://www.gnu.org/licenses/gpl-3.0.en.html
+
+** Guidelines
+
+To contribute to Pamhyr2, we expect a minimum of respect between
+contributors. We therefore ask you to respect the following rules
+regarding communication and contribution content:
++ No gender, racial, religious or social discrimination
++ No insults, personal attacks or potentially offensive remarks
++ Pamhyr2 is free software, and intended to remain so, so take care
+  with the licensing of libraries and external content you want to add
+  to the project
++ Humour or hidden easter eggs are welcome if they respect the
+  previous rules
+
+** Make a contribution
+
+There are several ways to contribute: you can report a bug by creating
+an issue on the project's gitlab page[fn:p2-gitlab], or you can create
+a merge request on the same page with the changes you have made to the
+code, translation or documentation.
+
+The Pamhyr2 copyright is owned by INRAE[fn:inrae], but we keep a
+record of each contributors. If you made a modification to pamhyr2
+software, please add your name at the end of {{{file(AUTHORS)}}} file
+and respect the Listing [[auth-format]] format. You can update this file
+information for following contribution.
+
+#+NAME: auth-format
+#+CAPTION: =AUTHORS= file format
+#+begin_src text
+<first name> <last name> [(optional) email], <organisation>, <years>
+#+end_src
+
+#+CAPTION: Current =AUTHORS= file
+#+INCLUDE: "../../AUTHORS" src text
+
+[fn:p2-gitlab] The Pamhyr2 Gitlab project page:
+https://gitlab.irstea.fr/theophile.terraz/pamhyr
+[fn:inrae] The INRAE web site: https://www.inrae.fr/
+
+** Translate
+
+You can improve or add translation for the project. To contribute to
+Pamhyr2 translate, you need to use Qt Linguist[fn:qt-linguist]. Open
+Qt-linguist and edite the translation ({{{file(.ts)}}}) file, finally,
+commit the new version of file and make a merge request.
+
+If you want add a new language, edit the script
+{{{file(src/lang/create_ts.sh)}}} like Listing [[ts-it]]. Run the script
+and open the new file with Qt-linguist, setup target language (Figure
+[[qt-linguist-setup]]) and complete translation. Finally, commit the new
+file and make a merge request.
+
+#+NAME: ts-it
+#+CAPTION: Example of modified {{{file(src/lang/create_ts.sh)}}} to add italian (it) translate for Pamhyr2
+#+begin_src shell
+  ...
+  LANG="fr it"
+  ...
+#+end_src
+
+#+NAME: qt-linguist-setup
+#+ATTR_HTML: :width 8cm
+#+ATTR_LATEX: :width 8cm
+#+CAPTION: Qt linguist lang setup example with italian.
+[[../images/Qt-linguist-setup-lang.png]]
+
+[fn:qt-linguist] The Qt linguist documentation web page:
+https://doc.qt.io/qt-5/qtlinguist-index.html (last access 2023-09-18)
+
+** Code contribution
+
+If you are developper you can improve and/or add features to
+Pamhyr2. Please, follow the architecture described in section
+[[Architecture]] as closely as possible. Keep the code simple, clear and
+efficient as possible. The master branch is reserved for the project
+maintainer; you can create a new branch or fork the project before the
+request.
+
+
+{{{biblio}}}

doc/paper/2023/simhydro/ajustement-for-paper.patch

@@ -0,0 +1,692 @@
+diff --git a/packages/linux.sh b/packages/linux.sh
+index 215047e..a19bfb6 100755
+--- a/packages/linux.sh
++++ b/packages/linux.sh
+@@ -76,11 +76,11 @@ mv ../pamhyr-src.tar.gz ./
+ 
+ echo " *** MAKE BIN PACKAGE"
+ 
+-OLD_PWD=$PWD
+-cd dist/
+-tar --xz -cf pamhyr-gnulinux-amd64.tar.xz pamhyr --checkpoint=.100
+-cd $OLD_PWD
++# OLD_PWD=$PWD
++# cd dist/
++# tar --xz -cf pamhyr-gnulinux-amd64.tar.xz pamhyr --checkpoint=.100
++# cd $OLD_PWD
+ 
+-mv dist/pamhyr-gnulinux-amd64.tar.xz ./
++# mv dist/pamhyr-gnulinux-amd64.tar.xz ./
+ 
+ echo " *** DONE"
+diff --git a/src/Model/Friction/Friction.py b/src/Model/Friction/Friction.py
+index 18167a8..5bf2302 100644
+--- a/src/Model/Friction/Friction.py
++++ b/src/Model/Friction/Friction.py
+@@ -170,6 +170,7 @@ class Friction(SQLSubModel):
+     @begin_strickler.setter
+     def begin_strickler(self, strickler):
+         self._begin_strickler = strickler
++        self._end_strickler = strickler # HACK: Only one Strickler
+         self._status.modified()
+ 
+     @property
+diff --git a/src/Model/Stricklers/Stricklers.py b/src/Model/Stricklers/Stricklers.py
+index baa484f..a5670d2 100644
+--- a/src/Model/Stricklers/Stricklers.py
++++ b/src/Model/Stricklers/Stricklers.py
+@@ -38,7 +38,7 @@ class Stricklers(SQLSubModel):
+         else:
+             self.id = id
+ 
+-        Stricklers._id_cnt = max(Stricklers._id_cnt + 1, self.id)
++        Stricklers._id_cnt = max(Stricklers._id_cnt + 1, id + 1)
+ 
+         self._name = name
+         self._comment = comment
+diff --git a/src/View/BoundaryCondition/translate.py b/src/View/BoundaryCondition/translate.py
+index c4ff8f8..172999a 100644
+--- a/src/View/BoundaryCondition/translate.py
++++ b/src/View/BoundaryCondition/translate.py
+@@ -29,9 +29,9 @@ _translate = QCoreApplication.translate
+ long_types = {
+     "ND": _translate("BoundaryCondition", "Not defined"),
+     "PC": _translate("BoundaryCondition", "Ponctual contribution"),
+-    "TZ": _translate("BoundaryCondition", "Time over Z"),
+-    "TD": _translate("BoundaryCondition", "Time over Discharge"),
+-    "ZD": _translate("BoundaryCondition", "Z over Discharge"),
++    "TZ": _translate("BoundaryCondition", "TimeZ"),
++    "TD": _translate("BoundaryCondition", "TimeDischarge"),
++    "ZD": _translate("BoundaryCondition", "ZDischarge"),
+     "SL": _translate("BoundaryCondition", "Solid"),
+ }
+ 
+diff --git a/src/View/Frictions/PlotStricklers.py b/src/View/Frictions/PlotStricklers.py
+index 71ff6e6..b9060bc 100644
+--- a/src/View/Frictions/PlotStricklers.py
++++ b/src/View/Frictions/PlotStricklers.py
+@@ -106,6 +106,8 @@ class PlotStricklers(APlot):
+                 )
+                 self.draw_frictions(lst, color="b")
+ 
++        self.canvas.axes.invert_xaxis()
++
+         self.canvas.figure.tight_layout()
+         self.canvas.figure.canvas.draw_idle()
+         if self.toolbar is not None:
+diff --git a/src/View/Frictions/Window.py b/src/View/Frictions/Window.py
+index 4af3b66..f46fae3 100644
+--- a/src/View/Frictions/Window.py
++++ b/src/View/Frictions/Window.py
+@@ -112,10 +112,10 @@ class FrictionsWindow(ASubMainWindow, ListedSubWindow):
+         )
+ 
+         table.setItemDelegateForColumn(
+-            3, self._delegate_stricklers
++            2, self._delegate_stricklers
+         )
+         table.setItemDelegateForColumn(
+-            4, self._delegate_stricklers
++            3, self._delegate_stricklers
+         )
+ 
+         table.setSelectionBehavior(QAbstractItemView.SelectRows)
+diff --git a/src/View/Frictions/translate.py b/src/View/Frictions/translate.py
+index 21981ab..82071e9 100644
+--- a/src/View/Frictions/translate.py
++++ b/src/View/Frictions/translate.py
+@@ -21,10 +21,10 @@ from PyQt5.QtCore import QCoreApplication
+ _translate = QCoreApplication.translate
+ 
+ table_headers = {
+-    "name": _translate("Frictions", "Name"),
++    # "name": _translate("Frictions", "Name"),
+     # "edge": _translate("Frictions", "Reach"),
+     "begin_kp": _translate("Frictions", "Begin kp (m)"),
+     "end_kp": _translate("Frictions", "End kp (m)"),
+-    "begin_strickler": _translate("Frictions", "Begin strickler"),
+-    "end_strickler": _translate("Frictions", "End strickler"),
++    "begin_strickler": _translate("Frictions", "Strickler"),
++    # "end_strickler": _translate("Frictions", "End strickler"),
+ }
+diff --git a/src/View/Geometry/PlotAC.py b/src/View/Geometry/PlotAC.py
+index fe6d933..e43f39e 100644
+--- a/src/View/Geometry/PlotAC.py
++++ b/src/View/Geometry/PlotAC.py
+@@ -67,18 +67,18 @@ class PlotAC(APlot):
+         gl = self.data.profile(selected_profile).names()
+ 
+         self.canvas.axes.set_xlabel(
+-            _translate("MainWindow_reach", "Abscisse en travers (m)"),
++            _translate("MainWindow_reach", "Transverse abscissa (m)"),
+             color='green', fontsize=12
+         )
+         self.canvas.axes.set_ylabel(
+-            _translate("MainWindow_reach", "Cote (m)"),
++            _translate("MainWindow_reach", "Height (m)"),
+             color='green', fontsize=12
+         )
+         self.canvas.figure.tight_layout()
+ 
+-        label_before_plot_selected = _translate("MainWindow_reach", "Profil précédent")
+-        label_plot_selected = _translate("MainWindow_reach", "Profil sélectionné")
+-        label_after_plot_selected = _translate("MainWindow_reach", "Profil suivant")
++        label_before_plot_selected = _translate("MainWindow_reach", "Previous cross-section")
++        label_plot_selected = _translate("MainWindow_reach", "Cross-section")
++        label_after_plot_selected = _translate("MainWindow_reach", "Next cross-section")
+         color_before_plot_selected = "k"  # 'grey'
+         color_plot_selected = 'b'
+         color_after_plot_selected = 'm'
+diff --git a/src/View/Geometry/PlotKPZ.py b/src/View/Geometry/PlotKPZ.py
+index eb7d5e2..64057aa 100644
+--- a/src/View/Geometry/PlotKPZ.py
++++ b/src/View/Geometry/PlotKPZ.py
+@@ -158,6 +158,8 @@ class PlotKPZ(APlot):
+             color='lightgrey'
+         )
+ 
++        self.canvas.axes.invert_xaxis()
++
+         self.canvas.figure.tight_layout()
+         self.canvas.figure.canvas.draw_idle()
+         if self.toolbar is not None:
+diff --git a/src/View/Geometry/Window.py b/src/View/Geometry/Window.py
+index 7776d73..1c8cb5d 100644
+--- a/src/View/Geometry/Window.py
++++ b/src/View/Geometry/Window.py
+@@ -111,7 +111,7 @@ class GeometryWindow(ASubMainWindow, ListedSubWindow):
+         self.ui.btn_sort_desc.triggered.connect(self.sort_descending)
+         self.ui.btn_move_up.triggered.connect(self.move_row_up)
+         self.ui.btn_move_down.triggered.connect(self.move_row_down)
+-        self.ui.btn_end_editing.triggered.connect(self.handleSave)
++        # self.ui.btn_end_editing.triggered.connect(self.handleSave)
+         self.ui.btn_add.triggered.connect(self.insert_row)
+         self.ui.btn_delete.triggered.connect(self.delete_rows)
+         self.ui.btn_edit.triggered.connect(self.edit_profile)
+diff --git a/src/View/Geometry/mainwindow_ui_reach.py b/src/View/Geometry/mainwindow_ui_reach.py
+index 897b698..60bcfa6 100644
+--- a/src/View/Geometry/mainwindow_ui_reach.py
++++ b/src/View/Geometry/mainwindow_ui_reach.py
+@@ -95,29 +95,29 @@ class Ui_MainWindow(object):
+         self.btn_edit.setIcon(icon_btn_edit)
+         self.toolBar.addAction(self.btn_edit)
+ 
+-        self.btn_copy_selected_profile = QtWidgets.QAction(self.widget)
+-        self.btn_copy_selected_profile.setObjectName("btn_copy_selected_profile")
+-        icon_btn_copy_selected_profile = QtGui.QIcon()
+-        icon_btn_copy_selected_profile.addPixmap(
+-            QtGui.QPixmap(f"{icon_path}gtk-copy.png"))
+-        self.btn_copy_selected_profile.setIcon(icon_btn_copy_selected_profile)
+-        self.toolBar.addAction(self.btn_copy_selected_profile)
+-
+-        self.btn_paste_selected_profile = QtWidgets.QAction(self.widget)
+-        self.btn_paste_selected_profile.setObjectName("btn_paste_selected_profile")
+-        icon_btn_paste_selected_profile = QtGui.QIcon()
+-        icon_btn_paste_selected_profile.addPixmap(
+-            QtGui.QPixmap(f"{icon_path}gtk-paste.png"))
+-        self.btn_paste_selected_profile.setIcon(icon_btn_paste_selected_profile)
+-        self.toolBar.addAction(self.btn_paste_selected_profile)
+-
+-        self.btn_duplicate_selected_profile = QtWidgets.QAction(self.widget)
+-        self.btn_duplicate_selected_profile.setObjectName("btn_duplicate_selected_profile")
+-        icon_btn_duplicate_selected_profile = QtGui.QIcon()
+-        icon_btn_duplicate_selected_profile.addPixmap(
+-            QtGui.QPixmap(f"{icon_path}gtk-dnd-multiple.png"))
+-        self.btn_duplicate_selected_profile.setIcon(icon_btn_duplicate_selected_profile)
+-        self.toolBar.addAction(self.btn_duplicate_selected_profile)
++        # self.btn_copy_selected_profile = QtWidgets.QAction(self.widget)
++        # self.btn_copy_selected_profile.setObjectName("btn_copy_selected_profile")
++        # icon_btn_copy_selected_profile = QtGui.QIcon()
++        # icon_btn_copy_selected_profile.addPixmap(
++        #     QtGui.QPixmap(f"{icon_path}gtk-copy.png"))
++        # self.btn_copy_selected_profile.setIcon(icon_btn_copy_selected_profile)
++        # self.toolBar.addAction(self.btn_copy_selected_profile)
++
++        # self.btn_paste_selected_profile = QtWidgets.QAction(self.widget)
++        # self.btn_paste_selected_profile.setObjectName("btn_paste_selected_profile")
++        # icon_btn_paste_selected_profile = QtGui.QIcon()
++        # icon_btn_paste_selected_profile.addPixmap(
++        #     QtGui.QPixmap(f"{icon_path}gtk-paste.png"))
++        # self.btn_paste_selected_profile.setIcon(icon_btn_paste_selected_profile)
++        # self.toolBar.addAction(self.btn_paste_selected_profile)
++
++        # self.btn_duplicate_selected_profile = QtWidgets.QAction(self.widget)
++        # self.btn_duplicate_selected_profile.setObjectName("btn_duplicate_selected_profile")
++        # icon_btn_duplicate_selected_profile = QtGui.QIcon()
++        # icon_btn_duplicate_selected_profile.addPixmap(
++        #     QtGui.QPixmap(f"{icon_path}gtk-dnd-multiple.png"))
++        # self.btn_duplicate_selected_profile.setIcon(icon_btn_duplicate_selected_profile)
++        # self.toolBar.addAction(self.btn_duplicate_selected_profile)
+ 
+         self.btn_sort_asc = QtWidgets.QAction(self.widget)
+         self.btn_sort_asc.setObjectName("btn_sort_asc")
+@@ -151,13 +151,13 @@ class Ui_MainWindow(object):
+         self.btn_move_down.setIcon(icon_btn_move_down)
+         self.toolBar.addAction(self.btn_move_down)
+ 
+-        self.btn_end_editing = QtWidgets.QAction(self.widget)
+-        self.btn_end_editing.setObjectName("btn_end_editing")
+-        icon_btn_save = QtGui.QIcon()
+-        icon_btn_save.addPixmap(
+-            QtGui.QPixmap(f"{icon_path}gtk-apply.png"))
+-        self.btn_end_editing.setIcon(icon_btn_save)
+-        self.toolBar.addAction(self.btn_end_editing)
++        # self.btn_end_editing = QtWidgets.QAction(self.widget)
++        # self.btn_end_editing.setObjectName("btn_end_editing")
++        # icon_btn_save = QtGui.QIcon()
++        # icon_btn_save.addPixmap(
++        #     QtGui.QPixmap(f"{icon_path}gtk-apply.png"))
++        # self.btn_end_editing.setIcon(icon_btn_save)
++        # self.toolBar.addAction(self.btn_end_editing)
+ 
+         self.verticalLayout.addLayout(self.horizontalLayout)
+         self.tableView = QtWidgets.QTableView(self.widget)
+@@ -314,30 +314,23 @@ class Ui_MainWindow(object):
+     def retranslateUi(self, MainWindow_reach):
+         _translate = QtCore.QCoreApplication.translate
+         self.name_reach = ""
+-        self.mainwindow_title = _translate("MainWindow_reach", "Jeu de sections du Bief")
++        self.mainwindow_title = _translate("MainWindow_reach", "Geometry reach edit")
+         MainWindow_reach.setWindowTitle(self.mainwindow_title + f"{self.name_reach}")
+-        self.text_section_order = _translate("MainWindow_reach", "\nOrdre des sections : "
+-                                                                 "Amont --> Aval")
+         # self.label_section_order.setText(f"{self.text_section_order}")
+         self.vertical_slider_value = ""
+-        self.vertical_slider_label.setText(_translate("MainWindow_reach", "Pk = ") + f"{self.vertical_slider_value}")
+-        self.btn_open.setToolTip(_translate("MainWindow_reach", " Importer une géométrie"))
+-        self.btn_add.setToolTip(_translate("MainWindow_reach", " Nouveau profil"))
+-        self.btn_delete.setToolTip(_translate("MainWindow_reach", " Supprimer le profil sélectionné"))
+-        self.btn_edit.setToolTip(_translate("MainWindow_reach", " Éditer le profil sélectionné"))
+-        self.btn_copy_selected_profile.setToolTip(_translate("MainWindow_reach", " Copier le profil sélectionné"))
+-        self.btn_paste_selected_profile.setToolTip(_translate("MainWindow_reach", "Coller le profil en fin de liste ("
+-                                                                                  "penser à modifier le Pk avant de "
+-                                                                                  "trier)"))
+-        self.btn_duplicate_selected_profile.setToolTip(
+-            _translate("MainWindow_reach", " Dupliquer la section sélectionnée"))
++        self.vertical_slider_label.setText(_translate("MainWindow_reach", "Kp = ") + f"{self.vertical_slider_value}")
++        self.btn_open.setToolTip(_translate("MainWindow_reach", " Import geometry"))
++        self.btn_add.setToolTip(_translate("MainWindow_reach", " New profile"))
++        self.btn_delete.setToolTip(_translate("MainWindow_reach", " Delete selected profile"))
++        self.btn_edit.setToolTip(_translate("MainWindow_reach", " Editer selected profile"))
++
+         self.btn_sort_asc.setToolTip(_translate("MainWindow_reach", " Trier les profils par ordre croissant des Pk"))
+         self.btn_sort_desc.setToolTip(_translate("MainWindow_reach", " Trier les profils par ordre décroissant des Pk"))
+         self.btn_move_up.setToolTip(_translate("MainWindow_reach", " Changer l'ordre des profils (en décalant le "
+                                                                    "profil sélectionné vers le haut)"))
+         self.btn_move_down.setToolTip(_translate("MainWindow_reach", " Changer l'ordre des profils (en décalant le "
+                                                                      "profil sélectionné vers le bas)"))
+-        self.btn_end_editing.setToolTip(_translate("MainWindow_reach", " Terminer l'édition"))
++        # self.btn_end_editing.setToolTip(_translate("MainWindow_reach", " Terminer l'édition"))
+         self.tableView_header = [_translate("MainWindow_reach", "Name"),
+                                  _translate("MainWindow_reach", "Kp (m)"),
+                                  _translate("MainWindow_reach", "Type")]
+diff --git a/src/View/InitialConditions/PlotDKP.py b/src/View/InitialConditions/PlotDKP.py
+index 58a8321..445cc08 100644
+--- a/src/View/InitialConditions/PlotDKP.py
++++ b/src/View/InitialConditions/PlotDKP.py
+@@ -71,6 +71,8 @@ class PlotDKP(APlot):
+                 color='b', marker='+', lw=1.
+             )
+ 
++        self.canvas.axes.invert_xaxis()
++
+         self.canvas.figure.tight_layout()
+         self.canvas.figure.canvas.draw_idle()
+         if self.toolbar is not None:
+diff --git a/src/View/InitialConditions/PlotDischarge.py b/src/View/InitialConditions/PlotDischarge.py
+index 1fe05de..0703c59 100644
+--- a/src/View/InitialConditions/PlotDischarge.py
++++ b/src/View/InitialConditions/PlotDischarge.py
+@@ -65,6 +65,8 @@ class PlotDischarge(APlot):
+                 color='r', lw=1.
+             )
+ 
++        self.canvas.axes.invert_xaxis()
++
+         self.canvas.figure.tight_layout()
+         self.canvas.figure.canvas.draw_idle()
+         if self.toolbar is not None:
+diff --git a/src/View/InitialConditions/translate.py b/src/View/InitialConditions/translate.py
+index 370d4a2..3fe45f7 100644
+--- a/src/View/InitialConditions/translate.py
++++ b/src/View/InitialConditions/translate.py
+@@ -21,11 +21,11 @@ from PyQt5.QtCore import QCoreApplication
+ _translate = QCoreApplication.translate
+ 
+ table_headers = {
+-    "name": _translate("LateralContribution", "Name"),
++    # "name": _translate("LateralContribution", "Name"),
+     "kp": _translate("LateralContribution", "KP (m)"),
+     # "speed": _translate("LateralContribution", "Speed (m/s)"),
+     "discharge": _translate("LateralContribution", "Discharge (m³/s)"),
+     "elevation": _translate("LateralContribution", "Elevation (m)"),
+     "height": _translate("LateralContribution", "Height (m)"),
+-    "comment": _translate("LateralContribution", "Comment"),
++    # "comment": _translate("LateralContribution", "Comment"),
+ }
+diff --git a/src/View/SedimentLayers/Edit/translate.py b/src/View/SedimentLayers/Edit/translate.py
+index fe9e44a..f9f02cb 100644
+--- a/src/View/SedimentLayers/Edit/translate.py
++++ b/src/View/SedimentLayers/Edit/translate.py
+@@ -6,7 +6,7 @@ _translate = QCoreApplication.translate
+ 
+ table_headers = {
+     "name": _translate("SedimentLayers", "Name"),
+-    "type": _translate("SedimentLayers", "Type"),
++    # "type": _translate("SedimentLayers", "Type"),
+     "height": _translate("Sedimentlayers", "Height"),
+     "d50": _translate("Sedimentlayers", "D50"),
+     "sigma": _translate("Sedimentlayers", "Sigma"),
+diff --git a/src/View/SedimentLayers/Reach/Plot.py b/src/View/SedimentLayers/Reach/Plot.py
+index a312a8c..5764abf 100644
+--- a/src/View/SedimentLayers/Reach/Plot.py
++++ b/src/View/SedimentLayers/Reach/Plot.py
+@@ -78,6 +78,8 @@ class Plot(APlot):
+                 color='grey' if i == len(z_sl) - 1 else None
+             )
+ 
++        self.canvas.axes.invert_xaxis()
++
+         self.canvas.figure.tight_layout()
+         self.canvas.figure.canvas.draw_idle()
+         if self.toolbar is not None:
+diff --git a/src/lang/fr.ts b/src/lang/fr.ts
+index 8034809..031f349 100644
+--- a/src/lang/fr.ts
++++ b/src/lang/fr.ts
+@@ -60,21 +60,6 @@
+         <source>Ponctual contribution</source>
+         <translation>Contributions ponctuelles</translation>
+     </message>
+-    <message>
+-        <location filename="../View/BoundaryCondition/translate.py" line="32"/>
+-        <source>Time over Z</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message>
+-        <location filename="../View/BoundaryCondition/translate.py" line="33"/>
+-        <source>Time over Discharge</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message>
+-        <location filename="../View/BoundaryCondition/translate.py" line="34"/>
+-        <source>Z over Discharge</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+     <message>
+         <location filename="../View/BoundaryCondition/translate.py" line="39"/>
+         <source>Name</source>
+@@ -110,6 +95,21 @@
+         <source>Solid</source>
+         <translation type="unfinished">Solide</translation>
+     </message>
++    <message>
++        <location filename="../View/BoundaryCondition/translate.py" line="32"/>
++        <source>TimeZ</source>
++        <translation type="unfinished"></translation>
++    </message>
++    <message>
++        <location filename="../View/BoundaryCondition/translate.py" line="33"/>
++        <source>TimeDischarge</source>
++        <translation type="unfinished"></translation>
++    </message>
++    <message>
++        <location filename="../View/BoundaryCondition/translate.py" line="34"/>
++        <source>ZDischarge</source>
++        <translation type="unfinished"></translation>
++    </message>
+ </context>
+ <context>
+     <name>Checker</name>
+@@ -533,7 +533,7 @@
+     <message>
+         <location filename="../View/Frictions/translate.py" line="24"/>
+         <source>Name</source>
+-        <translation>Nom</translation>
++        <translation type="obsolete">Nom</translation>
+     </message>
+     <message>
+         <location filename="../View/Frictions/translate.py" line="26"/>
+@@ -548,12 +548,17 @@
+     <message>
+         <location filename="../View/Frictions/translate.py" line="28"/>
+         <source>Begin strickler</source>
+-        <translation>strickler de départ</translation>
++        <translation type="obsolete">strickler de départ</translation>
+     </message>
+     <message>
+         <location filename="../View/Frictions/translate.py" line="29"/>
+         <source>End strickler</source>
+-        <translation>Strickler de fin</translation>
++        <translation type="obsolete">Strickler de fin</translation>
++    </message>
++    <message>
++        <location filename="../View/Frictions/translate.py" line="28"/>
++        <source>Strickler</source>
++        <translation type="unfinished"></translation>
+     </message>
+ </context>
+ <context>
+@@ -587,7 +592,7 @@
+ <context>
+     <name>LateralContribution</name>
+     <message>
+-        <location filename="../View/InitialConditions/translate.py" line="24"/>
++        <location filename="../View/LateralContribution/translate.py" line="35"/>
+         <source>Name</source>
+         <translation>Nom</translation>
+     </message>
+@@ -602,7 +607,7 @@
+         <translation>Lit moyen</translation>
+     </message>
+     <message>
+-        <location filename="../View/InitialConditions/translate.py" line="30"/>
++        <location filename="../View/Stricklers/translate.py" line="31"/>
+         <source>Comment</source>
+         <translation>Commentaire</translation>
+     </message>
+@@ -1715,7 +1720,7 @@ sur le plan défini par les deux points nommés extrêmes </source>
+         <translation type="unfinished"></translation>
+     </message>
+ </context>
+-<context>
++<context encoding="UTF-8">
+     <name>MainWindow_reach</name>
+     <message>
+         <location filename="../View/Frictions/PlotStricklers.py" line="76"/>
+@@ -1723,7 +1728,7 @@ sur le plan défini par les deux points nommés extrêmes </source>
+         <translation></translation>
+     </message>
+     <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="342"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="335"/>
+         <source>Kp (m)</source>
+         <translation>Pk (m)</translation>
+     </message>
+@@ -1778,68 +1783,47 @@ sur le plan défini par les deux points nommés extrêmes </source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="317"/>
+-        <source>Jeu de sections du Bief</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="319"/>
+-        <source>
+-Ordre des sections : Amont --&gt; Aval</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="323"/>
+-        <source>Pk = </source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="325"/>
+-        <source> Nouveau profil</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="334"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="327"/>
+         <source> Trier les profils par ordre croissant des Pk</source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="341"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="334"/>
+         <source>Name</source>
+         <translation>Nom</translation>
+     </message>
+     <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="343"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="336"/>
+         <source>Type</source>
+         <translation>Type</translation>
+     </message>
+     <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="345"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="338"/>
+         <source>Alt+Z</source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="348"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="341"/>
+         <source>Alt+E</source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="351"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="344"/>
+         <source>Alt+R</source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="355"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="348"/>
+         <source>Vue globale automatique (Alt+S)</source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="359"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="352"/>
+         <source>Vue globale automatique (Alt+D)</source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message>
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="363"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="356"/>
+         <source>Vue globale automatique (Alt+F)</source>
+         <translation type="unfinished"></translation>
+     </message>
+@@ -1889,68 +1873,33 @@ Seuls les 5 premiers seront édités.</source>
+         <source>Profil sélectionné</source>
+         <translation type="unfinished"></translation>
+     </message>
+-    <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="324"/>
+-        <source> Importer une géométrie</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="326"/>
+-        <source> Supprimer le profil sélectionné</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="327"/>
+-        <source> Éditer le profil sélectionné</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+     <message encoding="UTF-8">
+         <location filename="../View/Geometry/mainwindow_ui_reach.py" line="328"/>
+-        <source> Copier le profil sélectionné</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="329"/>
+-        <source>Coller le profil en fin de liste (penser à modifier le Pk avant de trier)</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="332"/>
+-        <source> Dupliquer la section sélectionnée</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="335"/>
+         <source> Trier les profils par ordre décroissant des Pk</source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="336"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="329"/>
+         <source> Changer l&apos;ordre des profils (en décalant le profil sélectionné vers le haut)</source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="338"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="331"/>
+         <source> Changer l&apos;ordre des profils (en décalant le profil sélectionné vers le bas)</source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="340"/>
+-        <source> Terminer l&apos;édition</source>
+-        <translation type="unfinished"></translation>
+-    </message>
+-    <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="346"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="339"/>
+         <source>Vue isométrique (Alt+Z)</source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="349"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="342"/>
+         <source>Vue isométrique (Alt+E)</source>
+         <translation type="unfinished"></translation>
+     </message>
+     <message encoding="UTF-8">
+-        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="352"/>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="345"/>
+         <source>Vue isométrique (Alt+R)</source>
+         <translation type="unfinished"></translation>
+     </message>
+@@ -1960,7 +1909,7 @@ Seuls les 5 premiers seront édités.</source>
+         <translation type="unfinished">Altitude (m)</translation>
+     </message>
+     <message>
+-        <location filename="../View/Geometry/PlotKPZ.py" line="66"/>
++        <location filename="../View/Geometry/PlotAC.py" line="73"/>
+         <source>Height (m)</source>
+         <translation type="unfinished">Hauteur (m)</translation>
+     </message>
+@@ -1969,6 +1918,56 @@ Seuls les 5 premiers seront édités.</source>
+         <source>Select destination file</source>
+         <translation type="unfinished"></translation>
+     </message>
++    <message>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="317"/>
++        <source>Geometry reach edit</source>
++        <translation type="unfinished"></translation>
++    </message>
++    <message>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="321"/>
++        <source>Kp = </source>
++        <translation type="unfinished"></translation>
++    </message>
++    <message>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="322"/>
++        <source> Import geometry</source>
++        <translation type="unfinished"></translation>
++    </message>
++    <message>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="323"/>
++        <source> New profile</source>
++        <translation type="unfinished"></translation>
++    </message>
++    <message>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="324"/>
++        <source> Delete selected profile</source>
++        <translation type="unfinished"></translation>
++    </message>
++    <message>
++        <location filename="../View/Geometry/mainwindow_ui_reach.py" line="325"/>
++        <source> Editer selected profile</source>
++        <translation type="unfinished"></translation>
++    </message>
++    <message>
++        <location filename="../View/Geometry/PlotAC.py" line="69"/>
++        <source>Transverse abscissa (m)</source>
++        <translation type="unfinished"></translation>
++    </message>
++    <message>
++        <location filename="../View/Geometry/PlotAC.py" line="79"/>
++        <source>Previous cross-section</source>
++        <translation type="unfinished"></translation>
++    </message>
++    <message>
++        <location filename="../View/Geometry/PlotAC.py" line="80"/>
++        <source>Cross-section</source>
++        <translation type="unfinished"></translation>
++    </message>
++    <message>
++        <location filename="../View/Geometry/PlotAC.py" line="81"/>
++        <source>Next cross-section</source>
++        <translation type="unfinished"></translation>
++    </message>
+ </context>
+ <context>
+     <name>Network</name>
+@@ -2071,7 +2070,7 @@ Seuls les 5 premiers seront édités.</source>
+     <message>
+         <location filename="../View/SedimentLayers/Edit/translate.py" line="9"/>
+         <source>Type</source>
+-        <translation type="unfinished">Type</translation>
++        <translation type="obsolete">Type</translation>
+     </message>
+     <message>
+         <location filename="../View/SedimentLayers/Reach/Profile/Window.py" line="68"/>

doc/paper/2023/simhydro/paper.txt

@@ -0,0 +1,131 @@
+Pamhyr2: a graphical user interface for 1D hydro-sedimentary modelling of rivers 
+Pierre-Antoine Rouby1, Benoît Camenen, Lionel Pénard, Léa Kieffer, Théophile Terraz1
+INRAE, UR RiverLy, centre de Lyon-Grenoble, 5 rue de la Doua, CS20244, 69625, Villeurbanne, France.
+pierre-antoine.rouby@inrae.fr, benoit.camenen@inrae.fr, lionel.penard@inrae.fr, 
+lea.kieffer@inrae.fr, theophile.terraz@inrae.fr
+KEY WORDS
+1D modelling, sedimentary modelling, user interface
+Abstract
+In numerical river simulation, one-dimensional models are commonly used to study water level and discharge for large domains or long time series. These models are less time-consuming than two- and three-dimensional numerical models, and require fewer input parameters and allow ensemble runs. To build a one-dimensional hydraulic model, a pre- and post-processing tool is needed for creating reach geometry, specifying initial and boundary conditions, friction coefficients and other numerical parameters. Such a tool needs to ensure the consistency of the model and provide a user-friendly graphical user interface. In this article, we present Pamhyr2, the fully rebuilt version of the PAMHYR modelling platform. It is developed using Python, PyQt and Matplotlib. Pamhyr2 is free and open-source, multilingual, cross-platform (Linux, Windows) and is generic enough to accept various one-dimensional solvers. Pamhyr2 includes and enhances all the features from the previous version: multiple-reach modelling, geometry definition from cross-section data, initial and boundary conditions, friction coefficients, hydraulic structures, lateral inflow, punctual intake and the results visualization window. In addition, this version includes new features such as pollutants modelling, bed-load and bed evolution. We describe several windows: the creation of sedimentary layers in the river bed, the sediment characteristics for each layer, the sediment, and pollutant boundary conditions and the lateral inputs. These functionalities are illustrated with simple examples. We finally show the visualization windows for bed evolution. 
+1.	introduction
+Numerous tools are available to simulate river flow. Among them, one-dimensional (1D) numerical hydraulic models have the advantage of being accurate while being inexpensive in terms of computation time, making it possible to simulate long rivers over very long timescales. A lot of 1D numerical solvers for river hydraulics are available, including HEC-RAS2, Mike 113 and Mascaret4 for example. However, in order to investigate specific research issues it is essential to have the ability to edit the source code of numeric solvers and to have an ergonomic Graphical User Interace (GUI). Most of the widespread 1D modeling tools are either written with a proprietary source code (such as HEC-RAS and Mike 11) or are research tools with irregularly maintained GUI. On its side, INRAE has been developing two 1D solvers for over 40 years: Mage5 [5]  and RubarBE6. In parallel, a GUI named PAMHyR has been developed to enable modelers to build river models, parameterize calculations, launch Mage and RubarBE solvers and post-process the results [3].
+
+	PAMHyR, developed in Java since the end of 90s (originally named X-NT), it is written with the following constraints in mind: 
+    • Multi-platform Windows, Unix-like,
+    • no dependence on an external graphics library with a paid license for user, 
+    • generate as much code as possible from a UML diagram.
+Now, PAMHyR suffers from complexity due to highly specific interface, as well as a significant technological debt [4]: old external library no longer maintained, over dated project architecture, and numerous modifications by multiple developers (and internship) over short periods of time, with no vision of long-term maintainability. What's more, the software is only in French and the code is closed (the sources are not available), which is a hindrance to the notoriety of the software as well as to the vitality of development. In addition, no significant changes have been made in the application since 2016, even if PAMHyR is still often used for river modeling [2].
+INRAE needs a new tool to keep up with today’s technologies and with the development of its numerical solvers. In this paper we introduce the current work on the new graphical user interface Pamhyr2.
+2.	Objective
+2.1	History and need 
+PAMHyR is available for two solvers, with two different types of studies: Mage and RubarBE. These two different kinds of studies have mostly similar interfaces, but have their own specific features. They can both be used to: define reach geometry, boundary conditions, lateral inflows, overflows, friction, and structures. In addition, you can run the solver on the study you are editing, and view the solver's results. Firstly, although both Mage and RubarBE are multiple-reach compatible, only the Mage interface enables multiple-reach study editing. Secondly, many small details change from one interface to the next. In addition, the latest Mage features, for example, are not available for modeling in the interface.
+
+	A total overhaul of the Java code would be too much work, given the technological debt accumulated since the beginning of its development. Given the workload involved in modernizing PAMHyR, the decision is made to develop a new version of the software from scratch. This version should be equivalent to the Java version in terms of functionality, and add new features, notably the new Mage functionalities. The aim is also to have an interface that is more modern, more user-friendly, more flexible, available on multiple platforms and in different languages. The interface must also have a certain degree of generality and be able to adapt to different solvers without major changes to the interface. It is also important that the code is open to contributions, to keep the project as lively as possible.
+2.2	Technical choices
+From a technical point of view, choices were made on the basis of the above-mentioned objectives (modernity, maintainability, multiplatform, multilingual, user-friendliness, flexibility), as well as: the ease of learning the technology to facilitate contribution to the software, the popularity of the technology in the scientific community and the interoperability of this technology.
+
+	We chose Python7 as our programming language. It is a popular language in the scientific world, but also in many other fields, with lot of very active software libraries. Moreover, it is a simple and interpreted language, often recommended for beginners, and practiced by beginners and experienced developers. In addition, there is the PyQT8 software library that lets you use the Qt9 GUI library, which is ideal for this project. It is very popular and active, and available on many platforms (GNU/Linux, Windows, MacOS, Unix, BSD, Android, iOS). Qt also includes a number of tools for translating interfaces into different languages. Here we use the community version of Qt, with no paid features. For graphics and visualization, the choice fell on Matplotlib10, which is very popular in the scientific community and can be used to draw a wide variety of graphics. For saving studies, the choice was made to use SQLite11, unlike the previous version which used a text file formatted in XML. The advantages of a SQLite database, in addition to being supplied with Python, include ease of use, interoperability, fast, formal structuring by table and consistency checking between tables. We have chosen to version the database to allow it to evolve with the software, and to make it possible to add columns to tables or make other minor modifications to the database without breaking the compatibility of studies from previous versions. However, PAMHyR and Pamhyr2 study files will not be compatible anymore, and previous studies have to be recreated for this version.
+
+	This new version Pamhyr2 remains closely linked to the Mage solver with which it is supplied. But the ambition of flexibility and generalization gives us the opportunity to develop an interface that can be adapted to other solvers. There can be two solutions for using a solver with Pamhyr2: either the solver adapts to Pamhyr2, or Pamhyr2 adapts to the solver. So Pamhyr2 must integrate a generic output containing the study data that can be used by other solvers, which will be launched via an external executable, and then a generic result input that must be supplied by the solver. In addition, Pamhyr2 must be able to handle the execution and reading of results for specific solvers, and the addition of a solver to Pamhyr2 must be reasonably easy for the solver developer. We have opted for the free and open source software (FOSS) [6] license GPLv312, which guarantees freedom to access, study, modify, redistribute and share. In addition, it guarantees that Pamhyr2's code will remain free and is not allowed to be added any closed-source modifications. However, it will still be possible to use closed-source solvers, because solver is not a part of Pamhyr2.
+3.	DEVELOPMENTS
+3.1	New version of PAMHYR 
+Pamhyr2 currently features the following interfaces: Editing of river network, geometry, boundary conditions, lateral contribution, friction, initial conditions, solver numerical parameters, as well as the interface for launching a solver and displaying hydraulic and sediment results.
+3.1.1	River network
+The river network corresponds to the topology of the hydraulic network of the river. River network editing (Figure 1) is performed via a window featuring a graphical network area and two tables, one with a list of nodes (zone 3) and the other with a list of reaches (zone 2). The graphic area (zone 1) lets you add, delete and move nodes. It can also be used to link nodes in the network by arrows that represent reaches. The color of the node in the graphic area depends on whether it is a downstream, upstream or internal network node. If there is more than one reach, you can use this interface to select a reach, which will then be selected for editing in the following reach-dependent windows. It is also possible to deactivate reaches, in which case they will not be used for the solver run. Tables can be used to edit additional information on nodes and reaches, such as giving them names.
+
+
+
+Figure 1: River network window for a dummy river with 8 reaches, 2 upstream node (yellow), 1 downstream node (green) and 5 internal node (blue).
+
+3.1.2	Geometry
+Geometry editing (Figure 2) is linked to a reach; the reach used will be the last reach selected in the river network window. This window contains a table with the list of the reach’s cross-sections (zone 1) and three graphics. A top view based on the XY positions of the cross-sections’ points and guidelines (zone 2). A longitudinal view of all cross-sections with the height of the guidelines’ points for each profile (zone 3). Finally, a cross-sectional view is displayed (zone 4)  showing the selected cross-section with the previous and the following one in dotted line. Each cross-section is defined by a Kilometer Point (KP) and an optional name. The KP represents the longitudinal coordinate of the cross-section in the numerical solver. It is possible to add cross-sections and then select them to access an editing window, or to import cross-sections from file. Cross-section editing window takes the form of a table showing the raw data of the points that make up the cross-section, and a graph showing the current cross-section.
+
+
+
+Figure 2: Geometry window for Saar river on the unique reach of the Saar river study.
+
+3.1.3	Boundary conditions
+Boundary conditions (BC) are associated with a network node. There are three different types: hydrograph, water level chart or rating curve, in three different categories: Liquid, solid and suspension. Boundary conditions can be added, deleted and edited. Editing a BC (Figure 3) takes the form of a two-column table (the type of column varies according to the type of BC) and a graph representing the data. Data can be copied and pasted from a spreadsheet, or edited by hand.
+
+
+
+Figure 3: Edition window of the boundary conditions hydrograph for the Saar river upstream node.
+
+
+3.1.3	Lateral contribution
+The user interface for lateral contributions is similar to that for boundary conditions, except that it applies to a part of a reach between two KP, and not to a node.
+3.1.4	Friction
+Friction editing is linked to a reach, a friction coefficient is defined and then applied to a portion of the reach between two KP (Figure 4). Friction is defined by two Strickler coefficients, one for the river's main channel and one for the flood plain. The list of available Strickler coefficients is defined at the study level.
+
+
+
+Figure 4: The reach portions for friction coefficient edition window for the Saar study.
+
+3.1.5	Initial conditions
+The initial conditions (Figure 5) are also linked to the reach, corresponding to the water level and flow rate at the initial time of the simulation. These two values can be defined separately for each cross-section of the reach. It is also possible to use predefined functions to generate a water line and flow rate.
+
+
+
+Figure 5:  Initial conditions edition window for the Saar study and the Saar reach. Initial conditions initialized with minimal height of 3 meter and a constant flow rate of 100 m3/s.
+
+3.1.6	Running a solver
+Solver numerical parameters are defined using a key/value table specific to each type of solver known to Pamhyr2. You can select and run the desired solver from the list of configured solvers. Then a window appears with a text area where you can follow the progress of file export and solver output in Figure 6 (zone 3). You can also pause, stop or restart the solver using buttons (zone 1). Once the simulation has finished, you can open a solver trace file, and if the simulation has finished correctly, you can open the results using the “results” button (zone 2).
+
+
+
+Figure 6:  Solver log window for “Elargissement” tests case with a mage solver named “src”.
+
+3.1.7	Results processing and visualization
+The post processing window is made up of tables allowing you to select the reaches, profiles, and timestamps for visualizing the results. The results are presented in different tabs according to the type of data displayed. The Figure 7 shows hydraulic data with water level at a given reach, profile, and time. The graphs shown here are the same as those in the geometry window, with the addition of water level results at a given time.
+
+
+
+Figure 7: Hydraulics results window for “Elargissement” tests case with a Mage solver at timestamps 2113200.0 sec for reach R1 and cross-section at KP 427.7778.
+
+3.2	New features
+Some new features are already implemented in Pamhyr2. First of all, the software is multilingual with English and French currently available. An other very useful enhancement, is the possibility to use the “Ctrl+Z” or “Ctrl+Shift+Z” shortcuts to undo or redo the user actions.
+
+	Another novelty lies in the initial conditions window. Previously the water line could be generated by constant height or altitude. It is now possible to give a minimum height with an increasing function from downstream to upstream, thus avoiding physically unrealistic water holes and getting a little closer to a real and physically coherent water line (Figure 5).
+
+	An important new feature in this version is the presence of sediment layers. They can be defined as shown in Figure 8. Sediment layers can be applied separately for each cross-section (Figure 9). As with friction, they must first be defined and then applied to a profile. They can also be applied separately to each point of a cross-section. Layers are defined by a thickness, a median diameter (D50), a grain size sorting coefficient (Sigma) and a critical shear stress as described in Pierre Balayn PhD thesis [1]. The visualization of sedimentary results is currently done using a plot of the river profile, which shows the evolution of the bed bottom and the thickness of layers as a function of time (Figure 10).
+
+
+
+Figure 8: Sediment layers definition window for Saar study. This sediment layers has three layers (L1, L2 and F).
+
+
+
+
+Figure 9: Sediment layers defined in Figure 8 applied to the reach Saar.
+
+
+
+
+Figure 10: Sedimentary results window for “Elargissement” tests case at timestamps 2113200.0 sec for reach R1 and cross-section at KP 427.7778. This example has only two sediment layers.
+
+4.	Future work
+Old PAMHyR functionalities and additional functionalities, allowed the definition of structures, traps, suspended solids or dissolved pollutants in water, have yet to be developed. Additional options for visualizing results also need to be added. One possibility would be a graph generation tool that would allow the user to select x-axis and y-axis values, as well as the time range, to create a customized graph. We also planned an advanced feature with code-editing window to enable the user to define his own Matplotlib drawing function, which would be saved in the study. Presentation work is also planned to harmonize the application's icons and visuals (images, graphics, etc.). The aim is to make the application clearer and more user-friendly.
+
+	Supporting new solvers and generic solvers will also be an important part of the remaining work. To use any solver, we need to define a generic output format for Pamhyr2 that is simple enough to be given as solver input without major code modification, or that can be processed by a simple Python script (for example). The same applies to the results file that Pamhyr2 will have to read. The main candidates for these generic formats are :
+    • SQLite, because it is easy to use and the library is compatible with most programming languages,
+    • HDF513, a hierarchical file format widely used in scientific circles and compatible with many programming languages.
+
+	As the project progresses, it becomes more and more complicated to test all the functionalities by hand. A suite of tests and unit tests is needed to prevent regression bugs and ensure that the code remains maintainable. For the moment, no such test has been defined. Unit tests will have to be developed on the model and each sub-module will have to include tests, but higher-level functionalities will also have to be tested, such as the recording and playback of a complete study. It will be difficult to test the graphical part automatically, so tests will continue to be carried out manually.
+
+	Full documentation must also be provided to give users an overview of the features offered by Pamhyr2. This documentation must include descriptions of the configuration options, as well as explanations and examples of how to use each feature. In addition, this documentation must include a section for contributors and developers, with a technical presentation of the application, as well as a presentation of the translation, development and debugging tools available.
+
+	For translations into other languages than English and French, we hope that external contributors will take on this work. But no languages other are currently planned by the development team. The application packages are currently available for GNU/Linux and Windows. Under GNU/Linux, a “.tar.xz” archive is available for download. For Windows, this is an installer for application and a compatible version of Mage. An installer for GNU/Linux systems is a work in progress. There are many choices for creating an installer, but their scope depends on the distribution for which they are intended. “.deb”, for example, is intended for Debian distributions and their derivatives (Ubuntu, Linux Mint, and so on). We do not exclude other, more generic options.
+5.	Conclusion
+After a great deal of preparatory work, experimentation, technical selection, proof-of-concept and several months of development, Pamhyr2 is currently in version “v0.0.0”14 and contains almost 20,000 lines of Python code. It allows to create a minimal 1D hydro-sedimentary study, run a Mage simulation and visualize some of the results obtained, all with a more user-friendly interface. But the software is still unstable and under development, and many modifications and improvements are still to be made before it can really be used by everyone. Furthermore, even though many of the features of the original version of PAMHyR are already available (see 3.1 page 3), there are still many to be added to make Pamhyr2 a complete, useful and fully functional piece of software.
+ACKNOWLEDGEMENTS
+The results presented here partially come from the Pamhyr PITI project funded by the Carnot Eau & Environnement institute and the French National Research Agency. Carnot Eau & Environnement was accredited in March 2020 by the French Ministry of Research and Higher Education on the proposal of the Agence Nationale de la Recherche. The signatory members are INRAE, CNRS, Aix Marseille University, IFTS, INSA Lyon, IRD, Université de Montpellier and Université de Savoie Mont-Blanc.
+
+The authors thank Sylvain Coulibaly for his preparatory work and his technical choices.
+REFERENCES and Citations
+[1]	Balayn, P. (2001). Contribution à la modélisation numérique de l'évolution morphologique des cours d'eau aménagés lors de crues. Université Claude Bernard - Lyon 1
+[2]	Camenen, B., Faure, J. Décanis, S. Diéval, L. (2021) A 1D numerical tool for real time modelling of a complex river network, SimHydro 2021
+[3]	Faure, J. (2013). An Integrated Simulation platform - PAMHyR. In John Wiley & Sons, Inc. eBooks (pp. 245–252). https://doi.org/10.1002/9781118557891.ch21
+[4]	Faure, J. (2017). Développement d’Adis-TS. OSR4 | Action V.1. https://hal.science/hal-03749339
+[5]	Faure, J-B., (2019). Mage, Résolution des équations de Barré de St Venant 1D en réseaux complexes, Documentation théorique et mode d’emploi. Techreport (in French)
+[6]	Walt, S. 2007. Free/open source software development. In Proceedings of the the 6th joint meeting of the European software engineering conference and the ACM SIGSOFT symposium on The foundations of software engineering (ESEC-FSE '07). Association for Computing Machinery, New York, NY, USA, 459–468. https://doi.org/10.1145/1287624.1287689
+

doc/tools/PamhyrDoc.cls

@@ -0,0 +1,163 @@
+\ProvidesClass{Pamhyr}[Pamhyr class]
+\NeedsTeXFormat{LaTeX2e}
+\ProcessOptions\relax
+\LoadClass[11pt,a4paper]{article}
+
+\usepackage[hmargin=2cm, vmargin=3cm]{geometry}
+\usepackage[british,UKenglish,USenglish,english,american]{babel}
+\usepackage[T1]{fontenc}
+\usepackage{times}
+\usepackage[usenames,dvipsnames]{xcolor}
+\usepackage[absolute]{textpos}
+\usepackage{titlesec}
+\usepackage{amssymb}
+\usepackage{tcolorbox}
+
+\usepackage{enumitem}
+\setlist[enumerate]{itemsep=-0.5em}
+\setlist[itemize]{itemsep=-0.5em}
+
+%%
+%% Bibtex
+%%
+
+\bibliographystyle{plain}
+
+%%
+%% Colors
+%%
+
+\definecolor{primary}{HTML}{212121}
+\definecolor{lightheadings}{HTML}{5c5c5c}
+\definecolor{headings}{HTML}{212121}
+\definecolor{links}{HTML}{e12121}
+
+%%
+%% Title
+%%
+
+\setlength{\columnsep}{20pt}
+\setlength{\parskip}{0pt}      %% inter paragraph space
+
+\renewcommand{\maketitle}{
+    {
+      \centering
+      {
+        \fontsize{35}{40}
+        \selectfont
+        \scshape
+        Pamhyr2 \\
+        \medskip
+      }
+      {
+        %\centering
+        \scshape
+        \fontsize{20pt}{24pt}
+        \selectfont
+        A graphical user interface for 1D hydro-sedimentary modelling
+        of rivers \\
+        \medskip
+      }
+      \rule{3cm}{0.4pt}\\
+      \vspace{0.3cm}
+      {
+        %\centering
+        \scshape
+        \fontsize{26pt}{24pt}
+        \selectfont
+        \@title \\
+        Document licence: GPLv3 \\
+        \medskip
+        \includegraphics[width=3cm]{../images/GPLv3_Logo.png} \\
+        \medskip
+        \fontsize{14pt}{16pt}
+        \selectfont
+        © \@author \\
+        \medskip
+      }
+      {
+        %\centering
+        %\scshape
+        \fontsize{12pt}{14pt}
+        \selectfont
+        \@date
+        %\medskip
+      }\\
+      \vspace{0.3cm}
+      \rule{3cm}{0.4pt}\\
+      \vspace{0.3cm}
+      {
+        \begin{minipage}[h]{0.24\textwidth}
+          \includegraphics[width=4cm]{../../src/View/ui/ressources/Logo-INRAE.png}
+        \end{minipage}
+        \begin{minipage}[h]{0.24\textwidth}
+          \includegraphics[width=4cm]{../images/riverly.png}
+        \end{minipage}
+      }\\
+      \vspace{0.3cm}
+      %%\rule{\textwidth}{0.4pt}
+      %%\newpage
+    }
+}
+
+
+%%
+%% Texttt
+%%
+
+\usepackage[htt]{hyphenat}      % Make texttt split at end of line
+
+%%
+%% Verbatim
+%%
+
+\AtBeginEnvironment{verbatim}{\footnotesize}
+
+%% Add box arrond of verbatim
+\BeforeBeginEnvironment{verbatim}{\begin{tcolorbox}[boxsep=0pt, left=0.1cm, right=0.1cm, arc=0pt, boxrule=0.5pt, colback=white]}%
+\AfterEndEnvironment{verbatim}{\end{tcolorbox}}%
+
+\makeatletter
+\AtBeginEnvironment{verbatim}{\dontdofcolorbox}
+
+\def\dontdofcolorbox{\renewcommand\fcolorbox[4][]{##4}}
+\makeatother
+
+%%
+%% Source code blocks
+%%
+
+%% Use minted
+\usepackage{minted}
+\usemintedstyle{tango}
+\setminted[c]{fontsize=\footnotesize,encoding=utf8,linenos}
+\setminted[c++]{breaklines,fontsize=\footnotesize,encoding=utf8,linenos}
+\setminted[python]{breaklines,fontsize=\footnotesize,encoding=utf8,linenos}
+\setminted[shell]{breaklines,fontsize=\footnotesize,encoding=utf8}
+\setminted[scheme]{breaklines,fontsize=\footnotesize,encoding=utf8}
+\setminted[commun-lisp]{breaklines,fontsize=\footnotesize,encoding=utf8}
+\setminted[text]{breaklines,fontsize=\footnotesize,encoding=utf8}
+\setminted[llvm]{breaklines,fontsize=\footnotesize,encoding=utf8}
+
+%% Add box arrond of minted
+\BeforeBeginEnvironment{minted}{\begin{tcolorbox}[boxsep=0pt, left=0.1cm, right=0.1cm, arc=0pt, boxrule=0.5pt, colback=white]}%
+\AfterEndEnvironment{minted}{\end{tcolorbox}}%
+
+\makeatletter
+\AtBeginEnvironment{minted}{\dontdofcolorbox}
+
+\def\dontdofcolorbox{\renewcommand\fcolorbox[4][]{##4}}
+\makeatother
+
+
+%%
+%% Water mark
+%%
+
+%% FIXME: Delete me for the first realease
+\usepackage{draftwatermark}
+\SetWatermarkLightness{0.8}
+\SetWatermarkAngle{25}
+\SetWatermarkScale{3}
+\SetWatermarkFontSize{1cm}
+\SetWatermarkText{Work in progress}

doc/tools/latex.org

@@ -0,0 +1,19 @@
+# latex.org -- Pamhyr documentation latex config
+# Copyright (C) 2023  INRAE
+#
+# This program 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.
+#
+# This program 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 this program.  If not, see <https://www.gnu.org/licenses/>.
+
+# -*- coding: utf-8 -*-
+
+#+LaTeX_CLASS: PamhyrDoc

doc/tools/macro.org

@@ -0,0 +1,51 @@
+# macro.org -- Pamhyr documentation macro
+# Copyright (C) 2023  INRAE
+#
+# This program 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.
+#
+# This program 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 this program.  If not, see <https://www.gnu.org/licenses/>.
+
+# -*- coding: utf-8 -*-
+
+# Constant name
+#+MACRO: oldPamhyr     PAMHyR
+#+MACRO: Pamhyr        Pamhyr2
+#+MACRO: Mage          Mage
+#+MACRO: Rubarbe       RubarBE
+#+MACRO: Inrae         INRAE
+#+MACRO: latex         \LaTeX
+
+# Information
+#+MACRO: version       (eval (pamhyr-version))
+#+MACRO: current-file  (eval (pamhyr-current-file))
+
+# Special format
+#+MACRO: file          =$1=
+
+# Biblio
+#+MACRO: cite          [cite:$1]
+#+MACRO: biblio        \bibliography{documentation}
+
+# Wrapper
+#+NAME: attr_wrap
+#+HEADER: :var width="\\textwidth"
+#+HEADER: :var caption=""
+#+HEADER: :var smallcaption=""
+#+HEADER: :var name=""
+#+HEADER: :var data=""
+#+HEADER: :var float="nil"
+#+BEGIN_SRC sh :results output :exports none
+  echo "#+CAPTION[$smallcaption]: $caption"
+  echo "#+NAME: $name"
+  echo "#+ATTR_LATEX: :width $width :float $float"
+  echo "$data"
+#+END_SRC

doc/tools/setup.el

@@ -0,0 +1,66 @@
+(require 'org)
+(require 'subr-x)
+
+;; LaTeX config
+
+(add-to-list
+ 'org-latex-classes
+ '("PamhyrDoc"
+   "\\documentclass{../tools/PamhyrDoc}"
+   ("\\section{%s}" . "\\section*{%s}")
+   ("\\subsection{%s}" . "\\subsection*{%s}")
+   ("\\subsubsection{%s}" . "\\subsubsection*{%s}")
+   ("\\paragraph{%s}" . "\\paragraph*{%s}")
+   ("\\subparagraph{%s}" . "\\subparagraph*{%s}")))
+
+(setq org-confirm-babel-evaluate nil)
+(setq org-latex-caption-above nil)
+
+(org-babel-do-load-languages
+   'org-babel-load-languages
+   (add-to-list 'org-babel-load-languages '(shell . t)))
+
+(org-babel-do-load-languages
+ 'org-babel-load-languages
+ (add-to-list 'org-babel-load-languages '(dot . t)))
+
+(org-babel-do-load-languages
+ 'org-babel-load-languages
+ (add-to-list 'org-babel-load-languages '(python . t)))
+
+(add-to-list 'org-latex-packages-alist '("" "minted"))
+(setq org-latex-listings 'minted)
+(setq org-src-fontify-natively t)
+
+(setq org-latex-pdf-process
+      '("pdflatex -shell-escape -interaction nonstopmode -output-directory %o %f"
+        "bibtex documentation"
+        "pdflatex -shell-escape -interaction nonstopmode -output-directory %o %f"
+        "pdflatex -shell-escape -interaction nonstopmode -output-directory %o %f"))
+
+;; Citations
+
+(org-add-link-type
+ "cite" 'ebib
+ (lambda (path desc format)
+   (cond
+    ((eq format 'html)
+     (format "(<cite>%s</cite>)" path))
+    ((eq format 'latex)
+     (if (or (not desc) (equal 0 (search "cite:" desc)))
+         (format "\\cite{%s}" path)
+       (format "\\cite[%s][%s]{%s}"
+               (cadr (split-string desc ";"))
+               (car (split-string desc ";"))  path))))))
+
+;; Functools
+
+(defun pamhyr-version ()
+  "Return the contents of the pamhyr version file."
+  (with-temp-buffer
+    (insert-file-contents "../../VERSION")
+    (buffer-string)))
+
+(defun pamhyr-current-file ()
+  (let ((lst (split-string (buffer-file-name) "/")))
+    (string-join (nthcdr (- (length lst) 3) lst) "/")))

doc/users/build.sh

@@ -0,0 +1,10 @@
+#! /bin/sh
+
+# Build the PDF
+
+emacs -Q -q -nw --batch                                 \
+      --eval "(require 'org)"                           \
+      --eval "(require 'ox-latex)"                      \
+      --eval "(load-file \"../tools/setup.el\")"        \
+      --file ./documentation.org                        \
+      -f org-latex-export-to-pdf

doc/users/documentation.bib

@@ -0,0 +1 @@
+../tools/ref.bib

doc/users/documentation.org

@@ -0,0 +1,20 @@
+#+STARTUP: indent
+
+#+INCLUDE: ../tools/macro.org
+#+INCLUDE: ../tools/latex.org
+
+#+TITLE: User documentation
+#+SUBTITLE: Version: {{{version}}}
+#+AUTHOR: {{{INRAE}}}
+
+#+OPTIONS: toc:t
+#+LANGUAGE: UKenglish
+
+* TODO Over view
+* TODO Get started
+** Installation
+*** Linux
+*** Windows
+* TODO Usage
+
+{{{biblio}}}

packages/linux.sh

@@ -44,8 +44,9 @@ cd $OLD_PWD
 echo " *** COPY DATA"
 
 mkdir -p dist/pamhyr/View/ui
+mkdir -p dist/pamhyr/View/ui/Widgets
 cp -r ../src/View/ui/ressources/ dist/pamhyr/View/ui/
-cp -r ../src/View/ui/Widgets/*.ui dist/pamhyr/View/ui/
+cp -r ../src/View/ui/Widgets/*.ui dist/pamhyr/View/ui/Widgets/
 cp -r ../src/View/ui/*.ui dist/pamhyr/View/ui/
 
 mkdir -p dist/pamhyr/lang

packages/make-windows-dir.bat

@@ -0,0 +1,54 @@
+rem windows.bat -- Pamhyr Windows batch for windows version building
+rem Copyright (C) 2023  INRAE
+rem
+rem This program is free software: you can redistribute it and/or modify
+rem it under the terms of the GNU General Public License as published by
+rem the Free Software Foundation, either version 3 of the License, or
+rem (at your option) any later version.
+rem
+rem This program is distributed in the hope that it will be useful,
+rem but WITHOUT ANY WARRANTY; without even the implied warranty of
+rem MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+rem GNU General Public License for more details.
+rem
+rem You should have received a copy of the GNU General Public License
+rem along with this program.  If not, see <https://www.gnu.org/licenses/>.
+
+@ECHO ON
+
+mkdir pamhyr\pamhyr
+xcopy /s /e /y dist\pamhyr pamhyr\pamhyr
+
+rem Copy data
+mkdir pamhyr\pamhyr\lang
+mkdir pamhyr\pamhyr\View\ui\ressources
+mkdir pamhyr\pamhyr\View\ui\Widgets
+
+rem UI
+copy /y ..\src\View\ui\ressources\ pamhyr\pamhyr\View\ui\ressources
+copy /y ..\src\View\ui\Widgets\*.ui pamhyr\pamhyr\View\ui\Widgets
+copy /y ..\src\View\ui\*.ui pamhyr\pamhyr\View\ui\
+
+rem Lang
+copy /y ..\src\lang\*.qm pamhyr\pamhyr\lang\
+
+rem Information
+copy /y ..\VERSION pamhyr\pamhyr
+copy /y ..\AUTHORS pamhyr\pamhyr
+copy /y ..\LICENSE pamhyr\pamhyr
+
+rem MAGE
+mkdir pamhyr\mage
+copy /y ..\mage-windows\mage.exe pamhyr\mage\
+copy /y ..\mage-windows\mage_extraire.exe pamhyr\mage\
+copy /y ..\mage-windows\mailleurPF.exe pamhyr\mage\
+
+rem Copy tests_cases
+mkdir pamhyr\tests_cases
+mkdir pamhyr\tests_cases\Saar
+copy /y ..\tests_cases\Saar\Saar.pamhyr pamhyr\tests_cases\Saar\
+
+rem Documentations
+mkdir pamhyr\doc
+copy /y ..\doc\users\documentation.pdf pamhyr\doc\Pamhyr2-users.pdf
+copy /y ..\doc\dev\documentation.pdf pamhyr\doc\Pamhyr2-dev.pdf

packages/make-windows-exe.bat

@@ -0,0 +1,20 @@
+rem windows.bat -- Pamhyr Windows batch for windows version building
+rem Copyright (C) 2023  INRAE
+rem
+rem This program is free software: you can redistribute it and/or modify
+rem it under the terms of the GNU General Public License as published by
+rem the Free Software Foundation, either version 3 of the License, or
+rem (at your option) any later version.
+rem
+rem This program is distributed in the hope that it will be useful,
+rem but WITHOUT ANY WARRANTY; without even the implied warranty of
+rem MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+rem GNU General Public License for more details.
+rem
+rem You should have received a copy of the GNU General Public License
+rem along with this program.  If not, see <https://www.gnu.org/licenses/>.
+
+@ECHO ON
+
+rem Make installer
+"C:\Program Files (x86)\NSIS\makensis.exe" pamhyr.nsi

packages/pamhyr.nsi

@@ -3,17 +3,17 @@
 !define LIC_NAME "LICENSE"
 !define APP_NAME "PAMHYR"
 
-Name "PAMHYR"
+Name "Pamhyr2"
-OutFile "pamhyr-win-amd64.exe"
+OutFile "pamhyr-windows.exe"
-LicenseData "..\LICENSE"
+LicenseData "pamhyr\pamhyr\LICENSE"
 LicenseText "I Agree"
 
 RequestExecutionLevel admin
 
 Unicode True
-InstallDir $PROGRAMFILES\PAMHYR
+InstallDir $PROGRAMFILES\Pamhyr2
 
-InstallDirRegKey HKLM "Software\PAMHYR" "Install_Dir"
+InstallDirRegKey HKLM "Software\Pamhyr2" "Install_Dir"
 
 Page license
 Page components
@@ -23,35 +23,35 @@ Page instfiles
 UninstPage uninstConfirm
 UninstPage instfiles
 
-Section "PAMHYR"
+Section "Pamhyr2"
   SectionIn RO
 
   SetOutPath $INSTDIR
 
-  File /r "dist\pamhyr\*"
+  File /r "pamhyr\*"
 
-  WriteRegStr HKLM "Software\Microsoft\Windows\CurrentVersion\Uninstall\PAMHYR" "DisplayName" "NSIS PAMHYR"
+  WriteRegStr HKLM "Software\Microsoft\Windows\CurrentVersion\Uninstall\Pamhyr2" "DisplayName" "NSIS Pamhyr2"
-  WriteRegStr HKLM "Software\Microsoft\Windows\CurrentVersion\Uninstall\PAMHYR" "UninstallString" "$INSTDIR\uninstall.exe"
+  WriteRegStr HKLM "Software\Microsoft\Windows\CurrentVersion\Uninstall\Pamhyr2" "UninstallString" "$INSTDIR\uninstall.exe"
-  WriteRegDWORD HKLM "Software\Microsoft\Windows\CurrentVersion\Uninstall\PAMHYR" "NoModify" 1
+  WriteRegDWORD HKLM "Software\Microsoft\Windows\CurrentVersion\Uninstall\Pamhyr2" "NoModify" 1
-  WriteRegDWORD HKLM "Software\Microsoft\Windows\CurrentVersion\Uninstall\PAMHYR" "NoRepair" 1
+  WriteRegDWORD HKLM "Software\Microsoft\Windows\CurrentVersion\Uninstall\Pamhyr2" "NoRepair" 1
   WriteUninstaller "$INSTDIR\uninstall.exe"
 
 SectionEnd
 
 Section "Start Menu Shortcuts"
 
-  CreateDirectory "$SMPROGRAMS\PAMHYR"
+  CreateDirectory "$SMPROGRAMS\Pamhyr2"
-  CreateShortcut "$SMPROGRAMS\PAMHYR\Uninstall.lnk" "$INSTDIR\uninstall.exe"
+  CreateShortcut "$SMPROGRAMS\Pamhyr2\Uninstall.lnk" "$INSTDIR\uninstall.exe"
-  CreateShortcut "$SMPROGRAMS\PAMHYR\PAMHYR.lnk" "$INSTDIR\pamhyr.exe"
+  CreateShortcut "$SMPROGRAMS\Pamhyr2\Pamhyr2.lnk" "$INSTDIR\pamhyr\pamhyr.exe"
 
 SectionEnd
 
 Section "Uninstall"
-  DeleteRegKey HKLM "Software\Microsoft\Windows\CurrentVersion\Uninstall\PAMHYR"
+  DeleteRegKey HKLM "Software\Microsoft\Windows\CurrentVersion\Uninstall\Pamhyr2"
-  DeleteRegKey HKLM SOFTWARE\NSIS_PAMHYR
+  DeleteRegKey HKLM SOFTWARE\NSIS_Pamhyr2
 
-  Delete "$SMPROGRAMS\PAMHYR\*.lnk"
+  Delete "$SMPROGRAMS\Pamhyr2\*.lnk"
 
-  RMDir "$SMPROGRAMS\PAMHYR"
+  RMDir "$SMPROGRAMS\Pamhyr2"
   RMDir /R "$INSTDIR"
 SectionEnd

src/Model/AModelComponent.py

@@ -1,10 +0,0 @@
-# -*- coding: utf-8 -*-
-
-class AbstractModelComponent(object):
-    def __init__(self):
-        super(AbstractModelComponent, self).__init__()
-
-    def new(self, filename):
-
-
-    def export(self, filename):

src/Model/BoundaryCondition/BoundaryCondition.py

@@ -20,7 +20,7 @@ import logging
 
 from tools import trace, timer, old_pamhyr_date_to_timestamp
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrDB import SQLSubModel
 from Model.Except import NotImplementedMethodeError
 
 logger = logging.getLogger()
@@ -83,7 +83,8 @@ class BoundaryCondition(SQLSubModel):
     def _get_ctor_from_type(cls, t):
         from Model.BoundaryCondition.BoundaryConditionTypes import (
             NotDefined, PonctualContribution,
-            TimeOverZ, TimeOverDischarge, ZOverDischarge
+            TimeOverZ, TimeOverDischarge, ZOverDischarge,
+            Solid,
         )
 
         res = NotDefined
@@ -95,6 +96,8 @@ class BoundaryCondition(SQLSubModel):
             res = TimeOverDischarge
         elif t == "ZD":
             res = ZOverDischarge
+        elif t == "SL":
+            res = Solid
         return res
 
     @classmethod

src/Model/BoundaryCondition/BoundaryConditionList.py

@@ -19,40 +19,23 @@
 from copy import copy
 from tools import trace, timer
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrList import PamhyrModelListWithTab
 from Model.Except import NotImplementedMethodeError
 
 from Model.BoundaryCondition.BoundaryCondition import BoundaryCondition
 from Model.BoundaryCondition.BoundaryConditionTypes import (
     NotDefined,
     PonctualContribution,
-    TimeOverZ, TimeOverDischarge, ZOverDischarge
+    TimeOverZ, TimeOverDischarge, ZOverDischarge,
+    Solid,
 )
 
-class BoundaryConditionList(SQLSubModel):
+class BoundaryConditionList(PamhyrModelListWithTab):
+    _tabs_list = ["liquid", "solid", "suspenssion"]
     _sub_classes = [
         BoundaryCondition,
     ]
 
-    def __init__(self, status = None):
-        super(BoundaryConditionList, self).__init__()
-
-        self._status = status
-
-        self._tabs = {
-            "liquid" : [],
-            "solid" : [],
-            "suspenssion" : []
-        }
-
-    @classmethod
-    def _sql_create(cls, execute):
-        return cls._create_submodel(execute)
-
-    @classmethod
-    def _sql_update(cls, execute, version):
-        return True
-
     @classmethod
     def _sql_load(cls, execute, data = None):
         new = cls(status = data['status'])
@@ -81,66 +64,12 @@ class BoundaryConditionList(SQLSubModel):
 
         return True
 
-    def len(self, lst):
-        return len(self._tabs[lst])
-
-    def get_tab(self, lst):
-        return self._tabs[lst].copy()
-
-    def get(self, lst, row):
-        return self._tabs[lst][row]
-
-    def set(self, lst, row, new):
-        self._tabs[lst][row] = new
-        self._status.modified()
-
     def new(self, lst, index):
         n = NotDefined(status=self._status)
         self._tabs[lst].insert(index, n)
         self._status.modified()
         return n
 
-    def insert(self, lst, index, new):
-        self._tabs[lst].insert(index, new)
-        self._status.modified()
-
-    def delete(self, lst, bcs):
-        for bc in bcs:
-            self._tabs[lst].remove(bc)
-        self._status.modified()
-
-    def delete_i(self, lst, indexes):
-        bcs = list(
-            map(
-                lambda x: x[1],
-                filter(
-                    lambda x: x[0] in indexes,
-                    enumerate(self._tabs[lst])
-                )
-            )
-        )
-        self.delete(lst, bcs)
-
-    def sort(self, lst, reverse=False, key=None):
-        self._tabs[lst].sort(reverse=reverse, key=key)
-        self._status.modified()
-
-    def move_up(self, lst, index):
-        if index < len(self._tabs[lst]):
-            next = index - 1
-
-            l = self._tabs[lst]
-            l[index], l[next] = l[next], l[index]
-            self._status.modified()
-
-    def move_down(self, lst, index):
-        if index >= 0:
-            prev = index + 1
-
-            l = self._tabs[lst]
-            l[index], l[prev] = l[prev], l[index]
-            self._status.modified()
-
     def __copy__(self):
         new = BoundaryConditionList()
 

src/Model/BoundaryCondition/BoundaryConditionTypes.py

@@ -83,3 +83,19 @@ class ZOverDischarge(BoundaryCondition):
     @property
     def _default_0(self):
         return 0.0
+
+
+class Solid(BoundaryCondition):
+    def __init__(self, id:int = -1, name:str = "", status = None):
+        super(Solid, self).__init__(id=id, name=name, status=status)
+
+        self.d50 = 0.002
+        self.sigma = 1
+
+        self._type = "SL"
+        self._header = ["time", "solid"]
+        self._types = [TimeOverDischarge.time_convert, float]
+
+    @classmethod
+    def compatibility(cls):
+        return ["solid"]

src/Model/Friction/Friction.py

@@ -18,7 +18,7 @@
 
 from tools import trace, timer
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrDB import SQLSubModel
 
 class Friction(SQLSubModel):
     def __init__(self, name:str = "", status = None):

src/Model/Friction/FrictionList.py

@@ -21,26 +21,18 @@ import logging
 from copy import copy
 from tools import trace, timer
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrDB import SQLSubModel
+from Model.Tools.PamhyrList import PamhyrModelList
+
 from Model.Friction.Friction import Friction
 
 logger = logging.getLogger()
 
-class FrictionList(SQLSubModel):
+class FrictionList(PamhyrModelList):
     _sub_classes = [
         Friction
     ]
 
-    def __init__(self, status = None):
-        super(FrictionList, self).__init__()
-
-        self._status = status
-        self._frictions = []
-
-    @classmethod
-    def _sql_create(cls, execute):
-        return cls._create_submodel(execute)
-
     @classmethod
     def _sql_update_0_0_1(cls, execute, version):
         execute("ALTER TABLE `section` RENAME TO `friction`")
@@ -57,7 +49,7 @@ class FrictionList(SQLSubModel):
     def _sql_load(cls, execute, data = None):
         new = cls(status = data['status'])
 
-        new._frictions = Friction._sql_load(
+        new._lst = Friction._sql_load(
             execute, data
         )
 
@@ -69,70 +61,19 @@ class FrictionList(SQLSubModel):
 
         ok = True
         ind = 0
-        for friction in self._frictions:
+        for friction in self._lst:
             data["ind"] = ind
             ok &= friction._sql_save(execute, data = data)
             ind += 1
 
         return ok
 
-    def __len__(self):
-        return len(self._frictions)
-
     @property
     def frictions(self):
-        return self._frictions.copy()
+        return self.lst
-
-    def get(self, row):
-        return self._frictions[row]
-
-    def set(self, row, new):
-        self._frictions[row] = new
-        self._status.modified()
 
     def new(self, index):
         n = Friction(status = self._status)
-        self._frictions.insert(index, n)
+        self._lst.insert(index, n)
         self._status.modified()
         return n
-
-    def insert(self, index, new):
-        self._frictions.insert(index, new)
-        self._status.modified()
-
-    def delete(self, frictions):
-        for friction in frictions:
-            self._frictions.remove(friction)
-        self._status.modified()
-
-    def delete_i(self, indexes):
-        frictions = list(
-            map(
-                lambda x: x[1],
-                filter(
-                    lambda x: x[0] in indexes,
-                    enumerate(self._frictions)
-                )
-            )
-        )
-        self.delete(frictions)
-
-    def sort(self, reverse=False, key=None):
-        self._frictions.sort(reverse=reverse, key=key)
-        self._status.modified()
-
-    def move_up(self, index):
-        if index < len(self._frictions):
-            next = index - 1
-
-            l = self._frictions
-            l[index], l[next] = l[next], l[index]
-            self._status.modified()
-
-    def move_down(self, index):
-        if index >= 0:
-            prev = index + 1
-
-            l = self._frictions
-            l[index], l[prev] = l[prev], l[index]
-            self._status.modified()

src/Model/Geometry/Point.py

@@ -19,12 +19,29 @@
 from Model.Except import NotImplementedMethodeError
 
 class Point(object):
-    def __init__(self, name:str = "", status = None):
+    def __init__(self, name:str = "", profile=None, status = None):
         super(Point, self).__init__()
 
         self._status = status
 
         self._name = name
+        self._profile = profile
+        self._sl = None
+
+
+    @property
+    def sl(self):
+        if self._sl is None:
+            return self._profile.sl
+
+        return self._sl
+
+    @sl.setter
+    def sl(self, sl):
+        if sl == self._profile.sl:
+            self._sl = None
+
+        self._sl = sl
 
     @property
     def name(self):

src/Model/Geometry/PointXYZ.py

@@ -19,15 +19,15 @@
 from math import dist
 import numpy as np
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrDB import SQLSubModel
 from Model.Geometry.Point import Point
 
 class PointXYZ(Point, SQLSubModel):
     _sub_classes = []
 
     def __init__(self, x:float = 0.0, y:float = 0.0, z:float = 0.0,
-                 name:str = "", status = None):
+                 name:str = "", profile = None, status = None):
-        super(PointXYZ, self).__init__(name=name, status=status)
+        super(PointXYZ, self).__init__(name=name, profile=profile, status=status)
 
         self._x = float(x)
         self._y = float(y)
@@ -44,7 +44,9 @@ class PointXYZ(Point, SQLSubModel):
             y INTEGER NOT NULL,
             z INTEGER NOT NULL,
             profile INTEGER NOT NULL,
-            FOREIGN KEY(profile) REFERENCES profileXYZ(id)
+            sl INTEGER,
+            FOREIGN KEY(profile) REFERENCES profileXYZ(id),
+            FOREIGN KEY(sl) REFERENCES sedimentary_layer(id)
           )
         """)
 
@@ -53,6 +55,17 @@ class PointXYZ(Point, SQLSubModel):
     @classmethod
     def _sql_update(cls, execute, version):
         cls._update_submodel(execute, version)
+
+        major, minor, release = version.strip().split(".")
+        if major == minor == "0":
+            if int(release) < 2:
+                execute(
+                    """
+                    ALTER TABLE geometry_pointXYZ
+                    ADD COLUMN sl INTEGER
+                    REFERENCES sedimentary_layer(id)
+                    """
+                )
         return True
 
     @classmethod
@@ -63,9 +76,9 @@ class PointXYZ(Point, SQLSubModel):
         profile = data["profile"]
 
         table = execute(
-            "SELECT ind, name, x, y, z " +
+            "SELECT ind, name, x, y, z, sl " +
             "FROM geometry_pointXYZ " +
-            f"WHERE profile = {profile}"
+            f"WHERE profile = {profile.id}"
         )
 
         # Create points list
@@ -79,13 +92,25 @@ class PointXYZ(Point, SQLSubModel):
             x = row[2]
             y = row[3]
             z = row[4]
+            sl = row[5]
 
             new = cls(
                 name = name,
                 x = x, y = y, z = z,
+                profile = profile,
                 status = status
             )
 
+            if sl == -1 or sl == None:
+                new._sl = None
+            else:
+                new._sl = next(
+                    filter(
+                        lambda s: s.id == sl,
+                        data["sediment_layers_list"].sediment_layers
+                    )
+                )
+
             points[ind] = new
 
         return points
@@ -94,13 +119,15 @@ class PointXYZ(Point, SQLSubModel):
         profile = data["profile"]
         ind = data["ind"]
 
+        sl = self._sl.id if self._sl is not None else -1
+
         sql = (
             "INSERT OR REPLACE INTO " +
-            "geometry_pointXYZ(ind, name, x, y, z, profile) "+
+            "geometry_pointXYZ(ind, name, x, y, z, profile, sl) "+
             "VALUES (" +
             f"{ind}, '{self._sql_format(self._name)}', " +
             f"{self.x}, {self.y}, {self.z}, " +
-            f"{profile}" +
+            f"{profile.id}, {sl}" +
             ")"
         )
         execute(sql)
@@ -117,7 +144,7 @@ class PointXYZ(Point, SQLSubModel):
                     *data
                 )
             else:
-                valid_header = {'name', 'x', 'y', 'z'}
+                valid_header = {'name', 'x', 'y', 'z', 'profile'}
                 d = {}
                 for i, v in enumerate(data):
                     h = header[i].strip().lower().split(' ')[0]

src/Model/Geometry/Profile.py

@@ -50,6 +50,7 @@ class Profile(object):
         self._kp = float(kp)
         self._name = str(name)
         self._reach = reach
+        self._sl = None
 
         self._points: List[Point] = []
 
@@ -63,6 +64,9 @@ class Profile(object):
     def points(self):
         return self._points.copy()
 
+    def point(self, index):
+        return self._points[index]
+
     @property
     def reach(self):
         return self._reach
@@ -136,6 +140,19 @@ class Profile(object):
         self._name = value.strip()
         self._status.modified()
 
+    @property
+    def sl(self):
+        """
+        Returns:
+            Profile sediment layers.
+        """
+        return self._sl
+
+    @sl.setter
+    def sl(self, value):
+        self._sl = value
+        self._status.modified()
+
     @property
     def profile_type(self):
         """
@@ -274,6 +291,39 @@ class Profile(object):
         )
         self._status.modified()
 
+    # Sediment Layers
+
+    def get_sl(self):
+        """Get sediment layer height of points
+
+        Get sediment layer of points (without spesific point sl)
+
+        Returns:
+            List of sediment layers height
+        """
+        res = []
+        psl = [point.sl for point in self.points]
+
+        # Compute max number of layers
+        sl_max = 0
+        for sl in psl:
+            n = 0 if sl is None else len(sl)
+            sl_max = max(n, sl_max)
+
+        # Create list of height for each sl and each layer
+        for i in range(0, sl_max):
+            cur = []
+            # Compute new layer line for each sl
+            for sl in psl:
+                if sl is not None and i < len(sl):
+                    cur.append(sl.get(i).height)
+                else:
+                    cur.append(0)
+            # Add layer line to result
+            res.append(cur)
+
+        return res
+
     # Abstract method, must be implemented for in non abstract class
     def get_station(self):
         raise NotImplementedMethodeError(self, self.get_station)

src/Model/Geometry/ProfileXYZ.py

@@ -21,7 +21,7 @@ from typing import List
 
 from tools import timer
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrDB import SQLSubModel
 from Model.Except import ClipboardFormatError
 from Model.Geometry.Profile import Profile
 from Model.Geometry.PointXYZ import PointXYZ
@@ -76,7 +76,9 @@ class ProfileXYZ(Profile, SQLSubModel):
             num INTEGER NOT NULL,
             code1 INTEGER NOT NULL,
             code2 INTEGER NOT NULL,
-            FOREIGN KEY(reach) REFERENCES river_reach(id)
+            sl INTEGER,
+            FOREIGN KEY(reach) REFERENCES river_reach(id),
+            FOREIGN KEY(sl) REFERENCES sedimentary_layer(id)
           )
         """)
 
@@ -84,6 +86,17 @@ class ProfileXYZ(Profile, SQLSubModel):
 
     @classmethod
     def _sql_update(cls, execute, version):
+        major, minor, release = version.strip().split(".")
+        if major == minor == "0":
+            if int(release) < 2:
+                execute(
+                    """
+                    ALTER TABLE geometry_profileXYZ
+                    ADD COLUMN sl INTEGER
+                    REFERENCES sedimentary_layer(id)
+                    """
+                )
+
         return cls._update_submodel(execute, version)
 
     @classmethod
@@ -93,7 +106,7 @@ class ProfileXYZ(Profile, SQLSubModel):
         reach = data["reach"]
 
         table = execute(
-            "SELECT id, ind, name, kp, num, code1, code2 " +
+            "SELECT id, ind, name, kp, num, code1, code2, sl " +
             "FROM geometry_profileXYZ " +
             f"WHERE reach = {reach}"
         )
@@ -109,6 +122,7 @@ class ProfileXYZ(Profile, SQLSubModel):
             num = row[5]
             code1 = row[5]
             code2 = row[6]
+            sl = row[7]
 
             new = cls(
                 id=id, num = num,
@@ -118,7 +132,17 @@ class ProfileXYZ(Profile, SQLSubModel):
                 status = status
             )
 
-            data["profile"] = id
+            if sl == -1 or sl == None:
+                new._sl = None
+            else:
+                new._sl = next(
+                    filter(
+                        lambda s: s.id == sl,
+                        data["sediment_layers_list"].sediment_layers
+                    )
+                )
+
+            data["profile"] = new
             new._points = PointXYZ._sql_load(execute, data)
 
             profiles[ind] = new
@@ -129,18 +153,20 @@ class ProfileXYZ(Profile, SQLSubModel):
         ok = True
         ind = data["ind"]
 
+        sl = self._sl.id if self._sl is not None else -1
+
         sql = (
             "INSERT OR REPLACE INTO " +
-            "geometry_profileXYZ(id, ind, name, reach, kp, num, code1, code2) "+
+            "geometry_profileXYZ(id, ind, name, reach, kp, num, code1, code2, sl) "+
             "VALUES (" +
             f"{self.id}, {ind}, '{self._sql_format(self._name)}', " +
             f"{self.reach.id}, {self.kp}, {self.num}, " +
-            f"{self.code1}, {self.code1}" +
+            f"{self.code1}, {self.code1}, {sl}" +
             ")"
         )
         execute(sql)
 
-        data["profile"] = self.id
+        data["profile"] = self
         execute(f"DELETE FROM geometry_pointXYZ WHERE profile = {self.id}")
 
         ind = 0
@@ -178,11 +204,11 @@ class ProfileXYZ(Profile, SQLSubModel):
         try:
             if len(header) == 0:
                 point = PointXYZ(
-                    *data, status=self._status
+                    *data, profile=self, status=self._status
                 )
             else:
                 valid_header = {'name', 'x', 'y', 'z'}
-                d = {"status": self._status}
+                d = {"status": self._status, "profile": self}
                 for i, v in enumerate(data):
                     h = header[i].strip().lower().split(' ')[0]
                     if h in valid_header:
@@ -235,7 +261,7 @@ class ProfileXYZ(Profile, SQLSubModel):
             Nothing.
         """
         for point in list_points:
-            pt = PointXYZ(*point, status=self._status)
+            pt = PointXYZ(*point, profile=self, status=self._status)
             self._points.append(pt)
         self._status.modified()
 
@@ -259,7 +285,7 @@ class ProfileXYZ(Profile, SQLSubModel):
         Returns:
             Nothing.
         """
-        point_xyz = PointXYZ(0., 0., 0., status=self._status)
+        point_xyz = PointXYZ(0., 0., 0., profile=self, status=self._status)
         self._points.append(point_xyz)
         self._status.modified()
 
@@ -272,7 +298,7 @@ class ProfileXYZ(Profile, SQLSubModel):
         Returns:
             The new point.
         """
-        point = PointXYZ(0., 0., 0., status=self._status)
+        point = PointXYZ(0., 0., 0., profile=self, status=self._status)
         self._points.insert(index, point)
         self._status.modified()
         return point

src/Model/Geometry/Reach.py

@@ -27,7 +27,7 @@ from functools import reduce
 
 from tools import flatten, timer, trace
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrDB import SQLSubModel
 
 from Model.Geometry.Profile import Profile
 from Model.Geometry.ProfileXYZ import ProfileXYZ
@@ -276,6 +276,39 @@ class Reach(SQLSubModel):
         return max(self.get_kp())
 
 
+    # Sediment Layers
+
+    def get_sl(self):
+        """Get sediment layer height of profile
+
+        Get sediment layer of profile (without spesific point sl)
+
+        Returns:
+            List of sediment layers height
+        """
+        res = []
+        psl = [profile.sl for profile in self.profiles]
+
+        # Compute max number of layers
+        sl_max = 0
+        for sl in psl:
+            n = 0 if sl is None else len(sl)
+            sl_max = max(n, sl_max)
+
+        # Create list of height for each sl and each layer
+        for i in range(0, sl_max):
+            cur = []
+            # Compute new layer line for each sl
+            for sl in psl:
+                if sl is not None and i < len(sl):
+                    cur.append(sl.get(i).height)
+                else:
+                    cur.append(0)
+            # Add layer line to result
+            res.append(cur)
+
+        return res
+
     # Guidelines
 
     @timer

src/Model/InitialConditions/InitialConditions.py

@@ -20,7 +20,7 @@ from copy import copy, deepcopy
 from tools import trace, timer
 from functools import reduce
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrDB import SQLSubModel
 
 class Data(SQLSubModel):
     def __init__(self, name:str = "",

src/Model/InitialConditions/InitialConditionsDict.py

@@ -19,29 +19,14 @@
 from copy import copy
 from tools import trace, timer
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrDict import PamhyrModelDict
 from Model.InitialConditions.InitialConditions import InitialConditions
 
-class InitialConditionsDict(SQLSubModel):
+class InitialConditionsDict(PamhyrModelDict):
     _sub_classes = [
         InitialConditions,
     ]
 
-    def __init__(self, status = None):
-        super(InitialConditionsDict, self).__init__()
-
-        self._status = status
-
-        self._reach = {}
-
-    @classmethod
-    def _sql_create(cls, execute):
-        return cls._create_submodel(execute)
-
-    @classmethod
-    def _sql_update(cls, execute, version):
-        return cls._update_submodel(execute, version)
-
     @classmethod
     def _sql_load(cls, execute, data = None):
         new = cls(status = data["status"])
@@ -55,7 +40,7 @@ class InitialConditionsDict(SQLSubModel):
             )
 
             if ic is not None:
-                new._reach[reach] = ic
+                new._dict[reach] = ic
 
         return new
 
@@ -66,30 +51,12 @@ class InitialConditionsDict(SQLSubModel):
 
         execute("DELETE FROM initial_conditions")
 
-        for reach in self._reach:
+        for reach in self._dict:
             data["reach"] = reach
-            ok &= self._reach[reach]._sql_save(execute, data)
+            ok &= self._dict[reach]._sql_save(execute, data)
 
         return ok
 
-    def __len__(self):
-        return len(self._reach)
-
-    def is_defined(self, reach):
-        return reach in self._reach
-
-    def get(self, reach):
-        if reach in self._reach:
-            return self._reach[reach]
-
-        new = self.new(reach)
-        self.set(reach, new)
-        return new
-
-    def set(self, reach, new):
-        self._reach[reach] = new
-        self._status.modified()
-
     def new(self, reach):
         new = InitialConditions(reach = reach, status = self._status)
         self.set(reach, new)

src/Model/LateralContribution/LateralContribution.py

@@ -20,7 +20,7 @@ import logging
 
 from tools import trace, timer, old_pamhyr_date_to_timestamp
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrDB import SQLSubModel
 from Model.Except import NotImplementedMethodeError
 
 logger = logging.getLogger()

src/Model/LateralContribution/LateralContributionList.py

@@ -19,7 +19,7 @@
 from copy import copy
 from tools import trace, timer
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrList import PamhyrModelListWithTab
 from Model.Except import NotImplementedMethodeError
 
 from Model.LateralContribution.LateralContribution import LateralContribution
@@ -27,30 +27,12 @@ from Model.LateralContribution.LateralContributionTypes import (
     NotDefined, LateralContrib, Rain, Evaporation,
 )
 
-class LateralContributionList(SQLSubModel):
+class LateralContributionList(PamhyrModelListWithTab):
+    _tabs_list = ["liquid", "solid", "suspenssion"]
     _sub_classes = [
         LateralContribution,
     ]
 
-    def __init__(self, status = None):
-        super(LateralContributionList, self).__init__()
-
-        self._status = status
-
-        self._tabs = {
-            "liquid" : [],
-            "solid" : [],
-            "suspenssion" : []
-        }
-
-    @classmethod
-    def _sql_create(cls, execute):
-        return cls._create_submodel(execute)
-
-    @classmethod
-    def _sql_update(cls, execute, version):
-        return True
-
     @classmethod
     def _sql_load(cls, execute, data = None):
         new = cls(status = data['status'])
@@ -79,66 +61,12 @@ class LateralContributionList(SQLSubModel):
 
         return True
 
-    def len(self, lst):
-        return len(self._tabs[lst])
-
-    def get_tab(self, lst):
-        return self._tabs[lst].copy()
-
-    def get(self, lst, row):
-        return self._tabs[lst][row]
-
-    def set(self, lst, row, new):
-        self._tabs[lst][row] = new
-        self._status.modified()
-
     def new(self, lst, index):
         n = NotDefined(status=self._status)
         self._tabs[lst].insert(index, n)
         self._status.modified()
         return n
 
-    def insert(self, lst, index, new):
-        self._tabs[lst].insert(index, new)
-        self._status.modified()
-
-    def delete(self, lst, bcs):
-        for bc in bcs:
-            self._tabs[lst].remove(bc)
-        self._status.modified()
-
-    def delete_i(self, lst, indexes):
-        bcs = list(
-            map(
-                lambda x: x[1],
-                filter(
-                    lambda x: x[0] in indexes,
-                    enumerate(self._tabs[lst])
-                )
-            )
-        )
-        self.delete(lst, bcs)
-
-    def sort(self, lst, reverse=False, key=None):
-        self._tabs[lst].sort(reverse=reverse, key=key)
-        self._status.modified()
-
-    def move_up(self, lst, index):
-        if index < len(self._tabs[lst]):
-            next = index - 1
-
-            l = self._tabs[lst]
-            l[index], l[next] = l[next], l[index]
-            self._status.modified()
-
-    def move_down(self, lst, index):
-        if index >= 0:
-            prev = index + 1
-
-            l = self._tabs[lst]
-            l[index], l[prev] = l[prev], l[index]
-            self._status.modified()
-
     def __copy__(self):
         new = LateralContributionList()
 

src/Model/Results/Results.py

@@ -0,0 +1,50 @@
+# Results.py -- Pamhyr
+# Copyright (C) 2023  INRAE
+#
+# This program 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.
+#
+# This program 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 this program.  If not, see <https://www.gnu.org/licenses/>.
+
+import logging
+import numpy as np
+
+from copy import deepcopy
+from datetime import datetime
+
+from Model.Results.River.River import River
+
+logger = logging.getLogger()
+
+class Results(object):
+    def __init__(self, study = None):
+        self._study = study
+        self._river = River(self._study)
+
+        self._meta_data = {
+            # Keep results creation date
+            "creation_date": datetime.now(),
+        }
+
+    @property
+    def date(self):
+        date = self._meta_data["creation_date"]
+        return f"{date.isoformat(sep=' ')}"
+
+    @property
+    def river(self):
+        return self._river
+
+    def set(self, key, value):
+        self._meta_data[key] = value
+
+    def get(self, key):
+        return self._meta_data[key]

src/Model/Results/River/River.py

@@ -0,0 +1,110 @@
+# River.py -- Pamhyr
+# Copyright (C) 2023  INRAE
+#
+# This program 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.
+#
+# This program 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 this program.  If not, see <https://www.gnu.org/licenses/>.
+
+import logging
+
+from datetime import datetime
+
+logger = logging.getLogger()
+
+class Profile(object):
+    def __init__(self, profile, study):
+        self._study = study
+        self._profile = profile # Source profile in the study
+        self._data = {} # Dict of dict {<ts>: {<key>: <value>, ...}, ...}
+
+    @property
+    def name(self):
+        return self._profile.name
+
+    @property
+    def kp(self):
+        return self._profile.kp
+
+    @property
+    def geometry(self):
+        return self._profile
+
+    def set(self, timestamp, key, data):
+        if timestamp not in self._data:
+            self._data[timestamp] = {}
+
+        self._data[timestamp][key] = data
+
+    def get_ts(self, timestamp):
+        return self._data[timestamp]
+
+    def get_key(self, key):
+        return list(
+            map(lambda ts: self._data[ts][key], self._data)
+        )
+
+    def get_ts_key(self, timestamp, key):
+        return self._data[timestamp][key]
+
+class Reach(object):
+    def __init__(self, reach, study):
+        self._study = study
+        self._reach = reach # Source reach in the study
+        self._profiles = list(
+            map(
+                lambda p: Profile(p, self._study),
+                reach.profiles
+            )
+        )
+
+    @property
+    def name(self):
+        return self._reach.name
+
+    @property
+    def geometry(self):
+        return self._reach
+
+    @property
+    def profiles(self):
+        return self._profiles.copy()
+
+    def profile(self, id):
+        return self._profiles[id]
+
+    def set(self, profile_id, timestamp, key, data):
+        self._profiles[profile_id].set(timestamp, key, data)
+
+class River(object):
+    def __init__(self, study):
+        self._study = study
+
+        # Dict with timestamps as key
+        self._reachs = []
+
+    @property
+    def reachs(self):
+        return self._reachs.copy()
+
+    def reach(self, id):
+        return self._reachs[id]
+
+    def has_reach(self, id):
+        return 0 <= id < len(self._reachs)
+
+    def add(self, reach_id):
+        reachs = self._study.river.enable_edges()
+
+        new = Reach(reachs[reach_id].reach, self._study)
+
+        self._reachs.append(new)
+        return new

src/Model/River.py

@@ -16,7 +16,9 @@
 
 # -*- coding: utf-8 -*-
 
-from Model.DB import SQLSubModel
+from tools import flatten
+
+from Model.Tools.PamhyrDB import SQLSubModel
 
 from Model.Network.Node import Node
 from Model.Network.Edge import Edge
@@ -31,6 +33,7 @@ from Model.InitialConditions.InitialConditionsDict import InitialConditionsDict
 from Model.Stricklers.StricklersList import StricklersList
 from Model.Friction.FrictionList import FrictionList
 from Model.SolverParameters.SolverParametersList import SolverParametersList
+from Model.SedimentLayer.SedimentLayerList import SedimentLayerList
 
 from Solver.Solvers import solver_type_list
 
@@ -208,6 +211,7 @@ class River(Graph, SQLSubModel):
         InitialConditionsDict,
         StricklersList,
         SolverParametersList,
+        SedimentLayerList,
     ]
 
     def __init__(self, status=None):
@@ -223,6 +227,7 @@ class River(Graph, SQLSubModel):
         self._initial_conditions = InitialConditionsDict(status=self._status)
         self._stricklers = StricklersList(status=self._status)
         self._parameters = {}
+        self._sediment_layers = SedimentLayerList(status=self._status)
 
     @classmethod
     def _sql_create(cls, execute):
@@ -246,6 +251,14 @@ class River(Graph, SQLSubModel):
         )
         data["stricklers"] = new._stricklers
 
+
+        # Initial conditions
+        new._sediment_layers = SedimentLayerList._sql_load(
+            execute,
+            data
+        )
+        data["sediment_layers_list"] = new._sediment_layers
+
         # Network
         new._nodes = RiverNode._sql_load(
             execute,
@@ -289,6 +302,7 @@ class River(Graph, SQLSubModel):
         objs.append(self._boundary_condition)
         objs.append(self._initial_conditions)
         objs.append(self._lateral_contribution)
+        objs.append(self._sediment_layers)
         objs.append(self._stricklers)
 
         for solver in self._parameters:
@@ -297,14 +311,6 @@ class River(Graph, SQLSubModel):
         self._save_submodel(execute, objs, data)
         return True
 
-    @property
-    def reach(self):
-        return self._reach
-
-    @property
-    def frictions(self):
-        return self._frictions
-
     @property
     def boundary_condition(self):
         return self._boundary_condition
@@ -317,6 +323,10 @@ class River(Graph, SQLSubModel):
     def initial_conditions(self):
         return self._initial_conditions
 
+    @property
+    def sediment_layers(self):
+        return self._sediment_layers
+
     @property
     def stricklers(self):
         return self._stricklers
@@ -366,3 +376,16 @@ class River(Graph, SQLSubModel):
 
     def set_current_reach(self, reach):
         self._current_reach = reach
+
+    def has_sediment(self):
+        has = len(self._sediment_layers) != 0
+        has &= any(
+            filter(
+                lambda p: p.sl != None,
+                flatten(
+                    map(lambda e: e.reach.profiles, self.edges())
+                )
+            )
+        )
+
+        return has

src/Model/SedimentLayer/SedimentLayer.py

@@ -0,0 +1,339 @@
+# -*- coding: utf-8 -*-
+
+from tools import trace, timer
+
+from Model.Tools.PamhyrDB import SQLSubModel
+from Model.Except import NotImplementedMethodeError
+
+class Layer(SQLSubModel):
+    _sub_classes = []
+    _id_cnt = 0
+
+    def __init__(self,
+                 id:int = -1, name:str = "",
+                 type = "",
+                 height = 0.0, d50 = 0.0, sigma = 0.0,
+                 critical_constraint = 0.0,
+                 sl = None, status = None):
+        super(Layer, self).__init__()
+
+        self._status = status
+
+        self._name = name
+        self._type = type
+
+        self._height = height
+        self._d50 = d50
+        self._sigma = sigma
+        self._critical_constraint = critical_constraint
+
+        if id == -1:
+            self.id = Layer._id_cnt
+        else:
+            self.id = id
+
+        Layer._id_cnt = max(id, Layer._id_cnt+1)
+
+    @property
+    def name(self):
+        return self._name
+
+    @name.setter
+    def name(self, name):
+        self._name = name
+
+    @property
+    def type(self):
+        return self._type
+
+    @type.setter
+    def type(self, type):
+        self._type = type
+
+    @property
+    def height(self):
+        return self._height
+
+    @height.setter
+    def height(self, height):
+        self._height = float(height)
+
+    @property
+    def d50(self):
+        return self._d50
+
+    @d50.setter
+    def d50(self, d50):
+        self._d50 = float(d50)
+
+    @property
+    def sigma(self):
+        return self._sigma
+
+    @sigma.setter
+    def sigma(self, sigma):
+        self._sigma = float(sigma)
+
+    @property
+    def critical_constraint(self):
+        return self._critical_constraint
+
+    @critical_constraint.setter
+    def critical_constraint(self, critical_constraint):
+        self._critical_constraint = float(critical_constraint)
+
+    @classmethod
+    def _sql_create(cls, execute):
+        execute("""
+          CREATE TABLE sedimentary_layer_layer(
+            id INTEGER NOT NULL PRIMARY KEY,
+            ind INTEGER NOT NULL,
+            name TEXT NOT NULL,
+            type TEXT NOT NULL,
+            height REAL NOT NULL,
+            d50 REAL NOT NULL,
+            sigma REAL NOT NULL,
+            critical_constraint REAL NOT NULL,
+            sl INTEGER,
+            FOREIGN KEY(sl) REFERENCES sedimentary_layer(id)
+          )
+        """)
+
+        return cls._create_submodel(execute)
+
+    @classmethod
+    def _sql_update(cls, execute, version):
+        major, minor, release = version.strip().split(".")
+        if major == minor == "0":
+            if int(release) < 2:
+                cls._sql_create(execute)
+
+        return True
+
+    @classmethod
+    def _sql_load(cls, execute, data = None):
+        new = []
+        sl = data["sl"]
+
+        table = execute(
+            "SELECT id, ind, name, type, height, d50, sigma, critical_constraint " +
+            "FROM sedimentary_layer_layer " +
+            f"WHERE sl = {sl}"
+        )
+
+        for _ in table:
+            new.append(None)
+
+        for row in table:
+            ind = row[1]
+
+            layer = cls(
+                id = row[0], name = row[2],
+                type = row[3], height = row[4],
+                d50 = row[5], sigma = row[6],
+                critical_constraint = row[7],
+                sl = sl, status = data['status']
+            )
+
+            new[ind] = layer
+
+        return new
+
+    def _sql_save(self, execute, data = None):
+        ind = data["ind"]
+        sl = data["sl"]
+
+        sql = (
+            "INSERT INTO " +
+            "sedimentary_layer_layer(id, ind, name, type, height, d50, sigma, critical_constraint, sl) "+
+            "VALUES (" +
+            f"{self.id}, {ind}, '{self._sql_format(self._name)}', " +
+            f"'{self._sql_format(self._type)}', {self._height}, " +
+            f"{self._d50}, {self._sigma}, {self._critical_constraint}, " +
+            f"{sl.id}" +
+            ")"
+        )
+        execute(sql)
+
+        return True
+
+
+class SedimentLayer(SQLSubModel):
+    _sub_classes = [Layer]
+    _id_cnt = 0
+
+    def __init__(self, id:int = -1,
+                 name:str = "", comment:str = "",
+                 status = None):
+        super(SedimentLayer, self).__init__()
+
+        self._status = status
+        self._name = name
+        self._comment = comment
+        self._layers = []
+
+        if id == -1:
+            self.id = SedimentLayer._id_cnt
+        else:
+            self.id = id
+
+        SedimentLayer._id_cnt = max(id, SedimentLayer._id_cnt+1)
+
+    def __str__(self):
+        s = f"{self.name} ({len(self)})"
+
+        if self.comment != "":
+            s += f" - {self.comment}"
+
+        return s
+
+    def __len__(self):
+        return len(self._layers)
+
+    @property
+    def layers(self):
+        return self._layers.copy()
+
+    def height(self):
+        return list(
+            map(lambda l: l.height, self._layers)
+        )
+
+    @property
+    def name(self):
+        return self._name
+
+    @name.setter
+    def name(self, name):
+        self._name = name
+
+    def names(self):
+        return list(
+            map(lambda l: l.name, self._layers)
+        )
+
+    @property
+    def comment(self):
+        return self._comment
+
+    @comment.setter
+    def comment(self, comment):
+        self._comment = comment
+
+    @classmethod
+    def _sql_create(cls, execute):
+        execute("""
+          CREATE TABLE sedimentary_layer(
+            id INTEGER NOT NULL PRIMARY KEY,
+            name TEXT NOT NULL,
+            comment TEXT NOT NULL
+          )
+        """)
+
+        return cls._create_submodel(execute)
+
+    @classmethod
+    def _sql_update(cls, execute, version):
+        major, minor, release = version.strip().split(".")
+        if major == minor == "0":
+            if int(release) < 2:
+                cls._sql_create(execute)
+
+        return True
+
+    @classmethod
+    def _sql_load(cls, execute, data = None):
+        new = []
+
+        table = execute(
+            "SELECT id, name, comment " +
+            "FROM sedimentary_layer "
+        )
+
+        for row in table:
+            sl = cls(
+                id = row[0],
+                name = row[1],
+                comment = row[2],
+                status = data['status']
+            )
+
+            data["sl"] = sl.id
+            sl._layers = Layer._sql_load(execute, data)
+
+            new.append(sl)
+
+        return new
+
+    def _sql_save(self, execute, data = None):
+        if data is None:
+            data = {}
+
+        sql = (
+            "INSERT INTO sedimentary_layer (id, name, comment) " +
+            f"VALUES ({self.id}, '{self._sql_format(self._name)}', " +
+            f"'{self._sql_format(self._comment)}')"
+        )
+        execute(sql)
+
+        data["sl"] = self
+
+        ind = 0
+        for l in self._layers:
+            data["ind"] = ind
+            l._sql_save(execute, data)
+            ind += 1
+
+        return True
+
+    def get(self, index):
+        return self._layers[index]
+
+    def set(self, index, new):
+        self._layers[index] = new
+        self._status.modified()
+
+    def insert(self, index, new):
+        self._layers.insert(index, new)
+        self._status.modified()
+
+    def new(self, index):
+        n = Layer(sl=self, status = self._status)
+        self.insert(index, n)
+        self._status.modified()
+        return n
+
+    def delete(self, els):
+        for el in els:
+            self._layers.remove(el)
+        self._status.modified()
+
+    def delete_i(self, indexes):
+        els = list(
+            map(
+                lambda x: x[1],
+                filter(
+                    lambda x: x[0] in indexes,
+                    enumerate(self._layers)
+                )
+            )
+        )
+        self.delete(els)
+
+    def move_up(self, index):
+        if index >= 0:
+            next = index - 1
+
+            l = self._layers
+            l[index], l[next] = l[next], l[index]
+
+            self._status.modified()
+
+    def move_down(self, index):
+        if index + 1 < len(self._layers):
+            prev = index + 1
+
+            l = self._layers
+            l[index], l[prev] = l[prev], l[index]
+
+            self._status.modified()

src/Model/SedimentLayer/SedimentLayerList.py

@@ -0,0 +1,40 @@
+# -*- coding: utf-8 -*-
+
+from tools import trace, timer
+
+from Model.Except import NotImplementedMethodeError
+from Model.Tools.PamhyrList import PamhyrModelList
+from Model.SedimentLayer.SedimentLayer import SedimentLayer
+
+class SedimentLayerList(PamhyrModelList):
+    _sub_classes = [SedimentLayer]
+
+    @classmethod
+    def _sql_load(cls, execute, data = None):
+        new = cls(status = data["status"])
+
+        new._lst = SedimentLayer._sql_load(execute, data)
+
+        return new
+
+    def _sql_save(self, execute, data = None):
+        ok = True
+
+        # Delete previous data
+        execute("DELETE FROM sedimentary_layer")
+        execute("DELETE FROM sedimentary_layer_layer")
+
+        for sl in self._lst:
+            ok &= sl._sql_save(execute, data)
+
+        return ok
+
+    @property
+    def sediment_layers(self):
+        return self.lst
+
+    def new(self, index):
+        n = SedimentLayer(status = self._status)
+        self.insert(index, n)
+        self._status.modified()
+        return n

src/Model/SolverParameters/SolverParametersList.py

@@ -20,7 +20,7 @@ from copy import copy
 from tools import trace, timer
 from Solver.Solvers import solver_type_list
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrList import PamhyrModelList
 
 class Parameter():
     def __init__(self,
@@ -63,13 +63,13 @@ class Parameter():
 
         return new
 
-class SolverParametersList(SQLSubModel):
+class SolverParametersList(PamhyrModelList):
     def __init__(self, solver_type = None, status = None):
         super(SolverParametersList, self).__init__()
 
         self._status = status
         self._solver = solver_type
-        self._parameters = list(
+        self._lst = list(
             map(
                 lambda t: Parameter.from_tuple(t, self._status),
                 solver_type.default_parameters()
@@ -93,6 +93,66 @@ class SolverParametersList(SQLSubModel):
 
     @classmethod
     def _sql_update(cls, execute, version):
+        major, minor, release = version.strip().split(".")
+
+        if major == minor == "0":
+            if int(release) < 3:
+                execute(f"UPDATE solver_parameter SET name='mage_implicitation' WHERE name='mage_implication'")
+                execute(f"UPDATE solver_parameter SET name='mage_iteration_type' WHERE name='mage_iter_type'")
+
+            if int(release) < 4:
+                solvers = execute("SELECT DISTINCT solver FROM solver_parameter")
+
+                for solver in solvers:
+                    if solver != "mage8":
+                        continue
+
+                    data = execute("SELECT ind, name, value " +
+                                  "FROM solver_parameter "+
+                                  f"WHERE solver = '{solver}' " +
+                                  "ORDER BY ind ACS")
+
+                    execute(
+                        "DELETE FROM solver_parameter " +
+                        f"WHERE solver = '{solver}'"
+                    )
+
+                    insert = lambda name, value, ind: execute(
+                            "INSERT INTO " +
+                            "solver_parameter(ind, name, value, solver) "+
+                            "VALUES (" +
+                            f"{ind}, " +
+                            f"'{cls._sql_format(name)}', " +
+                            f"'{cls._sql_format(value)}', " +
+                            f"'{cls._sql_format(solver)}'" +
+                            ")"
+                    )
+
+                    ind = 0
+                    lst = []
+                    for v in data:
+                        insert(v[0], v[1], ind)
+                        ind += 1
+
+                    new = [
+                        ("mage_sediment_masse_volumique", "2650.0"),
+                        ("mage_sediment_angle_repos", "40.0"),
+                        ("mage_sediment_porosity", "0.40"),
+                        ("mage_distance_han", "0.0"),
+                        ("mage_distance_chargement_d50", "100.0"),
+                        ("mage_distance_chargement_sigma", "100.0"),
+                        ("mage_methode_modification_geometrie", "1"),
+                        ("mage_shields_critique", "1"),
+                        ("mage_shields_correction", "1"),
+                        ("mage_capacite_solide", "1"),
+                        ("mage_pas_de_temps_charriage", "1"),
+                        ("mage_facteur_multiplicateur", "1.0"),
+                    ]
+
+                    for v in new:
+                        insert(v[0], v[1], ind)
+                        ind += 1
+
         return cls._update_submodel(execute, version)
 
     @classmethod
@@ -134,7 +194,7 @@ class SolverParametersList(SQLSubModel):
         )
 
         ind = 0
-        for param in self._parameters:
+        for param in self._lst:
             sql = (
                 "INSERT INTO " +
                 "solver_parameter(ind, name, value, solver) "+
@@ -150,28 +210,36 @@ class SolverParametersList(SQLSubModel):
 
         return True
 
-    def __len__(self):
-        return len(self._parameters)
-
     @property
     def parameters(self):
-        return self._parameters.copy()
+        return self.lst
 
     def get(self, index):
-        return self._parameters[index]
+        return self._lst[index]
+
+    def get_by_key(self, key):
+        try:
+            return next(
+                filter(
+                    lambda p: p["name"] == key,
+                    self._lst
+                )
+            )["value"]
+        except:
+            return None
 
     def set(self, index, new):
-        self._parameters[index] = new
+        self._lst[index] = new
         self._status.modified()
 
     def set_value(self, key, value):
-        for p in self._parameters:
+        for p in self._lst:
             if p["name"] == key:
                 p["value"] = value
                 self._status.modified()
                 return
 
-        self._parameters.append(
+        self._lst.append(
             Parameter(
                 name = key,
                 value = value,
@@ -182,47 +250,6 @@ class SolverParametersList(SQLSubModel):
 
     def new(self, index):
         n = Parameter(status = self._status)
-        self._parameters.insert(index, n)
+        self._lst.insert(index, n)
         self._status.modified()
         return n
-
-    def insert(self, index, new):
-        self._parameters.insert(index, new)
-        self._status.modified()
-
-    def delete(self, parameters):
-        for parameter in parameters:
-            self._parameters.remove(parameter)
-        self._status.modified()
-
-    def delete_i(self, indexes):
-        parameters = list(
-            map(
-                lambda x: x[1],
-                filter(
-                    lambda x: x[0] in indexes,
-                    enumerate(self._parameters)
-                )
-            )
-        )
-        self.delete(parameters)
-
-    def sort(self, reverse=False, key=None):
-        self._parameters.sort(reverse=reverse, key=key)
-        self._status.modified()
-
-    def move_up(self, index):
-        if index < len(self._parameters):
-            next = index - 1
-
-            l = self._parameters
-            l[index], l[next] = l[next], l[index]
-            self._status.modified()
-
-    def move_down(self, index):
-        if index >= 0:
-            prev = index + 1
-
-            l = self._parameters
-            l[index], l[prev] = l[prev], l[index]
-            self._status.modified()

src/Model/Stricklers/Stricklers.py

@@ -18,7 +18,7 @@
 
 from tools import trace, timer
 
-from Model.DB import SQLSubModel
+from Model.Tools.PamhyrDB import SQLSubModel
 
 class Stricklers(SQLSubModel):
     _id_cnt = 0
@@ -38,7 +38,7 @@ class Stricklers(SQLSubModel):
         else:
             self.id = id
 
-        Stricklers._id_cnt = max(Stricklers._id_cnt + 1, self.id)
+        Stricklers._id_cnt = max(Stricklers._id_cnt + 1, id + 1)
 
         self._name = name
         self._comment = comment

src/Model/Stricklers/StricklersList.py

@@ -18,22 +18,15 @@
 
 from tools import trace, timer
 
-from Model.DB import SQLSubModel
 from Model.Saved import SavedStatus
+from Model.Tools.PamhyrList import PamhyrModelList
 from Model.Stricklers.Stricklers import Stricklers
 
-class StricklersList(SQLSubModel):
+class StricklersList(PamhyrModelList):
     _sub_classes = [
         Stricklers,
     ]
 
-    def __init__(self, status = None):
-        if status is None:
-            status = SavedStatus()
-        self._status = status
-
-        self._stricks = []
-
     @classmethod
     def _sql_create(cls, execute):
         return cls._create_submodel(execute)
@@ -46,7 +39,7 @@ class StricklersList(SQLSubModel):
     def _sql_load(cls, execute, data = None):
         new = cls(status = data["status"])
 
-        new._stricks = Stricklers._sql_load(
+        new._lst = Stricklers._sql_load(
             execute,
             data = data
         )
@@ -56,60 +49,26 @@ class StricklersList(SQLSubModel):
     def _sql_save(self, execute, data = None):
         execute("DELETE FROM stricklers")
 
-        objs = self._stricks
+        objs = self._lst
         return self._save_submodel(execute, objs, data)
 
-    def __len__(self):
-        return len(self._stricks)
-
     @property
     def stricklers(self):
-        return self._stricks.copy()
+        return self.lst
 
-    def get(self, index):
+    def new(self, index):
-        if 0 <= index < len(self._stricks):
-            return self._stricks[index]
-
-        return None
-
-    def insert(self, index, strick):
-        self._stricks.insert(index, strick)
-        self._status.modified()
-
-    def add(self, index):
         s = Stricklers(status = self._status)
         self.insert(index, s)
         return s
 
-    def delete(self, stricks):
-        self._stricks = list(
-            filter(
-                lambda s: s not in stricks,
-                self._stricks
-            )
-        )
-        self._status.modified()
-
-    def delete_i(self, indexes):
-        stricks = set(
-            map(
-                lambda e: e[1],
-                filter(
-                    lambda e: e[0] in indexes,
-                    enumerate(self._stricks)
-                )
-            )
-        )
-        self.delete(stricks)
-
     @timer
     def sort(self, reverse:bool = False, key=None):
         f = lambda st: st.name
         if key is not None:
             f = key
 
-        self._stricks = sorted(
+        self._lst = sorted(
-            self._stricks,
+            self._lst,
             key = f,
             reverse = reverse,
         )

src/Model/Study.py

@@ -23,7 +23,7 @@ from datetime import datetime
 
 from tools import timer, timestamp
 
-from Model.DB import SQLModel
+from Model.Tools.PamhyrDB import SQLModel
 from Model.Saved import SavedStatus
 from Model.Serializable import Serializable
 from Model.Except import NotImplementedMethodeError
@@ -40,7 +40,7 @@ class Study(SQLModel):
 
     def __init__(self, filename = None, init_new = True):
         # Metadata
-        self._version = "0.0.1"
+        self._version = "0.0.4"
         self.creation_date = datetime.now()
         self.last_modification_date = datetime.now()
         self.last_save_date = datetime.now()

src/Model/Tools/PamhyrDB.py

@@ -1,4 +1,4 @@
-# DB.py -- Pamhyr abstract model database classes
+# PamhyrDB.py -- Pamhyr abstract model database classes
 # Copyright (C) 2023  INRAE
 #
 # This program is free software: you can redistribute it and/or modify

src/Model/Tools/PamhyrDict.py

@@ -0,0 +1,74 @@
+# PamhyrDict.py -- Pamhyr abstract dict for model classes
+# Copyright (C) 2023  INRAE
+#
+# This program 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.
+#
+# This program 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 this program.  If not, see <https://www.gnu.org/licenses/>.
+
+# -*- coding: utf-8 -*-
+
+import logging
+
+from copy import copy
+from tools import trace, timer
+
+from Model.Tools.PamhyrDB import SQLSubModel
+from Model.Except import NotImplementedMethodeError
+
+logger = logging.getLogger()
+
+
+class PamhyrModelDict(SQLSubModel):
+    _sub_classes = []
+
+    def __init__(self, status = None):
+        super(PamhyrModelDict, self).__init__()
+
+        self._status = status
+
+        self._dict = {}
+
+    @classmethod
+    def _sql_create(cls, execute):
+        return cls._create_submodel(execute)
+
+    @classmethod
+    def _sql_update(cls, execute, version):
+        return cls._update_submodel(execute, version)
+
+    @classmethod
+    def _sql_load(cls, execute, data = None):
+        raise NotImplementedMethodeError(cls, cls._sql_load)
+
+    def _sql_save(self, execute, data = None):
+        raise NotImplementedMethodeError(self, self._sql_save)
+
+    def __len__(self):
+        return len(self._dict)
+
+    def is_defined(self, key):
+        return key in self._dict
+
+    def set(self, key, new):
+        self._dict[key] = new
+        self._status.modified()
+
+    def get(self, key):
+        if key in self._dict:
+            return self._dict[key]
+
+        new = self.new(key)
+        self.set(key, new)
+        return new
+
+    def new(self, key):
+        raise NotImplementedMethodeError(self, self.new)

src/Model/Tools/PamhyrList.py

@@ -0,0 +1,229 @@
+# PamhyrList.py -- Pamhyr Abstract List object for the Model
+# Copyright (C) 2023  INRAE
+#
+# This program 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.
+#
+# This program 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 this program.  If not, see <https://www.gnu.org/licenses/>.
+
+# -*- coding: utf-8 -*-
+
+import logging
+
+from copy import copy
+from tools import trace, timer
+
+from Model.Tools.PamhyrDB import SQLSubModel
+from Model.Except import NotImplementedMethodeError
+
+logger = logging.getLogger()
+
+class PamhyrModelList(SQLSubModel):
+    _sub_classes = [
+    ]
+
+    def __init__(self, status = None):
+        super(PamhyrModelList, self).__init__()
+
+        self._status = status
+        self._lst = []
+
+    #######
+    # SQL #
+    #######
+
+    @classmethod
+    def _sql_create(cls, execute):
+        return cls._create_submodel(execute)
+
+    @classmethod
+    def _sql_update(cls, execute, version):
+        return cls._update_submodel(execute, version)
+
+    @classmethod
+    def _sql_load(cls, execute, data = None):
+        raise NotImplementedMethodeError(cls, cls._sql_load)
+
+    def _sql_save(self, execute, data = None):
+        raise NotImplementedMethodeError(self, self._sql_save)
+
+    ################
+    # MODEL METHOD #
+    ################
+
+    def __len__(self):
+        return len(self._lst)
+
+    @property
+    def lst(self):
+        return self._lst.copy()
+
+    def get(self, row):
+        return self._lst[row]
+
+    def set(self, row, new):
+        self._lst[row] = new
+        self._status.modified()
+
+    def new(self, index):
+        """Create new elements and add it to list
+
+        Args:
+            index: The index of new elements
+
+        Returns:
+            The new elements
+        """
+        raise NotImplementedMethodeError(self, self.new)
+
+    def insert(self, index, new):
+        self._lst.insert(index, new)
+        self._status.modified()
+
+    def delete(self, lst):
+        for el in lst:
+            self._lst.remove(el)
+        self._status.modified()
+
+    def delete_i(self, indexes):
+        lst = list(
+            map(
+                lambda x: x[1],
+                filter(
+                    lambda x: x[0] in indexes,
+                    enumerate(self._lst)
+                )
+            )
+        )
+        self.delete(lst)
+
+    def sort(self, reverse=False, key=None):
+        self._lst.sort(reverse=reverse, key=key)
+        self._status.modified()
+
+    def move_up(self, index):
+        if index < len(self._lst):
+            next = index - 1
+
+            l = self._lst
+            l[index], l[next] = l[next], l[index]
+            self._status.modified()
+
+    def move_down(self, index):
+        if index >= 0:
+            prev = index + 1
+
+            l = self._lst
+            l[index], l[prev] = l[prev], l[index]
+            self._status.modified()
+
+
+class PamhyrModelListWithTab(SQLSubModel):
+    _tabs_list = []
+    _sub_classes = [
+    ]
+
+    def __init__(self, status = None):
+        super(PamhyrModelListWithTab, self).__init__()
+
+        self._status = status
+        self._tabs = {}
+        for tab in self._tabs_list:
+            self._tabs[tab] = []
+
+    #######
+    # SQL #
+    #######
+
+    @classmethod
+    def _sql_create(cls, execute):
+        return cls._create_submodel(execute)
+
+    @classmethod
+    def _sql_update(cls, execute, version):
+        return cls._update_submodel(execute, version)
+
+    @classmethod
+    def _sql_load(cls, execute, data = None):
+        raise NotImplementedMethodeError(cls, cls._sql_load)
+
+    def _sql_save(self, execute, data = None):
+        raise NotImplementedMethodeError(self, self._sql_save)
+
+    ################
+    # MODEL METHOD #
+    ################
+
+    def len(self, lst):
+        return len(self._tabs[lst])
+
+    def get_tab(self, lst):
+        return self._tabs[lst].copy()
+
+    def get(self, lst, row):
+        return self._tabs[lst][row]
+
+    def set(self, lst, row, new):
+        self._tabs[lst][row] = new
+        self._status.modified()
+
+    def new(self, lst, index):
+        """Create new elements and add it to list
+
+        Args:
+            lst: The tab name
+            index: The index of new elements
+
+        Returns:
+            The new elements
+        """
+        raise NotImplementedMethodeError(self, self.new)
+
+    def insert(self, lst, index, new):
+        self._tabs[lst].insert(index, new)
+        self._status.modified()
+
+    def delete(self, lst, els):
+        for el in els:
+            self._tabs[lst].remove(el)
+        self._status.modified()
+
+    def delete_i(self, lst, indexes):
+        els = list(
+            map(
+                lambda x: x[1],
+                filter(
+                    lambda x: x[0] in indexes,
+                    enumerate(self._tabs[lst])
+                )
+            )
+        )
+        self.delete(lst, els)
+
+    def sort(self, lst, reverse=False, key=None):
+        self._tabs[lst].sort(reverse=reverse, key=key)
+        self._status.modified()
+
+    def move_up(self, lst, index):
+        if index < len(self._tabs[lst]):
+            next = index - 1
+
+            l = self._tabs[lst]
+            l[index], l[next] = l[next], l[index]
+            self._status.modified()
+
+    def move_down(self, lst, index):
+        if index >= 0:
+            prev = index + 1
+
+            l = self._tabs[lst]
+            l[index], l[prev] = l[prev], l[index]
+            self._status.modified()

src/Scripts/AScript.py

@@ -0,0 +1,36 @@
+# AScript.py -- Pamhyr abstract script class
+# Copyright (C) 2023  INRAE
+#
+# This program 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.
+#
+# This program 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 this program.  If not, see <https://www.gnu.org/licenses/>.
+
+# -*- coding: utf-8 -*-
+
+from Model.Except import NotImplementedMethodeError
+
+class AScript(object):
+    name = ""
+    description = ""
+
+    def __init__(self, app, conf, args):
+        super(AScript, self).__init__()
+
+        self._app = app
+        self._conf = conf
+        self._args = args.copy()
+
+    def usage(self):
+        raise NotImplementedMethodeError(self, self.usage)
+
+    def run(self):
+        raise NotImplementedMethodeError(self, self.run)

src/Scripts/Hello.py

@@ -0,0 +1,34 @@
+# help.py -- Pamhyr
+# Copyright (C) 2023  INRAE
+#
+# This program 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.
+#
+# This program 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 this program.  If not, see <https://www.gnu.org/licenses/>.
+
+# -*- coding: utf-8 -*-
+
+import logging
+from Scripts.AScript import AScript
+
+logger = logging.getLogger()
+
+
+class ScriptHello(AScript):
+    name = "Hello"
+    description = "Display hello"
+
+    def usage(self):
+        logger.info(f"Usage : {self._args[0]} hello")
+
+    def run(self):
+        logger.info("Hello from Pamhyr2 script !")
+        return 0

src/Solver/ASolver.py

@@ -22,6 +22,7 @@ import logging
 from tools import timer
 
 try:
+    # Installation allow Unix-like signal
     from signal import SIGTERM, SIGSTOP, SIGCONT
     _signal = True
 except:
@@ -67,6 +68,9 @@ class AbstractSolver(object):
         self._process = None
         self._output = None
 
+        # Last study running
+        self._study = None
+
     def __str__(self):
         return f"{self._name} : {self._type} : {self._description}"
 
@@ -86,6 +90,7 @@ class AbstractSolver(object):
             ("all_init_time", "000:00:00:00"),
             ("all_final_time", "999:99:00:00"),
             ("all_timestep", "300.0"),
+            ("all_command_line_arguments", ""),
         ]
 
         return lst
@@ -153,6 +158,17 @@ class AbstractSolver(object):
     def export(self, study, repertory, qlog = None):
         raise NotImplementedMethodeError(self, self.export)
 
+    def cmd_args(self, study):
+        """Return solver command line arguments list
+
+        Returns:
+            Command line arguments list
+        """
+        params = study.river.get_params(self.type)
+        args = params.get_by_key("all_command_line_arguments")
+
+        return args.split(" ")
+
     def input_param(self):
         """Return input command line parameter(s)
 
@@ -186,11 +202,11 @@ class AbstractSolver(object):
         return os.path.abspath(
             os.path.join(
                 os.path.dirname(__file__),
-                ".."
+                "..", ".."
             )
         )
 
-    def _format_command(self, cmd, path = ""):
+    def _format_command(self, study, cmd, path = ""):
         """Format command line
 
         Args:
@@ -201,11 +217,11 @@ class AbstractSolver(object):
             The executable and list of arguments
         """
         # HACK: Works in most case... Trust me i'm an engineer
-
         cmd = cmd.replace("@install_dir", self._install_dir())
         cmd = cmd.replace("@path", path.replace(" ", "\ "))
         cmd = cmd.replace("@input", self.input_param())
         cmd = cmd.replace("@dir", self._process.workingDirectory())
+        cmd = cmd.replace("@args", " ".join(self.cmd_args(study)))
 
         logger.debug(f"! {cmd}")
 
@@ -213,23 +229,33 @@ class AbstractSolver(object):
             # Command line executable path is between " char
             cmd = cmd.split("\"")
             exe = cmd[1].replace("\ ", " ")
-            args = "\"".join(cmd[2:]).split(" ")[1:]
+            args = list(
+                filter(
+                    lambda s: s != "",
+                    "\"".join(cmd[2:]).split(" ")[1:]
+                )
+            )
         else:
             # We suppose the command line executable path as no space char
             cmd = cmd.replace("\ ", "&_&").split(" ")
             exe = cmd[0].replace("&_&", " ")
-            args = list(map(lambda s: s.replace("&_&", "\ "), cmd[1:]))
+            args = list(
+                filter(
+                    lambda s: s != "",
+                    map(lambda s: s.replace("&_&", "\ "), cmd[1:])
+                )
+            )
 
         logger.info(f"! {exe} {args}")
         return exe, args
 
-    def run_input_data_fomater(self):
+    def run_input_data_fomater(self, study):
         if self._cmd_input == "":
-            self._run_next()
+            self._run_next(study)
             return True
 
         cmd = self._cmd_input
-        exe, args = self._format_command(cmd, self._path_input)
+        exe, args = self._format_command(study, cmd, self._path_input)
 
         if not os.path.exists(exe):
             error = f"[ERROR] Path {exe} do not exists"
@@ -242,13 +268,13 @@ class AbstractSolver(object):
 
         return True
 
-    def run_solver(self):
+    def run_solver(self, study):
         if self._cmd_solver == "":
-            self._run_next()
+            self._run_next(study)
             return True
 
         cmd = self._cmd_solver
-        exe, args = self._format_command(cmd, self._path_solver)
+        exe, args = self._format_command(study, cmd, self._path_solver)
 
         if not os.path.exists(exe):
             error = f"[ERROR] Path {exe} do not exists"
@@ -262,13 +288,13 @@ class AbstractSolver(object):
         self._status = STATUS.RUNNING
         return True
 
-    def run_output_data_fomater(self):
+    def run_output_data_fomater(self, study):
         if self._cmd_output == "":
-            self._run_next()
+            self._run_next(study)
             return True
 
         cmd = self._cmd_output
-        exe, args = self._format_command(cmd, self._path_output)
+        exe, args = self._format_command(study, cmd, self._path_output)
 
         if not os.path.exists(exe):
             error = f"[ERROR] Path {exe} do not exists"
@@ -287,29 +313,29 @@ class AbstractSolver(object):
             for x in s.split('\n'):
                 self._output.put(x)
 
-    def _run_next(self):
+    def _run_next(self, study):
         self._step += 1
         if self._step < len(self._runs):
-            res = self._runs[self._step]()
+            res = self._runs[self._step](study)
             if res is not True:
                 self._output.put(res)
         else:
             self._status = STATUS.STOPED
 
-    def _finished(self, exit_code, exit_status):
+    def _finished(self, study, exit_code, exit_status):
         if self._output is not None:
             self._output.put(exit_code)
 
-        self._run_next()
+        self._run_next(study)
 
-    def run(self, process = None, output_queue = None):
+    def run(self, study, process = None, output_queue = None):
         if process is not None:
             self._process = process
         if output_queue is not None:
             self._output = output_queue
 
         self._process.readyRead.connect(self._data_ready)
-        self._process.finished.connect(self._finished)
+        self._process.finished.connect(lambda c, s: self._finished(study, c, s))
 
         self._runs = [
             self.run_input_data_fomater,
@@ -318,7 +344,7 @@ class AbstractSolver(object):
         ]
         self._step = 0
         # Run first step
-        res = self._runs[0]()
+        res = self._runs[0](study)
         if res is not True:
             self._output.put(res)
 
@@ -330,14 +356,15 @@ class AbstractSolver(object):
         self._status = STATUS.STOPED
         return True
 
-    def start(self, process = None):
+    def start(self, study, process = None):
         if _signal:
+            # Solver is PAUSED, so continue execution
             if self._status == STATUS.PAUSED:
                 os.kill(self._process.pid(), SIGCONT)
                 self._status = STATUS.RUNNING
                 return True
 
-        self.run(process)
+        self.run(study, process)
         return True
 
     def pause(self):
@@ -345,6 +372,7 @@ class AbstractSolver(object):
             if self._process is None:
                 return False
 
+            # Send SIGSTOP to PAUSED solver
             os.kill(self._process.pid(), SIGSTOP)
             self._status = STATUS.PAUSED
             return True

src/Solver/Mage.py

@@ -61,11 +61,11 @@ class Mage(AbstractSolver):
             ("mage_timestep_tra", "3600"),
             ("mage_timestep_bin", "0"),
             # ("mage_timestep_melissa", "0"),
-            ("mage_implication", "0.70"),
+            ("mage_implicitation", "0.70"),
             ("mage_continuity_discretization", "S"),
             ("mage_qsj_discretization", "B"),
             ("mage_stop_criterion_iterations", "R"),
-            ("mage_iter_type", "0"),
+            ("mage_iteration_type", "0"),
             ("mage_smooth_coef", "0"),
             ("mage_cfl_max", "-1."),
             ("mage_min_height", "0.1"),
@@ -97,19 +97,31 @@ class Mage(AbstractSolver):
     # Export #
     ##########
 
+    def cmd_args(self, study):
+        l = super(Mage, self).cmd_args(study)
+
+        l.append("-r")
+
+        return l
+
     def input_param(self):
-        return "0.REP"
+        name = self._study.name
+        return f"{name}.REP"
 
     def log_file(self):
-        return "0.TRA"
+        name = self._study.name
+        return f"{name}.TRA"
 
     @timer
-    def _export_ST(self, study, repertory, qlog):
+    def _export_ST(self, study, repertory, qlog, name="0"):
         files = []
 
         if qlog is not None:
             qlog.put("Export ST file")
 
+        os.makedirs(os.path.join(repertory, "net"), exist_ok=True)
+        gra_file = f"{name}.GRA"
+
         # Write header
         edges = study.river.edges()
         edges = list(
@@ -124,33 +136,57 @@ class Mage(AbstractSolver):
             if edge._name == "":
                 name = f"Reach_{edge.id}"
 
-            with mage_file_open(os.path.join(repertory, f"{name}.ST"), "w+") as f:
+            with mage_file_open(os.path.join(repertory, "net", f"{name}.ST"), "w+") as f:
-                files.append(f"{name}.ST")
+                files.append(str(os.path.join("net", f"{name}.ST")))
 
+                cnt_num = 1
                 for profile in edge.reach.profiles:
-                    num = f"{profile.num:>6}"
+                    num = f"{cnt_num:>6}"
                     c1 = f"{profile.code1:>6}"
                     c2 = f"{profile.code2:>6}"
                     t = f"{len(profile.points):>6}"
                     kp = f"{profile.kp:>13.4f}"
-                    name = profile.name
+                    pname = profile.name
+                    if profile.name == "":
+                        pname = f"p{profile.id:>3}".replace(" ", "p")
+                    name = f"{pname}"
 
-                    f.write(f"{num}{c1}{c2}{t}{kp}  {name}\n")
+                    sediment = ""
+                    if profile.sl is not None:
+                        if not any(filter(lambda f: ".GRA" in f, files)):
+                            files.append(gra_file)
+
+                        nl = len(profile.sl)
+                        sediment = f" {nl:>3}"
+                        for l in profile.sl.layers:
+                            sediment += f" {l.height:>10} {l.d50:>10} {l.sigma:>10} {l.critical_constraint:>10}"
+
+                    f.write(f"{num}{c1}{c2}{t}{kp}  {name} {sediment}\n")
+                    cnt_num += 1
 
                     for point in profile.points:
                         x = f"{point.x:>13.4f}"
                         y = f"{point.y:>13.4f}"
                         z = f"{point.z:>13.4f}"
-                        n = point.name
+                        n = f"{point.name:<3}"
 
-                        f.write(f"{x}{y}{z} {n}\n")
+                        sediment = ""
+                        prev = point.z
+                        if point.sl is not None:
+                            nl = len(point.sl)
+                            sediment = f"{nl:>3}"
+                            for l in point.sl.layers:
+                                prev = round(prev - l.height, 5)
+                                sediment += f" {prev:>10} {l.d50:>10} {l.sigma:>10} {l.critical_constraint:>10}"
+
+                        f.write(f"{x}{y}{z} {n} {sediment}\n")
 
                     f.write(f"     999.9990     999.9990     999.9990\n")
 
         return files
 
     @timer
-    def _export_BC(self, t, bounds, repertory, qlog):
+    def _export_BC(self, t, bounds, repertory, qlog, name="0"):
         files = []
 
         if len(bounds) == 0:
@@ -159,8 +195,8 @@ class Mage(AbstractSolver):
         if qlog is not None:
             qlog.put(f"Export {t} file")
 
-        with mage_file_open(os.path.join(repertory, f"0.{t}"), "w+") as f:
+        with mage_file_open(os.path.join(repertory, f"{name}.{t}"), "w+") as f:
-            files.append(f"0.{t}")
+            files.append(f"{name}.{t}")
 
             for bound in bounds:
                 name = f"{bound.node.id:3}".replace(" ", "x")
@@ -170,18 +206,25 @@ class Mage(AbstractSolver):
                 f.write(f"*{header[0]:>9}|{header[1]:>10}\n")
 
                 for d in bound.data:
-                    f.write(f"{d[0]:10.3f}{d[1]:10.3f}\n")
+                    v0 = d[0]
+                    v1 = d[1]
+
+                    if t in ["HYD", "QSO", "LIM"]:
+                        v0 /= 60 # Convert first column to minute
+
+                    f.write(f"{v0:10}{v1:10}\n")
 
         return files
 
     @timer
-    def _export_bound_cond(self, study, repertory, qlog):
+    def _export_bound_cond(self, study, repertory, qlog, name="0"):
         files = []
         lst = study.river.boundary_condition
 
         AVA = []
         HYD = []
         LIM = []
+        QSO = []
 
         for tab in ["liquid", "solid", "suspenssion"]:
             for bound in lst.get_tab(tab):
@@ -191,22 +234,25 @@ class Mage(AbstractSolver):
                     HYD.append(bound)
                 elif bound.bctype == "TZ":
                     LIM.append(bound)
+                elif bound.bctype == "SL":
+                    QSO.append(bound)
 
-        files = files + self._export_BC("AVA", AVA, repertory,  qlog)
+        files = files + self._export_BC("AVA", AVA, repertory, qlog, name=name)
-        files = files + self._export_BC("HYD", HYD, repertory,  qlog)
+        files = files + self._export_BC("HYD", HYD, repertory, qlog, name=name)
-        files = files + self._export_BC("LIM", LIM, repertory,  qlog)
+        files = files + self._export_BC("LIM", LIM, repertory, qlog, name=name)
+        files = files + self._export_QSO(QSO, repertory, qlog, name=name)
 
         return files
 
     # @timer
-    # def _export_LC(self, lateral, repertory, qlog):
+    # def _export_LC(self, lateral, repertory, qlog, name="0"):
     #     files = []
 
     #     if qlog is not None:
     #         qlog.put(f"Export LAT file")
 
-    #     with mage_file_open(os.path.join(repertory, f"0.LAT"), "w+") as f:
+    #     with mage_file_open(os.path.join(repertory, f"{name}.LAT"), "w+") as f:
-    #         files.append(f"0.LAT")
+    #         files.append(f"{name}.LAT")
 
     #         name = f"{lateral.node.id:3}".replace(" ", "x")
     #         f.write(f"* {lateral.node.name} ({name}) {lateral.bctype}\n")
@@ -215,12 +261,12 @@ class Mage(AbstractSolver):
     #         f.write(f"*{header[0]:>9}|{header[1]:>10}\n")
 
     #         for d in lateral.data:
-    #             f.write(f"{d[0]:10.3f}{d[1]:10.3f}\n")
+    #             f.write(f"{d[0]:1{name}.3f}{d[1]:10.3f}\n")
 
     #     return files
 
     # @timer
-    # def _export_lateral_contrib(self, study, repertory, qlog):
+    # def _export_lateral_contrib(self, study, repertory, qlog, name="0"):
     #     files = []
     #     lst = study.river.lateral_contribution
 
@@ -231,15 +277,15 @@ class Mage(AbstractSolver):
     #     return files
 
     @timer
-    def _export_RUG(self, study, repertory, qlog):
+    def _export_RUG(self, study, repertory, qlog, name="0"):
         files = []
 
         if qlog is not None:
             qlog.put("Export RUG file")
 
         # Write header
-        with mage_file_open(os.path.join(repertory, "0.RUG"), "w+") as f:
+        with mage_file_open(os.path.join(repertory, f"{name}.RUG"), "w+") as f:
-            files.append("0.RUG")
+            files.append(f"{name}.RUG")
 
             edges = study.river.edges()
             edges = list(
@@ -272,14 +318,14 @@ class Mage(AbstractSolver):
         return files
 
     @timer
-    def _export_INI(self, study, repertory, qlog):
+    def _export_INI(self, study, repertory, qlog, name="0"):
         files = []
 
         if qlog is not None:
             qlog.put("Export INI file")
 
         # Write header
-        with mage_file_open(os.path.join(repertory, "0.INI"), "w+") as f:
+        with mage_file_open(os.path.join(repertory, f"{name}.INI"), "w+") as f:
             has_ini = False
             id = 1
             reachs = study.river.edges()
@@ -317,30 +363,41 @@ class Mage(AbstractSolver):
                 id += 1
 
         if has_ini:
-            files.append("0.INI")
+            files.append(f"{name}.INI")
         return files
 
     @timer
-    def _export_REP(self, study, repertory, files, qlog):
+    def _export_REP(self, study, repertory, files, qlog, name="0"):
         if qlog is not None:
             qlog.put("Export REP file")
 
         # Write header
-        with mage_file_open(os.path.join(repertory, f"0.REP"), "w+") as f:
+        with mage_file_open(os.path.join(repertory, f"{name}.REP"), "w+") as f:
             f.write("confirmation=non\n")
 
             for file in files:
                 EXT = file.split('.')[1]
 
-                f.write(f"{EXT} {file}\n")
+                if EXT not in ["ST", "GRA"]:
+                    f.write(f"{EXT} {file}\n")
 
             f.write("* OUTPUT\n")
-            f.write(f"TRA 0.TRA\n")
+            f.write(f"TRA {name}.TRA\n")
-            f.write(f"BIN 0.BIN\n")
+            f.write(f"BIN {name}.BIN\n")
+
+            for file in files:
+                EXT = file.split('.')[1]
+
+                if EXT in ["GRA"]:
+                    f.write(f"{EXT} {file}\n")
+
 
     @timer
     def export(self, study, repertory, qlog = None):
-        self._export_ST(study, repertory, qlog)
+        self._study = study
+        name = study.name.replace(" ", "_")
+
+        self._export_ST(study, repertory, qlog, name = name)
 
         return True
 
@@ -348,14 +405,14 @@ class Mage(AbstractSolver):
     # RESULTS #
     ###########
 
-    def read_bin(self, study, repertory, results, qlog = None):
+    def read_bin(self, study, repertory, results, qlog = None, name="0"):
         return
 
     @timer
-    def results(self, study, repertory, qlog = None):
+    def results(self, study, repertory, qlog = None, name="0"):
         results = Results(study = study)
 
-        self.read_bin(study, repertory, results, qlog)
+        self.read_bin(study, repertory, results, qlog, name=name)
 
         return results
 
@@ -398,21 +455,43 @@ class Mage8(Mage):
         i = names.index("mage_precision_reduction_factor_Q")
         lst.insert(i+1, ("mage_precision_reduction_factor_r", "1"))
 
+        # Mage parameter for sediment module (added in DB 0.0.4)
+        lst.append(("mage_sediment_masse_volumique", "2650.0"))
+        lst.append(("mage_sediment_angle_repos", "40.0"))
+        lst.append(("mage_sediment_porosity", "0.40"))
+        lst.append(("mage_distance_Han", "0.0"))
+        lst.append(("mage_distance_chargement_d50", "100.0"))
+        lst.append(("mage_distance_chargement_sigma", "100.0"))
+        lst.append(("mage_methode_modification_geometrie", "1"))
+        lst.append(("mage_shields_critique", "1"))
+        lst.append(("mage_shields_correction", "1"))
+        lst.append(("mage_capacite_solide", "1"))
+        lst.append(("mage_pas_de_temps_charriage", "1"))
+        lst.append(("mage_facteur_multiplicateur", "1.0"))
+
         return lst
 
     ##########
     # Export #
     ##########
 
+    def cmd_args(self, study):
+        l = super(Mage8, self).cmd_args(study)
+
+        if study.river.has_sediment():
+            l.append("-c=3")
+
+        return l
+
     @timer
-    def _export_PAR(self, study, repertory, qlog = None):
+    def _export_PAR(self, study, repertory, qlog = None, name = "0"):
         files = []
 
         if qlog is not None:
             qlog.put("Export PAR file")
 
-        with mage_file_open(os.path.join(repertory, "0.PAR"), "w+") as f:
+        with mage_file_open(os.path.join(repertory, f"{name}.PAR"), "w+") as f:
-            files.append("0.PAR")
+            files.append(f"{name}.PAR")
 
             params = study.river.get_params(self.type).parameters
             for p in params:
@@ -421,20 +500,25 @@ class Mage8(Mage):
                         .replace("mage_", "")
                 value = p.value
 
-                f.write(f"{name} {value}\n")
+                if name in ["command_line_arguments"]:
+                    continue
 
+                if name is "mage_compute_reach_volume_balance":
+                    value = "O" if value == "y" else "N"
+
+                f.write(f"{name} {value}\n")
 
         return files
 
     @timer
-    def _export_NET(self, study, repertory, qlog = None):
+    def _export_NET(self, study, repertory, qlog = None, name="0"):
         files = []
 
         if qlog is not None:
             qlog.put("Export NET file")
 
-        with mage_file_open(os.path.join(repertory, "0.NET"), "w+") as f:
+        with mage_file_open(os.path.join(repertory, f"{name}.NET"), "w+") as f:
-            files.append("0.NET")
+            files.append(f"{name}.NET")
 
             edges = study.river.edges()
             edges = list(
@@ -453,23 +537,63 @@ class Mage8(Mage):
 
                 n1 = f"{e.node1.id:3}".replace(" ", "x")
                 n2 = f"{e.node2.id:3}".replace(" ", "x")
-                file = name + ".ST"
+                file = os.path.join("net", name + ".ST")
 
                 f.write(f"{id} {n1} {n2} {file}\n")
 
         return files
 
     @timer
-    def export(self, study, repertory, qlog = None):
+    def _export_QSO(self, bounds, repertory, qlog, name="0"):
         files = []
 
-        self._export_ST(study, repertory, qlog)
+        if len(bounds) == 0:
-        files = files + self._export_PAR(study, repertory, qlog)
+            return files
-        files = files + self._export_NET(study, repertory, qlog)
+
-        files = files + self._export_bound_cond(study, repertory, qlog)
+        if qlog is not None:
-        files = files + self._export_RUG(study, repertory, qlog)
+            qlog.put(f"Export QSO file")
-        files = files + self._export_INI(study, repertory, qlog)
+
-        self._export_REP(study, repertory, files, qlog)
+        with mage_file_open(os.path.join(repertory, f"{name}.QSO"), "w+") as f:
+            files.append(f"{name}.QSO")
+
+            for bound in bounds:
+                # File header
+                name = f"{bound.node.id:3}".replace(" ", "x")
+                f.write(f"* {bound.node.name} ({name}) {bound.bctype}\n")
+
+                d50 = bound.d50
+                sigma = bound.sigma
+
+                if len(bound.data) == 0:
+                    f.write(f"${name} {d50} {sigma} default\n")
+                else:
+                    f.write(f"${name} {d50} {sigma}\n")
+
+                # Table header
+                header = bound.header
+                f.write(f"*{header[0]:>9}|{header[1]:>10}\n")
+
+                # Data
+                for d in bound.data:
+                    f.write(f"{d[0]:10.3f}{d[1]:10.3f}\n")
+
+        return files
+
+
+    @timer
+    def export(self, study, repertory, qlog = None, name="0"):
+        self._study = study
+        name = study.name.replace(" ", "_")
+
+        files = []
+
+        files = self._export_ST(study, repertory, qlog, name=name)
+        files = files + self._export_PAR(study, repertory, qlog, name=name)
+        files = files + self._export_NET(study, repertory, qlog, name=name)
+        files = files + self._export_bound_cond(study, repertory, qlog, name=name)
+        files = files + self._export_RUG(study, repertory, qlog, name=name)
+        files = files + self._export_INI(study, repertory, qlog, name=name)
+        self._export_REP(study, repertory, files, qlog, name=name)
 
         return True
 
@@ -477,10 +601,12 @@ class Mage8(Mage):
     # RESULTS #
     ###########
 
-    def read_bin(self, study, repertory, results, qlog = None):
+    @timer
-        with mage_file_open(os.path.join(repertory, f"0.BIN"), "r") as f:
+    def read_bin(self, study, repertory, results, qlog = None, name="0"):
-            logger.info("read_bin: Start ...")
+        fname = os.path.join(repertory, f"{name}.BIN")
+        logger.info(f"read_bin: Start reading '{fname}' ...")
 
+        with mage_file_open(fname, "r") as f:
             newline = lambda: np.fromfile(f, dtype=np.int32, count=1)
             endline = lambda: np.fromfile(f, dtype=np.int32, count=1)
 
@@ -525,16 +651,6 @@ class Mage8(Mage):
             iprofiles = {}
             reach_offset = {}
 
-            ip_to_r = lambda i: iprofiles[
-                next(
-                    filter(
-                        lambda k: k[0] <= i <= k[1],
-                        iprofiles
-                    )
-                )
-            ]
-            ip_to_ri = lambda r, i: i - reach_offset[r]
-
             data = read_int(2*nb_reach)
 
             for i in range(nb_reach):
@@ -570,6 +686,17 @@ class Mage8(Mage):
             # Data
             newline()
 
+            ip_to_r = lambda i: iprofiles[
+                next(
+                    filter(
+                        lambda k: k[0] <= i <= k[1],
+                        iprofiles
+                    )
+                )
+            ]
+            ip_to_ri = lambda r, i: i - reach_offset[r]
+
+            ts = set()
             end = False
             while not end:
                 n = read_int(1)[0]
@@ -577,9 +704,10 @@ class Mage8(Mage):
                 key = bytearray(np.fromfile(f, dtype=np.byte, count=1)).decode()
                 data = read_float(n)
 
-                if key in ['Z', 'Q']:
+                logger.debug(f"read_bin: timestamp = {timestamp} sec")
-                    logger.debug(f"read_bin: timestamp = {timestamp} sec")
+                ts.add(timestamp)
 
+                if key in ["Z", "Q"]:
                     for i, d in enumerate(data):
                         # Get reach corresponding to profile ID
                         reach = ip_to_r(i)
@@ -592,4 +720,169 @@ class Mage8(Mage):
                 endline()
                 end = newline().size <= 0
 
-            logger.info("read_bin: ... end")
+            logger.debug(reachs[0].profiles[0]._data)
+            results.set("timestamps", ts)
+            logger.info(f"read_bin: ... end with {len(ts)} timestamp read")
+
+    @timer
+    def read_gra(self, study, repertory, results, qlog = None, name="0"):
+        if not study.river.has_sediment():
+            return
+
+        fname = os.path.join(repertory, f"{name}.GRA")
+        logger.info(f"read_gra: Start reading '{fname}' ...")
+
+        with mage_file_open(fname, "r") as f:
+            newline = lambda: np.fromfile(f, dtype=np.int32, count=1)
+            endline = lambda: np.fromfile(f, dtype=np.int32, count=1)
+
+            read_int = lambda size: np.fromfile(f, dtype=np.int32, count=size)
+            read_float = lambda size: np.fromfile(f, dtype=np.float32, count=size)
+            read_float64 = lambda size: np.fromfile(f, dtype=np.float64, count=size)
+
+            # Meta data (1st line)
+            newline()
+
+            data = read_int(3)
+
+            nb_reach = data[0]
+            nb_profile = data[1]
+            mage_version = data[2]
+
+            logger.debug(f"read_gra: nb_reach = {nb_reach}")
+            logger.debug(f"read_gra: nb_profile = {nb_profile}")
+            logger.debug(f"read_gra: mage_version = {mage_version}")
+
+            if mage_version <= 80:
+                msg = (
+                    "Read GRA files: " +
+                    f"Possible incompatible mage version '{mage_version}', " +
+                    "please check your solver configuration..."
+                )
+                logger.warning(msg)
+
+                if qlog is not None:
+                    qlog.put("[WARNING] " + msg)
+
+            results.set("gra_solver_version", f"Mage8 ({mage_version})")
+            results.set("gra_nb_reach", f"{nb_reach}")
+            results.set("gra_nb_profile", f"{nb_profile}")
+
+            endline()
+
+            # Reach information (2nd line)
+            newline()
+
+            reachs = []
+            iprofiles = {}
+            reach_offset = {}
+
+            data = read_int(2*nb_reach)
+
+            for i in range(nb_reach):
+                # Get results reach to reach list
+                r = results.river.reach(i)
+                reachs.append(r)
+
+                # ID of first and last reach profiles
+                i1 = data[2*i] - 1
+                i2 = data[2*i+1] - 1
+
+                # Add profile id correspondance to reach
+                key = (i1, i2)
+                iprofiles[key] = r
+
+                # Profile ID offset
+                reach_offset[r] = i1
+
+            logger.debug(f"read_gra: iprofiles = {iprofiles}")
+
+            endline()
+
+            # X (3rd line)
+            newline()
+            _ = read_float(nb_profile)
+            endline()
+
+            # npts (4th line)
+            newline()
+            _ = read_int(nb_profile)
+            endline()
+
+            # Data
+            ip_to_r = lambda i: iprofiles[
+                next(
+                    filter(
+                        lambda k: k[0] <= i <= k[1],
+                        iprofiles
+                    )
+                )
+            ]
+            ip_to_ri = lambda r, i: i - reach_offset[r]
+
+            ts = set()
+            ind = 0
+            end = False
+
+            newline()
+            while not end:
+                n = read_int(1)[0]
+                timestamp = read_float64(1)[0]
+
+                logger.debug(f"read_gra: timestamp = {timestamp} sec")
+                ts.add(timestamp)
+
+                endline()
+
+                for i in range(n):
+                    newline()
+                    nsl = read_int(1)[0]
+                    endline()
+
+                    # Get current profile id
+                    reach = ip_to_r(i)
+                    ri = ip_to_ri(reach, i)
+
+                    if nsl > 1:
+                        logger.warning(
+                            "read_gra: " +
+                            "Multiple sediment layers for one profile " +
+                            "is not implemented yet..."
+                        )
+
+                    for j in range(nsl):
+                        newline()
+                        nl = read_int(1)[0]
+                        endline()
+
+                        sl = []
+
+                        for k in range(nl):
+                            newline()
+                            data = read_float64(3)
+                            endline()
+
+                            h = data[0]
+                            d50 = data[1]
+                            sigma = data[2]
+
+                            sl.append((h, d50, sigma))
+
+                        reach.set(ri, timestamp, "sl", sl)
+
+                ind += 1
+                end = newline().size <= 0
+
+            logger.debug(reachs[0].profiles[0]._data)
+            results.set("timestamps", ts)
+            logger.info(f"read_gra: ... end with {len(ts)} timestamp read")
+
+    @timer
+    def results(self, study, repertory, qlog = None):
+        self._study = study
+        name = study.name.replace(" ", "_")
+
+        results = super(Mage8, self).results(study, repertory, qlog, name = name)
+        self.read_gra(study, repertory, results, qlog, name = name)
+
+        return results

src/View/BoundaryCondition/Edit/Plot.py

@@ -16,6 +16,8 @@
 
 # -*- coding: utf-8 -*-
 
+import logging
+
 from datetime import datetime
 
 from tools import timer, trace
@@ -29,6 +31,8 @@ from View.BoundaryCondition.Edit.translate import *
 
 _translate = QCoreApplication.translate
 
+logger = logging.getLogger()
+
 class Plot(APlot):
     def __init__(self, canvas=None, data=None,
                  mode = "time", toolbar=None):
@@ -41,6 +45,9 @@ class Plot(APlot):
         self._mode = mode
 
     def custom_ticks(self):
+        if self.data.header[0] != "time":
+            return
+
         t0 = datetime.fromtimestamp(0)
         nb = len(self.data.data)
         mod = int(nb / 5)

src/View/BoundaryCondition/Edit/Table.py

@@ -25,6 +25,7 @@ from tools import trace, timer
 
 from View.ASubWindow import ASubMainWindow, AWidget
 from View.ListedSubWindow import ListedSubWindow
+from View.Tools.PamhyrTable import PamhyrTableModel
 
 from PyQt5.QtCore import (
     Qt, QVariant, QAbstractTableModel,
@@ -53,125 +54,7 @@ _translate = QCoreApplication.translate
 
 logger = logging.getLogger()
 
-class ExtendedTimeEdit(AWidget):
+class TableModel(PamhyrTableModel):
-    def __init__(self, parent=None):
-        super(ExtendedTimeEdit, self).__init__(
-            ui="extendedTimeEdit",
-            parent=parent
-        )
-        self.parent = parent
-
-        self.spinBox_days = self.find(QSpinBox, "spinBox_days")
-        self.timeEdit = self.find(QTimeEdit, "timeEdit")
-
-    def set_time(self, time):
-        days = 0
-        stime = time
-
-        # if ',' in time, time format is 'DD days, HH:MM:SS',
-        # otherelse is 'HH:MM:SS'
-        if "," in time:
-            s = time.strip().split(" ")
-            days = int(s[0])
-            stime = s[-1]
-
-        qtime = QTime.fromString(
-            stime,
-            "h:mm:ss"
-        )
-        self.spinBox_days.setValue(days)
-        self.timeEdit.setTime(qtime)
-
-    def get_time(self):
-        days = self.spinBox_days.value()
-        time = self.timeEdit.time().toPyTime()
-        secs = (
-            (time.hour * 3600) +
-            (time.minute * 60) +
-            time.second
-        )
-
-        return timedelta(days=days, seconds=secs)
-
-class ExtendedDateTimeEdit(AWidget):
-    def __init__(self, parent=None):
-        super(ExtendedDateTimeEdit, self).__init__(
-            ui="extendedDateTimeEdit",
-            parent=parent
-        )
-        self.parent = parent
-
-        self.dateTimeEdit = self.find(QDateTimeEdit, "dateTimeEdit")
-
-    def set_time(self, time):
-        qdatetime = QDateTime.fromString(
-            time,
-            "yyyy-MM-dd hh:mm:ss"
-        )
-        self.dateTimeEdit.setDateTime(qdatetime)
-
-    def get_time(self):
-        time = self.dateTimeEdit.dateTime().toPyDateTime()
-        return time
-
-class ExTimeDelegate(QItemDelegate):
-    def __init__(self, data=None, mode="time", parent=None):
-        super(ExTimeDelegate, self).__init__(parent)
-
-        self._data = data
-        self._mode = mode
-
-    def createEditor(self, parent, option, index):
-        if self._mode == "time":
-            self.editor = ExtendedTimeEdit(parent=parent)
-        else:
-            self.editor = ExtendedDateTimeEdit(parent=parent)
-        value = index.data(Qt.DisplayRole)
-        self.editor.set_time(value)
-        logger.debug(str(value))
-        return self.editor
-
-    def setModelData(self, editor, model, index):
-        time = editor.get_time()
-        if self._mode == "time":
-            model.setData(index, int(time.total_seconds()))
-        else:
-            logger.debug(str(time.timestamp()))
-            model.setData(index, int(time.timestamp()))
-        editor.close()
-        editor.deleteLater()
-
-    def updateEditorGeometry(self, editor, option, index):
-        r = QRect(option.rect)
-        if self.editor.windowFlags() & Qt.Popup and editor.parent() is not None:
-            r.setTopLeft(self.editor.parent().mapToGlobal(r.topLeft()))
-        editor.setGeometry(r)
-
-    @pyqtSlot()
-    def currentItemChanged(self):
-        self.commitData.emit(self.sender())
-
-
-class TableModel(QAbstractTableModel):
-    def __init__(self, data=None, mode="time", undo=None):
-        super(QAbstractTableModel, self).__init__()
-        self._headers = data.header
-        self._data = data
-        self._mode = mode
-        self._undo = undo
-
-    def flags(self, index):
-        options = Qt.ItemIsEnabled | Qt.ItemIsSelectable
-        options |= Qt.ItemIsEditable
-
-        return options
-
-    def rowCount(self, parent):
-        return len(self._data)
-
-    def columnCount(self, parent):
-        return len(self._headers)
-
     def data(self, index, role):
         if role == Qt.TextAlignmentRole:
             return Qt.AlignHCenter | Qt.AlignVCenter
@@ -189,7 +72,7 @@ class TableModel(QAbstractTableModel):
             if self._data.get_type_column(column) == float:
                 value = f"{v:.4f}"
             elif self._data.header[column] == "time":
-                if self._mode == "time":
+                if self._opt_data == "time":
                     t0 = datetime.fromtimestamp(0)
                     t = datetime.fromtimestamp(v)
                     value = str(t - t0)
@@ -200,12 +83,6 @@ class TableModel(QAbstractTableModel):
 
         return value
 
-    def headerData(self, section, orientation, role):
-        if role == Qt.ItemDataRole.DisplayRole and orientation == Qt.Orientation.Horizontal:
-            return table_headers[self._headers[section]]
-
-        return QVariant()
-
     def setData(self, index, value, role=Qt.EditRole):
         if not index.isValid() or role != Qt.EditRole:
             return False
@@ -315,12 +192,3 @@ class TableModel(QAbstractTableModel):
 
         self.layoutAboutToBeChanged.emit()
         self.layoutChanged.emit()
-
-
-    def undo(self):
-        self._undo.undo()
-        self.layoutChanged.emit()
-
-    def redo(self):
-        self._undo.redo()
-        self.layoutChanged.emit()

src/View/BoundaryCondition/Edit/UndoCommand.py

@@ -16,6 +16,8 @@
 
 # -*- coding: utf-8 -*-
 
+import logging
+
 from copy import deepcopy
 from tools import trace, timer
 
@@ -25,6 +27,8 @@ from PyQt5.QtWidgets import (
 
 from Model.BoundaryCondition.BoundaryCondition import BoundaryCondition
 
+logger = logging.getLogger()
+
 class SetDataCommand(QUndoCommand):
     def __init__(self, data, index, column, new_value):
         QUndoCommand.__init__(self)
@@ -42,6 +46,32 @@ class SetDataCommand(QUndoCommand):
     def redo(self):
         self._data._set_i_c_v(self._index, self._column, self._new)
 
+class SetMetaDataCommand(QUndoCommand):
+    def __init__(self, data, column, new_value):
+        QUndoCommand.__init__(self)
+
+        self._data = data
+        self._column = column
+        if self._column == "d50":
+            self._old = self._data.d50
+        elif self._column == "sigma":
+            self._old = self._data.sigma
+
+        self._new = float(new_value)
+
+    def undo(self):
+        if self._column == "d50":
+            self._data.d50 = self._old
+        elif self._column == "sigma":
+            self._data.sigma = self._old
+
+    def redo(self):
+        if self._column == "d50":
+            self._data.d50 = self._new
+        elif self._column == "sigma":
+            self._data.sigma = self._new
+
+
 class AddCommand(QUndoCommand):
     def __init__(self, data, index):
         QUndoCommand.__init__(self)

src/View/BoundaryCondition/Edit/Window.py

@@ -16,34 +16,85 @@
 
 # -*- coding: utf-8 -*-
 
+import logging
+
 from tools import timer, trace
 
-from View.ASubWindow import ASubMainWindow
+from View.ASubWindow import ASubMainWindow, AWidget
 from View.ListedSubWindow import ListedSubWindow
+from View.Tools.PamhyrDelegate import PamhyrExTimeDelegate
 
 from PyQt5.QtGui import (
     QKeySequence,
 )
 
+from PyQt5 import QtCore
 from PyQt5.QtCore import (
     Qt, QVariant, QAbstractTableModel, QCoreApplication,
+    pyqtSlot, pyqtSignal,
 )
 
 from PyQt5.QtWidgets import (
     QDialogButtonBox, QPushButton, QLineEdit,
     QFileDialog, QTableView, QAbstractItemView,
     QUndoStack, QShortcut, QAction, QItemDelegate,
-    QHeaderView,
+    QHeaderView, QDoubleSpinBox, QVBoxLayout,
 )
 
 from View.Plot.MplCanvas import MplCanvas
+from View.Plot.PamhyrToolbar import PamhyrPlotToolbar
 
 from View.BoundaryCondition.translate import long_types
-from View.BoundaryCondition.Edit.Table import TableModel, ExTimeDelegate
+from View.BoundaryCondition.Edit.translate import table_headers, retranslate
+from View.BoundaryCondition.Edit.UndoCommand import SetMetaDataCommand
+from View.BoundaryCondition.Edit.Table import TableModel
 from View.BoundaryCondition.Edit.Plot import Plot
 
 _translate = QCoreApplication.translate
 
+logger = logging.getLogger()
+
+class WD50Sigma(AWidget):
+    d50Changed = pyqtSignal(float)
+    sigmaChanged = pyqtSignal(float)
+
+    def __init__(self, parent=None):
+        super(WD50Sigma, self).__init__(
+            ui="d50sigma",
+            parent=parent
+        )
+        self.parent = parent
+
+        self.spinBox_d50 = self.find(QDoubleSpinBox, "doubleSpinBox_d50")
+        self.spinBox_sigma = self.find(QDoubleSpinBox, "doubleSpinBox_sigma")
+
+        self.spinBox_d50.valueChanged.connect(self.valueChangedD50)
+        self.spinBox_sigma.valueChanged.connect(self.valueChangedSigma)
+
+    def set_d50(self, d50):
+        self.spinBox_d50.valueChanged.disconnect(self.valueChangedD50)
+        self.spinBox_d50.setValue(float(d50))
+        self.spinBox_d50.valueChanged.connect(self.valueChangedD50)
+
+    def get_d50(self):
+        return float(self.spinBox_d50.value())
+
+    def set_sigma(self, sigma):
+        self.spinBox_sigma.valueChanged.disconnect(self.valueChangedSigma)
+        self.spinBox_sigma.setValue(float(sigma))
+        self.spinBox_sigma.valueChanged.connect(self.valueChangedSigma)
+
+    def get_sigma(self):
+        return float(self.spinBox_sigma.value())
+
+    @QtCore.pyqtSlot(float)
+    def valueChangedD50(self, value):
+        self.d50Changed.emit(value)
+
+    @QtCore.pyqtSlot(float)
+    def valueChangedSigma(self, value):
+        self.sigmaChanged.emit(value)
+
 class EditBoundaryConditionWindow(ASubMainWindow, ListedSubWindow):
     def __init__(self, title="Edit boundary condition",
                  data=None, study=None, parent=None):
@@ -62,6 +113,7 @@ class EditBoundaryConditionWindow(ASubMainWindow, ListedSubWindow):
         self.setup_sc()
         self.setup_table()
         self.setup_plot()
+        self.setup_data()
         self.setup_connections()
 
     def compute_title(self):
@@ -83,24 +135,41 @@ class EditBoundaryConditionWindow(ASubMainWindow, ListedSubWindow):
         self.copy_sc = QShortcut(QKeySequence.Copy, self)
         self.paste_sc = QShortcut(QKeySequence.Paste, self)
 
+    def setup_data(self):
+        self._is_solid = self._data.bctype == "SL"
+
+        if self._is_solid:
+            layout = self.find(QVBoxLayout, "verticalLayout_table")
+            self._d50sigma = WD50Sigma(parent = self)
+            layout.addWidget(self._d50sigma)
+
+            self._d50sigma.set_d50(self._data.d50)
+            self._d50sigma.set_sigma(self._data.sigma)
+
     def setup_table(self):
-        table = self.find(QTableView, "tableView")
+        retranslate()
-        self._table = TableModel(
+        headers = {}
+        for h in self._data.header:
+            headers[h] = table_headers[h]
+
+        self._delegate_time = PamhyrExTimeDelegate(
             data = self._data,
-            undo = self._undo_stack,
+            mode = self._study.time_system,
-            mode = self._study.time_system
+            parent = self
         )
 
-        if self._data.header[0] == "time":
+        table = self.find(QTableView, "tableView")
-            self._delegate_time = ExTimeDelegate(
+        self._table = TableModel(
-                data = self._data,
+            table_view = table,
-                mode = self._study.time_system,
+            table_headers = headers,
-                parent = self
+            editable_headers = self._data.header,
-            )
+            delegates = {
-
+                "time": self._delegate_time,
-            table.setItemDelegateForColumn(
+            },
-                0, self._delegate_time
+            data = self._data,
-            )
+            undo = self._undo_stack,
+            opt_data = self._study.time_system
+        )
 
         table.setModel(self._table)
         table.setSelectionBehavior(QAbstractItemView.SelectRows)
@@ -110,12 +179,17 @@ class EditBoundaryConditionWindow(ASubMainWindow, ListedSubWindow):
     def setup_plot(self):
         self.canvas = MplCanvas(width=5, height=4, dpi=100)
         self.canvas.setObjectName("canvas")
+        self.toolbar = PamhyrPlotToolbar(
+            self.canvas, self
+        )
+        self.verticalLayout.addWidget(self.toolbar)
         self.verticalLayout.addWidget(self.canvas)
 
         self.plot = Plot(
             canvas = self.canvas,
             data = self._data,
             mode = self._study.time_system,
+            toolbar = self.toolbar,
         )
         self.plot.draw()
 
@@ -132,6 +206,31 @@ class EditBoundaryConditionWindow(ASubMainWindow, ListedSubWindow):
 
         self._table.dataChanged.connect(self.update)
 
+        if self._is_solid:
+            self._d50sigma.d50Changed.connect(self.d50_changed)
+            self._d50sigma.sigmaChanged.connect(self.sigma_changed)
+
+    def d50_changed(self, value):
+        self._undo_stack.push(
+            SetMetaDataCommand(
+                self._data,
+                "d50", value
+            )
+        )
+
+    def sigma_changed(self, value):
+        self._undo_stack.push(
+            SetMetaDataCommand(
+                self._data,
+                "sigma", value
+            )
+        )
+
+    def widget_update(self):
+        if self._is_solid:
+            self._d50sigma.set_d50(self._data.d50)
+            self._d50sigma.set_sigma(self._data.sigma)
+
     def update(self):
         self.plot.update()
 
@@ -215,7 +314,9 @@ class EditBoundaryConditionWindow(ASubMainWindow, ListedSubWindow):
     def undo(self):
         self._table.undo()
         self.plot.update()
+        self.widget_update()
 
     def redo(self):
         self._table.redo()
         self.plot.update()
+        self.widget_update()

src/View/BoundaryCondition/Edit/translate.py

@@ -26,5 +26,15 @@ table_headers = {
     "time": _translate("BoundaryCondition", "Time"),
     "date": _translate("BoundaryCondition", "Date"),
     "discharge": _translate("BoundaryCondition", "Discharge (m³/s)"),
-    "z": _translate("BoundaryCondition", "Z (m)")
+    "z": _translate("BoundaryCondition", "Z (m)"),
+    "solid": _translate("BoundaryCondition", "Solid (kg/s)"),
 }
+
+def retranslate():
+    table_headers["x"] = _translate("BoundaryCondition", "X")
+    table_headers["y"] = _translate("BoundaryCondition", "Y")
+    table_headers["time"] = _translate("BoundaryCondition", "Time")
+    table_headers["date"] = _translate("BoundaryCondition", "Date")
+    table_headers["discharge"] = _translate("BoundaryCondition", "Discharge (m³/s)")
+    table_headers["z"] = _translate("BoundaryCondition", "Z (m)")
+    table_headers["solid"] = _translate("BoundaryCondition", "Solid (kg/s)")

src/View/BoundaryCondition/Table.py

@@ -39,6 +39,8 @@ from Model.BoundaryCondition.BoundaryConditionTypes import (
     TimeOverZ, TimeOverDischarge, ZOverDischarge
 )
 
+from View.Tools.PamhyrTable import PamhyrTableModel
+
 from View.BoundaryCondition.UndoCommand import (
     SetNameCommand, SetNodeCommand, SetTypeCommand,
     AddCommand, DelCommand, SortCommand,
@@ -104,26 +106,13 @@ class ComboBoxDelegate(QItemDelegate):
         self.commitData.emit(self.sender())
 
 
-class TableModel(QAbstractTableModel):
+class TableModel(PamhyrTableModel):
-    def __init__(self, data=None, undo=None, tab=""):
+    def _setup_lst(self):
-        super(QAbstractTableModel, self).__init__()
+        self._lst = self._data.boundary_condition
-        self._headers = list(table_headers.keys())
+        self._tab = self._opt_data
-        self._data = data
-        self._undo = undo
-        self._tab = tab
-        self._bcs = self._data.boundary_condition
-
-    def flags(self, index):
-        options = Qt.ItemIsEnabled | Qt.ItemIsSelectable
-        options |= Qt.ItemIsEditable
-
-        return options
 
     def rowCount(self, parent):
-        return self._bcs.len(self._tab)
+        return self._lst.len(self._tab)
-
-    def columnCount(self, parent):
-        return len(self._headers)
 
     def data(self, index, role):
         if role != Qt.ItemDataRole.DisplayRole:
@@ -133,24 +122,18 @@ class TableModel(QAbstractTableModel):
         column = index.column()
 
         if self._headers[column] == "name":
-            return self._bcs.get(self._tab, row).name
+            return self._lst.get(self._tab, row).name
         elif self._headers[column] == "type":
-            t = self._bcs.get(self._tab, row).bctype
+            t = self._lst.get(self._tab, row).bctype
             return long_types[t]
         elif self._headers[column] == "node":
-            n = self._bcs.get(self._tab, row).node
+            n = self._lst.get(self._tab, row).node
             if n is None:
                 return _translate("BoundaryCondition", "Not associate")
             return n.name
 
         return QVariant()
 
-    def headerData(self, section, orientation, role):
-        if role == Qt.ItemDataRole.DisplayRole and orientation == Qt.Orientation.Horizontal:
-            return table_headers[self._headers[section]]
-
-        return QVariant()
-
     def setData(self, index, value, role=Qt.EditRole):
         if not index.isValid() or role != Qt.EditRole:
             return False
@@ -162,20 +145,20 @@ class TableModel(QAbstractTableModel):
             if self._headers[column] == "name":
                 self._undo.push(
                     SetNameCommand(
-                        self._bcs, self._tab,row, value
+                        self._lst, self._tab,row, value
                     )
                 )
             elif self._headers[column] == "type":
                 key = next(k for k, v in long_types.items() if v == value)
                 self._undo.push(
                     SetTypeCommand(
-                        self._bcs, self._tab,row, BC_types[key]
+                        self._lst, self._tab,row, BC_types[key]
                     )
                 )
             elif self._headers[column] == "node":
                 self._undo.push(
                     SetNodeCommand(
-                        self._bcs, self._tab,row, self._data.node(value)
+                        self._lst, self._tab,row, self._data.node(value)
                     )
                 )
         except Exception as e:
@@ -190,7 +173,7 @@ class TableModel(QAbstractTableModel):
 
         self._undo.push(
             AddCommand(
-                self._bcs, self._tab,row
+                self._lst, self._tab,row
             )
         )
 
@@ -202,7 +185,7 @@ class TableModel(QAbstractTableModel):
 
         self._undo.push(
             DelCommand(
-                self._bcs, self._tab,rows
+                self._lst, self._tab, rows
             )
         )
 
@@ -214,7 +197,7 @@ class TableModel(QAbstractTableModel):
 
         self._undo.push(
             SortCommand(
-                self._bcs, self._tab,False
+                self._lst, self._tab, False
             )
         )
 
@@ -231,7 +214,7 @@ class TableModel(QAbstractTableModel):
 
         self._undo_stack.push(
             MoveCommand(
-                self._bcs, self._tab,"up", row
+                self._lst, self._tab, "up", row
             )
         )
 
@@ -239,7 +222,7 @@ class TableModel(QAbstractTableModel):
         self.layoutChanged.emit()
 
     def move_down(self, index, parent=QModelIndex()):
-        if row > len(self._bcs):
+        if row > len(self._lst):
             return
 
         target = row
@@ -248,7 +231,7 @@ class TableModel(QAbstractTableModel):
 
         self._undo_stack.push(
             MoveCommand(
-                self._bcs, self._tab,"down", row
+                self._lst, self._tab, "down", row
             )
         )
 

src/View/BoundaryCondition/Window.py

@@ -56,7 +56,7 @@ from View.BoundaryCondition.Table import (
 )
 
 from View.Network.GraphWidget import GraphWidget
-from View.BoundaryCondition.translate import *
+from View.BoundaryCondition.translate import table_headers, retranslate
 from View.BoundaryCondition.Edit.Window import EditBoundaryConditionWindow
 
 _translate = QCoreApplication.translate
@@ -90,17 +90,10 @@ class BoundaryConditionWindow(ASubMainWindow, ListedSubWindow):
         self.paste_sc = QShortcut(QKeySequence.Paste, self)
 
     def setup_table(self):
+        retranslate()
         self._table = {}
 
         for t in ["liquid", "solid", "suspenssion"]:
-            table = self.find(QTableView, f"tableView_{t}")
-            self._table[t] = TableModel(
-                data = self._study.river,
-                undo = self._undo_stack,
-                tab = t,
-            )
-            table.setModel(self._table[t])
-
             self._delegate_type = ComboBoxDelegate(
                 data = self._study.river,
                 mode = "type",
@@ -114,13 +107,20 @@ class BoundaryConditionWindow(ASubMainWindow, ListedSubWindow):
                 parent=self
             )
 
-            table.setItemDelegateForColumn(
+            table = self.find(QTableView, f"tableView_{t}")
-                1, self._delegate_type
+            self._table[t] = TableModel(
+                table_view = table,
+                table_headers = table_headers,
+                editable_headers = ["name", "type", "node"],
+                delegates = {
+                    "type": self._delegate_type,
+                    "node": self._delegate_node,
+                },
+                data = self._study.river,
+                undo = self._undo_stack,
+                opt_data = t,
             )
-            table.setItemDelegateForColumn(
+            table.setModel(self._table[t])
-                2, self._delegate_node
-            )
-
             table.setSelectionBehavior(QAbstractItemView.SelectRows)
             table.horizontalHeader().setSectionResizeMode(QHeaderView.Stretch)
             table.setAlternatingRowColors(True)

src/View/BoundaryCondition/translate.py

@@ -20,7 +20,8 @@ from PyQt5.QtCore import QCoreApplication
 
 from Model.BoundaryCondition.BoundaryConditionTypes import (
     NotDefined, PonctualContribution,
-    TimeOverZ, TimeOverDischarge, ZOverDischarge
+    TimeOverZ, TimeOverDischarge, ZOverDischarge,
+    Solid
 )
 
 _translate = QCoreApplication.translate
@@ -28,9 +29,10 @@ _translate = QCoreApplication.translate
 long_types = {
     "ND": _translate("BoundaryCondition", "Not defined"),
     "PC": _translate("BoundaryCondition", "Ponctual contribution"),
-    "TZ": _translate("BoundaryCondition", "Time over Z"),
+    "TZ": _translate("BoundaryCondition", "TimeZ"),
-    "TD": _translate("BoundaryCondition", "Time over Discharge"),
+    "TD": _translate("BoundaryCondition", "TimeDischarge"),
-    "ZD": _translate("BoundaryCondition", "Z over Discharge"),
+    "ZD": _translate("BoundaryCondition", "ZDischarge"),
+    "SL": _translate("BoundaryCondition", "Solid"),
 }
 
 table_headers = {
@@ -44,5 +46,11 @@ BC_types = {
     "PC": PonctualContribution,
     "TZ": TimeOverZ,
     "TD": TimeOverDischarge,
-    "ZD": ZOverDischarge
+    "ZD": ZOverDischarge,
+    "SL": Solid,
 }
+
+def retranslate():
+    table_headers["name"] = _translate("BoundaryCondition", "Name")
+    table_headers["type"] = _translate("BoundaryCondition", "Type")
+    table_headers["node"] = _translate("BoundaryCondition", "Node")

src/View/CheckList/Worker.py

@@ -37,10 +37,11 @@ class Worker(QObject):
         self.signalStatus.emit('end')
 
     def _compute(self):
+        # The computation loop
         for checker in self._checker_list:
             self.signalStatus.emit(checker.name)
 
-            # time.sleep(1)
+            # Run checker
             checker.run(self._study)
 
             self.signalStatus.emit("progress")

src/View/Debug/Window.py

@@ -92,7 +92,9 @@ class ReplWindow(ASubMainWindow, ListedSubWindow):
         self.set_plaintext_edit_text("plainTextEdit", "")
         self._history_ind = 0
 
+        # Code to rich_code
         rich_code = code.strip().split("\n")
+        # Add return variable for results display at last row
         rich_code[-1] = "self.__debug_exec_result__ = " + rich_code[-1]
         rich_code = "\n".join(rich_code)
         logger.debug(f"User debug command : {rich_code}")
@@ -102,6 +104,8 @@ class ReplWindow(ASubMainWindow, ListedSubWindow):
         except Exception as e:
             value = f"<font color=\"red\">" + str(e) + "</font>"
 
+        # Display code
         msg = f"<font color=\"grey\"> # " + code + " #</font>"
         self.find(QTextEdit, "textEdit").append(msg)
+        # Display results
         self.find(QTextEdit, "textEdit").append(str(value))

src/View/Frictions/Window.py

@@ -49,9 +49,9 @@ from View.Frictions.Table import (
 )
 
 from View.Plot.MplCanvas import MplCanvas
-from View.Geometry.PlotKPC import PlotKPC
+from View.Geometry.PlotKPZ import PlotKPZ
 from View.Frictions.PlotStricklers import PlotStricklers
-from View.Frictions.translate import *
+from View.Frictions.translate import table_headers, retranslate
 
 from View.Stricklers.Window import StricklersWindow
 
@@ -94,6 +94,7 @@ class FrictionsWindow(ASubMainWindow, ListedSubWindow):
         self.paste_sc = QShortcut(QKeySequence.Paste, self)
 
     def setup_table(self):
+        retranslate()
         self._table = {}
 
         table = self.find(QTableView, f"tableView")
@@ -128,7 +129,7 @@ class FrictionsWindow(ASubMainWindow, ListedSubWindow):
         self.plot_layout = self.find(QVBoxLayout, "verticalLayout")
         self.plot_layout.addWidget(self.canvas)
 
-        self.plot = PlotKPC(
+        self.plot = PlotKPZ(
             canvas = self.canvas,
             data = self._reach.reach,
             toolbar = None,
@@ -193,7 +194,7 @@ class FrictionsWindow(ASubMainWindow, ListedSubWindow):
             sec = self._frictions.get(rows[0])
             highlight = (sec.begin_kp, sec.end_kp)
 
-        self.plot = PlotKPC(
+        self.plot = PlotKPZ(
             canvas = self.canvas,
             data = reach,
             toolbar = None,

src/View/Frictions/translate.py

@@ -28,3 +28,11 @@ table_headers = {
     "begin_strickler": _translate("Frictions", "Begin strickler"),
     "end_strickler": _translate("Frictions", "End strickler"),
 }
+
+def retranslate():
+    table_headers["name"] = _translate("Frictions", "Name")
+    # table_headers["edge"] = _translate("Frictions", "Reach")
+    table_headers["begin_kp"] = _translate("Frictions", "Begin kp (m)")
+    table_headers["end_kp"] = _translate("Frictions", "End kp (m)")
+    table_headers["begin_strickler"] = _translate("Frictions", "Begin strickler")
+    table_headers["end_strickler"] = _translate("Frictions", "End strickler")

src/View/Geometry/PlotAC.py

@@ -67,18 +67,18 @@ class PlotAC(APlot):
         gl = self.data.profile(selected_profile).names()
 
         self.canvas.axes.set_xlabel(
-            _translate("MainWindow_reach", "Abscisse en travers (m)"),
+            _translate("MainWindow_reach", "Transverse abscissa (m)"),
             color='green', fontsize=12
         )
         self.canvas.axes.set_ylabel(
-            _translate("MainWindow_reach", "Cote (m)"),
+            _translate("MainWindow_reach", "Height (m)"),
             color='green', fontsize=12
         )
         self.canvas.figure.tight_layout()
 
-        label_before_plot_selected = _translate("MainWindow_reach", "Profil précédent")
+        label_before_plot_selected = _translate("MainWindow_reach", "Previous cross-section")
-        label_plot_selected = _translate("MainWindow_reach", "Profil sélectionné")
+        label_plot_selected = _translate("MainWindow_reach", "Cross-section")
-        label_after_plot_selected = _translate("MainWindow_reach", "Profil suivant")
+        label_after_plot_selected = _translate("MainWindow_reach", "Next cross-section")
         color_before_plot_selected = "k"  # 'grey'
         color_plot_selected = 'b'
         color_after_plot_selected = 'm'

src/View/Geometry/PlotKPZ.py

@@ -29,10 +29,10 @@ _translate = QCoreApplication.translate
 
 logger = logging.getLogger()
 
-class PlotKPC(APlot):
+class PlotKPZ(APlot):
     def __init__(self, canvas=None, data=None, toolbar=None,
                  display_current=True):
-        super(PlotKPC, self).__init__(
+        super(PlotKPZ, self).__init__(
             canvas=canvas,
             data=data,
             toolbar=toolbar
@@ -64,7 +64,7 @@ class PlotKPC(APlot):
             color='green', fontsize=12
         )
         self.canvas.axes.set_ylabel(
-            _translate("MainWindow_reach", "Cote (m)"),
+            _translate("MainWindow_reach", "Height (m)"),
             color='green', fontsize=12
         )
 
@@ -234,7 +234,7 @@ class PlotKPC(APlot):
                         kp, z_complete[i]
                     )
             except Exception as e:
-                logger.warning(f"Failed to update graphic KPC: {e}")
+                logger.warning(f"Failed to update graphic KPZ: {e}")
 
             self.canvas.axes.autoscale_view(True, True, True)
             self.canvas.figure.canvas.draw_idle()

src/View/Geometry/Profile/mainwindow_ui_profile.py

@@ -25,7 +25,8 @@ from PyQt5.QtWidgets import (
     QToolBar,
 )
 
-from View.Plot import MplCanvas, navigation_toolbar_2qt
+from View.Plot import MplCanvas
+from View.Plot.PamhyrToolbar import PamhyrPlotToolbar
 
 file_path = os.path.abspath(os.path.dirname(__file__))
 icon_path = f"{os.path.dirname(os.path.dirname(file_path))}/../View/ui/ressources/"
@@ -197,7 +198,7 @@ class Ui_MainWindow(object):
         self.canvas = MplCanvas.MplCanvas(width=5, height=4, dpi=100)
         self.canvas.setObjectName("canvas")
 
-        self.my_toolbar = navigation_toolbar_2qt.PamHyrNavigationToolbar2QT(
+        self.my_toolbar = PamhyrPlotToolbar(
             self.canvas, self.widget1
         )
 

src/View/Geometry/Window.py

@@ -38,7 +38,7 @@ from PyQt5.QtWidgets import (
 )
 
 from View.Geometry.PlotXY import PlotXY
-from View.Geometry.PlotKPC import PlotKPC
+from View.Geometry.PlotKPZ import PlotKPZ
 from View.Geometry.PlotAC import PlotAC
 
 from View.ASubWindow import ASubMainWindow, WindowToolKit
@@ -111,7 +111,7 @@ class GeometryWindow(ASubMainWindow, ListedSubWindow):
         self.ui.btn_sort_desc.triggered.connect(self.sort_descending)
         self.ui.btn_move_up.triggered.connect(self.move_row_up)
         self.ui.btn_move_down.triggered.connect(self.move_row_down)
-        self.ui.btn_end_editing.triggered.connect(self.handleSave)
+        # self.ui.btn_end_editing.triggered.connect(self.handleSave)
         self.ui.btn_add.triggered.connect(self.insert_row)
         self.ui.btn_delete.triggered.connect(self.delete_rows)
         self.ui.btn_edit.triggered.connect(self.edit_profile)
@@ -230,7 +230,7 @@ class GeometryWindow(ASubMainWindow, ListedSubWindow):
     def plot_kpc(self):
         self.tableView.model().blockSignals(True)
 
-        self._plot_kpc = PlotKPC(
+        self._plot_kpc = PlotKPZ(
             canvas = self.ui.canvas_2,
             data = self._reach,
             toolbar = self.ui.toolbar_2

src/View/Geometry/mainwindow_ui_reach.py

@@ -25,9 +25,9 @@ from PyQt5.QtWidgets import (
     QAbstractItemView, QHeaderView, QToolBar,
 )
 
-from View.Plot import (
+from View.Plot import MplCanvas
-    MplCanvas, navigation_toolbar_2qt
+from View.Plot.PamhyrToolbar import PamhyrPlotToolbar
-)
+
 
 file_path = os.path.abspath(os.path.dirname(__file__))
 icon_path = f"{os.path.dirname(os.path.dirname(file_path))}/View/ui/ressources/"
@@ -95,29 +95,29 @@ class Ui_MainWindow(object):
         self.btn_edit.setIcon(icon_btn_edit)
         self.toolBar.addAction(self.btn_edit)
 
-        self.btn_copy_selected_profile = QtWidgets.QAction(self.widget)
+        # self.btn_copy_selected_profile = QtWidgets.QAction(self.widget)
-        self.btn_copy_selected_profile.setObjectName("btn_copy_selected_profile")
+        # self.btn_copy_selected_profile.setObjectName("btn_copy_selected_profile")
-        icon_btn_copy_selected_profile = QtGui.QIcon()
+        # icon_btn_copy_selected_profile = QtGui.QIcon()
-        icon_btn_copy_selected_profile.addPixmap(
+        # icon_btn_copy_selected_profile.addPixmap(
-            QtGui.QPixmap(f"{icon_path}gtk-copy.png"))
+        #     QtGui.QPixmap(f"{icon_path}gtk-copy.png"))
-        self.btn_copy_selected_profile.setIcon(icon_btn_copy_selected_profile)
+        # self.btn_copy_selected_profile.setIcon(icon_btn_copy_selected_profile)
-        self.toolBar.addAction(self.btn_copy_selected_profile)
+        # self.toolBar.addAction(self.btn_copy_selected_profile)
 
-        self.btn_paste_selected_profile = QtWidgets.QAction(self.widget)
+        # self.btn_paste_selected_profile = QtWidgets.QAction(self.widget)
-        self.btn_paste_selected_profile.setObjectName("btn_paste_selected_profile")
+        # self.btn_paste_selected_profile.setObjectName("btn_paste_selected_profile")
-        icon_btn_paste_selected_profile = QtGui.QIcon()
+        # icon_btn_paste_selected_profile = QtGui.QIcon()
-        icon_btn_paste_selected_profile.addPixmap(
+        # icon_btn_paste_selected_profile.addPixmap(
-            QtGui.QPixmap(f"{icon_path}gtk-paste.png"))
+        #     QtGui.QPixmap(f"{icon_path}gtk-paste.png"))
-        self.btn_paste_selected_profile.setIcon(icon_btn_paste_selected_profile)
+        # self.btn_paste_selected_profile.setIcon(icon_btn_paste_selected_profile)
-        self.toolBar.addAction(self.btn_paste_selected_profile)
+        # self.toolBar.addAction(self.btn_paste_selected_profile)
 
-        self.btn_duplicate_selected_profile = QtWidgets.QAction(self.widget)
+        # self.btn_duplicate_selected_profile = QtWidgets.QAction(self.widget)
-        self.btn_duplicate_selected_profile.setObjectName("btn_duplicate_selected_profile")
+        # self.btn_duplicate_selected_profile.setObjectName("btn_duplicate_selected_profile")
-        icon_btn_duplicate_selected_profile = QtGui.QIcon()
+        # icon_btn_duplicate_selected_profile = QtGui.QIcon()
-        icon_btn_duplicate_selected_profile.addPixmap(
+        # icon_btn_duplicate_selected_profile.addPixmap(
-            QtGui.QPixmap(f"{icon_path}gtk-dnd-multiple.png"))
+        #     QtGui.QPixmap(f"{icon_path}gtk-dnd-multiple.png"))
-        self.btn_duplicate_selected_profile.setIcon(icon_btn_duplicate_selected_profile)
+        # self.btn_duplicate_selected_profile.setIcon(icon_btn_duplicate_selected_profile)
-        self.toolBar.addAction(self.btn_duplicate_selected_profile)
+        # self.toolBar.addAction(self.btn_duplicate_selected_profile)
 
         self.btn_sort_asc = QtWidgets.QAction(self.widget)
         self.btn_sort_asc.setObjectName("btn_sort_asc")
@@ -151,13 +151,13 @@ class Ui_MainWindow(object):
         self.btn_move_down.setIcon(icon_btn_move_down)
         self.toolBar.addAction(self.btn_move_down)
 
-        self.btn_end_editing = QtWidgets.QAction(self.widget)
+        # self.btn_end_editing = QtWidgets.QAction(self.widget)
-        self.btn_end_editing.setObjectName("btn_end_editing")
+        # self.btn_end_editing.setObjectName("btn_end_editing")
-        icon_btn_save = QtGui.QIcon()
+        # icon_btn_save = QtGui.QIcon()
-        icon_btn_save.addPixmap(
+        # icon_btn_save.addPixmap(
-            QtGui.QPixmap(f"{icon_path}gtk-apply.png"))
+        #     QtGui.QPixmap(f"{icon_path}gtk-apply.png"))
-        self.btn_end_editing.setIcon(icon_btn_save)
+        # self.btn_end_editing.setIcon(icon_btn_save)
-        self.toolBar.addAction(self.btn_end_editing)
+        # self.toolBar.addAction(self.btn_end_editing)
 
         self.verticalLayout.addLayout(self.horizontalLayout)
         self.tableView = QtWidgets.QTableView(self.widget)
@@ -190,8 +190,9 @@ class Ui_MainWindow(object):
         self.canvas_1 = MplCanvas.MplCanvas(width=5, height=4, dpi=100)
         self.canvas_1.setObjectName("canvas_1")
 
-        self.toolbar_1 = navigation_toolbar_2qt.PamHyrNavigationToolbar2QT(
+        self.toolbar_1 = PamhyrPlotToolbar(
-            self.canvas_1, self.centralwidget
+            self.canvas_1, self.centralwidget,
+            items = ["home", "zoom", "save", "iso", "back/forward", "move"]
         )
 
         self.verticalLayout_canvas_1.addWidget(self.toolbar_1)
@@ -206,8 +207,9 @@ class Ui_MainWindow(object):
         self.canvas_2 = MplCanvas.MplCanvas(width=5, height=4, dpi=100)
         self.canvas_2.setObjectName("canvas_2")
 
-        self.toolbar_2 = navigation_toolbar_2qt.PamHyrNavigationToolbar2QT(
+        self.toolbar_2 = PamhyrPlotToolbar(
-            self.canvas_2, self.centralwidget
+            self.canvas_2, self.centralwidget,
+            items = ["home", "zoom", "save", "iso", "back/forward", "move"]
         )
 
         self.verticalLayout_canvas_2.addWidget(self.toolbar_2)
@@ -221,8 +223,9 @@ class Ui_MainWindow(object):
         self.canvas_3 = MplCanvas.MplCanvas(width=5, height=4, dpi=100)
         self.canvas_3.setObjectName("canvas_3")
 
-        self.toolbar_3 = navigation_toolbar_2qt.PamHyrNavigationToolbar2QT(
+        self.toolbar_3 = PamhyrPlotToolbar(
-            self.canvas_3, self.centralwidget
+            self.canvas_3, self.centralwidget,
+            items = ["home", "zoom", "save", "iso", "back/forward", "move"]
         )
 
         self.verticalLayout_canvas_3.addWidget(self.toolbar_3)
@@ -311,30 +314,23 @@ class Ui_MainWindow(object):
     def retranslateUi(self, MainWindow_reach):
         _translate = QtCore.QCoreApplication.translate
         self.name_reach = ""
-        self.mainwindow_title = _translate("MainWindow_reach", "Jeu de sections du Bief")
+        self.mainwindow_title = _translate("MainWindow_reach", "Geometry reach edit")
         MainWindow_reach.setWindowTitle(self.mainwindow_title + f"{self.name_reach}")
-        self.text_section_order = _translate("MainWindow_reach", "\nOrdre des sections : "
-                                                                 "Amont --> Aval")
         # self.label_section_order.setText(f"{self.text_section_order}")
         self.vertical_slider_value = ""
-        self.vertical_slider_label.setText(_translate("MainWindow_reach", "Pk = ") + f"{self.vertical_slider_value}")
+        self.vertical_slider_label.setText(_translate("MainWindow_reach", "Kp = ") + f"{self.vertical_slider_value}")
-        self.btn_open.setToolTip(_translate("MainWindow_reach", " Importer une géométrie"))
+        self.btn_open.setToolTip(_translate("MainWindow_reach", " Import geometry"))
-        self.btn_add.setToolTip(_translate("MainWindow_reach", " Nouveau profil"))
+        self.btn_add.setToolTip(_translate("MainWindow_reach", " New profile"))
-        self.btn_delete.setToolTip(_translate("MainWindow_reach", " Supprimer le profil sélectionné"))
+        self.btn_delete.setToolTip(_translate("MainWindow_reach", " Delete selected profile"))
-        self.btn_edit.setToolTip(_translate("MainWindow_reach", " Éditer le profil sélectionné"))
+        self.btn_edit.setToolTip(_translate("MainWindow_reach", " Editer selected profile"))
-        self.btn_copy_selected_profile.setToolTip(_translate("MainWindow_reach", " Copier le profil sélectionné"))
+
-        self.btn_paste_selected_profile.setToolTip(_translate("MainWindow_reach", "Coller le profil en fin de liste ("
-                                                                                  "penser à modifier le Pk avant de "
-                                                                                  "trier)"))
-        self.btn_duplicate_selected_profile.setToolTip(
-            _translate("MainWindow_reach", " Dupliquer la section sélectionnée"))
         self.btn_sort_asc.setToolTip(_translate("MainWindow_reach", " Trier les profils par ordre croissant des Pk"))
         self.btn_sort_desc.setToolTip(_translate("MainWindow_reach", " Trier les profils par ordre décroissant des Pk"))
         self.btn_move_up.setToolTip(_translate("MainWindow_reach", " Changer l'ordre des profils (en décalant le "
                                                                    "profil sélectionné vers le haut)"))
         self.btn_move_down.setToolTip(_translate("MainWindow_reach", " Changer l'ordre des profils (en décalant le "
                                                                      "profil sélectionné vers le bas)"))
-        self.btn_end_editing.setToolTip(_translate("MainWindow_reach", " Terminer l'édition"))
+        # self.btn_end_editing.setToolTip(_translate("MainWindow_reach", " Terminer l'édition"))
         self.tableView_header = [_translate("MainWindow_reach", "Name"),
                                  _translate("MainWindow_reach", "Kp (m)"),
                                  _translate("MainWindow_reach", "Type")]

src/View/InitialConditions/PlotDKP.py

@@ -42,7 +42,7 @@ class PlotDKP(APlot):
             return
 
         self.canvas.axes.set_ylabel(
-            _translate("MainWindow_reach", "Draft (m)"),
+            _translate("MainWindow_reach", "Elevation (m)"),
             color='green', fontsize=11
         )
         self.canvas.axes.set_xlabel(

src/View/InitialConditions/Window.py

@@ -49,9 +49,11 @@ from View.InitialConditions.UndoCommand import (
 from View.InitialConditions.Table import TableModel, ComboBoxDelegate
 
 from View.Plot.MplCanvas import MplCanvas
+from View.Plot.PamhyrToolbar import PamhyrPlotToolbar
+
 from View.InitialConditions.PlotDKP import PlotDKP
 from View.InitialConditions.PlotDischarge import PlotDischarge
-from View.InitialConditions.translate import *
+from View.InitialConditions.translate import table_headers, retranslate
 from View.InitialConditions.DialogHeight import HeightDialog
 from View.InitialConditions.DialogDischarge import DischargeDialog
 
@@ -95,6 +97,8 @@ class InitialConditionsWindow(ASubMainWindow, ListedSubWindow):
         self.paste_sc = QShortcut(QKeySequence.Paste, self)
 
     def setup_table(self):
+        retranslate()
+
         table = self.find(QTableView, f"tableView")
         self._table = TableModel(
             river = self._study.river,
@@ -120,25 +124,33 @@ class InitialConditionsWindow(ASubMainWindow, ListedSubWindow):
     def setup_graph(self):
         self.canvas_1 = MplCanvas(width=5, height=4, dpi=100)
         self.canvas_1.setObjectName("canvas_1")
+        self.toolbar_1 = PamhyrPlotToolbar(
+            self.canvas_1, self
+        )
         self.plot_layout_1 = self.find(QVBoxLayout, "verticalLayout_1")
+        self.plot_layout_1.addWidget(self.toolbar_1)
         self.plot_layout_1.addWidget(self.canvas_1)
 
         self.plot_1 = PlotDKP(
             canvas = self.canvas_1,
             data = self._ics,
-            toolbar = None,
+            toolbar = self.toolbar_1,
         )
         self.plot_1.draw()
 
         self.canvas_2 = MplCanvas(width=5, height=4, dpi=100)
         self.canvas_2.setObjectName("canvas_2")
+        self.toolbar_2 = PamhyrPlotToolbar(
+            self.canvas_2, self
+        )
         self.plot_layout_2 = self.find(QVBoxLayout, "verticalLayout_2")
+        self.plot_layout_2.addWidget(self.toolbar_2)
         self.plot_layout_2.addWidget(self.canvas_2)
 
         self.plot_2 = PlotDischarge(
             canvas = self.canvas_2,
             data = self._ics,
-            toolbar = None,
+            toolbar = self.toolbar_2,
         )
         self.plot_2.draw()
 

src/View/InitialConditions/translate.py

@@ -21,11 +21,20 @@ from PyQt5.QtCore import QCoreApplication
 _translate = QCoreApplication.translate
 
 table_headers = {
-    "name": _translate("LateralContribution", "Name"),
+    # "name": _translate("LateralContribution", "Name"),
     "kp": _translate("LateralContribution", "KP (m)"),
     # "speed": _translate("LateralContribution", "Speed (m/s)"),
     "discharge": _translate("LateralContribution", "Discharge (m³/s)"),
     "elevation": _translate("LateralContribution", "Elevation (m)"),
     "height": _translate("LateralContribution", "Height (m)"),
-    "comment": _translate("LateralContribution", "Comment"),
+    # "comment": _translate("LateralContribution", "Comment"),
 }
+
+def retranslate():
+    # table_headers["name"] = _translate("LateralContribution", "Name")
+    table_headers["kp"] = _translate("LateralContribution", "KP (m)")
+    # table_headers["speed"] = _translate("LateralContribution", "Speed (m/s)")
+    table_headers["discharge"] = _translate("LateralContribution", "Discharge (m³/s)")
+    table_headers["elevation"] = _translate("LateralContribution", "Elevation (m)")
+    table_headers["height"] = _translate("LateralContribution", "Height (m)")
+    # table_headers["comment"] = _translate("LateralContribution", "Comment")

src/View/LateralContribution/Edit/Plot.py

@@ -41,6 +41,9 @@ class Plot(APlot):
         self._mode = mode
 
     def custom_ticks(self):
+        if self.data.header[0] != "time":
+            return
+
         t0 = datetime.fromtimestamp(0)
         nb = len(self.data.data)
         mod = int(nb / 5)

src/View/LateralContribution/Edit/Table.py

@@ -24,6 +24,7 @@ from tools import trace, timer
 
 from View.ASubWindow import ASubMainWindow, AWidget
 from View.ListedSubWindow import ListedSubWindow
+from View.Tools.PamhyrTable import PamhyrTableModel
 
 from PyQt5.QtCore import (
     Qt, QVariant, QAbstractTableModel,
@@ -51,125 +52,8 @@ _translate = QCoreApplication.translate
 
 logger = logging.getLogger()
 
-class ExtendedTimeEdit(AWidget):
-    def __init__(self, parent=None):
-        super(ExtendedTimeEdit, self).__init__(
-            ui="extendedTimeEdit",
-            parent=parent
-        )
-        self.parent = parent
-
-        self.spinBox_days = self.find(QSpinBox, "spinBox_days")
-        self.timeEdit = self.find(QTimeEdit, "timeEdit")
-
-    def set_time(self, time):
-        days = 0
-        stime = time
-
-        # if ',' in time, time format is 'DD days, HH:MM:SS',
-        # otherelse is 'HH:MM:SS'
-        if "," in time:
-            s = time.strip().split(" ")
-            days = int(s[0])
-            stime = s[-1]
-
-        qtime = QTime.fromString(
-            stime,
-            "h:mm:ss"
-        )
-        self.spinBox_days.setValue(days)
-        self.timeEdit.setTime(qtime)
-
-    def get_time(self):
-        days = self.spinBox_days.value()
-        time = self.timeEdit.time().toPyTime()
-        secs = (
-            (time.hour * 3600) +
-            (time.minute * 60) +
-            time.second
-        )
-
-        return timedelta(days=days, seconds=secs)
-
-class ExtendedDateTimeEdit(AWidget):
-    def __init__(self, parent=None):
-        super(ExtendedDateTimeEdit, self).__init__(
-            ui="extendedDateTimeEdit",
-            parent=parent
-        )
-        self.parent = parent
-
-        self.dateTimeEdit = self.find(QDateTimeEdit, "dateTimeEdit")
-
-    def set_time(self, time):
-        qdatetime = QDateTime.fromString(
-            time,
-            "yyyy-MM-dd hh:mm:ss"
-        )
-        self.dateTimeEdit.setDateTime(qdatetime)
-
-    def get_time(self):
-        time = self.dateTimeEdit.dateTime().toPyDateTime()
-        return time
-
-class ExTimeDelegate(QItemDelegate):
-    def __init__(self, data=None, mode="time", parent=None):
-        super(ExTimeDelegate, self).__init__(parent)
-
-        self._data = data
-        self._mode = mode
-
-    def createEditor(self, parent, option, index):
-        if self._mode == "time":
-            self.editor = ExtendedTimeEdit(parent=parent)
-        else:
-            self.editor = ExtendedDateTimeEdit(parent=parent)
-        value = index.data(Qt.DisplayRole)
-        self.editor.set_time(value)
-        logger.debug(str(value))
-        return self.editor
-
-    def setModelData(self, editor, model, index):
-        time = editor.get_time()
-        if self._mode == "time":
-            model.setData(index, int(time.total_seconds()))
-        else:
-            logger.debug(str(time.timestamp()))
-            model.setData(index, int(time.timestamp()))
-        editor.close()
-        editor.deleteLater()
-
-    def updateEditorGeometry(self, editor, option, index):
-        r = QRect(option.rect)
-        if self.editor.windowFlags() & Qt.Popup and editor.parent() is not None:
-            r.setTopLeft(self.editor.parent().mapToGlobal(r.topLeft()))
-        editor.setGeometry(r)
-
-    @pyqtSlot()
-    def currentItemChanged(self):
-        self.commitData.emit(self.sender())
-
-
-class TableModel(QAbstractTableModel):
-    def __init__(self, data=None, mode="time", undo=None):
-        super(QAbstractTableModel, self).__init__()
-        self._headers = data.header
-        self._data = data
-        self._mode = mode
-        self._undo = undo
-
-    def flags(self, index):
-        options = Qt.ItemIsEnabled | Qt.ItemIsSelectable
-        options |= Qt.ItemIsEditable
-
-        return options
-
-    def rowCount(self, parent):
-        return len(self._data)
-
-    def columnCount(self, parent):
-        return len(self._headers)
 
+class TableModel(PamhyrTableModel):
     def data(self, index, role):
         if role == Qt.TextAlignmentRole:
             return Qt.AlignHCenter | Qt.AlignVCenter
@@ -187,7 +71,7 @@ class TableModel(QAbstractTableModel):
             if self._data.get_type_column(column) == float:
                 value = f"{v:.4f}"
             elif self._data.header[column] == "time":
-                if self._mode == "time":
+                if self._opt_data == "time":
                     t0 = datetime.fromtimestamp(0)
                     t = datetime.fromtimestamp(v)
                     value = str(t - t0)
@@ -198,12 +82,6 @@ class TableModel(QAbstractTableModel):
 
         return value
 
-    def headerData(self, section, orientation, role):
-        if role == Qt.ItemDataRole.DisplayRole and orientation == Qt.Orientation.Horizontal:
-            return table_headers[self._headers[section]]
-
-        return QVariant()
-
     def setData(self, index, value, role=Qt.EditRole):
         if not index.isValid() or role != Qt.EditRole:
             return False
@@ -313,12 +191,3 @@ class TableModel(QAbstractTableModel):
 
         self.layoutAboutToBeChanged.emit()
         self.layoutChanged.emit()
-
-
-    def undo(self):
-        self._undo.undo()
-        self.layoutChanged.emit()
-
-    def redo(self):
-        self._undo.redo()
-        self.layoutChanged.emit()

src/View/LateralContribution/Edit/Window.py

@@ -20,6 +20,7 @@ from tools import timer, trace
 
 from View.ASubWindow import ASubMainWindow
 from View.ListedSubWindow import ListedSubWindow
+from View.Tools.PamhyrDelegate import PamhyrExTimeDelegate
 
 from PyQt5.QtGui import (
     QKeySequence,
@@ -37,9 +38,11 @@ from PyQt5.QtWidgets import (
 )
 
 from View.Plot.MplCanvas import MplCanvas
+from View.Plot.PamhyrToolbar import PamhyrPlotToolbar
 
 from View.LateralContribution.translate import long_types
-from View.LateralContribution.Edit.Table import TableModel, ExTimeDelegate
+from View.LateralContribution.Edit.translate import table_headers, retranslate
+from View.LateralContribution.Edit.Table import TableModel
 from View.LateralContribution.Edit.Plot import Plot
 
 _translate = QCoreApplication.translate
@@ -84,23 +87,29 @@ class EditLateralContributionWindow(ASubMainWindow, ListedSubWindow):
         self.paste_sc = QShortcut(QKeySequence.Paste, self)
 
     def setup_table(self):
-        table = self.find(QTableView, "tableView")
+        retranslate()
-        self._table = TableModel(
+        headers = {}
+        for h in self._data.header:
+            headers[h] = table_headers[h]
+
+        self._delegate_time = PamhyrExTimeDelegate(
             data = self._data,
-            undo = self._undo_stack,
+            mode = self._study.time_system,
-            mode = self._study.time_system
+            parent = self
         )
 
-        if self._data.header[0] == "time":
+        table = self.find(QTableView, "tableView")
-            self._delegate_time = ExTimeDelegate(
+        self._table = TableModel(
-                data = self._data,
+            table_view = table,
-                mode = self._study.time_system,
+            table_headers = headers,
-                parent = self
+            editable_headers = self._data.header,
-            )
+            delegates = {
-
+                "time": self._delegate_time,
-            table.setItemDelegateForColumn(
+            },
-                0, self._delegate_time
+            data = self._data,
-            )
+            undo = self._undo_stack,
+            opt_data = self._study.time_system
+        )
 
         table.setModel(self._table)
         table.setSelectionBehavior(QAbstractItemView.SelectRows)
@@ -110,12 +119,17 @@ class EditLateralContributionWindow(ASubMainWindow, ListedSubWindow):
     def setup_plot(self):
         self.canvas = MplCanvas(width=5, height=4, dpi=100)
         self.canvas.setObjectName("canvas")
+        self.toolbar = PamhyrPlotToolbar(
+            self.canvas, self
+        )
+        self.verticalLayout.addWidget(self.toolbar)
         self.verticalLayout.addWidget(self.canvas)
 
         self.plot = Plot(
             canvas = self.canvas,
             data = self._data,
             mode = self._study.time_system,
+            toolbar = self.toolbar,
         )
         self.plot.draw()
 

src/View/LateralContribution/Edit/translate.py

@@ -28,3 +28,11 @@ table_headers = {
     "discharge": _translate("LateralContribution", "Discharge (m³/s)"),
     "z": _translate("LateralContribution", "Z (m)")
 }
+
+def retranslate():
+    table_headers["x"] = _translate("LateralContribution", "X")
+    table_headers["y"] = _translate("LateralContribution", "Y")
+    table_headers["time"] = _translate("LateralContribution", "Time")
+    table_headers["date"] = _translate("LateralContribution", "Date")
+    table_headers["discharge"] = _translate("LateralContribution", "Discharge (m³/s)")
+    table_headers["z"] = _translate("LateralContribution", "Z (m)")

src/View/LateralContribution/Table.py

@@ -44,6 +44,8 @@ from View.LateralContribution.UndoCommand import (
 from Model.LateralContribution.LateralContributionTypes import (
     NotDefined, LateralContrib, Rain, Evaporation,
 )
+
+from View.Tools.PamhyrTable import PamhyrTableModel
 from View.LateralContribution.translate import *
 
 logger = logging.getLogger()
@@ -102,26 +104,13 @@ class ComboBoxDelegate(QItemDelegate):
         self.commitData.emit(self.sender())
 
 
-class TableModel(QAbstractTableModel):
+class TableModel(PamhyrTableModel):
-    def __init__(self, data=None, undo=None, tab=""):
+    def _setup_lst(self):
-        super(QAbstractTableModel, self).__init__()
+        self._lst = self._data.lateral_contribution
-        self._headers = list(table_headers.keys())
+        self._tab = self._opt_data
-        self._data = data
-        self._undo = undo
-        self._tab = tab
-        self._lcs = self._data.lateral_contribution
-
-    def flags(self, index):
-        options = Qt.ItemIsEnabled | Qt.ItemIsSelectable
-        options |= Qt.ItemIsEditable
-
-        return options
 
     def rowCount(self, parent):
-        return self._lcs.len(self._tab)
+        return self._lst.len(self._tab)
-
-    def columnCount(self, parent):
-        return len(self._headers)
 
     def data(self, index, role):
         if role != Qt.ItemDataRole.DisplayRole:
@@ -131,26 +120,19 @@ class TableModel(QAbstractTableModel):
         column = index.column()
 
         if self._headers[column] == "name":
-            return self._lcs.get(self._tab, row).name
+            return self._lst.get(self._tab, row).name
         elif self._headers[column] == "type":
-            t = self._lcs.get(self._tab, row).lctype
+            t = self._lst.get(self._tab, row).lctype
             return long_types[t]
         elif self._headers[column] == "edge":
-            n = self._lcs.get(self._tab, row).edge
+            n = self._lst.get(self._tab, row).edge
             if n is None:
                 return _translate("LateralContribution", "Not associate")
             return n.name
         elif self._headers[column] == "begin_kp":
-            return self._lcs.get(self._tab, row).begin_kp
+            return self._lst.get(self._tab, row).begin_kp
         elif self._headers[column] == "end_kp":
-            return self._lcs.get(self._tab, row).end_kp
+            return self._lst.get(self._tab, row).end_kp
-
-
-        return QVariant()
-
-    def headerData(self, section, orientation, role):
-        if role == Qt.ItemDataRole.DisplayRole and orientation == Qt.Orientation.Horizontal:
-            return table_headers[self._headers[section]]
 
         return QVariant()
 
@@ -165,32 +147,32 @@ class TableModel(QAbstractTableModel):
             if self._headers[column] == "name":
                 self._undo.push(
                     SetNameCommand(
-                        self._lcs, self._tab, row, value
+                        self._lst, self._tab, row, value
                     )
                 )
             elif self._headers[column] == "type":
                 key = next(k for k, v in long_types.items() if v == value)
                 self._undo.push(
                     SetTypeCommand(
-                        self._lcs, self._tab, row, LC_types[key]
+                        self._lst, self._tab, row, LC_types[key]
                     )
                 )
             elif self._headers[column] == "edge":
                 self._undo.push(
                     SetEdgeCommand(
-                        self._lcs, self._tab, row, self._data.edge(value)
+                        self._lst, self._tab, row, self._data.edge(value)
                     )
                 )
             elif self._headers[column] == "begin_kp":
                 self._undo.push(
                     SetBeginCommand(
-                        self._lcs, self._tab, row, value
+                        self._lst, self._tab, row, value
                     )
                 )
             elif self._headers[column] == "end_kp":
                 self._undo.push(
                     SetEndCommand(
-                        self._lcs, self._tab, row, value
+                        self._lst, self._tab, row, value
                     )
                 )
         except Exception as e:
@@ -205,7 +187,7 @@ class TableModel(QAbstractTableModel):
 
         self._undo.push(
             AddCommand(
-                self._lcs, self._tab, row
+                self._lst, self._tab, row
             )
         )
 
@@ -217,7 +199,7 @@ class TableModel(QAbstractTableModel):
 
         self._undo.push(
             DelCommand(
-                self._lcs, self._tab, rows
+                self._lst, self._tab, rows
             )
         )
 
@@ -229,7 +211,7 @@ class TableModel(QAbstractTableModel):
 
         self._undo.push(
             SortCommand(
-                self._lcs, self._tab, False
+                self._lst, self._tab, False
             )
         )
 
@@ -246,7 +228,7 @@ class TableModel(QAbstractTableModel):
 
         self._undo_stack.push(
             MoveCommand(
-                self._lcs, self._tab,"up", row
+                self._lst, self._tab, "up", row
             )
         )
 
@@ -254,7 +236,7 @@ class TableModel(QAbstractTableModel):
         self.layoutChanged.emit()
 
     def move_down(self, index, parent=QModelIndex()):
-        if row > len(self._lcs):
+        if row > len(self._lst):
             return
 
         target = row
@@ -263,17 +245,9 @@ class TableModel(QAbstractTableModel):
 
         self._undo_stack.push(
             MoveCommand(
-                self._lcs, self._tab,"down", row
+                self._lst, self._tab, "down", row
             )
         )
 
         self.endMoveRows()
         self.layoutChanged.emit()
-
-    def undo(self):
-        self._undo.undo()
-        self.layoutChanged.emit()
-
-    def redo(self):
-        self._undo.redo()
-        self.layoutChanged.emit()

src/View/LateralContribution/Window.py

@@ -56,7 +56,10 @@ from View.LateralContribution.Table import (
 
 from View.Plot.MplCanvas import MplCanvas
 from View.Geometry.PlotXY import PlotXY
-from View.LateralContribution.translate import *
+from View.LateralContribution.translate import (
+    long_types, table_headers, LC_types,
+    retranslate,
+)
 from View.LateralContribution.Edit.Window import EditLateralContributionWindow
 
 _translate = QCoreApplication.translate
@@ -90,17 +93,10 @@ class LateralContributionWindow(ASubMainWindow, ListedSubWindow):
         self.paste_sc = QShortcut(QKeySequence.Paste, self)
 
     def setup_table(self):
+        retranslate()
         self._table = {}
 
         for t in ["liquid", "solid", "suspenssion"]:
-            table = self.find(QTableView, f"tableView_{t}")
-            self._table[t] = TableModel(
-                data = self._study.river,
-                undo = self._undo_stack,
-                tab = t,
-            )
-            table.setModel(self._table[t])
-
             self._delegate_type = ComboBoxDelegate(
                 data = self._study.river,
                 mode = "type",
@@ -114,13 +110,20 @@ class LateralContributionWindow(ASubMainWindow, ListedSubWindow):
                 parent=self
             )
 
-            table.setItemDelegateForColumn(
+            table = self.find(QTableView, f"tableView_{t}")
-                1, self._delegate_type
+            self._table[t] = TableModel(
+                table_view = table,
+                table_headers = table_headers,
+                editable_headers = True,
+                delegates = {
+                    "type": self._delegate_type,
+                    "edge": self._delegate_edge,
+                },
+                data = self._study.river,
+                undo = self._undo_stack,
+                opt_data = t,
             )
-            table.setItemDelegateForColumn(
+            table.setModel(self._table[t])
-                2, self._delegate_edge
-            )
-
             table.setSelectionBehavior(QAbstractItemView.SelectRows)
             table.horizontalHeader().setSectionResizeMode(QHeaderView.Stretch)
             table.setAlternatingRowColors(True)

src/View/LateralContribution/translate.py

@@ -45,3 +45,15 @@ LC_types = {
     "RA": Rain,
     "EV": Evaporation,
 }
+
+def retranslate():
+    table_headers["name"] = _translate("LateralContribution", "Name")
+    table_headers["type"] = _translate("LateralContribution", "Type")
+    table_headers["edge"] = _translate("LateralContribution", "Reach")
+    table_headers["begin_kp"] = _translate("LateralContribution", "Begin kp (m)")
+    table_headers["end_kp"] = _translate("LateralContribution", "End kp (m)")
+
+    long_types["ND"] = _translate("LateralContribution", "Not defined")
+    long_types["LC"] = _translate("LateralContribution", "Lateral contribution")
+    long_types["RA"] = _translate("LateralContribution", "Rain")
+    long_types["EV"] = _translate("LateralContribution", "Evaporation")

src/View/MainWindow.py

@@ -18,6 +18,7 @@
 
 import os
 import logging
+import subprocess
 from queue import Queue
 
 from PyQt5 import QtGui
@@ -49,9 +50,12 @@ from View.LateralContribution.Window import LateralContributionWindow
 from View.InitialConditions.Window import InitialConditionsWindow
 from View.Stricklers.Window import StricklersWindow
 from View.Frictions.Window import FrictionsWindow
+from View.SedimentLayers.Window import SedimentLayersWindow
+from View.SedimentLayers.Reach.Window import ReachSedimentLayersWindow
 from View.SolverParameters.Window import SolverParametersWindow
 from View.RunSolver.Window import SelectSolverWindow, SolverLogWindow
 from View.CheckList.Window import CheckListWindow
+from View.Results.Window import ResultsWindow
 from View.Debug.Window import ReplWindow
 
 from Model.Study import Study
@@ -86,7 +90,8 @@ define_model_action = [
     "action_menu_edit_network", "action_menu_edit_geometry",
     "action_menu_boundary_conditions", "action_menu_initial_conditions",
     "action_menu_edit_friction", "action_menu_edit_lateral_contribution",
-    "action_menu_run_solver",
+    "action_menu_run_solver", "action_menu_sediment_layers",
+    "action_menu_edit_reach_sediment_layers"
 ]
 
 action = (
@@ -105,6 +110,9 @@ class ApplicationWindow(QMainWindow, ListedSubWindow, WindowToolKit):
         # Model
         self.model = None
 
+        # Results
+        self._last_results = None
+
         # UI
         self.ui = loadUi(
             os.path.join(os.path.dirname(__file__), "ui", "MainWindow.ui"),
@@ -123,10 +131,14 @@ class ApplicationWindow(QMainWindow, ListedSubWindow, WindowToolKit):
             self.dialog_reopen_study()
 
     def set_title(self):
+        title = "(dbg) " if self.conf.debug else ""
+
         if self.model is not None:
-            self.setWindowTitle(f"PAMHYR - {self.model.name}")
+            title += f"Pamhyr2 - {self.model.name}"
+            self.setWindowTitle(title)
         else:
-            self.setWindowTitle("PAMHYR")
+            title += "Pamhyr2"
+            self.setWindowTitle(title)
 
     def enable_actions(self, action:str, enable:bool):
         """Enable of disable an action componant
@@ -166,6 +178,8 @@ class ApplicationWindow(QMainWindow, ListedSubWindow, WindowToolKit):
             "action_menu_edit_friction": self.open_frictions,
             "action_menu_edit_lateral_contribution": self.open_lateral_contrib,
             "action_menu_run_solver": self.run_solver,
+            "action_menu_sediment_layers": self.open_sediment_layers,
+            "action_menu_edit_reach_sediment_layers": self.open_reach_sediment_layers,
             ## Help
             "action_menu_about": self.open_about,
             # ToolBar action
@@ -245,24 +259,31 @@ class ApplicationWindow(QMainWindow, ListedSubWindow, WindowToolKit):
             logger.setLevel(logging.INFO)
             logger.info("Set logging level to INFO")
 
-
     def setup_debug_mode(self, init = False):
         menu = self.findChild(QMenu, "menu_help")
+        self.set_title()
 
         if init:
             self.debug_action = QAction("Debug", self)
             self.debug_action.setStatusTip("Debug")
             self.debug_action.triggered.connect(self.open_debug)
 
+            self.debug_sqlite_action = QAction("Debug SQLite", self)
+            self.debug_sqlite_action.setStatusTip("Open SQLite debuging tool (sqlitebrowser)")
+            self.debug_sqlite_action.triggered.connect(self.open_sqlite)
+
             if self.conf.debug:
                 menu.addAction(self.debug_action)
+                menu.addAction(self.debug_sqlite_action)
                 self.set_debug_lvl(debug = True)
         else:
             if self.conf.debug:
                 menu.addAction(self.debug_action)
+                menu.addAction(self.debug_sqlite_action)
                 self.set_debug_lvl(debug = True)
             else:
                 menu.removeAction(self.debug_action)
+                menu.removeAction(self.debug_sqlite_action)
                 self.set_debug_lvl(debug = False)
 
     #########
@@ -304,6 +325,9 @@ class ApplicationWindow(QMainWindow, ListedSubWindow, WindowToolKit):
         for action in model_action:
             self.enable_actions(action, not no_model)
 
+    def set_results(self, results):
+        self._last_results = results
+
     ############
     # FEATURES #
     ############
@@ -610,6 +634,20 @@ class ApplicationWindow(QMainWindow, ListedSubWindow, WindowToolKit):
         else:
             params.activateWindow()
 
+    def open_sediment_layers(self):
+        sl = SedimentLayersWindow(
+            study = self.model,
+            parent = self
+        )
+        sl.show()
+
+    def open_reach_sediment_layers(self):
+        sl = ReachSedimentLayersWindow(
+            study = self.model,
+            parent = self
+        )
+        sl.show()
+
     def run_solver(self):
         if self.model is None:
             return
@@ -638,6 +676,32 @@ class ApplicationWindow(QMainWindow, ListedSubWindow, WindowToolKit):
         )
         sol.show()
 
+    def open_solver_results(self, solver, results = None):
+        # If no specific results, get last results
+        if results is None:
+            results = self._last_results
+
+        # No results available
+        if results is None:
+            return
+
+        # Windows already opened
+        res = self.sub_win_filter_first(
+            "Results",
+            contain = [solver.name, results.date]
+        )
+
+        if res is None:
+            res = ResultsWindow(
+                study = self.model,
+                solver = solver,
+                results = results,
+                parent = self
+            )
+            res.show()
+        else:
+            res.activateWindow()
+
     #########
     # DEBUG #
     #########
@@ -650,6 +714,17 @@ class ApplicationWindow(QMainWindow, ListedSubWindow, WindowToolKit):
         )
         repl.show()
 
+    def open_sqlite(self):
+        if self.model is None:
+            logger.debug("No study open for sql debuging...")
+            return
+
+        file = self.model.filename
+        _ = subprocess.Popen(
+            f"sqlitebrowser {file}",
+            shell=True
+        )
+
     # TODO: Delete me !
     ###############
     # DUMMY STUFF #

src/View/Network/GraphWidget.py

@@ -584,7 +584,7 @@ class GraphWidget(QGraphicsView):
                 if item:
                     item.update()
         except Exception as e:
-            logger.debug(str(e))
+            logger.warning(str(e))
 
     def selected_new_edge_src_node(self):
         """The current node item selected to add new edge
@@ -626,12 +626,15 @@ class GraphWidget(QGraphicsView):
         Returns:
             Nothing
         """
-        previous_edge = self._current_edge
+        try:
-        self._current_edge = edge
+            previous_edge = self._current_edge
-        self.graph.set_current_reach(edge.edge)
+            self._current_edge = edge
+            self.graph.set_current_reach(edge.edge)
 
-        if previous_edge:
+            if previous_edge:
-            previous_edge.update()
+                previous_edge.update()
+        except Exception as e:
+            logger.warning(str(e))
 
     def reset_selection(self):
         """Reset all the selected items

src/View/Network/Table.py

@@ -25,6 +25,7 @@ from Model.Network.Graph import Graph
 from View.ASubWindow import ASubWindow
 from View.Network.GraphWidget import GraphWidget
 from View.Network.UndoCommand import *
+from View.Tools.PamhyrTable import PamhyrTableModel
 
 from PyQt5.QtCore import (
     Qt, QRect, QVariant, QAbstractTableModel, pyqtSlot, pyqtSignal,
@@ -69,151 +70,105 @@ class ComboBoxDelegate(QItemDelegate):
     def currentItemChanged(self):
         self.commitData.emit(self.sender())
 
-class TrueFalseComboBoxDelegate(QItemDelegate):
+class NodeTableModel(PamhyrTableModel):
-    def __init__(self, parent=None):
+    def _setup_lst(self):
-        super(TrueFalseComboBoxDelegate, self).__init__(parent)
+        self._lst = self._data.nodes()
-
-    def createEditor(self, parent, option, index):
-        self.editor = QComboBox(parent)
-        self.editor.addItems(["true", "false"])
-        self.editor.setCurrentText("true" if index.data(Qt.DisplayRole) else "false")
-        return self.editor
-
-    def setEditorData(self, editor, index):
-        value = str(index.data(Qt.DisplayRole))
-        self.editor.currentTextChanged.connect(self.currentItemChanged)
-
-    def setModelData(self, editor, model, index):
-        value = str(editor.currentText()) == "true"
-        model.setData(index, value)
-        editor.close()
-        editor.deleteLater()
-
-    def updateEditorGeometry(self, editor, option, index):
-        r = QRect(option.rect)
-        if self.editor.windowFlags() & Qt.Popup and editor.parent() is not None:
-            r.setTopLeft(self.editor.parent().mapToGlobal(r.topLeft()))
-        editor.setGeometry(r)
-
-    @pyqtSlot()
-    def currentItemChanged(self):
-        self.commitData.emit(self.sender())
-
-class GraphTableModel(QAbstractTableModel):
-    def __init__(self, headers=[], graph=None, undo = None, rows_type="nodes"):
-        super(QAbstractTableModel, self).__init__()
-        self.headers = headers
-        self.graph = graph
-        self._type = rows_type
-        self._undo = undo
-
-        if self._type == "nodes":
-            self.rows = graph.nodes()
-        elif self._type == "edges":
-            self.rows = graph.edges()
-
-    def flags(self, index):
-        options = Qt.ItemIsEnabled | Qt.ItemIsSelectable
-
-        if self.headers[index.column()] != "type":
-            options |= Qt.ItemIsEditable
-
-        return options
-
-    def rowCount(self, parent):
-        return len(self.rows)
-
-    def columnCount(self, parent):
-        return len(self.headers)
 
     def data(self, index, role):
         if role != Qt.ItemDataRole.DisplayRole:
             return QVariant()
 
-        if self.headers[index.column()] == "type":
+        if self._headers[index.column()] == "type":
-            node = self.rows[index.row()]
+            node = self._lst[index.row()]
             ret = "internal"
 
-            if not self.graph.is_enable_node(node):
+            if not self._data.is_enable_node(node):
                 ret = "disable"
-            elif self.graph.is_upstream_node(node):
+            elif self._data.is_upstream_node(node):
                 ret = "upstream"
-            elif self.graph.is_downstream_node(node):
+            elif self._data.is_downstream_node(node):
                 ret = "downstream"
 
             return ret
 
-        return self.rows[index.row()][self.headers[index.column()]]
+        return self._lst[index.row()][self._headers[index.column()]]
-
-    def headerData(self, section, orientation, role):
-        if (role == Qt.ItemDataRole.DisplayRole and
-            orientation == Qt.Orientation.Horizontal):
-            return self.headers[section].capitalize()
-
-        return QVariant()
 
     @pyqtSlot()
     def setData(self, index, value, role=Qt.EditRole):
-        if index.isValid():
+        if not index.isValid():
-            if role == Qt.EditRole:
+            return False
-                try:
+
-                    if (self.headers[index.column()] == "node1" or
+        if role == Qt.EditRole:
-                        self.headers[index.column()] == "node2"):
+            try:
-                        node = self.graph.node(value)
+                self._undo.push(
-                        self._undo.push(
+                    SetCommand(
-                            SetNodeCommand(
+                        self._lst[index.row()],
-                                self.graph,
+                        self._headers[index.column()],
-                                self.rows[index.row()],
+                        value
-                                self.headers[index.column()],
+                    )
-                                node
+                )
-                            )
+            except Exception as e:
-                        )
+                logger.info(e)
-                        # elif self.headers[index.column()] == "enable":
+                logger.debug(traceback.format_exc())
-                        #     self._undo.push(
-                        #         EnableEdgeCommand(
-                        #             self.rows[index.row()], value
-                        #         )
-                        #     )
-                    else:
-                        self._undo.push(
-                            SetCommand(
-                                self.rows[index.row()],
-                                self.headers[index.column()],
-                                value
-                            )
-                        )
-                except Exception as e:
-                    logger.info(e)
-                    logger.debug(traceback.format_exc())
 
                 self.dataChanged.emit(index, index, [Qt.DisplayRole])
                 self.layoutChanged.emit()
                 return True
 
             self.dataChanged.emit(index, index)
-        else:
-            return False
 
     def update(self):
-        if self._type == "nodes":
+        self._lst = self._data.nodes()
-            self.rows = self.graph.nodes()
-        elif self._type == "edges":
-            self.rows = self.graph.edges()
-
         self.layoutChanged.emit()
 
-    def reverse_edge(self, index):
+class EdgeTableModel(PamhyrTableModel):
-        if self._type == "edges":
+    def _setup_lst(self):
-            tmp = self.rows[index.row()].node1
+        self._lst = self._data.edges()
-            self.rows[index.row()].node1 = self.rows[index.row()].node2
+
-            self.rows[index.row()].node2 = tmp
+    def data(self, index, role):
+        if role != Qt.ItemDataRole.DisplayRole:
+            return QVariant()
+
+        return self._lst[index.row()][self._headers[index.column()]]
+
+    @pyqtSlot()
+    def setData(self, index, value, role=Qt.EditRole):
+        if not index.isValid():
+            return False
+
+        if role != Qt.EditRole:
+            return QVariant()
+
+        try:
+            if (self._headers[index.column()] == "node1" or
+                self._headers[index.column()] == "node2"):
+                node = self.graph.node(value)
+                self._undo.push(
+                    SetNodeCommand(
+                        self._data,
+                        self._lst[index.row()],
+                        self._headers[index.column()],
+                        node
+                    )
+                )
+            else:
+                self._undo.push(
+                    SetCommand(
+                        self._lst[index.row()],
+                        self._headers[index.column()],
+                        value
+                    )
+                )
+        except Exception as e:
+            logger.info(e)
+            logger.debug(traceback.format_exc())
+
             self.dataChanged.emit(index, index, [Qt.DisplayRole])
             self.layoutChanged.emit()
+            return True
 
-    def undo(self):
+        self.dataChanged.emit(index, index)
-        self._undo.undo()
-        self.layoutChanged.emit()
 
-    def redo(self):
+    def update(self):
-        self._undo.redo()
+        self._lst = self._data.edges()
         self.layoutChanged.emit()

src/View/Network/Window.py

@@ -16,6 +16,8 @@
 
 # -*- coding: utf-8 -*-
 
+from PyQt5.QtCore import QCoreApplication
+
 from PyQt5.QtGui import (
     QKeySequence,
 )
@@ -37,9 +39,16 @@ from Model.River import RiverNode, RiverReach, River
 from View.ASubWindow import ASubMainWindow
 from View.Network.GraphWidget import GraphWidget
 from View.Network.UndoCommand import *
-from View.Network.Table import (
+from View.Network.translate import (
-    GraphTableModel, ComboBoxDelegate, TrueFalseComboBoxDelegate,
+    table_headers_node, table_headers_edge,
+    retranslate,
 )
+from View.Network.Table import (
+    ComboBoxDelegate, NodeTableModel, EdgeTableModel,
+)
+
+_translate = QCoreApplication.translate
+
 
 class NetworkWindow(ASubMainWindow):
     def __init__(self, model=None, title="River network", parent=None):
@@ -69,43 +78,41 @@ class NetworkWindow(ASubMainWindow):
         self.redo_sc = QShortcut(QKeySequence.Redo, self)
 
     def setup_table(self):
-        # Nodes table
+        retranslate()
 
-        self._nodes_model = GraphTableModel(
+        # Nodes table
-            headers = ["name", "type"],
+        table = self.find(QTableView, "tableView_nodes")
-            graph = self._graph,
+        self._nodes_model = NodeTableModel(
-            rows_type = "nodes",
+            table_view = table,
+            table_headers = table_headers_node,
+            editable_headers = ["name"],
+            data = self._graph,
             undo = self._undo_stack,
         )
-        table = self.find(QTableView, "tableView_nodes")
         table.setModel(self._nodes_model)
-        #table.resizeColumnsToContents()
 
         table.setSelectionBehavior(QAbstractItemView.SelectRows)
         table.horizontalHeader().setSectionResizeMode(QHeaderView.Stretch)
 
         # Edges table
-
+        table = self.find(QTableView, "tableView_reachs")
-        self._reachs_model = GraphTableModel(
-            headers = ["name", # "enable",
-                       "node1", "node2"],
-            graph = self._graph,
-            rows_type = "edges",
-            undo = self._undo_stack,
-        )
         self.delegate_combobox = ComboBoxDelegate(
             graph = self._graph,
             parent = self,
         )
-        self.delegate_true_false_combobox = TrueFalseComboBoxDelegate(
-            parent = self,
-        )
 
-        table = self.find(QTableView, "tableView_reachs")
+        self._reachs_model = EdgeTableModel(
+            table_view = table,
+            table_headers = table_headers_edge,
+            editable_headers = ["name", "node1", "node2"],
+            delegates = {
+                "node1": self.delegate_combobox,
+                "node2": self.delegate_combobox,
+            },
+            data = self._graph,
+            undo = self._undo_stack,
+        )
         table.setModel(self._reachs_model)
-        # table.setItemDelegateForColumn(1, self.delegate_true_false_combobox)
-        table.setItemDelegateForColumn(1, self.delegate_combobox)
-        table.setItemDelegateForColumn(2, self.delegate_combobox)
         table.setSelectionBehavior(QAbstractItemView.SelectRows)
         table.horizontalHeader().setSectionResizeMode(QHeaderView.Stretch)
         #table.resizeColumnsToContents()