VGEST.PAS 49.6 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
PROGRAM VGEST;

{--JC Bader, nov 2016--}

uses Utilit,Interfas,DECLARA,ECRAN,PARAM,CALCULS1,CALCULS2,CALCULS3,CALCULS4;

{--------------------------------------------------------------------------}
{-Effacement des fichiers texte de rsultats prcdents et incrmentation--}
{-du rang de calcul.                                                     --}
{--------------------------------------------------------------------------}

procedure Ici_Effacement;
 begin
  assign(YFitext,YRepR+YNomFicResult6); rewrite(YFitext); close(YFitext);
  assign(YFitext,YRepR+YNomFicSortie); rewrite(YFitext); close(YFitext);
  assign(YFitext,YRepR+YNomFicResult1); rewrite(YFitext); close(YFitext);
  assign(YFitext,YrepR+YNomFicResult2); rewrite(YFitext); close(YFitext);
  assign(YFitext,YRepR+YNomFicResult4); rewrite(YFitext); close(YFitext);
  assign(YFitext,YRepR+YNomFicResult5[1]); rewrite(YFitext); close(YFitext);
  assign(YFitext,YRepR+YNomFicResult5[2]); rewrite(YFitext); close(YFitext);
  assign(YFitext,YRepR+YNomFicResChr1); rewrite(YFitext); close(YFitext);
  assign(YFitext,YRepR+YNomFicResChr2); rewrite(YFitext); close(YFitext);
  assign(YFitext,YRepR+YNomFicQuant); rewrite(YFitext); close(YFitext);
  for YRangLac:=1 to YNmaxLacs do
   begin
    str(YRanglac,YNum);
    assign(YFitext,YRepR+YNomFicResultVobj+YNum+YExIsofreq); rewrite(YFitext); close(YFitext);
    assign(YFitext,YRepR+YNomFicResultVop+YNum+YExIsofreq); rewrite(YFitext); close(YFitext);
   end;
  assign(YFitext,YRepR+YNomFicResDef+'0B'+YExIsofreq); rewrite(YFitext); close(YFitext);
  assign(YFitext,YRepR+YNomFicResDef+'1B'+YExIsofreq); rewrite(YFitext); close(YFitext);
  assign(YFitext,YRepR+YNomFicResultM+'A'+YExIsofreq); rewrite(YFitext); close(YFitext);
  assign(YFitext,YRepR+YNomFicResultM+'B'+YExIsofreq); rewrite(YFitext); close(YFitext);
  assign(YFitext,YNomFicRangCal);
  {$i-}   {dsactive le contrle des entres-sorties}
  reset(YFitext);
  if ioresult<>0 then
   begin
    rewrite(YFitext);
    YCompteur:=1;
    writeln(YFitext,YCompteur);
   end
  else
   begin
    readln(YFitext,YCompteur);
    if ioresult=0 then
     YCompteur:=YCompteur+1
    else
     YCompteur:=1;
    rewrite(YFitext);
    writeln(YFitext,YCompteur);
   end;
  close(YFitext);
  {$i+}   {active le contrle des entres-sorties}
 end;

 {------------------------------------------------------------------------}
 {--OUVERTURE DU FICHIER DES SORTIES dcrivant les paramtres du calcul --}
 {--et donnant certaines statistiques sur les rsultats obtenus.        --}
 {--Premires critures sur ce fichier.                                 --}
 {------------------------------------------------------------------------}

 Procedure Ici_Ecriture_EnTete_CondStat;
  begin
   assign(YFicSortie,YRepR+YNomFicSortie);
   rewrite(YFicSortie);
   writeln(YFicSortie,'Rang du calcul : ',YCompteur);
   B_Datation(YFicSortie);
   writeln(YFicSortie);
   writeln(YFicSortie,'Ce fichier contient les informations suivantes :');
   writeln(YFicSortie);
   writeln(YFicSortie,'RESUME DES CONDITIONS DU CALCUL : ');
   writeln(YFicSortie,'# noms des fichiers de parametres et de donnees choisis');
   writeln(YFicSortie,'# caracteristiques des differents reservoirs (contraintes et consignes de gestion, ',
                      'capacite, temps de propagation jusqu''a la station aval du systeme)');
   writeln(YFicSortie,'# mode de repartition entre les reservoirs du debit a stocker pour satisfaire l''objectif aval');
   writeln(YFicSortie,'# stations, superficies de bassin versant et temps de propagation utilises pour le ',
                      'calcul des debits naturels aux prises et restitutions');
   writeln(YFicSortie,'# hydrogramme objectif annuel : nature (minimal ou maximal), valeurs de debit et ',
                      'localisation a la station aval du systeme');
   writeln(yFicSortie);
   writeln(YFicSortie);
   writeln(YFicSortie,'STATISTIQUES DE RESULTATS concernant la repartition des debits a stocker : ');
   writeln(YFicSortie,'# possibilites de satisfaction de l''objectif');
   writeln(YFicSortie,'# valeurs optimales de debits stockes, preleves et restitues');
   writeln(YFicSortie);
   writeln(YFicSortie);
   writeln(yFicSortie,'===================================');
   writeln(YFicSortie,'! RESUME DES CONDITIONS DU CALCUL !');
   writeln(yFicSortie,'===================================');
   writeln(YFicSortie);
  end;
93

