SIC^2 : Logiciel de Simulation Intégrée des Canaux et de leur Contrôle

Import de géométrie à partir d'un MNT ou d'un relevé de géomètre

David Dorchies — 2019-10-22T11:27:22Z

Nous proposons un outil en ligne permettant d'extraire des profils en travers et de créer un fichier d'importation au format texte de SIC à partir d'un MNT ou d'un relevé de géomètre :

http://sic.g-eau.fr/sic_tools/import_xyz.html

Deux méthodes d'extractions des profils en travers sont proposées : l'une est plus adaptée à des relevés homogènes de points tels que les MNT, l'autre est plus adaptée à des relevés de géomètres où des profils en travers ont été relevés.

Format du fichier d'importation des points à coordonnées (X, Y, Z)

Le fichier d'importation est un fichier au format texte contenant 3 colonnes représentant l'abscisse, l'ordonnée et l'altitude dans un repère orthonormé et métrique des points de géométrie relevés sur le terrain (Par exemple des points référencés dans le système de coordonnées Lambert pour la France). Les séparateurs acceptés pour les colonnes sont la tabulation, le point-virgule et la virgule. Le caractère utilisé pour les décimales est le point.

Exemple de fichier de points en coordonnées Lambert avec séparateur tabulation :

Il est conseillé de mettre l'outil directement en plein écran avant d'effectuer l'importation du fichier. L'importation consiste à choisir le fichier sur son disque dur à l'aide du bouton "Parcourir..." et de cliquer sur le bouton "Process file".

Tracé de la ligne centrale d'écoulement

Une fois le fichier importé, il faut indiquer où se trouve le cours d'eau à numériser en dessinant à la souris la ligne centrale d'écoulement. Cela se fait en posant, par un clic de souris, des noeuds qui vont former un ensemble de segments.

Le point kilométrique (pK) est indiqué à chaque noeud tracé. Il est possible de recommencer son tracé en cliquant sur le bouton "Reset route".

Extraction de profils en travers à partir d'un MNT

Cet outil permet d'extraire des profils en travers à intervalles réguliers en transformant le relevé de point en un modèle numérique de terrain par la méthode de triangulation de Delaunay.

L'accès à cet outil se fait par le bouton "Export for SIC²". Les paramètres de l'outil sont les suivants :

Space step (m) : le pas d'espace désiré entre chaque section en travers (en mètres) ;
Section width (m) : la largeur maximale de la section (en mètres). Au delà de cette distance par rapport à la ligne centrale, les intersections de la section avec les arêtes des triangles de Delaunay seront ignorées ;
Max length of triangle edges (m) : Longueur maximale des arêtes des triangles de Delaunay.

Un appui sur le bouton "Generate section profiles" lance l'outil avec les paramètres choisis. Le résultat s'affiche graphiquement et produit un texte à copier-coller dans l'outil d'import de sections et de profils en travers au format texte d'EdiSIC. L'opération peut prendre du temps et il faut patienter même si le navigateur indique qu'un script bloque le fonctionnement de l'application.

Extraction de profils en travers à partir d'un relevé de profils en travers de géomètre

Cet outil permet de créer automatiquement les sections pour SIC à l'emplacement où ont été effectués des relevés de profils en travers. L'algorithme convertit les coordonnées des points dans un repère correspondant au tracé de la ligne centrale d'écoulement (axe pK). Il regroupe ensuite les points proches par section.

L'accès à l'outil se fait par le bouton "Export pK coords for SIC²". Les paramètres de l'outil sont les suivants :

Max Length between points (m) : distance maximale selon l'axe pK entre deux points pour un regroupement dans une même section (en mètres) ;
Max width of cross profiles (m) : largeur maximale des sections en travers (en mètres). Au delà de cette largeur, les points sont ignorés.
Min number of points in section : nombre minimum de points pour constituer une section en travers. Si un regroupement de points n'atteint pas cette limite minimale, la section est ignorée.

Le résultat s'affiche graphiquement et produit un texte à copier-coller dans l'outil d'import de sections et de profils en travers au format texte d'EdiSIC. L'opération peut prendre du temps et il faut patienter même si le navigateur indique qu'un script bloque le fonctionnement de l'application.

Import de fichier de géométrie au format Mike 11

David Dorchies — 2019-10-22T09:27:34Z

