Merge branch 'dev' into 'master'

Upgrade VGEST until V8

See merge request !1

.gitignore

@@ -9,3 +9,5 @@
 # Program outputs
 /compteur
 /*.TXT
+/FiniFini
+/*Optim

README.md

@@ -16,13 +16,15 @@ At the root folder of this repository type the following command line (Adapt the
 
 ```
 C:\FPC\bin\i386-win32\fpc.exe src\VGEST.PAS && move /Y src\VGEST.exe .
+C:\FPC\bin\i386-win32\fpc.exe src\ENDE.PAS && move /Y src\ENDE.exe .
+C:\FPC\bin\i386-win32\fpc.exe src\REGLAGE.PAS && move /Y src\REGLAGE.exe .
 ```
 
 ## Get started with the example
 
 The `example` directory contains all the configuration for calculating statistics about respecting a maximum flow of 110 m<sup>3</sup>/s at Paris. These statistics are calculated from the Seine River Basin naturalised flow database between 1993/01/01 and 2009/12/31 [^Hydratec2011].
 
-Copy the files contained in the `example` directory to the root directory and run `VGEST.exe`.
+Copy the files contained in the `example` directory to the root directory and run `SURVGEST.bat`.
 
 ## Known issues
 

SUIVOPTI/.gitignore

+# Ignore everything in this directory
+*
+# Except this file
+!.gitignore

SURVGEST.bat

+if exist FiniFini del FiniFini
+echo off
+if exist PasOptim del PasOptim
+if exist FinOptim del FinOptim
+if exist PARAMETR\reglage.txt del PARAMETR\reglage.txt
+if exist PARAMETR\OUITRI.TXT del PARAMETR\OUITRI.TXT
+if exist Temporai.re\etape del Temporai.re\etape
+if exist Temporai.re\result0 del Temporai.re\result0;
+if exist Suivopti\suivi del Suivopti\suivi
+if exist Temporai.re\dichot1 del Temporai.re\dichot1
+if exist Temporai.re\dichot2 del Temporai.re\dichot2
+if exist Temporai.re\dichot3 del Temporai.re\dichot3
+if exist Temporai.re\dichot4 del Temporai.re\dichot4
+:NouveauCalcul
+REGLAGE
+if exist PasOptim goto FindesCalculs
+VGEST
+if exist FinOptim goto FinDesCalculs
+goto NouveauCalcul
+:FinDesCalculs
+del PARAMETR\reglage.txt
+ENDE
\ No newline at end of file

example/CHOIX.TXT

-PARIS_05
-Q_OBS_OBS.txt
-BATCH
-1
-01/01/1994
-31/12/2009

example/MODE.TXT

-1

example/PARAMETR/AUBE.1 → example/PARAMETR/BARRAGES/1/AUBE.1

-paramètres des contraintes et consignes du réservoir Aube (scénario : réalité d'après rapport Florine Dehay, ENGEES 2012)
+contraintes et consignes du réservoir Aube (règles actuelles d'après F.Dehay, ENGEES 2012). Débits réservé et de référence considérés (chacun de façon identique) à la prise et à la restitution 
 AUBE									
-0			Qemin	(m3/s)							
-1.35E+02		Qemax	(m3/s)							
-0			Qsmin	(m3/s)							
-35			Qsmax	(m3/s)							
+1	0	0		Qemin	(m3/s)	fonction de V (hm3)							
+1	0	1.35E+02	Qemax	(m3/s)	fonction de V (hm3)						
+1	0	0		Qsmin	(m3/s)	fonction de V (hm3)						
+1	0	35		Qsmax	(m3/s)	fonction de V (hm3)						
 181.2			Vtot	(million de m3)							
 0	3		QMres						
 1	91	305							
@@ -25,21 +25,25 @@ AUBE
 1	-1	130							
 									
 FORMAT DU FICHIER :									
-ligne 0 : en-tête obligatoire									
-ligne 1 : nom du réservoir (chaîne de 8 caractères)									
-lignes 2 à 6 : Qemin, Qemax, Qsmin, Qsmax, V, 									
+ligne 1 : en-tête obligatoire (version)
+ligne 2	: nom du réservoir								
+ligne 3 : barème V-> Qemin (effectif (entre 1 et 10), valeurs de V (hm3) croissantes, valeurs de Qemin (m3/s))
+ligne 4 : barème V-> Qemax (effectif, valeurs de V (hm3) croissantes, valeurs de Qemax (m3/s))
+ligne 5 : barème V-> Qsmin (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmin (m3/s))
+ligne 6 : barème V-> Qsmax (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmax (m3/s))
+ligne 7 : volume Vtot (hm3, valeur max possible de V) 									
 4 groupes de lignes (pour QMres, QVres, QMref et QVref respectivement) composé chacun de :									
 	ligne 1 : Type (entre 0 et 2), Nombre N de dates pivot (entre 1 et 10) et Borne (si Type <>0). Le débit vaut alors Q défini dans les lignes suivantes si Type=0, max(Q,QV(t-1)-borne) si Type=1 et min(Q,QV(t-1)+borne) si Type = 2								
-	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, inférieurs à 366)								
+	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, compris entre 1 et 365)								
 	ligne 3 : N valeurs entières comprises entre 1 et M(<=N), correspondant aux catégories de valeurs de Q								
 	lignes 4 à 3+M : nombre X (entier entre 1 et 5), puis X couples de valeurs Ai (décroissant),Bi interprétés comme suit : si QMnat>A1 alors Q=B1, sinon, si QMnat>A2 alors Q=B2, sinon…,si QMnat>Ax alors Q=Bx, sinon Q=QMnat								
-notations :	
-	Qemin et Qemax : débitance min et max du canal de prise								
-	Qsmin et Qsmax : débitance min et max du canal de restitution								
-	QMres et QMref : débit réservé et débit de référence, juste à l'aval de la prise								
-	QVres et QVref : débit réservé et débit de référence juste à l'aval de la restitution								
-	Vtot : volume total du réservoir								
-	QMnat : débit observé juste à l'amont de la prise								
-	QV : débit observé juste à l'aval de la restitution										
+NOTATIONS :
+	V : volume d'eau exploitable présent dans le réservoir	
+	Qemin et Qemax : débitance min et max de la prise (limites physiques)								
+	Qsmin et Qsmax : débitance min et max de la restitution	(limites physiques)							
+	QMres et QMref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la prise								
+	QVres et QVref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la restitution								
+	Vtot : volume exploitable total du réservoir, compris entre le sommet de la tranche morte (ou le sommet du volume de garde) et le niveau maximal de sécurité								
+	QMnat : débit du cours d'eau à l'amont immédiat de la prise								
+	QV : débit du cours d'eau à l'aval immédiat de la restitution										
 									
-

example/PARAMETR/MARNE.1 → example/PARAMETR/BARRAGES/1/MARNE.1

-paramètres des contraintes et consignes du réservoir Marne (scénario : réalité d'après rapport Florine Dehay, ENGEES 2012)
+contraintes et consignes du réservoir Marne (règles actuelles d'après F.Dehay, ENGEES 2012). Débits réservé et de référence considérés (de façon identique pour le premier) à la prise et à la restitution 
 MARNE										
-0			Qemin	(m3/s)								
-4.08E+02		Qemax	(m3/s)								
-0			Qsmin	(m3/s)								
-50			Qsmax	(m3/s)								
+1	0	0		Qemin	(m3/s)	fonction de V (hm3)							
+1	0	4.08E+02	Qemax	(m3/s)	fonction de V (hm3)							
+1	0	0		Qsmin	(m3/s)	fonction de V (hm3)							
+1	0	50		Qsmax	(m3/s)	fonction de V (hm3)							
 354.5			Vtot	(million de m3)								
 0	3		QMres							
 1	91	305								
@@ -31,20 +31,25 @@ MARNE
 1	-1	120								
 										
 FORMAT DU FICHIER :									
-ligne 0 : en-tête obligatoire									
-ligne 1 : nom du réservoir (chaîne de 8 caractères)									
-lignes 2 à 6 : Qemin, Qemax, Qsmin, Qsmax, V, 									
+ligne 1 : en-tête obligatoire (version)
+ligne 2	: nom du réservoir								
+ligne 3 : barème V-> Qemin (effectif (entre 1 et 10), valeurs de V (hm3) croissantes, valeurs de Qemin (m3/s))
+ligne 4 : barème V-> Qemax (effectif, valeurs de V (hm3) croissantes, valeurs de Qemax (m3/s))
+ligne 5 : barème V-> Qsmin (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmin (m3/s))
+ligne 6 : barème V-> Qsmax (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmax (m3/s))
+ligne 7 : volume Vtot (hm3, valeur max possible de V) 									
 4 groupes de lignes (pour QMres, QVres, QMref et QVref respectivement) composé chacun de :									
 	ligne 1 : Type (entre 0 et 2), Nombre N de dates pivot (entre 1 et 10) et Borne (si Type <>0). Le débit vaut alors Q défini dans les lignes suivantes si Type=0, max(Q,QV(t-1)-borne) si Type=1 et min(Q,QV(t-1)+borne) si Type = 2								
-	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, inférieurs à 366)								
+	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, compris entre 1 et 365)								
 	ligne 3 : N valeurs entières comprises entre 1 et M(<=N), correspondant aux catégories de valeurs de Q								
 	lignes 4 à 3+M : nombre X (entier entre 1 et 5), puis X couples de valeurs Ai (décroissant),Bi interprétés comme suit : si QMnat>A1 alors Q=B1, sinon, si QMnat>A2 alors Q=B2, sinon…,si QMnat>Ax alors Q=Bx, sinon Q=QMnat								
-notations :	
-	Qemin et Qemax : débitance min et max du canal de prise								
-	Qsmin et Qsmax : débitance min et max du canal de restitution								
-	QMres et QMref : débit réservé et débit de référence, juste à l'aval de la prise								
-	QVres et QVref : débit réservé et débit de référence juste à l'aval de la restitution								
-	Vtot : volume total du réservoir								
-	QMnat : débit observé juste à l'amont de la prise								
-	QV : débit observé juste à l'aval de la restitution										
+NOTATIONS :
+	V : volume d'eau exploitable présent dans le réservoir	
+	Qemin et Qemax : débitance min et max de la prise (limites physiques)								
+	Qsmin et Qsmax : débitance min et max de la restitution	(limites physiques)							
+	QMres et QMref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la prise								
+	QVres et QVref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la restitution								
+	Vtot : volume exploitable total du réservoir, compris entre le sommet de la tranche morte (ou le sommet du volume de garde) et le niveau maximal de sécurité								
+	QMnat : débit du cours d'eau à l'amont immédiat de la prise								
+	QV : débit du cours d'eau à l'aval immédiat de la restitution										
 									

example/PARAMETR/SEINE.1 → example/PARAMETR/BARRAGES/1/SEINE.1

-paramètres des contraintes et consignes du réservoir Seine (scénario : réalité d'après rapport Florine Dehay, ENGEES 2012)
+contraintes et consignes du réservoir Seine (règles actuelles d'après F.Dehay, ENGEES 2012). Débits réservé et de référence considérés (chacun de façon identique) à la prise et à la restitution 
 SEINE									
-0			Qemin	(m3/s)							
-2.00E+02		Qemax	(m3/s)							
-0.8			Qsmin	(m3/s)							
-35			Qsmax	(m3/s)							
+1	0	0		Qemin	(m3/s) fonction de V (hm3)							
+1	0	2.00E+02	Qemax	(m3/s) fonction de V (hm3)							
+1	0	0.8		Qsmin	(m3/s) fonction de V (hm3)							
+1	0	35		Qsmax	(m3/s) fonction de V (hm3)							
 212.9			Vtot	(million de m3)							
 0	4		QMres						
 1	91	274	335						
@@ -29,20 +29,25 @@ SEINE
 1	-1	40							
 									
 FORMAT DU FICHIER :									
-ligne 0 : en-tête obligatoire									
-ligne 1 : nom du réservoir (chaîne de 8 caractères)									
-lignes 2 à 6 : Qemin, Qemax, Qsmin, Qsmax, V, 									
+ligne 1 : en-tête obligatoire (version)
+ligne 2	: nom du réservoir								
+ligne 3 : barème V-> Qemin (effectif (entre 1 et 10), valeurs de V (hm3) croissantes, valeurs de Qemin (m3/s))
+ligne 4 : barème V-> Qemax (effectif, valeurs de V (hm3) croissantes, valeurs de Qemax (m3/s))
+ligne 5 : barème V-> Qsmin (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmin (m3/s))
+ligne 6 : barème V-> Qsmax (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmax (m3/s))
+ligne 7 : volume Vtot (hm3, valeur max possible de V) 									
 4 groupes de lignes (pour QMres, QVres, QMref et QVref respectivement) composé chacun de :									
 	ligne 1 : Type (entre 0 et 2), Nombre N de dates pivot (entre 1 et 10) et Borne (si Type <>0). Le débit vaut alors Q défini dans les lignes suivantes si Type=0, max(Q,QV(t-1)-borne) si Type=1 et min(Q,QV(t-1)+borne) si Type = 2								
-	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, inférieurs à 366)								
+	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, compris entre 1 et 365)								
 	ligne 3 : N valeurs entières comprises entre 1 et M(<=N), correspondant aux catégories de valeurs de Q								
 	lignes 4 à 3+M : nombre X (entier entre 1 et 5), puis X couples de valeurs Ai (décroissant),Bi interprétés comme suit : si QMnat>A1 alors Q=B1, sinon, si QMnat>A2 alors Q=B2, sinon…,si QMnat>Ax alors Q=Bx, sinon Q=QMnat								
-notations :	
-	Qemin et Qemax : débitance min et max du canal de prise								
-	Qsmin et Qsmax : débitance min et max du canal de restitution								
-	QMres et QMref : débit réservé et débit de référence, juste à l'aval de la prise								
-	QVres et QVref : débit réservé et débit de référence juste à l'aval de la restitution								
-	Vtot : volume total du réservoir								
-	QMnat : débit observé juste à l'amont de la prise								
-	QV : débit observé juste à l'aval de la restitution										
+NOTATIONS :
+	V : volume d'eau exploitable présent dans le réservoir	
+	Qemin et Qemax : débitance min et max de la prise (limites physiques)								
+	Qsmin et Qsmax : débitance min et max de la restitution	(limites physiques)							
+	QMres et QMref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la prise								
+	QVres et QVref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la restitution								
+	Vtot : volume exploitable total du réservoir, compris entre le sommet de la tranche morte (ou le sommet du volume de garde) et le niveau maximal de sécurité								
+	QMnat : débit du cours d'eau à l'amont immédiat de la prise								
+	QV : débit du cours d'eau à l'aval immédiat de la restitution										
 									

example/PARAMETR/YONNE.1 → example/PARAMETR/BARRAGES/1/YONNE.1

-paramètres des contraintes et consignes du réservoir Yonne (scénario : réalité d'après rapport Florine Dehay, ENGEES 2012)
+contraintes et consignes du réservoir Yonne (règles actuelles d'après F.Dehay, ENGEES 2012). Débits réservés et de référence considérés à la restitution et nuls à la prise. Qsmax=Qsmax1=16 réglementaire  
 YONNE									
-0			Qemin	(m3/s)							
-1.00E+20		Qemax	(m3/s)							
-0			Qsmin	(m3/s)							
-16			Qsmax	(m3/s)							
-74.37			Vtot	(million de m3)							
+1	0	1.00E+20	Qemin	(m3/s) fonction de V (hm3)							
+1	0	1.00E+20	Qemax	(m3/s) fonction de V (hm3)							
+1	0	0		Qsmin	(m3/s) fonction de V (hm3)							
+1	0	16		Qsmax	(m3/s) fonction de V (hm3)							
+74			Vtot	(million de m3)							
 0	1		QMres						
 1									
 1									
@@ -29,20 +29,25 @@ YONNE
 1	-1	12							
 									
 FORMAT DU FICHIER :									
-ligne 0 : en-tête obligatoire									
-ligne 1 : nom du réservoir (chaîne de 8 caractères)									
-lignes 2 à 6 : Qemin, Qemax, Qsmin, Qsmax, V, 									
+ligne 1 : en-tête obligatoire (version)
+ligne 2	: nom du réservoir								
+ligne 3 : barème V-> Qemin (effectif (entre 1 et 10), valeurs de V (hm3) croissantes, valeurs de Qemin (m3/s))
+ligne 4 : barème V-> Qemax (effectif, valeurs de V (hm3) croissantes, valeurs de Qemax (m3/s))
+ligne 5 : barème V-> Qsmin (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmin (m3/s))
+ligne 6 : barème V-> Qsmax (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmax (m3/s))
+ligne 7 : volume Vtot (hm3, valeur max possible de V) 									
 4 groupes de lignes (pour QMres, QVres, QMref et QVref respectivement) composé chacun de :									
 	ligne 1 : Type (entre 0 et 2), Nombre N de dates pivot (entre 1 et 10) et Borne (si Type <>0). Le débit vaut alors Q défini dans les lignes suivantes si Type=0, max(Q,QV(t-1)-borne) si Type=1 et min(Q,QV(t-1)+borne) si Type = 2								
-	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, inférieurs à 366)								
+	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, compris entre 1 et 365)								
 	ligne 3 : N valeurs entières comprises entre 1 et M(<=N), correspondant aux catégories de valeurs de Q								
 	lignes 4 à 3+M : nombre X (entier entre 1 et 5), puis X couples de valeurs Ai (décroissant),Bi interprétés comme suit : si QMnat>A1 alors Q=B1, sinon, si QMnat>A2 alors Q=B2, sinon…,si QMnat>Ax alors Q=Bx, sinon Q=QMnat								
-notations :	
-	Qemin et Qemax : débitance min et max du canal de prise								
-	Qsmin et Qsmax : débitance min et max du canal de restitution								
-	QMres et QMref : débit réservé et débit de référence, juste à l'aval de la prise								
-	QVres et QVref : débit réservé et débit de référence juste à l'aval de la restitution								
-	Vtot : volume total du réservoir								
-	QMnat : débit observé juste à l'amont de la prise								
-	QV : débit observé juste à l'aval de la restitution										
+NOTATIONS :
+	V : volume d'eau exploitable présent dans le réservoir	
+	Qemin et Qemax : débitance min et max de la prise (limites physiques)								
+	Qsmin et Qsmax : débitance min et max de la restitution	(limites physiques)							
+	QMres et QMref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la prise								
+	QVres et QVref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la restitution								
+	Vtot : volume exploitable total du réservoir, compris entre le sommet de la tranche morte (ou le sommet du volume de garde) et le niveau maximal de sécurité								
+	QMnat : débit du cours d'eau à l'amont immédiat de la prise								
+	QV : débit du cours d'eau à l'aval immédiat de la restitution										
 									

example/PARAMETR/BARRAGES/1/version.txt

+2012, Qres et Qref prise et restit, QSmax(Yonne)=QSmax1 
\ No newline at end of file

example/PARAMETR/BARRAGES/2/AUBE.1

+Contraintes et consignes du réservoir Aube (règles actuelles d'après F.Dehay, ENGEES 2012). Débit réservé considéré (de façon identique) à la prise et à la restitution. Débit de référence non considéré. 
+AUBE									
+1	0	0			Qemin	(m3/s) fonction de V (hm3)							
+1	0	1.35E+02		Qemax	(m3/s) fonction de V (hm3)							
+1	0	0			Qsmin	(m3/s) fonction de V (hm3)							
+1	0	35			Qsmax	(m3/s) fonction de V (hm3)							
+181.2			Vtot	(million de m3)							
+0	3		QMres						
+1	91	305							
+1	2	1							
+1	-1	2							
+3	8	2	6	3	4	4			
+0	3		QVres						
+1	91	305							
+1	2	1							
+1	-1	2							
+3	8	2	6	3	4	4			
+0	1		QMref						
+1									
+1									
+1	-1	9999							
+0	1		QVref						
+1									
+1									
+1	-1	9999							
+									
+FORMAT DU FICHIER :									
+ligne 1 : en-tête obligatoire (version)
+ligne 2	: nom du réservoir								
+ligne 3 : barème V-> Qemin (effectif (entre 1 et 10), valeurs de V (hm3) croissantes, valeurs de Qemin (m3/s))
+ligne 4 : barème V-> Qemax (effectif, valeurs de V (hm3) croissantes, valeurs de Qemax (m3/s))
+ligne 5 : barème V-> Qsmin (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmin (m3/s))
+ligne 6 : barème V-> Qsmax (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmax (m3/s))
+ligne 7 : volume Vtot (hm3, valeur max possible de V) 									
+4 groupes de lignes (pour QMres, QVres, QMref et QVref respectivement) composé chacun de :									
+	ligne 1 : Type (entre 0 et 2), Nombre N de dates pivot (entre 1 et 10) et Borne (si Type <>0). Le débit vaut alors Q défini dans les lignes suivantes si Type=0, max(Q,QV(t-1)-borne) si Type=1 et min(Q,QV(t-1)+borne) si Type = 2								
+	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, compris entre 1 et 365)								
+	ligne 3 : N valeurs entières comprises entre 1 et M(<=N), correspondant aux catégories de valeurs de Q								
+	lignes 4 à 3+M : nombre X (entier entre 1 et 5), puis X couples de valeurs Ai (décroissant),Bi interprétés comme suit : si QMnat>A1 alors Q=B1, sinon, si QMnat>A2 alors Q=B2, sinon…,si QMnat>Ax alors Q=Bx, sinon Q=QMnat								
+NOTATIONS :
+	V : volume d'eau exploitable présent dans le réservoir	
+	Qemin et Qemax : débitance min et max de la prise (limites physiques)								
+	Qsmin et Qsmax : débitance min et max de la restitution	(limites physiques)							
+	QMres et QMref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la prise								
+	QVres et QVref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la restitution								
+	Vtot : volume exploitable total du réservoir, compris entre le sommet de la tranche morte (ou le sommet du volume de garde) et le niveau maximal de sécurité								
+	QMnat : débit du cours d'eau à l'amont immédiat de la prise								
+	QV : débit du cours d'eau à l'aval immédiat de la restitution										
+									

example/PARAMETR/BARRAGES/2/MARNE.1

+Contraintes et consignes du réservoir Marne (règles actuelles d'après F.Dehay, ENGEES 2012). Débit réservé considéré (de façon identique) à la prise et à la restitution. Débit de référence non considéré. 
+MARNE										
+1	0	0			Qemin	(m3/s) fonction de V (hm3)								
+1	0	4.08E+02		Qemax	(m3/s) fonction de V (hm3)								
+1	0	0			Qsmin	(m3/s) fonction de V (hm3)								
+1	0	50			Qsmax	(m3/s) fonction de V (hm3)								
+354.5			Vtot	(million de m3)								
+0	3		QMres							
+1	91	305								
+1	2	1								
+1	-1	7								
+1	-1	11								
+0	3		QVres							
+1	91	305								
+1	2	1								
+1	-1	7								
+1	-1	11								
+0	1		QMref						
+1									
+1									
+1	-1	9999							
+0	1		QVref						
+1									
+1									
+1	-1	9999							
+										
+FORMAT DU FICHIER :									
+ligne 1 : en-tête obligatoire (version)
+ligne 2	: nom du réservoir								
+ligne 3 : barème V-> Qemin (effectif (entre 1 et 10), valeurs de V (hm3) croissantes, valeurs de Qemin (m3/s))
+ligne 4 : barème V-> Qemax (effectif, valeurs de V (hm3) croissantes, valeurs de Qemax (m3/s))
+ligne 5 : barème V-> Qsmin (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmin (m3/s))
+ligne 6 : barème V-> Qsmax (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmax (m3/s))
+ligne 7 : volume Vtot (hm3, valeur max possible de V) 									
+4 groupes de lignes (pour QMres, QVres, QMref et QVref respectivement) composé chacun de :									
+	ligne 1 : Type (entre 0 et 2), Nombre N de dates pivot (entre 1 et 10) et Borne (si Type <>0). Le débit vaut alors Q défini dans les lignes suivantes si Type=0, max(Q,QV(t-1)-borne) si Type=1 et min(Q,QV(t-1)+borne) si Type = 2								
+	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, compris entre 1 et 365)								
+	ligne 3 : N valeurs entières comprises entre 1 et M(<=N), correspondant aux catégories de valeurs de Q								
+	lignes 4 à 3+M : nombre X (entier entre 1 et 5), puis X couples de valeurs Ai (décroissant),Bi interprétés comme suit : si QMnat>A1 alors Q=B1, sinon, si QMnat>A2 alors Q=B2, sinon…,si QMnat>Ax alors Q=Bx, sinon Q=QMnat								
+NOTATIONS :
+	V : volume d'eau exploitable présent dans le réservoir	
+	Qemin et Qemax : débitance min et max de la prise (limites physiques)								
+	Qsmin et Qsmax : débitance min et max de la restitution	(limites physiques)							
+	QMres et QMref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la prise								
+	QVres et QVref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la restitution								
+	Vtot : volume exploitable total du réservoir, compris entre le sommet de la tranche morte (ou le sommet du volume de garde) et le niveau maximal de sécurité								
+	QMnat : débit du cours d'eau à l'amont immédiat de la prise								
+	QV : débit du cours d'eau à l'aval immédiat de la restitution										
+									

example/PARAMETR/BARRAGES/2/SEINE.1

+Contraintes et consignes du réservoir Seine (règles actuelles d'après F.Dehay, ENGEES 2012). Débit réservé considéré (de façon identique) à la prise et à la restitution. Débit de référence non considéré. 
+SEINE									
+1	0	0			Qemin	(m3/s) fonction de V (hm3)							
+1	0	2.00E+02		Qemax	(m3/s) fonction de V (hm3)							
+1	0	0.8			Qsmin	(m3/s) fonction de V (hm3)							
+1	0	35			Qsmax	(m3/s) fonction de V (hm3)							
+212.9			Vtot	(million de m3)							
+0	4		QMres						
+1	91	274	335						
+1	2	1	2						
+1	-1	3							
+3	10	3	8	4	6	5			
+0	4		QVres						
+1	91	274	335						
+1	2	1	2						
+1	-1	3							
+3	10	3	8	4	6	5			
+0	1		QMref						
+1									
+1									
+1	-1	9999							
+0	1		QVref						
+1									
+1									
+1	-1	9999							
+									
+FORMAT DU FICHIER :									
+ligne 1 : en-tête obligatoire (version)
+ligne 2	: nom du réservoir								
+ligne 3 : barème V-> Qemin (effectif (entre 1 et 10), valeurs de V (hm3) croissantes, valeurs de Qemin (m3/s))
+ligne 4 : barème V-> Qemax (effectif, valeurs de V (hm3) croissantes, valeurs de Qemax (m3/s))
+ligne 5 : barème V-> Qsmin (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmin (m3/s))
+ligne 6 : barème V-> Qsmax (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmax (m3/s))
+ligne 7 : volume Vtot (hm3, valeur max possible de V) 									
+4 groupes de lignes (pour QMres, QVres, QMref et QVref respectivement) composé chacun de :									
+	ligne 1 : Type (entre 0 et 2), Nombre N de dates pivot (entre 1 et 10) et Borne (si Type <>0). Le débit vaut alors Q défini dans les lignes suivantes si Type=0, max(Q,QV(t-1)-borne) si Type=1 et min(Q,QV(t-1)+borne) si Type = 2								
+	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, compris entre 1 et 365)								
+	ligne 3 : N valeurs entières comprises entre 1 et M(<=N), correspondant aux catégories de valeurs de Q								
+	lignes 4 à 3+M : nombre X (entier entre 1 et 5), puis X couples de valeurs Ai (décroissant),Bi interprétés comme suit : si QMnat>A1 alors Q=B1, sinon, si QMnat>A2 alors Q=B2, sinon…,si QMnat>Ax alors Q=Bx, sinon Q=QMnat								
+NOTATIONS :
+	V : volume d'eau exploitable présent dans le réservoir	
+	Qemin et Qemax : débitance min et max de la prise (limites physiques)								
+	Qsmin et Qsmax : débitance min et max de la restitution	(limites physiques)							
+	QMres et QMref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la prise								
+	QVres et QVref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la restitution								
+	Vtot : volume exploitable total du réservoir, compris entre le sommet de la tranche morte (ou le sommet du volume de garde) et le niveau maximal de sécurité								
+	QMnat : débit du cours d'eau à l'amont immédiat de la prise								
+	QV : débit du cours d'eau à l'aval immédiat de la restitution										
+									

example/PARAMETR/BARRAGES/2/YONNE.1

+Contraintes et consignes du réservoir Yonne (règles actuelles d'après F.Dehay, ENGEES 2012). Débit réservé considéré à la restitution et nul à la prise. Débit de référence non considéré à la restitution et nul à la prise. Qsmax=Qsmax1=16 réglementaire  
+YONNE									
+1	0	1.00E+20		Qemin	(m3/s) fonction de V (hm3)							
+1	0	1.00E+20		Qemax	(m3/s) fonction de V (hm3)							
+1	0	0			Qsmin	(m3/s) fonction de V (hm3)							
+1	0	16			Qsmax	(m3/s) fonction de V (hm3)							
+74			Vtot	(million de m3)							
+0	1		QMres						
+1									
+1									
+1	-1	0							
+1	10	2	QVres						
+1	60	91	121	152	182	244	274	305	335
+1	2	3	4	5	6	5	3	2	1
+1	-1	0.6							
+1	-1	1.2							
+1	-1	1.7							
+1	-1	1.8							
+1	-1	1.9							
+1	-1	2							
+0	1		QMref						
+1									
+1									
+1	-1	0							
+2	1	2	QVref						
+1									
+1									
+1	-1	9999						
+										
+FORMAT DU FICHIER :									
+ligne 1 : en-tête obligatoire (version)
+ligne 2	: nom du réservoir								
+ligne 3 : barème V-> Qemin (effectif (entre 1 et 10), valeurs de V (hm3) croissantes, valeurs de Qemin (m3/s))
+ligne 4 : barème V-> Qemax (effectif, valeurs de V (hm3) croissantes, valeurs de Qemax (m3/s))
+ligne 5 : barème V-> Qsmin (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmin (m3/s))
+ligne 6 : barème V-> Qsmax (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmax (m3/s))
+ligne 7 : volume Vtot (hm3, valeur max possible de V) 									
+4 groupes de lignes (pour QMres, QVres, QMref et QVref respectivement) composé chacun de :									
+	ligne 1 : Type (entre 0 et 2), Nombre N de dates pivot (entre 1 et 10) et Borne (si Type <>0). Le débit vaut alors Q défini dans les lignes suivantes si Type=0, max(Q,QV(t-1)-borne) si Type=1 et min(Q,QV(t-1)+borne) si Type = 2								
+	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, compris entre 1 et 365)								
+	ligne 3 : N valeurs entières comprises entre 1 et M(<=N), correspondant aux catégories de valeurs de Q								
+	lignes 4 à 3+M : nombre X (entier entre 1 et 5), puis X couples de valeurs Ai (décroissant),Bi interprétés comme suit : si QMnat>A1 alors Q=B1, sinon, si QMnat>A2 alors Q=B2, sinon…,si QMnat>Ax alors Q=Bx, sinon Q=QMnat								
+NOTATIONS :
+	V : volume d'eau exploitable présent dans le réservoir	
+	Qemin et Qemax : débitance min et max de la prise (limites physiques)								
+	Qsmin et Qsmax : débitance min et max de la restitution	(limites physiques)							
+	QMres et QMref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la prise								
+	QVres et QVref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la restitution								
+	Vtot : volume exploitable total du réservoir, compris entre le sommet de la tranche morte (ou le sommet du volume de garde) et le niveau maximal de sécurité								
+	QMnat : débit du cours d'eau à l'amont immédiat de la prise								
+	QV : débit du cours d'eau à l'aval immédiat de la restitution										
+									

example/PARAMETR/BARRAGES/2/version.txt

+2012, Qres prise et restit, pas de Qref, QSmax(Yonne)=QSmax1
\ No newline at end of file

example/PARAMETR/BARRAGES/3/AUBE.1

+contraintes et consignes du réservoir Aube (règles actuelles d'après F.Dehay, ENGEES 2012). Débit réservé considéré seulement à la prise et débit de référence considéré de façon identique à la prise et à la restitution. 
+AUBE									
+1	0	0			Qemin	(m3/s) en fonction de V (hm3)							
+1	0	1.35E+02		Qemax	(m3/s) en fonction de V (hm3)							
+1	0	0			Qsmin	(m3/s) en fonction de V (hm3)							
+1	0	35			Qsmax	(m3/s) en fonction de V (hm3)							
+181.2			Vtot	(million de m3)							
+0	3		QMres						
+1	91	305							
+1	2	1							
+1	-1	2							
+3	8	2	6	3	4	4			
+0	1		QVres						
+1							
+1							
+1	-1	0							
+0	1		QMref						
+1									
+1									
+1	-1	130							
+0	1		QVref						
+1									
+1									
+1	-1	130							
+									
+FORMAT DU FICHIER :									
+ligne 1 : en-tête obligatoire (version)
+ligne 2	: nom du réservoir								
+ligne 3 : barème V-> Qemin (effectif (entre 1 et 10), valeurs de V (hm3) croissantes, valeurs de Qemin (m3/s))
+ligne 4 : barème V-> Qemax (effectif, valeurs de V (hm3) croissantes, valeurs de Qemax (m3/s))
+ligne 5 : barème V-> Qsmin (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmin (m3/s))
+ligne 6 : barème V-> Qsmax (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmax (m3/s))
+ligne 7 : volume Vtot (hm3, valeur max possible de V) 									
+4 groupes de lignes (pour QMres, QVres, QMref et QVref respectivement) composé chacun de :									
+	ligne 1 : Type (entre 0 et 2), Nombre N de dates pivot (entre 1 et 10) et Borne (si Type <>0). Le débit vaut alors Q défini dans les lignes suivantes si Type=0, max(Q,QV(t-1)-borne) si Type=1 et min(Q,QV(t-1)+borne) si Type = 2								
+	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, compris entre 1 et 365)								
+	ligne 3 : N valeurs entières comprises entre 1 et M(<=N), correspondant aux catégories de valeurs de Q								
+	lignes 4 à 3+M : nombre X (entier entre 1 et 5), puis X couples de valeurs Ai (décroissant),Bi interprétés comme suit : si QMnat>A1 alors Q=B1, sinon, si QMnat>A2 alors Q=B2, sinon…,si QMnat>Ax alors Q=Bx, sinon Q=QMnat								
+NOTATIONS :
+	V : volume d'eau exploitable présent dans le réservoir	
+	Qemin et Qemax : débitance min et max de la prise (limites physiques)								
+	Qsmin et Qsmax : débitance min et max de la restitution	(limites physiques)							
+	QMres et QMref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la prise								
+	QVres et QVref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la restitution								
+	Vtot : volume exploitable total du réservoir, compris entre le sommet de la tranche morte (ou le sommet du volume de garde) et le niveau maximal de sécurité								
+	QMnat : débit du cours d'eau à l'amont immédiat de la prise								
+	QV : débit du cours d'eau à l'aval immédiat de la restitution										
+									

example/PARAMETR/BARRAGES/3/MARNE.1

+contraintes et consignes du réservoir Marne (règles actuelles d'après F.Dehay, ENGEES 2012). Débit réservé considéré seulement à la prise et débit de référence considéré à la prise et à la restitution. 
+MARNE										
+1	0	0			Qemin	(m3/s) en fonction de V (hm3)								
+1	0	4.08E+02		Qemax	(m3/s) en fonction de V (hm3)								
+1	0	0			Qsmin	(m3/s) en fonction de V (hm3)								
+1	0	50			Qsmax	(m3/s) en fonction de V (hm3)								
+354.5			Vtot	(million de m3)								
+0	3		QMres							
+1	91	305								
+1	2	1								
+1	-1	7								
+1	-1	11								
+0	1		QVres							
+1								
+1								
+1	-1	0																
+0	6		QMref							
+1	91	121	182	274	305					
+1	2	3	4	3	1					
+1	-1	160								
+1	-1	140								
+1	-1	120								
+1	-1	90								
+0	5		QVref							
+1	182	274	289	305						
+1	2	3	4	1						
+1	-1	140								
+1	-1	75								
+1	-1	100								
+1	-1	120								
+										
+FORMAT DU FICHIER :									
+ligne 1 : en-tête obligatoire (version)
+ligne 2	: nom du réservoir								
+ligne 3 : barème V-> Qemin (effectif (entre 1 et 10), valeurs de V (hm3) croissantes, valeurs de Qemin (m3/s))
+ligne 4 : barème V-> Qemax (effectif, valeurs de V (hm3) croissantes, valeurs de Qemax (m3/s))
+ligne 5 : barème V-> Qsmin (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmin (m3/s))
+ligne 6 : barème V-> Qsmax (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmax (m3/s))
+ligne 7 : volume Vtot (hm3, valeur max possible de V) 									
+4 groupes de lignes (pour QMres, QVres, QMref et QVref respectivement) composé chacun de :									
+	ligne 1 : Type (entre 0 et 2), Nombre N de dates pivot (entre 1 et 10) et Borne (si Type <>0). Le débit vaut alors Q défini dans les lignes suivantes si Type=0, max(Q,QV(t-1)-borne) si Type=1 et min(Q,QV(t-1)+borne) si Type = 2								
+	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, compris entre 1 et 365)								
+	ligne 3 : N valeurs entières comprises entre 1 et M(<=N), correspondant aux catégories de valeurs de Q								
+	lignes 4 à 3+M : nombre X (entier entre 1 et 5), puis X couples de valeurs Ai (décroissant),Bi interprétés comme suit : si QMnat>A1 alors Q=B1, sinon, si QMnat>A2 alors Q=B2, sinon…,si QMnat>Ax alors Q=Bx, sinon Q=QMnat								
+NOTATIONS :
+	V : volume d'eau exploitable présent dans le réservoir	
+	Qemin et Qemax : débitance min et max de la prise (limites physiques)								
+	Qsmin et Qsmax : débitance min et max de la restitution	(limites physiques)							
+	QMres et QMref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la prise								
+	QVres et QVref : débit réservé (valeur qu'on souhaite voir dépassée) et débit de référence (valeur qu'on souhaite ne pas voir dépassée) dans le cours d'eau à l'aval immédiat de la restitution								
+	Vtot : volume exploitable total du réservoir, compris entre le sommet de la tranche morte (ou le sommet du volume de garde) et le niveau maximal de sécurité								
+	QMnat : débit du cours d'eau à l'amont immédiat de la prise								
+	QV : débit du cours d'eau à l'aval immédiat de la restitution										
+									

example/PARAMETR/BARRAGES/3/SEINE.1

+contraintes et consignes du réservoir Seine (règles actuelles d'après F.Dehay, ENGEES 2012). Débit réservé considéré seulement à la prise et débit de référence considéré de façon identique à la prise et à la restitution. 
+SEINE									
+1	0	0			Qemin	(m3/s) en fonction de V (hm3)							
+1	0	2.00E+02		Qemax	(m3/s) en fonction de V (hm3)							
+1	0	0.8			Qsmin	(m3/s) en fonction de V (hm3)							
+1	0	35			Qsmax	(m3/s) en fonction de V (hm3)							
+212.9			Vtot	(million de m3)							
+0	4		QMres						
+1	91	274	335						
+1	2	1	2						
+1	-1	3							
+3	10	3	8	4	6	5			
+0	1		QVres						
+1						
+1						
+1	-1	0										
+0	4		QMref						
+1	91	182	305						
+1	2	3	1						
+1	-1	120							
+1	-1	90							
+1	-1	40							
+0	4		QVref						
+1	91	182	305						
+1	2	3	1						
+1	-1	120							
+1	-1	90							
+1	-1	40							
+									
+FORMAT DU FICHIER :									
+ligne 1 : en-tête obligatoire (version)
+ligne 2	: nom du réservoir								
+ligne 3 : barème V-> Qemin (effectif (entre 1 et 10), valeurs de V (hm3) croissantes, valeurs de Qemin (m3/s))
+ligne 4 : barème V-> Qemax (effectif, valeurs de V (hm3) croissantes, valeurs de Qemax (m3/s))
+ligne 5 : barème V-> Qsmin (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmin (m3/s))
+ligne 6 : barème V-> Qsmax (effectif, valeurs de V (hm3) croissantes, valeurs de Qsmax (m3/s))
+ligne 7 : volume Vtot (hm3, valeur max possible de V) 									
+4 groupes de lignes (pour QMres, QVres, QMref et QVref respectivement) composé chacun de :									
+	ligne 1 : Type (entre 0 et 2), Nombre N de dates pivot (entre 1 et 10) et Borne (si Type <>0). Le débit vaut alors Q défini dans les lignes suivantes si Type=0, max(Q,QV(t-1)-borne) si Type=1 et min(Q,QV(t-1)+borne) si Type = 2								
+	ligne 2 : N dates pivot entières croissantes (quantièmes de changement de valeurs de Q, compris entre 1 et 365)