94
95
96
97
98
99
 {--------------------------------------------------------------------------}
 {--Si saisie clavier, CHOIX DU TYPE D'OBJECTIF : DEBIT MIN OU MAX A      --}
 {--RESPECTER.                                                            --}
 {--Description de l'hyrogramme objectif et de son type dans le fichier   --}
 {--des paramtres et options du calcul.                                  --}
 {--------------------------------------------------------------------------}
100

101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
 Procedure Ici_ChoixTypeObjectif;
  begin
   if YSaisie then
    E_Saisie2;
   write(YFicSortie,'Hydrogramme objectif (debit QXobj) que l''on souhaite ');
   if YTypeObjectif=0 then
    begin
     writeln(YFicSortie,'VOIR TOUJOURS DEPASSE a la station aval du systeme : ');
     YResume:=YResume+' (Qmin) ; ';
    end
   else
    begin
     writeln(YFicSortie,'NE JAMAIS VOIR DEPASSE a la station aval du systeme : ');
     YResume:=Yresume+' (Qmax) ; ';
    end;
   writeln(YFicSortie);
   for Yi:=1 to trunc((YNbQXobj-1)/4)+1 do
    begin
     write(YFicSortie,'date   : ');
     for Yj:=(Yi-1)*4+1 to B_minin(4*Yi,YNbQXObj-1) do
121
      write(YFicSortie,'              ',B_DateClaire(YDateObj[Yj]));
122
123
124
125
126
127
128
129
130
131
     writeln(YFicSortie);
     write(YFicSortie,'debit  : ');
     for Yj:=(Yi-1)*4+1 to B_minin(4*Yi,YNbQXObj-1) do
      write(YFicSortie,'     ',YDebitObj[Yj]:14);
     writeln(YFicSortie);
     writeln(YFicSortie);
    end;
   writeln(YFicSortie);
   flush(YFicSortie)
  end;
132

133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
 {-------------------------------------------------------------------}
 {--OUVERTURE DU FICHIER TEXTE DE DEBITS NATURELS SYNCHRONISES, et --}
 {--criture d'en-tte sur ce fichier                              --}
 {-------------------------------------------------------------------}

 Procedure Ici_EcritureEntete_QSynchro;
  begin
   assign(YFicResult2,YRepR+YNomFicResult2);
   rewrite(YFicResult2);
   writeln(YFicResult2,'Rang du calcul : ',YCompteur,' ;  ',YResume);
   B_Datation(YFicResult2);
   writeln(YFicResult2);
   writeln(YFicResult2,'Debits naturels synchronises a la station aval et aux points de prise et de restitution des reservoirs');
   writeln(YFicResult2,'Les valeurs donnees en haut des colonnes de dates sont les parties ',
                       'fractionnaires de celles-ci (exprimees en jour)');
   writeln(YFicResult2);
   writeln(YFicResult2,'Variables editees :');
   writeln(YFicResult2);
   writeln(YFicResult2,'# date aval : date a la station situee a l''aval du systeme (lieu de l''objectif)');
   writeln(YFicResult2,'# code : 1 si debits sans lacune et 0 sinon');
   writeln(YFicResult2,'# QXnat : debit naturel a la station aval, a 12h (date aval)');
   writeln(YFicResult2,'# QXobj : objectif de debit defini a la station aval, a 12h (date aval)');
   writeln(YFicResult2,'# date amont : date aval diminuee du temps de propagation entre le point de ',
                       'restitution et la station aval du systeme');
   writeln(YFicResult2,'# QMnat : debit naturel au point de prise (date amont)');
   writeln(YFicResult2,'# QVnat : debit naturel au point de restitution (date amont)');
   writeln(YFicResult2);
   write(YFicResult2,'   0.50000                                                            ');
   for YRanglac:=1 to YNbLacs do
    write(YFicResult2,'  ',YDecalDateQ[YRanglac]:8:5,
                      '                                                ');
   writeln(YFicResult2);
   write(YFicResult2,' date aval   code                       QXnat                QXobj   ');
   for YRangLac:=1 to YNbLacs do
    write(YFicResult2,'date amont',YRangLac,
                      '                QMnat',YRanglac,
                      '                QVnat',YRanglac,'   ');
   writeln(YFicResult2);
171

172
  end;
173
174


175
176
177
178
 {---------------------------------------------------------------------}
 {--criture d'en-tte sur fichier texte de rsultats des chroniques --}
 {--de volume                                                        --}
 {---------------------------------------------------------------------}
179

180
 Procedure Ici_EcritureEntete_Volumes;
181
  begin
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
   assign(YFicResult4,YRepR+YNomFicResult4);
   rewrite(YFicResult4);
   writeln(YFicResult4,'Rang du calcul : ',YCompteur,' ;  ',Yresume);
   B_Datation(YFicResult4);
   writeln(YFicResult4);
   writeln(YFicResult4,'Chroniques de volumes presents dans les reservoirs,',
                       ' les plus compatibles avec les chroniques de debits stockes objectifs');
   writeln(YFicResult4,'Les valeurs situees en haut des colonnes de date representent leur partie fractionnaire');
   writeln(YFicResult4);
   writeln(YFicResult4,'Variables editees :');
   writeln(YFicResult4);
   writeln(YFicResult4,'# date aval : date T a laquelle on vise l''objectif de debit a la station situee a l''aval du systeme');
   writeln(YFicResult4,'# date(i) : date aval diminuee du temps de propagation entre le point de ',
                       'restitution du reservoir de rang i et la station aval du systeme, plus 12 heures (= T-D(i)+0.5)');
   write(YFicResult4,'# Vobj(i) : ');
   if YTypeObjectif=0 then
    writeln(YFicResult4,'plus petit volume d''eau present dans le reservoir i, permettant la meilleure satisfaction posterieure possible de l''objectif aval ',
                        '(un volume superieur ou egal a Vobj permet cette meilleure satisfaction possible dans le futur, mais un volume plus petit implique plus de defaillances futures)')
   else
    writeln(YFicResult4,'plus grand volume d''eau present dans le reservoir i, permettant la meilleure satisfaction posterieure possible de l''objectif aval ',
                        '(un volume inferieur ou egal a Vobj permet cette meilleure satisfaction possible dans le futur, mais un volume plus grand implique plus de defaillances futures)');
   write(YFicResult4,'# Vges(i) : ');
   if YTypeObjectif=0 then
    writeln(YFicResult4,'plus grand volume d''eau pouvant etre present dans le reservoir i, apres la meilleure satisfaction anterieure possible de l''objectif aval ',
                        '(un volume superieur a Vges, plus favorable a la satisfaction future, implique plus de defaillances anterieures)')
   else
    writeln(YFicResult4,'plus petit volume d''eau pouvant etre present dans le reservoir i, apres la meilleure satisfaction anterieure possible de l''objectif aval ',
                        '(un volume inferieur a Vges, plus favorable a la satisfaction future, implique plus de defaillances anterieures)');
   write(YFicResult4,'# Vop(i) : ');
   if YTypeObjectif=0 then
    writeln(YFicResult4,'plus petit volume d''eau present dans le reservoir i, compatible avec la meilleure satisfaction possible de l''objectif aval,',
                        ' passee (premiere priorite) et future (seconde priorite) : Vop = min(Vobj ; Vges)')
   else
    writeln(YFicResult4,'plus grand volume d''eau present dans le reservoir i, compatible avec la meilleure satisfaction possible de l''objectif aval,',
                        ' passee (premiere priorite) et future (seconde priorite) : Vop = max(Vobj ; Vges)');
   writeln(YFicResult4);
   write(YFicResult4,'   0.50000 ');
   for YRangLac:=1 to YNbLacs do
    write(YFicResult4,' ',YDecalDateQ[YRangLac]:11:8,
                      '                                   ',
                      '                                    ');
   writeln(YFicResult4);
   write(YFicResult4,' date aval   ');
   if YAbsRel=1 then
    for YRangLac:=1 to YNbLacs do
     write(YFicResult4,'     date',YRanglac,'       Vobj',YRangLac,'(m3)',
                         '   Vobj',YRangLac,'(m3) a 24H   ',
                         '      Vges',YRangLac,'(m3)   ','       Vop',YRangLac,'(m3)   ')
   else
    for YRangLac:=1 to YNbLacs do
     write(YFicResult4,'     date',YRanglac,'     Vobj',YRangLac,'/Vtot',YRangLac,
                         '  Vobj',YRangLac,'/Vtot',YRangLac,'(24H)   ',
                         '    Vges',YRangLac,'/VTot',YRangLac,'   ','     Vop',YRangLac,'/Vtot',YRangLac,'   ');

   writeln(YFicResult4);
237
  end;
238
239
240
241
242
243
244
245
246
247

 {-----------------------------------------------------------------------}
 {Ecriture sur fichier texte de rsultats concernant les dfaillances de }
 {gestion (volumes manquants Vdef dans les rservoirs et dfaut de dbit }
 { la station objectif aval) pour une exploitation visant  optimiser la}
 {rpartition des taches entre rservoirs                                }
 {-----------------------------------------------------------------------}

 procedure Ici_EditionBilan;

248
249
250
  const
    XNul : double = 0;

251
252
253
254
  var
   Xi : integer;
   Xj : integer;
   XQXdef : array[1..2,1..2] of double;
255
   nbEof : integer;
256

257
  begin
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289

   assign(YFicResult4,YRepR+YNomFicResult6);
   rewrite(YFicResult4);
   writeln(YFicResult4,'Rang du calcul : ',YCompteur,'   ',Yresume);
   B_Datation(YFicResult4);
   writeln(YFicResult4);
   writeln(YFicResult4,YNbLacs,'  (nombre de reservoirs)');
   writeln(YFicResult4);
   writeln(YFicResult4,YOptionRepart,'  (cle de repartition des stockages entre reservoirs : 1 = figee ; autre = variable)');
   writeln(YFicResult4);
   writeln(YficResult4,YtypeObjectif,'  (type d''objectif : 0 = maintien du debit ',
                                     'au dessus de Qmin ; autre = maintien du debit ',
                                     'en dessous de Qmax)');
   writeln(YficResult4);
   for YRanglac:=1 to YNbLacs do
    write(YFicResult4,YVtot[YRangLac],'  ');
   writeln(YFicResult4,YVtotCumul,'  (volumes exploitables reels Vtot et leur somme (m3))');
   if YOptionRepart=1 then
    begin
     for YRanglac:=1 to YNbLacs do
      write(YFicResult4,YVRtot[YRangLac],'  ');
     writeln(YFicResult4,YVRtotCumul,'  (volumes exploitables fictifs VRtot et leur somme)');
     for YRanglac:=1 to YNbLacs do
      write(YFicResult4,YPartTheorique[YRangLac],'  ');
     writeln(YFicResult4, '  (repartition theorique des taches)');
     for YRanglac:=1 to YNbLacs do
      write(YFicResult4,YPart[YRangLac],'  ');
     writeln(YFicResult4, '  (repartition des taches apres reglage)');
     for YRanglac:=1 to YNblacs do
      write(YFicResult4,YReglage[Yranglac],'  ');
     writeln(YFicResult4,'  (coefficients de reglage)');
    end
290
   else
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
    begin
     writeln(yFicResult4);
     writeln(YFicResult4);
     writeln(yFicResult4);
     writeln(yFicResult4);
    end;
   writeln(YFicResult4);
   writeln(yFicResult4,'---------------------------------------------------------------------------');
   writeln(yFicResult4);
   writeln(YFicResult4,'Valeurs relatives moyennes du volume d''eau present dans chaque reservoir sur la periode de calcul ');
   for YRangLac:=1 to YNblacs do
    write(YFicResult4,YVobjMoyRel[1,YRangLac],' ');
   writeln(YFicResult4,'  (sens chronologique inverse : Vobj(24H)/Vtot)');
   for YRangLac:=1 to YNblacs do
    write(YFicResult4,YVobjMoyRel[2,YRangLac],' ');
   writeln(YFicResult4,'  (sens chronologique : Vges/Vtot)');
   writeln(YFicResult4);
    writeln(yFicResult4,'---------------------------------------------------------------------------');
   writeln(yFicResult4);
   writeln(YFicResult4,'taux moyens quadratiques de sollicitation des reservoirs sur la periode de calcul');
   for YRangLac:=1 to YNblacs do
    write(YFicResult4,YSollicitation[1,YRangLac],' ');
   write(YFicResult4,'  (sens chronologique inverse : ');
   if YTypeObjectif=0 then
    writeln(YFicResult4,'Vobj/Vtot)')
   else
    writeln(YFicResult4,'(Vtot-Vobj)/Vtot)');
   for YRangLac:=1 to YNblacs do
    write(YFicResult4,YSollicitation[2,YRangLac],' ');
   write(YFicResult4,'  (sens chronologique : ');
   if YTypeObjectif=1 then
    writeln(YFicResult4,'Vges/Vtot)')
   else
    writeln(YFicResult4,'(Vtot-Vges)/Vtot)');
   for YRangLac:=1 to YNblacs do
    write(YFicResult4,(YSollicitation[2,YRangLac]+YSollicitation[1,YRangLac])/2,' ');
   writeln(YFicResult4,'  (moyenne des taux quadratiques moyens des deux sens de calcul)');
   writeln(YFicResult4);
   writeln(yFicResult4,'---------------------------------------------------------------------------');
   writeln(yFicResult4);
   writeln(YFicResult4,'Valeurs maximales sur la periode de calcul, du volume manquant global Vdef rapporte a la capacite de l''ensemble des reservoirs : ');
   writeln(YFicResult4,YVManqueMax[1]/YVTotCumul:26:11,'  (sens chronologique inverse)');
   writeln(YFicResult4,YVManqueMax[2]/YVtotCumul:26:11,'  (sens chronologique)');
   writeln(YFicResult4);
   writeln(yFicResult4,'---------------------------------------------------------------------------');
   writeln(yFicResult4);
   write(YFicResult4,'Defaillances correspondant a un ');
   if YTypeObjectif=0 then
    write(YFicResult4,'manque')
   else
    write(YFicResult4,'exces');
   writeln(YFicResult4,' de debit (m3/s) par rapport a l''objectif defini a la station aval du systeme :');
   writeln(YFicResult4,'  QXdef0 : defaillances inevitables absolues dues aux contraintes physiques et consignes aux prises et restitutions des reservoirs');
   writeln(YFicResult4,'  QXdef1 : defaillances dues a la fois aux contraintes physiques et consignes aux prises et restitutions, aux volumes d''eau V presents et aux capacites Vtot des reservoirs');
   writeln(YFicResult4,'  QXdef2 : defaillance s''ajoutant a QXdef0, du fait des volumes d''eau V presents et des capacites Vtot des reservoirs (QXdef2=QXdef1-QXdef0)');
   writeln(YFicResult4);
   write(YFicResult4,'Effectifs des defaillances selon differentes valeurs du code Succes (cf fichier ',YNomficSortie,'). Le code ');
   if YTypeObjectif=0 then
    write(YFicResult4,'-1')
   else
    write(YFicResult4,'1');
   writeln(YFicResult4,' correspond a des defaillances qui auraient pu etre evitees avec une repartition differente du volume global present');
   if YTypeObjectif=0 then
    begin
     for Xi:=-1 downto -3 do
      write(YFicResult4,Xi:9,' ');
     writeln(YFicResult4,'    (code Succes)');
     for Xi:=-1 downto -3 do
      write(YFicResult4,YNSuccesA[Xi]:9,' ');
     writeln(YFicResult4,'    (sens chronologique inverse)');
     for Xi:=-1 downto -3 do
      write(YFicResult4,YNSuccesB[Xi]:9,' ');
     writeln(YFicResult4,'    (sens chronologique)');
    end
   else
    begin
     for Xi:=1 to 4 do
      write(YFicResult4,Xi:9,' ');
     writeln(YFicResult4,'    (code Succes)');
     for Xi:=1 to 4 do
      write(YFicResult4,YNSuccesA[Xi]:9,' ');
     writeln(YFicResult4,'    (sens chronologique inverse)');
     for Xi:=1 to 4 do
      write(YFicResult4,YNSuccesB[Xi]:9,' ');
     writeln(YFicResult4,'    (sens chronologique)');
    end;
   writeln(yFicResult4);
   writeln(YFicResult4,'Valeur moyenne de fonction puissance ((moyenne(QXdef^N))^(1/N) avec N = ',YPuiss,') : ');
   writeln(YFicResult4,'                QXdef2                 QXdef1                 QXdef0');
   writeln(YFicResult4,YQXdefmoyGlop[2,1]:22,' ',YQXdefmoyGlop[1,1]:22,' ',YQXdefmoyGlop[0,1]:22,' (sens chronologique inverse)');
   writeln(YFicResult4,YQXdefmoyGlop[2,2]:22,' ',YQXdefmoyGlop[1,2]:22,' ',YQXdefmoyGlop[0,2]:22,' (sens chronologique)');
   writeln(YFicResult4,((YQXdefmoyGlop[2,2]+YQXdefmoyGlop[2,1])/2):22,' ',((YQXdefmoyGlop[1,2]+YQXdefmoyGlop[1,1])/2):22,' ',
                       ((YQXdefmoyGlop[0,2]+YQXdefmoyglop[0,1])/2):22,' (moyenne des valeurs obtenues dans les deux sens de calcul)');

   {--cas o on dite les valeurs tries de dfaillances de dbit--}
   if YTestCalcTri then
    begin

     for Xi:=1 to 2 do  {Xi : 1 pour QXdef1 et 2 pour QXdef2}
      for Xj:=1 to 2 do  {Xj : 1 pour sens pas chrono et 2 pour sens chrono}
       begin
        assign(YFicQXdef [Xi,Xj],YRepT+YNomFicQXdef [Xj]+YExt[Xi]);
        assign(YFicQXdefT[Xi,Xj],YRepT+YNomFicQXdefT[Xj]+YExt[Xi]);
394
        C4_TriChronique(YFicQXdef[Xi,Xj],YFicQXdefT[Xi,Xj],YTouComNul,YQXdefmoy[Xi,Xj],YQXdefmed[Xi,Xj],YQXdefmax[Xi,Xj],YQXdefmin[Xi,Xj],YQXdefrare[Xi,Xj],YQXdeffrequent[Xi,Xj]);
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
       end;
     writeln(YFicResult4);
     writeln(YFicResult4,'Fonctions de repartition (frequence = (n-YA)/(Effectif+YB) avec n=rang, YA=',YA,' et YB=',YB,') :');
     for Xi:= 1 to 2 do
      begin
       writeln(YFicResult4,'  Repartition de QXdef',Xi);
       write(YFicResult4,'    Frequence-->               Moyenne    Maximum ');
       for Yi:=1 to YNbFreq do
        begin
         YFreq:=1/YTpsFreq[Yi];
         write(YFicResult4,'   ',Yfreq:7:5,' ');
        end;
       Yfreq:=0.5;
       write(YFicResult4,'   ',YFreq:7:5,' ');
       for Yi:=1 to YNbFreq do
        begin
         YFreq:=1-1/YTpsFreq[YNbFreq+1-Yi];
         write(YFicResult4,'   ',Yfreq:7:5,' ');
        end;
       writeln(YFicResult4,'   Minimum');
       for Xj:=1 to 2 do
        begin
         write(YFicResult4,'                ',YQXdefmoy[Xi,Xj]:22,' ',YQXdefmax[Xi,Xj]:10:3,' ');
         for Yi:=1 to YNbfreq do
          write(YFicResult4,YQXdefrare[Xi,Xj,Yi]:10:3,' ');
         write(YFicResult4,YQXdefmed[Xi,Xj]:10:3,' ');
         for Yi:=1 to YNbfreq do
          write(YFicResult4,YQXdeffrequent[Xi,Xj,Yi]:10:3,' ');
         write(YFicResult4,YQXdefmin[Xi,Xj]:10:3);
         if Xj=1 then
          writeln(YFicResult4,'  (sens chronologique inverse)')
         else
          writeln(YFicResult4,'  (sens chronologique)');
        end;
      end;
     writeln(YFicResult4);
     writeln(YFicResult4,'valeurs classees :');
     for Xi:=1 to 2 do
      for Xj:=1 to 2 do
       reset(YFicQXdefT[Xi,Xj]);
     writeln(YFicResult4,'                  QXdef1                  QXdef1                  QXdef2                  QXdef2');
     writeln(YFicResult4,'   (chronologie inverse)    (sens chronologique)   (chronologie inverse)    (sens chronologique)');
437
438
     nbEof:=0;
     while nbEof < 4 do
439
      begin
440
       nbEof:=0;
441
442
443
       for Xi:=1 to 2 do
        for Xj:=1 to 2 do
         begin
444
445
446
447
448
449
450
451
          if not eof(YFicQXdefT[Xi,Xj]) Then
            read(YFicQXdefT[Xi,Xj],XQXdef[Xi,Xj])
          else
            begin
              XQXdef[Xi,Xj]:=XNul;
              nbEof:=nbEof+1;
            end;

452
453
454
455
456
457
458
459
460
461
462
463
          write(YFicResult4,'  ',XQXdef[Xi,Xj]:22);
         end;
       writeln(YFicResult4);
      end;
     for Xi:=1 to 2 do
      for Xj:=1 to 2 do
      close(YFicQXdefT[Xi,Xj]);

    end; {fin du cas o on dite les valeurs tries de dfaillances de dbit}

   close(YficResult4);

464
  end;
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480

{===========================================================================}
{-- PROGRAMME PRINCIPAL                                                   --}
{===========================================================================}

BEGIN

 {--CHOIX DU MODE DE PARAMETRAGE DU CALCUL, PAR SAISIE CLAVIER OU LECTURE --}
 {--DES PARAMETRES SUR FICHIER TEXTE (et lecture des paramtres pour le   --}
 {--second cas).                                                          --}
 P_ChoixEntreeOptions;

 {--EFFACEMENT DE TOUS LES FICHIERS texte de rsultats prcdants,--}
 {--incrmentation du rang du calcul                              --}
 Ici_Effacement;

481
 if Ysaisie then
482
  E_Message1;
483
484
485

 YCalculPossibl:=true;

486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
 {--OUVERTURE DU FICHIER DES SORTIES (YNomFicSortie = 'Condstat') dcrivant--}
 {--les paramtres et options du calcul et donnant certaines statistiques  --}
 {--de rsultats obtenus. Premires critures sur ce fichier.              --}
 Ici_Ecriture_EnTete_Condstat;

 {--Si saisie clavier : saisie du nom de station aval                   --}
 {--Lecture du fichier de paramtres de station aval                    --}
 {--Lecture des rglages (ne sert que pour mode de rpartition fig)    --}
 {--Si saisie clavier : saisie du rang de paramtres lacs               --}
 {--si saisie clavier : saisie du rang de paramtres rseaux            --}
 {--criture sur fichier de sortie et rsum : station aval             --}
 {--criture sur fichier de sortie : rservoirs                         --}
 P_StationAval;

 {--Si saisie clavier : choix de la mthode de rpartition et du type--}
 {--d'dition des volumes                                            --}
 if Ysaisie then
  E_saisie7;

 {--Dans le cas de l'option 5 concernant la rpartition des objectifs --}
 {--entre rservoirs, lecture de la frquence de non dpassement      --}
 {--souhaite pour les valeurs de dbits QMN aux droit des prises    --}
 {--utiliser dans les calculs d'quilibrage des dures potentielles   --}
 {--de remplissage.                                                   --}
 if YOptionRepart=5 then
  C2_LectureFreqNonDep;

 {--LECTURE DES PARAMETRES DES CONTRAINTES ET CONSIGNES, calcul des YVrtot --}
 {--(utiles seulement pour mode de rpartition 1) et criture des rgles de--}
 {--gestion interprte sur fichier de sortie, criture des noms de        --}
 {--rservoirs et des rangs de paramtres rservoirs et rseaux sur rsum,--}
 {--criture du mode de rpartition des stockages sur fichier de sortie et --}
 {--rsum.                                                                --}
 P_ContraintesEtConsignesLacs;
520
521
522

 {--LECTURE DES PARAMETRES DE RESEAUX et criture de ces paramtres sur   --}
 {--fichier de sortie                                                     --}
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
 P_ParametresReseau;

 {--calcul des norme d'dition des volumes--}
 for YRangLac:=1 to YnbLacs do
  if YAbsRel=1 then
   YNorme[YRangLac]:=1
  else
   YNorme[YRangLac]:=YVtot[YRangLac];
 if YAbsRel=1 then
  YNorme[0]:=1
 else
  YNorme[0]:=YVtotCumul;

 {--Si saisie clavier, choix du fichier texte des dbits naturels          --}
 {--LECTURE DES DEBITS NATURELS AUX STATIONS, calcul des dbits correspon  --}
 {--dants  la station aval du systme et aux prises et restitutions, et   --}
 {--Ecriture de ceux-ci sur fichier texte (YNomFicResult1 = 'Qobs-P-R.txt')--}
 {--criture du nom de fichier des apports naturel sur fichier de sortie et--}
 {--rsum                                                                 --}
542
543
 if YCalculPossibl then
  begin
544
   C1_CalculeQPriseQrestitution;
545
546
   YCalculPossibl:=YGlopGlopDebit;
   if YCalculPossibl then
547
    P_EcritureSeparation
548
549
550
551
   else
    close(YFicSortie);
  end;

552
553
554
555
 {--Si saisie clavier, CHOIX DU FICHIER DES DEBITS OBJECIFS.         --}
 {--Lecture des dbits objectifs pivots et calcul des dbits objectif--}
 {--pour chaque jour de l'anne                                      --}
 {--criture du nom de fichier d'objectif sur fichier texte et rsum--}
556
557
 if YCalculPossibl then
  begin
558
559
   C1_LectureObjectif;
   if YNbQXobj<=0 then
560
561
562
563
564
565
    begin
     YCalculPossibl:=false;
     close(YFicSortie);
    end;
  end;

566
567
568
569
 {--Si saisie clavier CHOIX DU TYPE D'OBJECTIF A RESPECTER (Qmin ou Qmax  --}
 {--Description de l'hyrogramme objectif et de son type dans le fichier de--}
 {--sortie des paramtres et options du calcul. Description du type de    --}
 {--l'objectif dans le resum                                             --}
570
 if YCalculPossibl then
571
  Ici_ChoixTypeObjectif;
572
573
574
575
576
577
578
579
580
581

 {--CAS OU PARAMETRES ET DONNEES RENDENT LES CALCULS POSSIBLES--}
 if YCalculPossibl then

  begin

   {--OUVERTURE DU FICHIER DES DEBITS NATURELS AUX PRISES ET RESTITUTIONS--}
   assign(YFicResult1,YRepR+YNomFicResult1);
   reset(YFicResult1);

582
   {$R+} {active le contrle de validit d'intervalle}
583

584
585
586
   {--REMPLISSAGE INITIAL DU TABLEAU DES DATES (YDate) ET DEBITS NATURELS, --}
   {--(YQ) aval, prises et restitutions, d'indices YDelaimin  Ydelaimax-1,--}
   {--lus dans fichier texte des dbits aux prises et restitutions         --}
587
588
589
590
   for YRangDelai:=1 to 7 do
    readln(YFicResult1);
   for YRangDelai:=YDelaimin to YDelaimax-1 do
    {--lecture de YDate[YRangDelai] et YQ[YRangdelai]--}
591
    C1_LectureQAvalPriseRestit(YrangDelai);
592
593
594
595

   {--CALCUL DES DATES CRUCIALES : date associe  l'indice 0 du tableau --}
   {--des dbits pour la station aval, et dates correspondantes pour les --}
   {--restitutions (dcales en arrire selon les temps de propagation   --}
596
597
598
599
   {--La date du dbit  une restitution correspond  la date du dbit  --}
   {--la station aval (affecte  midi), dcale du temps de propagation.--}
   {--Elle est constitue d'une partie entire (YDateEnJour) et d'une    --}
   {--partie fractionnaire (YDecalDateQ).                                --}
600
   val(copy(YDate[0],1,4),YAn,YCode);
601
602
603
   val(copy(YDate[0],5,2),YMois[0],YCode);
   val(copy(YDate[0],7,2),YJour[0],YCode);
   YDeb2.YDAteEnJour[0]:=B_Tjour(YAn,YMois[0],YJour[0]);
604
605
   for YrangLac:=1 to YNbLacs do
    begin
606
607
     YDeb2.YDateEnJour[YrangLac]:=trunc(YDeb2.YDateEnJour[0]+0.5-YD[YRangLac]);
     YDecalDateQ[YRanglac]:=frac(YDeb2.YDateEnJour[0]+0.5-YD[YRangLac]);
608
609
    end;

610
611
612
613
614
615
616
   {--Ouverture du fichier typ temporaire d'enregistrements YDeb2        --}
   {--contenant l'objectif global de dbit  stocker (YQXsous) et les     --}
   {--variables suivantes aux diffrents rservoirs: date (YDateEnJour),  --}
   {--dbit naturel  la prise (YQMnat) et dbit naturel  la restitution --}
   {(YQVnat).                                                             --}
   assign(YFicDebTemp,YRepT+YNomFicDebTemp);
   rewrite(YFicDebTemp);
617

618
619
620
621
622
623
   {--Ouverture du fichier texte destin aux dbits synchroniss, et --}
   {--criture d'en tte sur ce fichier                              --}
   Ici_EcritureEntete_QSynchro;

   {--BOUCLE SUR TOUTES LES LIGNES A LIRE dans le fichier texte des--}
   {--dbits naturels (station aval, prises et restitutions)       --}
624
625
   YNbPbBornes:=0;
   for YrangLu:=YDelaimax-YDelaimin+1 to YNbEcrit do
626
    with YDeb2 do
627

628
     BEGIN
629

630
631
632
      YJJ_MM_AAAA[0]:=B_Date(YDateEnJour[0]);
      val(copy(YJJ_MM_AAAA[0],4,2),YMois[0],YCode);
      val(copy(YJJ_MM_AAAA[0],1,2),YJour[0],YCode);
633
      YQuantieme:=B_Quantieme(YJour[0],YMois[0]);
634

635
636
      {--LECTURE DES DEBITS NATURELS pour l'indice YDelaimax--}
      C1_LectureQavalPriseRestit(Ydelaimax);
637

638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
      {--POUR CHAQUE LAC, CALCUL DES DEBITS NATURELS  la prise et  la    --}
      {--restitution, dcals en arrire par rapport au dbit aval d'indice--}
      {--nul, selon les temps de propagation entre restitution et aval.    --}
      YQAntMin:=1e20;
      for YRangLac:=1 to YnbLacs do
       begin
        YQVnat[YrangLac]:=C1_Qinter(2*YrangLac,-YD[Yranglac]);
        YQMnat[YrangLac]:=C1_Qinter(2*YrangLac-1,-YD[Yranglac]);
        YQAntMin:=B_min(YQAntMin,YQVnat[YRangLac]);
        YQAntMin:=B_min(YQAntMin,YQMnat[YRangLac]);
       end; {fin du calc des Q nat. dcals.--}

      {--Calcul de l'objectif de dbit  stocker --}
      if (YQAntMin>=0) and (YQ[0,0]>=0) then
       begin
        YComplet:=1;
        YQXsous:=YQ[0,0]-YQXobj[YQuantieme];
       end
      else
       begin
        YComplet:=0;
        YQXsous:=YLacune;
       end;

      {--ECRITURE SUR FICHIER DES DEBITS SYNCHRONISES--}
      write(YFicResult2,B_Date(YdateEnJour[0]),'   ',YComplet:4,'    ',YQ[0,0]:24:15,'  ',YQXObj[YQuantieme]:19:10,'    ');
      for YRangLac:=1 to YNbLacs do
       write(YFicResult2,B_Date(YDateEnJour[YRangLac]),' ',YQMnat[YrangLac]:21:12,' ',
                         YQVnat[YrangLac]:21:12,'    ');
      writeln(YFicResult2);

      {--Ecriture sur fichier temporaire d'objectif de dbit stock et des--}
      {--dbits nat, res et ref aux prises et restitutions                --}
      write(YFicDebTemp,YDeb2);

      {--DECALAGE D'UN PAS DE TEMPS de toutes les valeurs dans le tableau --}
      {--des dbits naturels aval, prises et restitutions                 --}
      for YRangDelai:=Ydelaimin to Ydelaimax-1 do
       begin
        for YrangLac:=0 to 2*YNbLacs do
         YQ[YRangLac,YRangDelai]:=YQ[Yranglac,YRangDelai+1];
        YDate[YRangDelai]:=YDate[YRangDelai+1];
       end;

      {--DECALAGE D'UN PAS DE TEMPS de toutes les dates cruciales--}
      YDateEnJour[0]:=YDateEnJour[0]+1;
      for YRangLac:=1 to YNbLacs do
       begin
        YDateEnJour[YrangLac]:=YDateEnJour[YrangLac]+1;
        YQVVeille[YrangLac]:=YQV[YRanglac];
       end;
689

690
     END; {fin de boucle sur les lignes  lire dans le fichier des dbits naturels}
691

692
693
   {--FERMETURE DES FICHIERS texte des dbits naturels aux prises et      --}
   {--restitutions, et des dbits synchroniss                            --}
694
695
   close(YFicResult1);
   close(YFicResult2);
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
   close(yFicDebTemp);

   {--RECHERCHE DE PLUS LONGUE PERIODE DE DEBITS SANS LACUNE--}

   reset(YFicResult2);
   {--Achtung : lecture d'en-tte}
   for Yi:=1 to 19 do
    readln(YFicResult2);
   readln(YFicResult2,YJJ_MM_AAAA[0],YComplet);
   val(copy(YJJ_MM_AAAA[0],7,4),YAn,YCode);
   val(copy(YJJ_MM_AAAA[0],4,2),YMois[0],YCode);
   val(copy(YJJ_MM_AAAA[0],1,2),YJour[0],YCode);
   YDeb2.YDAteEnJour[0]:=B_Tjour(YAn,YMois[0],YJour[0]);
   YDateInitiale:=YDeb2.YDateEnJour[0];
   if YComplet=0 then
711
712
    begin
     YDernierLuGlop:=false;
713
     YDateFin:=YDeb2.YDateEnJour[0]-1;
714
715
716
717
    end
   else
    begin
     YdernierLuGlop:=true;
718
     YDateDebut:=YDeb2.YDateEnJour[0];
719
720
721
    end;
   YPlusLongDuree:=0;
   repeat
722
723
    readln(YFicResult2,YJJ_MM_AAAA[0],YComplet);
    YDeb2.YDateEnJour[0]:=YDeb2.YdateEnJour[0]+1;
724
725
    if not YDernierLuGlop then
     begin
726
      if YComplet=1 then
727
       begin
728
        YDateDebut:=YDeb2.YdateEnJour[0];
729
730
731
732
        YDernierLuGlop:=true;
       end
     end
    else
733
     if YComplet=0 then
734
      begin
735
       YDateFin:=YDeb2.YDateEnJour[0]-1;
736
737
738
739
740
741
742
743
       YdernierLuGlop:=false;
       if (YDateFin-YDateDebut+1)>=YPlusLongDuree then
        begin
         YPlusLongDuree:=YDateFin-YDateDebut+1;
         YMeilleurDebut:=YDateDebut;
         YMeilleurFin:=YDateFin;
        end;
      end;
744
   until eof(YFicResult2);
745
746
   if YDernierLuGlop then
    begin
747
     YDateFin:=YDeb2.YDateEnJour[0];
748
749
750
751
752
753
754
755
     if (YDateFin-YDateDebut+1)>=YPlusLongDuree then
      begin
       YPlusLongDuree:=YDateFin-YDateDebut+1;
       YMeilleurDebut:=YDateDebut;
       YMeilleurFin:=YDateFin;
      end;
    end;

756
757
758
759
760
761
   {--CAS OU LES RESULTATS DE DEBIT DEPASSENT 3 ANS D'AFFILE SANS LACUNE. --}
   {--Si saisie clavier, choix de la priode  analyser au sein de la plus--}
   {--longue priode de donnes sans lacunes, sinon utilisation des dates --}
   {--lues sur fichier pour dfinir la priode d'analyse au sein de cette --}
   {--priode sans lacune.                                                --}
   {--Ecriture de la priode d'analyse sur le rsum                      --}
762
763
764
   if YPlusLongDuree>(3*365) then
    begin
     if YSaisie then
765
      E_Saisie6
766
     else
767
768
769
      P_DatesLuesSurFichier;
     YResume:=YResume+B_Date(YDateDebut)+'-'+B_Date(YDateFin);
    end
770
771
772
773
774
775

   {--cas o la priode de rsultats sans lacune est trop courte}
   else
    begin
     YcalculPossibl:=false;
     if YSaisie then
776
      E_Message7;
777
778
779
780
    end;

  end; {fin du cas o paramtres et donnes rendent les calculs possibles}

781

782
783
784
785
786
 {--CAS OU ON FAIT LES CALCULS--}
 if YCalculPossibl then

  begin

787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
   {==CONSTITUTION D'UN FICHIER TYPE DES OBJECTIFS DE DEBIT STOCKE ET DES==}
   {==DEBITS NAT, RES ET REF AUX PRISES ET RESTITUTIONS, SUR LA PERIODE  ==}
   {==CHOISIE POUR LES TRAITEMENTS                                       ==}

   {--remplissage du fichier avec des lacunes--}
   with YDeb2 do
    begin
     YQXsous:=YLacune;
     for YRangLac:=1 to YNbLacs do
      begin
       YDateEnJour[YRangLac]:=0;
       YQMnat[YRangLac]:=YLacune;
      end;
    end;
   assign(YFicDeb,YRepT+YNomFicDeb);
   rewrite(YFicDeb);
   for Yi:=YDateDebut to YdateFin do
    write(YFicDeb,YDeb2);

   {--copie des enregistrements sur le fichier typ tendu sur toute la  --}
   {--priode, et collage sur le fichier typ limit  la priode choisie--}
   {--pour les traitements. Calcul des dbits moyens des cours d'eau aux --}
   {--prises.                                                            --}
   for YrangLac:=1 to YNbLacs do
    YQMmoyen[YRanglac]:=0;
   reset(YFicDebTemp);
   seek(YFicDebTemp,YDateDebut-YDateInitiale);
814
815
   for Yi:=YDateDebut to YDateFin do
    begin
816
817
818
819
820
821
822
     read(YFicDebTemp,YDeb2);
     seek(YFicDeb,filepos(YFicDeb)-1);
     write(YFicDeb,YDeb2);
     seek(YFicDeb,filepos(YFicDeb)-1);
     if Yi>(-YDelaimin) then
      for YRanglac:=1 to YNbLacs do
       YQMmoyen[YRangLac]:=YQMmoyen[YRangLac]+YDeb2.YQMnat[YRangLac];
823
    end;
824
825
826
827
   close(YFicDebTemp);
   close(YFicDeb);
   for YrangLac:=1 to YNbLacs do
    YQMmoyen[YRangLac]:=YQMmoyen[YRangLac]/(YDateFin-YDateDebut+1);
828

829
   {==ASSIGNATION DES FICHIERS TYPES DE VOLUME ET MANQUE DE CAPACITE==}
830
831
832
833
834
835
   for YRangLac:=1 to YNbLacs  do
    begin
     str(YRangLac,Yextension);
     assign(YFicVobj[YRangLac],YRepT+YNomFicV+YExtension);
     assign(YFicVges[YRangLac],YRepT+YNomFicVges+YExtension);
     assign(YFicVop[YRangLac],YRepT+YNomFicVop+YExtension);
836
837
838
    end;
   assign(YFicVdef[1],YRepT+YNomFicVdef+'1');
   assign(YFicVdef[2],YRepT+YNomFicVdef+'2');
839

840
841
842
843
   {--Lecture des paramtres YA et YB de formule de frquence, ainsi que--}
   {--des temps de retour sur fichier idoine, ou valeurs par dfaut le  --}
   {--cas chant.                                                      --}
   C2_LectureParametresStat;
844

845
846
847
848
   {--Lecture du facteur puissance qui sera utilis pour le calcul de la--}
   {--dfaillance moyenne QXdef, ou valeur 2 par dfaut (pour moyenne    --}
   {--quadratique                                                       --}
   C2_LecturePuissanceQXdef_DecideTri;
849
850


851
852
853
854
855
   {--COEUR DU PROGRAMME : CALCULS EN SENS CHRONO INVERS PUIS EN SENS  --}
   {--CHRONO, DE LA REPARTITION DES DEBITS A STOCKER ET DE L'EVOLUTION --}
   {--DU NIVEAU DANS LES RESERVOIRS                                    --}
   C3_CalculChroniquesVolume;
   close(YFicSortie);
856

857
858
859
860
861
862
863
864
865
866
   if YChoixSubProg=0 then begin

    {--CALCUL DE LA PARTIE FRACTIONNAIRE DES DATES DE VOLUME, pour chaque--}
    {--lac. Par rapport aux dates des dbits aval, ces dates de volume   --}
    {--dcales de YdecalDateQ+0.5. A cet instant dans le programme, les --}
    {--volumes prsents dans les fichiers typs de rsultats (stocks    --}
    {--juste avant, par C3_CalculChroniquesVolume), correspondent  ces  --}
    {--dates fractionnaires allant de YDateDebut+0.5+YdecalDateQ+0.5    --}
    {--YdateFin+0.5+YdecalDateQ+0.5.                                     --}
    for YRangLac:=1 to YNbLacs do
867
      begin
868
869
870
871
872
873
      YDecalDateQ[YRanglac]:=YDecalDateQ[YrangLac]+0.5;
      if YdecalDateQ[YrangLac] >=1 then
        begin
        YdateFinLac[YRangLac]:=YdateFinLac[YRanglac]+1;
        YDecalDateQ[YRanglac]:=YDecaldateQ[YRangLac]-1;
        end;
874
875
      end;

876
877
878
    {==A PARTIR DES FICHIERS TYPES DE VOLUME A DATE FRACTIONNAIRE, CALCUL==}
    {==DES VOLUMES A 24H ET ECRITURE DE CEUX-CI SUR LES MEMES FICHIERS   ==}
    {==TYPES, ET ECRITURE DES RESULTATS SUR FICHIER TEXTE DES VOLUMES    ==}
879

880
881
882
    {--criture d'en-tte sur fichier texte de rsultats des chroniques --}
    {--de volume                                                        --}
    Ici_EcritureEntete_Volumes;
883

884
885
886
887
888
889
890
891
892
893
894
    {--ouverture des fichiers typs de volume et de manque de capacit pour --}
    {--chaque lac, et  premire lecture de volume Vobj  date fractionnaire --}
    for YRangLac:=1 to YNbLacs do
      begin
      reset(YFicVobj[YRangLac]);
      reset(YFicVges[YRangLac]);
      reset(YFicVop[YRangLac]);
      read(YFicVobj[YRanglac],YVobjAnt[YrangLac]);
      YVobjMoyrel[1,YRangLac]:=0;
      YVobjMoyrel[2,YRangLac]:=0;
      end;
895

896
897
    {--boucle sur les temps, sur la priode choisie--}
    for Yi:=YDateDebut+1 to YDateFin+1 do
898
899
900

      begin

901
902
      {--criture sur fichier texte de date  la station aval--}
      write(YFicResult4,B_Date(Yi-1),'   ');
903

904
905
      {--boucle sur les lacs--}
      for YRangLac:=1 to YNbLacs do
906
907
908

        begin

909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
        {--cas o on n'est pas arriv  la fin du fichier typ de volume--}
        if Yi<=YdateFin then

          begin

          {--lecture de volume  heure fractionnaire sur fichier typ, --}
          {--CALCUL DU VOLUME A 24h par interpolation linaire et      --}
          {--criture de celui-ci sur mme fichier typ                --}
          read(YFicVobj[YRangLac],YVobjNouv[Yranglac]);
          YVobj[YRangLac]:=YVobjAnt[YrangLac]+(YVobjNouv[YRanglac]-YVobjAnt[Yranglac])*(1-YDecalDateQ[YRangLac]);
          seek(YFicVobj[YRangLac],filepos(YFicVobj[YRangLac])-2);
          write(YFicVobj[YRanglac],YVobj[YRangLac]);
          seek(YficVobj[YRangLac],filepos(YficVobj[Yranglac])+1);
          YVobjMoyRel[1,YRangLac]:=YVobjMoyRel[1,YRangLac]+YVobj[YRangLac];
          end

        {--cas o on est arriv  la fin du fichier typ de volume--}
        else
          begin
          YVobjNouv[YRanglac]:=Ylacune;
          YVobj[YRangLac]:=YLacune;
          end;

        {--lecture de volume de gestion Vobj sur fichier typ--}
        read(YFicVges[YRangLac],YVges[YRangLac]);
        read(YFicVop[YRangLac],YVop[YRangLac]);
        YVobjMoyRel[2,YRangLac]:=YVobjMoyRel[2,YRangLac]+YVges[YRangLac];

        {--criture sur fichier texte de la date lac, du volume Vobj  date--}
        {--fractionnaire et  24h, des volumes Vges et Vop  date          --}
        {--fractionnaire                                                   --}
        if YAbsRel=1 then
          begin
          write(YFicResult4,B_Date(Yi-1+YDateFinLac[YRanglac]-YDateFin),' ',
                          YVobjAnt[YRanglac]:15:2,'   ');
          if YVobj[YRangLac]<>YLacune then
            write(YFicResult4,YVobj[YRangLac]:15:2)
          else
            write(YFicResult4,'               ');
          write(YFicResult4,'   ',YVges[YRangLac]:15:2,'   ',YVop[YRangLac]:15:2,'   ');
          end
        else
          begin
          write(YFicResult4,B_Date(Yi-1+YDateFinLac[YRanglac]-YDateFin),' ',
                          (YVobjAnt[YRanglac]/YVtot[YRangLac]):15:13,'   ');
          if YVobj[YRangLac]<>YLacune then
            write(YficResult4,(YVobj[YRangLac]/YVtot[YRangLac]):15:13)
          else
            write(YFicResult4,'               ');
          write(YFicResult4,'   ',(YVges[YRangLac]/YVtot[YRangLac]):15:13,'   ',(YVop[YRangLac]/YVtot[YRangLac]):15:13,'   ');
          end;
        {--incrmentation du volume  date fractionnaire--}
        YVobjAnt[YRanglac]:=YVobjNouv[YRangLac];
962

963
        end; {fin de boucle sur les lacs--}
964

965
966
      {--passage  la ligne suivante dans le fichier texte de rsultats--}
      writeln(YFicResult4);
967

968
      end; {fin de boucle sur les temps}
969

970
971
972
973
974
975
976
977
978
979
    {--fermeture des fichiers typs de volume et de manque de capacit, --}
    {--ainsi que du fichier texte de rsultats de chronique de volume --}
    for Yranglac:=1 to yNbLacs do
      begin
      close(YFicVobj[YRangLac]);
      close(YFicVges[YRanglac]);
      YVobjMoyRel[1,YRangLac]:=YVobjMoyRel[1,YRangLac]/(YVtot[YrangLac]*(YdateFin-YDateDebut));
      YVobjMoyRel[2,YRangLac]:=YVobjMoyRel[2,YRangLac]/(YVtot[YrangLac]*(YdateFin-YDateDebut+1));
      end;
    close(YFicResult4);
980

981
982
983
    {--remplacement des valeurs  date fractionnaire par les valeurs  24H--}
    {--dans le fichier typ des volumes Vop                               --}
    for YRanglac:=1 to YNbLacs do
984
      begin
985
986
987
988
989
990
991
992
993
994
995
996
      reset(YficVop[YRangLac]);
      read(YFicVop[YRanglac],YVopAnt);
      for Yi:=YDateDebut+1 to YDateFin do
        begin
        read(YFicVop[YRangLac],YVopNouv);
        YVop[YRangLac]:=YVopAnt+(YVopNouv-YVopAnt)*(1-YDecalDateQ[YRangLac]);
        seek(YFicVop[YRangLac],filepos(YFicVop[YRangLac])-2);
        write(YFicVop[YRanglac],YVop[YRangLac]);
        seek(YficVop[YRangLac],filepos(YficVop[Yranglac])+1);
        YVopAnt:=YVopNouv;
        end;
      close(YFicVop[Yranglac]);
997
998
      end;

999
1000
1001
1002
1003
1004
    {Ecriture sur fichier texte de rsultats concernant les dfaillances de }
    {gestion (volumes manquants Vdef dans les rservoirs et dfaut de dbit }
    { la station objectif aval) pour une exploitation visant  optimiser la}
    {rpartition des taches entre rservoirs                                }
    E_Message12;
    Ici_EditionBilan;
1005

1006
1007
    {--message cran--}
    E_Message9;
1008

1009
1010
    {==EDITION DES RESULTATS CLASSES ET ISOFREQUENCE DE VOLUME ET DE ==}
    {==DEBIT DEFAILLANT                                              ==}
1011

1012
1013
1014
    {--boucle sur les rservoirs--}
    for YrangLac:=1 to YNbLacs do
      begin
1015

1016
1017
      {--calcul de la date de dbut de srie des volumes  24h--}
      YDateDebutV:=YDateDebut+YDateFinLac[YRanglac]-YDateFin;
1018

1019
1020
1021
1022
      {--date de fin de srie des volumes  24h : on perd une valeur par --}
      {--rapport  la srie des volumes aux heures brutes du rservoir,  --}
      {-- cause des interpolations linaires                            --}
      YDateFinV:=YDateFinLac[YRangLac]-1;
1023

1024
1025
      {--extension destine au nom des fichiers texte de rsultats--}
      str(YRangLac,Yextension);
1026

1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
1048
1049
1050
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065
      {--calcul et criture sur fichier texte des rsultats isofrquence   --}
      {--des volumes limites Vobj ( 24 H)  respecter dans le rservoir   --}
      assign(YFicResult4,YRepR+YNomFicResultVobj+YExtension+YExIsofreq);
      if YTypeObjectif=0 then
        begin
        YChaine:='Volume d''eau minimal Vobj';
        if YAbsRel=1 then
          YChaine:=YChaine+' (m3)'
        else
          YChaine:=YChaine+'/Vtot';
        YChaine:=YChaine+' en dessous duquel il ne faut pas descendre'
        end
      else
        begin
        YChaine:='Volume d''eau maximal Vobj';
        if YAbsRel=1 then
          YChaine:=YChaine+' (m3)'
        else
          YChaine:=YChaine+'/Vtot';
        YChaine:=YChaine+' qu''il ne faut pas depasser';
        end;
      YChaine:=YChaine+' (a 24h) dans le reservoir '+YNomLac[YRangLac]
                +' pour pouvoir satisfaire du mieux possible l''objectif aval par la suite'
                +' (en evitant au maximum les defaillances causees par un stock ';
      if YTypeObjectif=0 then
        YChaine:=YChaine+'trop vide)'
      else
        YChaine:=YChaine+'trop plein)';
      C2_ISOFREQ(YdateDebutV,YDateDebutV,YDateFinV,YFicVobj[YRangLac],Ychaine,YNorme[YRangLac],YAbsRel);

      {--calcul et criture sur fichier texte des rsultats isofrquence   --}
      {--des volumes limites Vop ( 24H)  respecter dans le rservoir     --}
      assign(YFicResult4,YRepR+YNomFicResultVop+YExtension+YExIsofreq);
      if YTypeObjectif=0 then
        YChaine:='minimale'
      else
        YChaine:='maximale';
      YChaine:='Valeur '+YChaine+' Vop';
      if YAbsRel=1 then
1066
        YChaine:=YChaine+' (m3)'
1067
      else
1068
        YChaine:=YChaine+'/Vtot';
1069
1070
1071
1072
1073
1074
1075
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
1086
1087
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
1103
1104
1105
1106
1107
1108
1109
1110
1111
1112
1113
1114
1115
1116
1117
      YChaine:=YChaine+' du volume dans le reservoir '+YNomLac[YRangLac]+' (a 24h), '
                +'pouvant etre produit quand on vise d''abord la meilleure satisfaction possible de '
                +'l''objectif commun dans l''immediat, puis sa meilleure satisfaction possible '
                +'par la suite';
      C2_ISOFREQ(YdateDebutV,YDateDebutV,YDateFinV,YFicVop[YRangLac],Ychaine,YNorme[YRangLac],YAbsRel);

      end; {fin de boucle sur les rservoirs}

    {--calcul et criture sur fichier texte des rsultats isofrquence   --}
    {--des manques globaux de capacit dans les rservoirs, calculs en  --}
    {--chrono inverse                                                    --}
    assign(YFicResult4,YRepR+YNomFicResultM+'A'+YExIsofreq);
    YChaine:='Manque global Vdef';
    if YAbsRel=2 then
      YChaine:=YChaine+'/somme(Vtot)';
    YChaine:=YChaine+' de capacite des reservoirs (';
    if YabsRel=1 then
      YChaine:=YChaine+'m3, ';
    yChaine:=YChaine+'date aval, calcul en chronologie inverse), responsable '
    +'de defaillances QXdef2 futures s''ajoutant aux inevitables absolues QXdef0 : capacite supplementaire necessaire ';
    if YTypeObjectif=0 then
      YChaine:=YChaine+'au dessus de la crete'
    else
      YChaine:=YChaine+'en dessous du fond';
    if YAbsRel=1 then
      C2_ISOFREQ(YdateDebut,YDateDebut,YDateFin,YFicVdef[1],YChaine,YNorme[0],1)
    else
      C2_ISOFREQ(YdateDebut,YDateDebut,YDateFin,YFicVdef[1],YChaine,YNorme[0],4);

    {--calcul et criture sur fichier texte des rsultats isofrquence   --}
    {--des manques globaux de capacit dans les rservoirs, calculs en  --}
    {--sens chronologique                                                --}
    assign(YFicResult4,YRepR+YNomFicResultM+'B'+YExIsofreq);
    YChaine:='Manque global Vdef';
    if YAbsRel=2 then
      YChaine:=YChaine+'/somme(Vtot)';
    YChaine:=YChaine+'de capacite des reservoirs (';
    if YAbsRel=1 then
      YChaine:=YChaine+'m3, ';
    YChaine:=YChaine+'date aval, calcul en sens chronologique), responsable '
    +'de defaillances QXdef2 passees s''ajoutant aux inevitables absolues QXdef0 : capacite supplementaire necessaire ';
    if YTypeObjectif=0 then
      YChaine:=YChaine+'en dessous du fond'
    else
      YChaine:=YChaine+'au dessus de la crete';
    if YAbsRel=1 then
      C2_ISOFREQ(YdateDebut,YDateDebut,YDateFin,YFicVdef[2],YChaine,YNorme[0],1)
    else
      C2_ISOFREQ(YdateDebut,YDateDebut,YDateFin,YFicVdef[2],YChaine,YNorme[0],4);
1118

1119
1120
1121
1122
1123
1124
1125
1126
1127
1128
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
1139
1140
1141
1142
1143
1144
1145
1146
1147
1148
    {--rsultats de dfaillances de dbit invitables absolus (mme avec des--}
    {--capacits de rservoirs sans limite                                  --}
    assign(YFicResult4,YRepR+YNomFicResDef+'0B'+YExIsofreq);
    if YTypeObjectif=0 then
      YChaine:='Manque'
    else
      YChaine:='Exces';
    YChaine:=YChaine+' de debit inevitable absolu QXdef0 a la station aval du systeme (m3/s, date aval, calcul en sens chronologique) par rapport a '
              +'l''hydrogramme objectif (du aux contraintes physiques et '
              +'consignes de gestion aux prises et restitutions des reservoirs)';
    C2_ISOFREQ(YDateDebut,YDateDebut,YDateFin,YficQXdef[0,2],YChaine,1,3);

    {--rsultats de dfaillances de dbit invitables --}
    assign(YFicResult4,YRepR+YNomFicResDef+'1B'+YExIsofreq);
    if YTypeObjectif=0 then
      YChaine:='Manque'
    else
      YChaine:='Exces';
    YChaine:=YChaine+' de debit QXdef1 a la station aval du systeme (m3/s, date aval, calcul en sens chronologique) par rapport a l''objectif '
              +'(du a la fois aux contraintes physiques et consignes aux prises et restitutions, et au manque de capacite des reservoir)';
    C2_ISOFREQ(YDateDebut,YDateDebut,YDateFin,YficQXdef[1,2],YChaine,1,3);

    {--message cran et stat du calcul--}
    assign(YFicSortie, YRepR+YNomFicSortie);
    append(YFicSortie);
    writeln(YFicSortie);
    P_EcritureSeparation;
    writeln(YFicSortie,'FIN DU CALCUL (LES ETAPES SUCCESSIVES SE SONT DEROULEES NORMALEMENT)');
    close(YFicSortie);
    E_Message10;
1149

1150
    end {if YChoixSubProg=0}
Dorchies David's avatar
Dorchies David committed
1151
1152
1153
  end {fin du cas o on fait les calculs}
 else {cas o les calculs n'ont pas pu tre faits}
  if YSaisie then
1154
   E_Message11;
Dorchies David's avatar
Dorchies David committed
1155
1156

 {cration du fichier marquant la fin programme}
1157
1158
1159
1160
1161
 assign(YFitext,YNomFicFin);
 rewrite(YFitext);
 close(YFitext);

END.