Le logiciel Mike 11 (https://www.mikepoweredbydhi.com/products/mike-11) est un code de modélisation hydraulique monodimensionnelle à surface libre, basé sur les équations de Saint-Venant développé par DHI (Dansk Hydraulisk Institut).

Pour importer un fichier au format Mike 11, il faut utiliser l'outil en ligne de conversion de fichier au format Mike 11 vers le format d'import de sections et de profils en travers au format texte de SIC : http://sic.g-eau.fr/sic_tools/import_mike11.html

Une autre option est, sous Mike 11, d'exporter les sections dans un fichier txt en mode raw data, puis de lancer une routine sous Matlab (script mike2sic.m) qui va convertir ce fichier dans un autre fichier txt au format d'importation de SIC, et enfin d'importer les données de ce fichier dans un bief. En option on peut convertir le format des sections de abscisse-cote en largeur-cote (routine ac2lc), ainsi que de réduire le nombre de point de description (routine F_Approx_Sections_v6). Ces routines sont disponibles sur simple demande.

Les sections en travers de MIKE11 peuvent être exportées en suivant les étapes ci-dessous :

Ouvrir MIKE Zéro
Ouvrez le .xns11 via Ouvrir (ou glisser-déposer)
Sélectionnez Fichier >> Exporter toutes les sections >> Exporter les données brutes…
Enregistrez les sections transversales sous forme de fichier .txt

Format du fichier SIR (SIC version 4.33)

Pierre-Olivier Malaterre — 2019-09-20T10:57:31Z

Ce fichier ASCII est créé par EDISIR ou au moyen d'un éditeur de fichier ASCII quelconque. Il y a certaines rubriques qui sont gérées par SIRENE mais pas entièrement par EDISIR. EDISIR peut les lire et écrire mais pas les saisir aux interfaces. Ces rubriques doivent donc être saisies par un éditeur de fichier ASCII. Il contient les données hydrauliques pour le programme de calcul en régime transitoire SIRENE. C'est directement ce fichier qui sera destiné à SIRENE, il n'y a pas d'équivalent binaire intermédiaire contrairement au calcul en régime permanent (DON équivalent binaire du fichier ASCII FLU).

Ce fichier est destiné aux anciennes versions de SIC 4.33 ou antérieures. Ces fichiers peuvent également être importés à partir de SIC versions 5 et ultérieures (en mode interface, mais pas en mode batch). Par ailleurs, si vous avez des anciens projets SIC version 4 vous avez peut être envie de comprendre ce format pour saisir les données correspondantes dans un projet SIC version 5. C'est pour cette raison que ce format, devenu obsolète est néanmoins décrit. Des versions préliminaires de la version 5 ont utilisé à une époque des fichiers à ce format, pour accompagner la transition entre la version 4 et la version 5, mais ce n'est plus le cas actuellement.

Aucun ordre spécial n'est suivi pour écrire ce fichier de données. Les différents types de lignes possibles sont décrits ci-dessous, chaque ligne commençant avec un emblème particulier. Par contre, pour importer des données à un objet il faut que cet objet soit déjà défini (ex : une vanne doit être définie avant d'y rentrer une loi d'ouverture fonction du temps).

Ce fichier contient toutes les informations de nature hydraulique nécessaires pour le calcul en régime transitoire du programme SIRENE. Les données géométriques sont dans des fichiers séparés, les mêmes que pour le régime permanent. Lorsque l'unité n'est pas précisée ce sont les unités SI qui sont utilisées.

Pour les fichiers de la version 4 toutes les lois fonctions du temps sont comprises en relatif, c'est à dire que s'il y a un écart entre la valeur initiale dans la ligne d'eau initiale et la valeur au temps initial de cette loi, alors les valeurs de la loi seront toutes décalées de cet écart pour avoir une continuité entre le régime initial et le transitoire que l'on veut simuler.

Pour les fichiers de la version 5 cette manière de gérer les lois dans le calcul a changé et ce décalage des lois n'est plus fait, les lois étant alors considérées comme en absolu. Cela est vrai pour toutes les lois, sauf les lois des débits latéraux en fonction du temps qui sont toujours décalées de l'écart initial éventuel, s'il existe. Pour la version 5 ce sont les éditeurs de données qui gèrent les recalages automatiques éventuels des lois par rapport aux valeurs initiales.

Titre_etude

TITM, VERFIC, IFLAGX, IFLAGQ
FORMAT (' !', A40, 1X, F5.3, 1X, I1, 1X, I1)

TITM titre de l'étude,
VERFIC version du fichier (3.0, 3.1, 4.0, 4.32, 4.33, etc.),
IFLAGX mode grandes abscisses (0 si non, 1 si oui ie les abscisses en long seront multipliées par 1000),
IFLAGQ mode grands débits (0 si non, 1 si oui ie les débit seront multipliées par 1000)

Biefs et Noeuds

IB, N1, N2, TITB
FORMAT ('N', 3I3, A40)

IB numéro de bief,
N1 numéro noeud amont du bief IB,
N2 numéro noeud aval du bief IB,
TITB titre du bief

A partir de la version 5 le format est modifié comme suit :
FORMAT ('N', 3I3, A40, 2I5, 2F10.5)
avec les mêmes données que si dessus et 4 nouvelles données :
IFBPC(I,1), IFBPC(I,2), PCAM(I), PCAV(I)
IFBPC(I,1)=1 si égalité des charges a l'amont, =0 si égalité des cotes,
IFBPC(I,2)=1 si égalité des charges a l'aval, =0 si égalité des cotes,
PCAM(I) perte de charge a l'entrée du bief,
PCAV(I) perte de charge a la sortie du bief

Paramètres du calcul 1

IFLTI, IREP, TREP, IDHA, NEWT, NQUAS, IDISCR
FORMAT ('B', 1X, I1, 1X, I1, 1X, F10.0, 1X, I1,1X,I1,1X,I3,1X,I1)

IFLTI flag_terme_inertie (1 : avec terme d'inertie classique, 2 : torrentiel local, 0 torrentiel global),
IREP flag_continuation_fichiers_res/rcs (0 : non),
TREP temps_de_reprise (démarrage a chaud) en mn,
IDHA indice de perte de charge automatique aux élargissements (0 : non, 1 : oui), ce flag est transmis à partir de celui indiqué dans le fichier .flu,
NEWT = 0 si point fixe, 1 si Newton, 2 si Point fixe + Newton,
NQUAS = 0 si Newton pur, N si Quasi-Newton : fréquence calcul de la dérivée,
IDISCR indice de discrétisation (0 : classique, 1 : homogène au permanent)

A partir de la version 5 le format est modifié comme suit :

IFLTI, IREP, TREP, IDHA, INISV, NLSV, IDISCR

INISV=0 initialisation classique valeur du pas précédent,
INISV=1 initialisation variation du pas précédent,
NLSV=0 calcul en linéaire,
NLSV>0 calcul en non-linéaire donne la mise à jour de la jacobienne tous les NLSV itérations

Paramètres du calcul 2

IPS, IFR, IRCS, DT, TDEB, TFIN, XX, WW, XLAMB, ITER1, ERZ0, ITER2, EPSI
FORMAT ('C', I3, 1X, 2(I1, 1X), F10.2, 1X, 2(F10.0, 1X), 2(F4.2, 1X), F8.6, 1X, 2(I4, 1X, F8.6))

IPS pas_de_stockage_resultat (tous les 1 ou n pas de calcul, pour limiter la taille des fichiers .RES),
IFR flag_de_regulation (0 : sans utilisation de fichier .reg de régulation, 1 : avec),
IRCS flag de fichier complet .RCS (0 : sans création de ce fichier, 1 : avec création de ce fichier qui est nécessaire pour une reprise à chaud ultérieure, ie un calcul en transitoire à partir d'un autre calcul en transitoire),
DT pas_de_temps (en minutes),
TDEB temps_debut du calcul (en minutes),
TFIN temps_fin du calcul (en minutes),
XX limite_superieure du debit_effectif_aux_prises,
WW limite_inferieure du debit_effectif_aux_prises,
XLAMB coefficient lambda de contractance (0 ITER1 nombre d'itérations (sections singulières, casiers),
ERZ0 valeur numérique de la précision du système non linéaire,
ITER2 nombre d'itérations (prises),
EPSI erreur max. sur les itérations locales pour les prises

Ligne d'eau initiale (temps_debut)

IB, IS, IQL, Q, Z, QL, SK, SKM
FORMAT ('I', 3I3, 2F10.5, F10.7, 2F10.5) en version 4.32
FORMAT ('I', I3, I4, I3, 2F10.5, F10.7, 2F10.5) en version 4.33

IB bief,
IS numero_section_relatif_bief,
IQL indicateur_debit_lateral (0 : sans débit latéral, 1 : avec débit latéral), ce flag provient de ce qui a été rentré dans le fichier .flu, c'est lui qui permettra de repérer les sections sur lesquelles répartir éventuellement les débits latéraux variables au cours du temps (cf rubrique ci-dessous),
Q débit (m3/s),
Z cote (m),
QL debit_lateral (m3/s/m, ie avec un facteur 10-6 par rapport à ce qui a été rentré dans le fichier .flu),
SK strickler_lit_moyen,
SKM Strickler_lit_mineur

Déversoirs ou vannes

IB, IS, KS, (CSG(NSG, I), I=1, 5), IFOL, ISF
FORMAT ('D' ou 'V', 3I3, 5F13.6, 2I2) en version 4.32
FORMAT ('D' ou 'V', I3, I4, I3, 5F13.6, 2I2) en version 4.33

IB numéro du bief,
IS numero_section_relatif,
KS numero_singularite_relatif,
(CSG(NSG, I), I=1, 5) caractéristiques de la singularité NSG (1 larg._seuil, 2 cote_seuil, 3 ouverture, 4 coef._debit, 5 haut._surverse),
IFOL = équation (0 : Cemagref 88 seuil-vanne pour V ou seuil-orifice pour D, 1 : Cunge 80, 2 : Classique noyée pour V ou dénoyé pour D, 3 : Cemagref 88 seuil-orifice pour V ou seuil-vanne pour D, 4 : MatLab, 5 : Fortran),
ISF = mode de calcul (0 : non-linéaire, 1 : linéaire)

Vannes Gec-Alsthom

IB, IS, KS, (CSG(NSG, I), I=1, 8), IDNS, IFOL, IFOC, IFLO, ISF
FORMAT ('G', 3I3, 7F13.6, F7.1, I3, 2I2, I3, I2) en version 4.32
FORMAT ('G', I3, I4, I3, 7F13.6, F7.1, I3, 2I2, I3, I2) en version 4.33

IB numéro du bief,
IS numero_section_relatif,
KS numero_singularite_relatif,
(CSG(NSG, I), I=1, 8) caractéristiques de la singularité NSG (1 rayon, 2 cote_seuil, 3 ouverture max, 4 Zamont_max, 5 décrément, 6 cote_axe cadrée par rapport a ZMIN, 7 ouverture réelle, 8 ratio S1/S2 surface libre / surface fente pour bac vannes AVIS, AVIO et Mixtes pour Versions >= 4.17),
IDNS=décallage de la section pour Z2 (utile pour vannes AVIS, AVIO et Mixtes) en nombre de sections de calcul (Versions >= 4.17),
IFOL = équation (0 : Goussard, 1 : Goussard pour position et Cemagref pour débit, 2 : W donné par l'utilisateur et Cemagref pour débit),
IFOC = type (0 : Amil, 1 : Avio Haute Chute, 2 : Avio BC, 3 : Avis HC, 4 : Avis BC),
IFLO = modèle,
ISF = mode de calcul (0 : non-linéaire, 1 : linéaire).

Noms des noeuds, somme débits et cotes aux noeuds

TC(NOMAX), ITN(NOMAX), QAE(NOMAX), ZNO(NOMAX), I
FORMAT ('S', A6, I3, 2F10.5, 7X, I3)

TC nom_noeud (alphanumérique),
ITN type_noeud = 1 si amont_modele, 2 si noeud_intermédiaire, -2 si aval_modele,
QAE somme_debits = (Somme_Q_partant - Somme_Q_arrivant) au noeud non aval_modele ou débit de fuite au noeud aval_modele,
ZNO somme_cotes,
I numéro du noeud

Prises définies par leurs caractéristiques géométriques

I, KP(I), KVP(I), QPO(I), QPR(I), X1P(I), X2P(I), X3P(I), X4P(I), X5P(I), X6P(I), X7P(I), X8P(I), X9P(I), IFPL(I)
FORMAT ('P', 3I3, 2F7.3, F6.2, 3F10.5, 3F6.2, F7.3, F6.2, I2)

I noeud (non amont_modele et non_aval_modele)
(KP) cas_prise =
1 Q(t) vanne calcul ouverture W
2 Q(t) déversoir calcul Z_seuil
3 Q(t) déversoir calcul largeur_seuil
4 vanne calcul Q
5 déversoir Z(t) calcul Q
6 apm calcul Q
7 tube calcul Q
8 Q(z) calcul Q
9 déversoir L (t) calcul Q
(KVP) cas_aval_prise =
1 Z_aval_prise constante
2 deversoir_aval_prise
3 Q(z)_aval_prise

les lois Q(z) aux prises sont introduites sous forme analytique :
Q=Q_base*( (Z-Z_seuil)/(Z_base-Z_seuil) )**exposant

(QPO) debit_prise_objectif
(QPR) debit_prise_reel
(X1P) coefficient_debit (KP=1, 2, 3, 4, 5, 6, 7, 9) ou exposant (KP=8)
(X2P) Z_seuil
(X3P) larg._seuil (KP=1, 2, 3, 4, 5, 6, 9) ou diamètre (KP=1, 4, 7) ou Q_base (KP=8)
(X4P) ouverture (KP=1, 4, 6) ou Z_seuil_mini (KP=2) ou Z_base (KP=8)
(X5P) ouverture_maxi (KP=1, 4 et X5P>0)
ou indique vanne circulaire (KP=1, 4 et X5P=0) ou Z_seuil_maxi (KP=2)
ou largeur_seuil_maxi (KP=3) ou switch (KP=6, 7, 8)
(X6P) coefficient_debit_aval_prise (KVP=2) ou exposant_aval_prise (KVP=3)
(X7P) Z_aval_prise (KVP=1) ou Z_seuil_aval_prise (KVP=2, 3)
(X8P) largeur_seuil_aval_prise (KVP=2) ou Q_base_aval_prise (KVP=3)
(X9P) Z_base_aval_prise (KVP=3)
IFPL = équation (0 : Cemagref ou 1 : Cunge)

Lois Q(z) au noeud aval de modèle

I, Z, Q
FORMAT ('A', I3, 6X, 2F10.3)

I noeud_aval,
Z cote_aval (m),
Q debit_aval (m3/s)

Paramètres display (affichage)

CODEM, BRANM, NAMEM, NBM, NSM, XLONG, ZREF
FORMAT ('Y', A1, A8, A6, 2I3, F8.0, F8.2) en version 4.32
FORMAT ('Y', A1, A8, A6, I3, I4, F8.0, F8.2) en version 4.33

CODEM code,
BRANM nom_branche,
NOEUDM nom_noeud,
NBM bief,
NSM num._relatif_section,
XLONG abscisse,
ZREF cote_refer

Variations ouverture des vannes

IB, IS, KS, IO, T, X, K, L, ITYPOUV
FORMAT ('U', 4I3, 2F10.3,3I3)

IB bief,
IS num._section_relatif,
KS num._singular._relatif,
IO type d'ouvrage : 1 (Déversoir), 2 (Vanne rectangulaire) et 3 (Vanne Gec-Alsthom) à partir de la version 4.07 incluse (avant la version 4.07 on ne gérait que les vannes en ouvrage mobile)
T temps (mn),
X ouverture,
K numéro de la loi,
L numéro du point de la loi,
ITYPOUV (0 : échelon, 1 : rampe)

Hydrogramme ou loi a la prise au noeud non aval de modèle

CH1, J, T, X, I, L, IOD, NOE, ITYPOFF
FORMAT ('H', A1, 2X, I3, 3X, 2F10.3, I3, I5, 3I3)

CH1 caractere_loi_prise : 'Q' débit_prise, 'W' ouverture_prise, 'S' switch_prise, 'L' largeur_seuil_prise, 'Z' Z_seuil_prise,
J noeud_non_aval,
T temps (mn),
X valeur_loi,
I numéro de la loi,
L numéro du point de la loi,
IOD(I) type de prise,
NOE(I,3) mode de calcul à la prise,
ITYPOFF (0 : échelon, 1 : rampe)

Débits objectifs

J, T, Q, I, L
FORMAT ('O', I3, 6X, 2F10.3, 2I3)

J noeud,
T temps (mn),
Q debit_objectif,
I numéro de la loi,
L numéro du point de la loi

Casier (loi S(z) au noeud)

I, Z, S, X
FORMAT ('K', I3, 6X, 3F10.3)

I noeud,
Z cote,
S surface (ha),
X vitesse_evaporation (cm/jour)

Hydrogramme latéral (géré par Sirene mais pas par Edisir)

Cette rubrique permet de faire évoluer le terme d'infiltration - évaporation (si le terme est négatif) ou d'apport (si le terme est positif) en fonction du temps. En régime permanent, ce terme était une constante éventuellement différente pour chaque section de calcul. La valeur entrée pour le régime permanent était exprimée en l/s/km. En régime transitoire, ce terme peut être défini comme une fonction du temps. La valeur entrée est un débit total (en m3/s) pour le bief considéré. Ce débit sera réparti sur les zones à l'intérieur du bief pour lesquelles il existe un débit latéral, au prorata de leurs longueurs. A priori, par défaut, ces zones correspondent à celles ainsi définies en régime permanent. Mais il est possible de les modifier en mettant des 0 ou des 1 dans la colonne correspondante des lignes "I" du fichier .SIR. Le débit n'est pas modulé proportionnellement aux débits latéraux initiaux. Il sera réparti de manière homogène sur toutes les zones ayant un débit latéral. On fait donc l'hypothèse que dans le bief il y a des zones avec et des zones sans débit latéral, et que, dans les zones où il y a un débit latéral, le taux d'apports latéraux ou de pertes latérales est identique sur toutes ces zones (en débit/unité de longueur). Si ce n'est pas le cas il faut découper ce bief en plusieurs biefs.

Si le débit latéral total initial indiqué dans cet hydrogramme latéral pour ce bief n'est pas égal au débit latéral total initial provenant du calcul en régime permanent, l'hydrogramme sera translaté pour assurer la continuité au démarrage entre le régime permanent et le régime transitoire. Cela est vrai pour les versions 4 de SIC mais également pour les version 5.

Par contre, il n'est pas possible, pour l'instant, de rentrer un terme d'infiltration - évaporation comme une fonction du tirant d'eau. Cela peut éventuellement être modélisé de manière approchée par une succession de prises en mode calcul de débit. Cette limitation nous semble peut contraignante dans la mesure où une estimation du terme d'infiltration - évaporation est très délicate avec un facteur d'incertitude très important et il est illusoire de considérer d'une loi q=f(z) donnera un modèle plus réaliste.

Cette rubrique "hydrogramme latéral fonction du temps" sera prochainement gérée entièrement par Edisir (saisie aux interfaces graphiques) dans une version ultérieure. Edisir permet déjà de le lire et de l'écrire sur fichier, mais pas de le garnir par les interfaces (depuis la version 3.8). Sirene le gère parfaitement au niveau des calculs.

IB, K, T, Q
FORMAT ('L', 2I3, 3X, 2F10.3)

IB numéro de bief,
K flag échelon (0) ou rampe (1)
T temps (mn),
Q debit_lateral (m3/s)

Limnigramme au noeud aval de modèle (géré par Sirene mais pas par Edisir)

Cette rubrique permet de remplacer la loi Q(z) rentrée pour le calcul en régime permanent par une loi z(t) permettant par exemple de représenter l'influence aval d'une marée.

Prochainement géré par Edisir dans une version ultérieure. Edisir permet déjà de le lire et de l'écrire sur fichier, mais pas de le garnir par les interfaces.

Les valeurs de z sont traitées en relatif par rapport à la première valeur. Par exemple, si vous voulez augmenter la cote z de 1m à l'instant 2h=120mn, il suffit de rentrer :

Z 2 0 0.000 0.000
Z 2 0 120.000 1.000

Attention : si vous rentrez une loi z(t) il faut penser à enlever la loi Q(z) car sinon elle sera conservée car elle est prioritaire par rapport à la z(t). Une loi Q(z) est également prioritaire par rapport à un casier qui fermerait le canal. Si ce casier est rentré à partir des interfaces, alors la Q(z) sera correctement supprimée du fichier .SIR. Par contre, si ce casier est rentré à la main directement dans le fichier .SIR il faut alors penser à supprimer la Q(z) pour qu'elle ne soit plus utilisée.

Le calcul avec une z(t) est beaucoup plus rapide qu'avec une Q(z) car il n'y a pas d'itérations non-linéaires pour le calcul au noeud aval. Par contre, cette condition aval z(t) est parfois considérée comme "dure" hydrauliquement et peut générer des oscillations sur les débits dans cette zone, par exemple dans le cas où un hydrogramme amont arrive sur un noeud aval avec une z(t) constante. Cette condition aval z(t) n'est en effet pas très physique car la "pincette à surface" n'existe pas. La marée est physique mais n'est pas une simple imposition d'une cote, mais un procédé qui relie également les débits et les cotes via des processus physiques complexes. Dans le cas d'une z(t) aval de type sinusoïdale mimant une marée il n'y a pas, a priori, d'oscillation parasite.

Même si cela est possible au niveau du fichier .SIR, il n'est pas conseillé d'avoir à la fois un casier et une z(t) pour un même noeud aval. Dans ce cas on imposerait le débit ainsi que la cote ce qui pose des problèmes numériques. Dans une version ultérieure ce mode sera explicitement interdit et dans ce cas c'est la z(t) qui sera effective et le casier inopérant.

I, K, T, Z

FORMAT ('Z', 2I3, 3X, 2F10.3)

I numéro du noeud_aval,
K flag échelon (0) ou rampe (1) [option depuis la version 3.7]
T temps (mn),
Z cote_aval

Format du fichier TAL (SIC version 4.33)

Pierre-Olivier Malaterre — 2018-09-27T13:35:02Z

Le fichier .TAL est un fichier texte décrivant l'ensemble des informations de nature bathymétrique nécessaires à EDITAL et à TALWEG pour les anciennes versions de SIC 4.33 ou antérieures. Ces fichiers peuvent être importés et/ou exportés à partir de SIC versions 5 et ultérieures (Editeur de projet EDISIC). C'est pour cette raison que ce format, devenu obsolète est néanmoins décrit dans cette documentation.

Sa structure est hiérarchisée en ce qui concerne l'ordre et la nature des informations, et comprend les parties suivantes, obligatoirement dans cet ordre :

Description du réseau
Description des noeuds du réseau
Description des biefs du réseau & Description des sections du bief
Description des branches.

Chaque partie comprend un certain nombre d'éléments, dont la structure de stockage est organisée en lignes. Le premier caractère de la ligne identifie le type de ligne, donc l'élément auquel la ligne se rapporte. Dans cette version du logiciel, les caractères suivants sont "réservés" :

! identifie un réseau
$ identifie un noeud
# identifie un bief
N identifie une section
% identifie une branche
* identifie une ligne commentaire
+ identifie une ligne suite.

Attention !
Dans les descriptions suivantes, chaque virgule est traduite en un caractère espace dans le fichier, et chaque ligne est terminée par un retour chariot (caractère '\r\n en c).

Description du réseau (1 ligne)

' !', titre, version, cohérence, iflagx, xlm, ylm
au format : A1, A40, F5.3, I1, I1, F10.4, F10.4
où :
version est la version d'EDITAL (ou EDISIC) utilisée pour générer ce fichier,
cohérence vaut 1 si le réseau est cohérent et 0 sinon,
iflagx : si vaut 1 les abscisses longitudinales sont à multiplier par 1000
xlm : facteur d'homothétie en largeurs pour tout le réseau
ylm : facteur d'homothétie en hauteur pour tout le réseau

Nota Bene : TALWEG pourrait tester "cohérence" pour analyser ou non le réseau ! Pour l'instant il re-vérifie systématiquement la cohérence.

Description d'un noeud

'$', nom, x, y, cohérence
au format : A1, A6, I5, I5, I1
où :
nom est le nom du noeud
x et y sont les coordonnées en pixels du noeud,
cohérence vaut 1 si le noeud est cohérent et 0 sinon (toujours 1 pour cette version).

Description d'un bief (2 lignes + les sections du bief)

1ère ligne

'#', nom branche, nom noeud amont, nom noeud aval, pas, sin, aff, nom, xlm, ylm
au format : A1, A8, A6, A6, I5, F5.3, I1, A40, F10.4, F10.4
où :
pas : est le pas de calcul
sin : est la sinuosité du bief = rapport de la longueur du lit moyen du bief / longueur du lit mineur du bief (entre 0 et 1).
aff : est un paramètre qui vaut 0 par défaut et 1 si l'on veut que toutes les sections de calcul de ce bief soient stockées sur le fichier résultat en régime transitoire (.RES). Cela permet d'avoir une plus grande précision dans la présentation des résultats en régime transitoire. C'est une nouveauté des versions à partir de la 3.3.
nom : le nom du bief
xlm : facteur d'homothétie en largeurs pour tout le bief
ylm : facteur d'homothétie en hauteur pour tout le bief

Nota Bene : dans l'ancienne version SIC 2.1 un signe - peut apparaître en 2ème colonne si le bief a été inversé. Cette particularité a disparu dans les fichiers des versions postérieures à 3.0. Maintenant le bief inversé est stocké de la même manière qu'un bief normal, avec ses vrais noeuds amont et aval.

2ème ligne

'+', numéro noeud amont, numéro noeud aval, x, y, cohérence
au format : A1, I5, I5, I5, I5, I1
où : numéro noeud est l'indice d'apparition du noeud dans la liste des noeuds qui précède,
x et y sont les coordonnées du "point milieu" du bief en pixels,
cohérence vaut 1 ou 0 si le bief est cohérent et 0 sinon.

Description d'une section (nombre de lignes variables)

1ère ligne

'N', nom, abscisse, type, singulière, display, cohérence, xlm, ylm
au format : A1, A14, F8.1, A1, A1, A1, I1, F10.4, F10.4
où :
abscisse est l'abscisse de la section,
type vaut :
C si la section est Circulaire
D si la section est Dalot
R si la section est Rectangulaire
P si la section est Puissance
T si la section est Triangulaire ou Trapézoïdale
L si la section est Largeur-cote
A si la section est Abscisse-cote

singulière : vaut S si la section est singulière
display : vaut D si la section est d'affichage
cohérence vaut 1 si la section est cohérente et 0 sinon
xlm : facteur d'homothétie en largeurs pour la section
ylm : facteur d'homothétie en hauteur pour la section

Le nombre de lignes dépend du type de section :

pour les sections paramétrées, le nombre de lignes est de 2,

pour les sections largeur-cote et abscisse-cote, le nombre de ligne est :
4+2*(nombre de points-1)/8
(sauf si aucun point auquel cas il n'y a que 2 lignes).

2ème ligne et suivantes

La deuxième ligne (et les lignes suivantes éventuelles) dépendent du type de section.

Section Circulaire
'+', cote fond, cote berge, rayon, cote de débordement des lits mineur - moyen
au format : A1, F7.2, F7.2, F7.2, 7X, F7.2

Section Rectangulaire
'+', cote fond, cote berge, largeur, cote de débordement des lits mineur - moyen
au format : A1, F7.2, F7.2, F7.2, 7X, F7.2

Section Dalot ou Trapézoïdale (ou Triangulaire)
'+', cote fond, cote berge, largeur, fruit, cote de débordement des lits mineur - moyen
au format : A1, F7.2, F7.2, F7.2, F7.2, F7.2

Section Puissance
'+', cote fond, cote berge, largeur, coefficient, cote de débordement des lits mineur - moyen
au format : A1, F7.2, F7.2, F7.2, F7.2, F7.2

Section Abscisse-cote et Largeur-cote
2ème ligne
'+', référence élévation, cote de débordement des lits mineur - moyen

au format : A1, F7.2, F7.2

lignes suivantes (par paires) (tant qu'il y a des points)
'+', (abscisse ou largeur [I] I = 1 à 8)
'+', (cote [I] I = 1 à 8)
au format : 1A, (F7.2), x 8

Description d'une branche

1ère ligne

'%', nom, couleur, cohérence
au format : A1, A8, I10, I1
où :
couleur est la couleur RGB de la branche
cohérence vaut 1 si la branche est cohérente et 0 sinon.

Le nombre de lignes dépend du nombre de biefs que contient la branche. Il vaut :
2 + (nombre de biefs - 1)/12
si la branche comporte au moins un bief et 1 sinon.

Chaque ligne suivante est composée de :
'+' (numéro du bief I = 1 à 12)
au format : (A1, (I5), x 12)
où :
numéro du bief est l'ordre d'apparition du bief dans le fichier EDITAL.

Format du fichier FLU (SIC version 4.33)

Pierre-Olivier Malaterre — 2018-09-27T13:34:57Z

Ce fichier ASCII est créé par EDIFLU ou au moyen d'un éditeur de fichier ASCII quelconque. Il contient les données hydrauliques pour le programme de calcul en régime permanent FLUVIA. Il sera vérifié au moyen du programme EDIFLU et un fichier (.DON) destiné à FLUVIA sera créé.

Ce fichier est destiné aux anciennes versions de SIC 4.33 ou antérieures. Ces fichiers peuvent être importés à partir de SIC versions 5 et ultérieures. C'est pour cette raison que ce format, devenu obsolète est néanmoins décrit. Ils peuvent être importés en mode interface, pour créer un nouveau scénario, ou également en mode batch (un ou plusieurs fichiers FLU), en complément d'un fichier de géométrie TAL.

Aucun ordre spécial n'est suivi pour écrire ce fichier de données. Les différents types de lignes possibles sont décrits ci-dessous, chaque ligne commençant avec un emblème particulier.

La ligne du commentaire

Contient un caractère "*" dans la première colonne et un commentaire à la suite.
Ligne titre_etude
' !', TITM, VERNUM au format (A1,1X,A40,1X,F5.3,I1)
' !' en colonne 1,
TITM titre de l'étude,
VERNUM numéro de version d'EDIFLU utilisée pour créer ce fichier
IFLAGQ mode grands débits (0 si non, 1 si oui ie les débit seront multipliées par 1000)

Ligne pour les paramètres de calcul

Ces données étaient stockées dans le fichier .LDA dans les versions antérieures à 3.5. Ce fichier n'est plus utilisé maintenant, ces informations étant directement dans le fichier .FLU.
'X', XLAMBQ, XLAMBZ, ERQ, ERZ, NMI, NMAX, IDHA, IWRIT, IPAS au format (A1, 4F10.5, 2I5, 3(1X, I1))
XLAMBQ coefficient_Lambda de norme_erreur_Q_apport (prises)
XLAMBZ coefficient_Lambda de norme_erreur_Z1 (diffluences)
ERQ erreur maxi sur la continuité aux noeuds avec prises
ERZ demi-erreur maxi aux diffluences= (Z_amont - Z_moyenne)
NMI nombre maxi d'itérations
NMAX nombre maxi d'itérations de calage
IDHA indice de perte de charge automatique aux élargissements (0 : non, 1 : oui)
IWRIT pour avoir plus de messages dans le fichier .LST créé par Fluvia (mode Debug)
IPAS pour avoir le petit pas automatique activé quand les pentes de l'eau sont fortes.

Ligne Seuil

'D', NB, NS, NUM, XL, ZD, COF, IFOL, IFOC au format (A1, I3, I4, I3, 2F10.4, 10X, F10.4, 30X, 2I2)
D in colonne 1
NB = nombre du bief
NS = nombre de la section dans le bief
NUM = numéro de l'ouvrage
XL = largeur du radier
ZD = élévation du radier
COF = coefficient de débit
IFOL = équation (0 : Cemagref déversoir/orifice, 1 : Cunge, 2 : Dénoyée, 3 : Cemagref déversoir/vanne de fond, 4 : Matlab, 5 : Fortran). Les modes 4 et 5 ne sont pas encore gérés
IFOC = 0 (pas utilisé pour l'instant)

Ligne Vanne

'V', NB, NS, NUM, XL, ZD, W, COF, Z1, HS, IFOL, IFOC au format (A1, I3, I4, I3, 6F10.4, 10X, 2I2)
V en colonne 1
NB, NS, NUM, XL, ZD même signification que ci-dessus
W = Ouverture si IFOC=0 et ouverture maximum si IFOC=1 (bouge ouverture) ou incrément cote radier maximum si IFOC=2 (bouge cote radier)
COF = coefficient de débit
Z1 = élévation objectif si IFOC=1 ou 2
HS = hauteur de la vanne (utile en cas de débordement)
IFOL = équation (0 : Cemagref déversoir/vanne de fond, 1 : Cunge, 2 : Noyée, 3 : Cemagref déversoir/orifice, 4 : Matlab, 5 : Fortran). Les modes 4 et 5 ne sont pas encore gérés
IFOC = 0 aux vannes et 1 aux vannes ajustables en régulation amont

Ligne Vanne Gec-Alsthom

'G', NB, NS, NUM, XL, ZD, W, COF, Z1, HS, S1S2, IDNS, IFOL, IFOC, NOT au format (A1, I3, I4, I3, 6F10.4, F7.1, I3, 3I2)
G en colonne 1
NB, NS, NUM = même signification que ci-dessus
XL = Rayon du tablier de la vanne (R)
ZD = Cote du radier
W = Ouverture (utilisé si IFOL = 2) ou sinon ouverture verticale max
COF = Hauteur d'eau amont max (ZMAX pour AMIL) par rapport à la cote de l'axe (entre 0 et d)
Z1 = Cote de l'axe
HS = Décrément (d)
S1S2 = Ratio S1/S2 surface libre / surface fente pour bac vannes AVIS, AVIO et Mixtes
IDNS = Décalage de la section pour Z2 (utile pour vannes AVIS, AVIO et Mixtes)
IFOL = Type de loi (0 Goussard, 1 Cemagref non bloqué, 2 Cemagref bloqué)
IFOC = type de vanne. AMIL, AVIS, AVIO, haute chute, basse chute (0 Amil, 1 Avio Haute Chute, 2 Avio BC, 3 Avis HC, 4 Avis BC, 5 Mixte 45, 6 Mixte 55, 7 Mixto 45)
NOT = Numéro du type (divers modèles AVIO, AVIS ET AMIL)

Ligne Prise

'P', NOEUD, TYPEP, IFPD, TYPEC, QP, XL, ZD, WMAX, W, COF, Z2, XMV, Q0V, Z0V, IFPL
au format (A1, A6, A1, I1, A1, F10.3, F9.3, F7.2, F8.3, F8.3, F4.2, F7.2, F4.2, F6.2, F7.2, I2)
P en colonne 1
NOEUD = nom du noeud sur 6 caractères des colonnes 2 à 7
TYPEP = type de prise sur 1 caractère en colonne 8
Q = débit imposé
V = vanne
D = déversoir
A = APM
P = tube
L = courbe de tarage
IFPD = condition aval sur 1 caractère en colonne 9
0 = rien
1 = cote constante
2 = seuil
3 = courbe de tarage
TYPEC = mode de calcul sur 1 caractère en colonne 10
F = off
N = on
Q = calcul de débit
W = calcul d'ouverture
Z = calcul de la cote de seuil
L = calcul de la largeur de seuil
I = débit imposé
QP = débit à la prise sur 10 caractères en colonnes 11 à 20
XL = largeur ou diamètre sur 9 caractères en colonnes 21 à 29
ZD = cote de seuil sur 7 caractères en colonnes 30 à 36
WMAX = ouverture maximum sur 8 caractères en colonnes 37 à 44
W = ouverture sur 8 caractères en colonnes 45 à 52
C0F = coefficient de débit sur 4 caractères en colonnes 53 à 56
Z2 = cote aval ou cote du seuil aval sur 7 caractères en colonnes 57 à 63
XMV = coefficient du seuil aval ou puissance sur 4 caractères en colonnes 64 à 67
Q0V = débit ou largeur sur 6 caractères en colonnes 68 à 73
Z0V = cote correspondant à Q0V sur 7 caractères en colonnes 74 à 80.
IFPL = équation (0 : Cemagref ou 1 : Cunge)

Ligne pour les coefficients de Strickler

'K', NB, NI, NS, SK au format (A1, I3, 2I4, F10.4)
K en colonne 1
NB = numéro du bief sur 3 caractères en colonnes 2 à 4
NI = numéro de la première section de calcul sur 4 caractères en colonnes 5 à 8
NS = numéro de la dernière section de calcul sur 4 caractères en colonnes 9 à 12
SK = valeur du coefficient de Strickler de la section NI à NS sur 10 caractères en colonnes 13 à 22

Ligne pour les débits latéraux

'L', NB, NI, NS, QL au format (A1, I3, 2I4, F10.4)
L en colonne 1
NB = numéro du bief sur 3 caractères en colonnes 2 à 4
NI = numéro de la première section de calcul sur 4 caractères en colonnes 5 à 8
NS = numéro de la dernière section de calcul sur 4 caractères en colonnes 9 à 12
QL = valeur du débit latéral de la section NI à NS sur 10 caractères en colonnes 13 à 22

Ligne pour les conditions aval

'A', NOEUD, QAV, ZAV au format (A1, A6, 3X, 2F10.3)
A en colonne 1
NOEUD = nom du noeud aval de modèle sur 6 caractères en colonnes 2 à 7
QAV = débit sur 10 caractères en colonnes 11 à 20
ZAV = cote sur 10 caractères en colonnes 21 à 30

Ligne pour les diffluences

'B', NB, X au format (I3, F10.3)
B en colonne 1
NB = numéro du bief
X = débit correspondant (en m3/s, % ou surface des biefs aval. Seules les valeurs relatives sont importantes)

Ligne pour cotes de calage

'C', NB, NS, COT au format (2I3, F8.2)
C en colonne 1
NB = numéro du bief
NS = numéro de la section
COT = cote de calage correspondante

Ligne pour le mode de calcul aux noeuds

Par défaut, Fluvia fait le calcul en écrivant l'égalité des cotes aux noeuds. L'option suivante propose d'utiliser à la place l'égalité des charges aux noeuds, avec éventuellement introduction d'une perte de charge en entrée et en sortie de bief. Cette perte de charge est de la forme DH=K.V*V/(2g).
'H', I, PCAM(I), PCAV(I) au format (A1,I3,6X,2F10.5)
H en colonne 1
I = numéro du bief
PCAM(I) = coefficient de perte de charge à l'amont du bief
PCAV(I) = coefficient de perte de charge à l'aval du bief

Changement de la bathymétrie avec un fichier externe

Pierre-Olivier Malaterre — 2017-11-24T16:42:46Z

On veut parfois changer la bathymétrie d'un canal ou d'un fleuve à partir d'un programme externe (ex. : script Matlab, Scilab, R, etc) pour faire de l'optimisation de profils par exemple, ou pour calculer la sensibilité de certaines variables à des paramètres de bathymétrie.

On peut faire cela avec un fichier nommé changebathy.txt, qui doit être localisé avec le fichier du projet xml. L'exemple \dat\ex32_lacguiers\LacGuiers_v1.xml utilise cette fonctionnalité.

Ce fichier a la structure suivante :

 27 2 0.00000 0.00000 0.10000 
 28 2 0.00000 0.00000 0.10000 
 28 3 0.00000 0.00000 0.15000

On peut avoir autant de lignes que l'on souhaite. Le format est libre ('*' en format Fortran), c'est à dire que les nombres peuvent être séparés par autant d'espaces que souhaité, ou a priori aussi des tabulations.

Sur chaque ligne :

le premier nombre est le numéro du bief
le second nombre est le numéro de la section dans ce bief
le troisième est le facteur d'homothétie en Y (latérale) que l'on veut appliquer
le quatrième est le facteur d'homothétie en Z (verticale) que l'on veut appliquer
le cinquième est le décalage vertical de la cote de fond, et donc de toute la section (en m)

Lors de l'exécution du programme de calcul (Fluvia ou Sirene) les sections ainsi modifiées seront indiquées dans la fenêtre de contrôle.

Les conditions aux limites

David Dorchies — 2017-09-24T10:20:43Z

Les prises sont obligatoirement pourvues d'une condition aux limites. Les conditions aux limites disponibles sont :

Débit imposé (Pompe) : Q(t)
Cote fixe : Z(t)
Débit fonction de la cote (tabulée) : Q(Z) (ou plutôt Z(Q))
Fonction Q puissance alpha : $Q(t)=Q_{ref} [(Z(t)-Z_{0})/(Z_{ref}-Z_{0})]^{\alpha}$

Les nœuds amont et aval sont considérés également comme des endroits où il faut définir des "prises", qui seront utilisées comme conditions aux limites dans les algorithmes de calcul hydraulique. La terminologie peut sembler abusive mais il faut bien choisir un terme générique. Cela est défendable, dans le sens où ce sont les points qui permettent d'injecter ou retirer du débit dans le réseau hydraulique, aussi bien aux nœuds amont et aval, qu'aux nœuds intermédiaires.

Les conditions aux nœuds amont et aval peuvent être choisis dans cette liste. Les détails donnés pour les noeuds aval sont donc applicables aussi ici pour un noeud amont ou intermédiaire.

Toutes les combinaisons sont possibles en régime transitoire.

Quelques combinaisons sont impossibles en régime permanent : Q(t) à l'amont et à l'aval. En effet dans ce cas soit il n'y a pas de solution, soit il y a une infinité de solutions. La configuration la plus classique est Q(t) à l'amont et Q(Z) à l'aval. Les autres options qui n'indiquent pas le débit à l'amont de manière explicite impliqueront des itérations du programme de calcul Fluvia afin de trouver le débit qui satisfait la conditions amont choisie.

Rappel : les débits entrants ou sortants, dans SIC, suivent une convention de signe qui est que les débits entrants sont positifs et les débits sortants sont négatifs. Cela est valable généralement aux conditions aux limites externes (amont et aval), internes (prises) et les infiltrations ou apports répartis le long des sections. Mais, pour éviter les erreurs classiques qu'on rencontrait souvent avec des utilisateurs, on fait des exceptions à cette convention quand il n'y a pas d'ambiguïté. Par exemple une condition limite aval de type Q(Z) ne peut être que sortante, et on attend donc des débits positifs aux interfaces de saisie des données alors que selon la convention il faudrait les choisir négatifs. La calculette les choisira positifs, mais ce seront bien des débits sortants.

Initialisation du calcul aux diffluences

David Dorchies — 2017-09-24T09:57:24Z

Sur un réseau pourvu de diffluences (cas courant sur les canaux d'irrigation contrairement aux rivières), le calcul de la ligne d'eau en régime permanent effectue des itérations pour distribuer le débit provenant de l'Amont dans chaque Bief Aval diffluent.

Par défaut, la répartition est initialisée à partir de la surface des biefs situés à l'aval de la diffluence. Dans le cas où cette répartition par défaut provoquerait des débits initiaux négatifs dans certaines Branches du réseau, l'utilisateur peut modifier ces paramètres en double-cliquant sur le Nœud posant problème et modifier la répartition du débit entre les Biefs à l'Aval du Nœud.

Définition du mode de calcul au noeud

David Dorchies — 2017-09-24T09:46:46Z

La résolution du calcul de la ligne d'eau aux nœuds peut se faire selon deux modes :

égalité des cotes (mode par défaut)
égalité des charges

Dans le cas de l'égalité des charges, il est possible de définir une perte de charge à l'entrée du nœud (pour les biefs Amont arrivant à ce nœud) et à la sortie du nœud (pour les biefs Aval quittant ce nœud). Il est possible d'avoir pour le même nœud des biefs arrivant à ou partant de ce nœud avec des modes différents (cote pour certains et charge pour d'autres). Dans ce cas d'une perte de charge une équation de type Borda (proportionnelle à la vitesse au carré) est utilisée.

Remarque : Dans le cas où l'écoulement à un nœud est en régime torrentiel sur un ou plusieurs des biefs arrivant ou partant de ce nœud, le mode de calcul sera basculé automatiquement en mode "égalité des charges", même si c'était le mode "égalité des cotes" qui avait été choisi ici. C'est en effet plus réaliste physiquement. Pour que les calculs en régime torrentiel soient autorisés il faut activer l'option correspondante dans les paramètres de calcul permanent.

Géolocalisation

Pierre-Olivier Malaterre — 2017-04-28T10:05:44Z

Cet outil permet de gérer la géolocalisation des objets du réseau (nœuds, points de préemption des biefs aussi appelés points milieu, sections en travers).

La partie supérieure de l'onglet "Géolocalisation" permet de traiter les nœuds et les points-milieu des biefs (appelés aussi points de préemption puisqu'ils permettent de sélectionner un bief, pour en éditer les paramètres ou en visualiser les résultats).

L'option "Calculer pxl -> geo" va utiliser les données des pixels (position x-y des objets graphiques sur le dessin en plan 2D) pour calculer la géolocalisation (lat-lon) des objets traités (nœuds et points de préemption ici). Cela se fera en utilisant la géolocalisation des nœuds amont et aval de la zone sélectionnée. Il faut donc renseigner cela avant d'utiliser cette option. Cela se fait par de simples relations de proportionnalité entre x-y et lat-lon (dite règle de 3). Ces données de géolocalisation seront ensuite stockées dans les champs correspondants des objets correspondants.

L'option "Calculer geo -> pxl" va faire le contraire, c'est à dire positionner (x-y) les objets sur le dessin en plan 2D en tenant compte des données de géolocalisation (lat-lon), en conservant les nœuds amont et aval comme références (x1-y1 et x2-y2) et en utilisant de simples règles de 3 pour calculer leurs positions. Il faut donc positionner ces points amont et aval où on le souhaite, par exemple sur un fond de carte importé auparavant (voir la page dédiée). Ces données seront stockées et les objets nœuds et points milieu seront ensuite dessinés à ces endroits sur le plan 2D.

La partie inférieure de l'onglet "Géolocalisation" permet de traiter les sections en travers. Cette option est proche de la précédente dans les idées derrière "Calculer pxl -> geo" et "Calculer geo -> pxl". Cependant les points de référence seront pris pour chaque bief avec son nœud amont et son nœud aval (et non pas comme précédemment le nœud amont de la sélection et le nœud aval de la sélection). Il faut donc les avoir traité précédemment avec les options de la partie supérieure de l'onglet, détaillées ci-dessus. En outre les positions (x-y) des sections en travers seront calculées pour l'option "Calculer geo -> pxl", et dessinées ainsi juste après avoir lancé ce calcul (cliqué sur le bouton correspondant), mais comme de base les biefs sont dessinés sur le plan 2D comme 2 segments de droite entre le nœud amont et le point milieu puis le point milieu et le point aval, et que les sections sont dessinées et réparties sur ces segments en fonction de leurs abscisses longitudinales au sein de leur bief (les coordonnées des pixels x-y seront donc recalculées à chaque affichage), ce dessin en plan détaillé des sections ne sera pas dessiné ainsi aux affichages ultérieurs.

Une évolution ultérieure de SIC autorisant plus de segments et potentiellement autant que les sections en travers permettrait de gérer cela plus esthétiquement.

Dans cette gestion de la géolocalisation des sections en travers, on peut/doit aussi donner d'autres informations, puisque la section entière est géolocalisée (1 seule donnée lat-lon pour chaque section) et non pas les divers points de la section. Cependant une section vue de près n'est pas ponctuelle (elle peut être de plusieurs centaines de m ou km), mais constituée par un ensemble de points. Il faut donc indiquer à quel point correspond la géolocalisation donnée, ainsi que comment se situe la section par rapport à ce point de géolocalisation, latéralement en angulairement. L'angle est donné par rapport à une référence indiquée dans les options globales de géolocalisation (Edition/Propriétés du Réseau/Paramètres de Géolocalisation), avec les choix du Nord Géographique, de la ligne centrale du bief ou du segment dessiné au sein du bief. L'"abscisse latérale Y" est donnée dans le cas du choix "pour une abscisse latérale donnée". Si on coche les boites pour ces valeurs elles seront ainsi renseignées pour les objets correspondants lors du lancement de calcul "Calculer pxl -> geo". Actuellement ces données sont stockées mais ne sont pas prises en compte lors du dessin des sections en travers (à faire dans une version ultérieure). Pour l'instant l'angle utilisé pour le dessin est de 90° par rapport au segment, et le point de géolocalisation supposé est celui central de la section.

En bas l'onglet on indique une échelle en X et en Y, faisant la correspondance entre les lat-lon et les x-y. Ces champs permettent de visualiser ces échelles après avoir utilisé une des 2 options de modification ci-dessus, mais pas de forcer un dessin à une échelle donnée (ce sont des variables en sortie pas en entrée de calcul). Comme indiqué par la bulle d'aide quand on positionne la souris sur les champs correspondants, les valeurs indiquées sont les nombres de pixels pour 1 degré d'angle.

Un exemple de projet géolocalisé est donné à :
\dat\ex30_garonne\Garonne_benchmark_MinMoy_v4.xml