VGEST.PAS 49.1 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
121
122
123
124
125
126
127
128
129
130
131
 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
      write(YFicSortie,'              ',B_DateClaire(2001,YDateObj[Yj]));
     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
248
249
250
251
252

 {-----------------------------------------------------------------------}
 {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;

  var
   Xi : integer;
   Xj : integer;
   XQXdef : array[1..2,1..2] of double;

253
  begin
254
255
256
257
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

   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
286
   else
287
288
289
290
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
394
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
437
438
439
440
441
442
443
444
445
446
447
448
449
450
    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]);
        C4_TriChronique(YFicQXdef[Xi,Xj],YFicQXdefT[Xi,Xj],YNbQXdefNul[Xi,Xj],YTouComNul,YQXdefmoy[Xi,Xj],YQXdefmed[Xi,Xj],YQXdefmax[Xi,Xj],YQXdefmin[Xi,Xj],YQXdefrare[Xi,Xj],YQXdeffrequent[Xi,Xj]);
       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)');
     while not eof(YFicQXdefT[1,1]) do
      begin
       for Xi:=1 to 2 do
        for Xj:=1 to 2 do
         begin
          read(YFicQXdefT[Xi,Xj],XQXdef[Xi,Xj]);
          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);

451
  end;
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467

{===========================================================================}
{-- 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;

468
 if Ysaisie then
469
  E_Message1;
470
471
472

 YCalculPossibl:=true;

473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
 {--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;
507
508
509

 {--LECTURE DES PARAMETRES DE RESEAUX et criture de ces paramtres sur   --}
 {--fichier de sortie                                                     --}
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
 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                                                                 --}
529
530
 if YCalculPossibl then
  begin
531
   C1_CalculeQPriseQrestitution;
532
533
   YCalculPossibl:=YGlopGlopDebit;
   if YCalculPossibl then
534
    P_EcritureSeparation
535
536
537
538
   else
    close(YFicSortie);
  end;

539
540
541
542
 {--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--}
543
544
 if YCalculPossibl then
  begin
545
546
   C1_LectureObjectif;
   if YNbQXobj<=0 then
547
548
549
550
551
552
    begin
     YCalculPossibl:=false;
     close(YFicSortie);
    end;
  end;

553
554
555
556
 {--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                                             --}
557
 if YCalculPossibl then
558
  Ici_ChoixTypeObjectif;
559
560
561
562
563
564
565
566
567
568

 {--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);

569
   {$R+} {active le contrle de validit d'intervalle}
570

571
572
573
   {--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         --}
574
575
576
577
   for YRangDelai:=1 to 7 do
    readln(YFicResult1);
   for YRangDelai:=YDelaimin to YDelaimax-1 do
    {--lecture de YDate[YRangDelai] et YQ[YRangdelai]--}
578
    C1_LectureQAvalPriseRestit(YrangDelai);
579
580
581
582

   {--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   --}
583
584
585
586
   {--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).                                --}
587
   val(copy(YDate[0],1,4),YAn,YCode);
588
589
590
   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]);
591
592
   for YrangLac:=1 to YNbLacs do
    begin
593
594
     YDeb2.YDateEnJour[YrangLac]:=trunc(YDeb2.YDateEnJour[0]+0.5-YD[YRangLac]);
     YDecalDateQ[YRanglac]:=frac(YDeb2.YDateEnJour[0]+0.5-YD[YRangLac]);
595
596
    end;

597
598
599
600
601
602
603
   {--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);
604

605
606
607
608
609
610
   {--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)       --}
611
612
   YNbPbBornes:=0;
   for YrangLu:=YDelaimax-YDelaimin+1 to YNbEcrit do
613
    with YDeb2 do
614

615
     BEGIN
616

617
618
619
620
      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);
      YQuantieme:=B_Quantieme(YJour[0],YMois[0],2001);
621

622
623
      {--LECTURE DES DEBITS NATURELS pour l'indice YDelaimax--}
      C1_LectureQavalPriseRestit(Ydelaimax);
624

625
626
627
628
629
630
631
632
633
634
635
636
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
      {--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;
676

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

679
680
   {--FERMETURE DES FICHIERS texte des dbits naturels aux prises et      --}
   {--restitutions, et des dbits synchroniss                            --}
681
682
   close(YFicResult1);
   close(YFicResult2);
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
   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
698
699
    begin
     YDernierLuGlop:=false;
700
     YDateFin:=YDeb2.YDateEnJour[0]-1;
701
702
703
704
    end
   else
    begin
     YdernierLuGlop:=true;
705
     YDateDebut:=YDeb2.YDateEnJour[0];
706
707
708
    end;
   YPlusLongDuree:=0;
   repeat
709
710
    readln(YFicResult2,YJJ_MM_AAAA[0],YComplet);
    YDeb2.YDateEnJour[0]:=YDeb2.YdateEnJour[0]+1;
711
712
    if not YDernierLuGlop then
     begin
713
      if YComplet=1 then
714
       begin
715
        YDateDebut:=YDeb2.YdateEnJour[0];
716
717
718
719
        YDernierLuGlop:=true;
       end
     end
    else
720
     if YComplet=0 then
721
      begin
722
       YDateFin:=YDeb2.YDateEnJour[0]-1;
723
724
725
726
727
728
729
730
       YdernierLuGlop:=false;
       if (YDateFin-YDateDebut+1)>=YPlusLongDuree then
        begin
         YPlusLongDuree:=YDateFin-YDateDebut+1;
         YMeilleurDebut:=YDateDebut;
         YMeilleurFin:=YDateFin;
        end;
      end;
731
   until eof(YFicResult2);
732
733
   if YDernierLuGlop then
    begin
734
     YDateFin:=YDeb2.YDateEnJour[0];
735
736
737
738
739
740
741
742
     if (YDateFin-YDateDebut+1)>=YPlusLongDuree then
      begin
       YPlusLongDuree:=YDateFin-YDateDebut+1;
       YMeilleurDebut:=YDateDebut;
       YMeilleurFin:=YDateFin;
      end;
    end;

743
744
745
746
747
748
   {--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                      --}
749
750
751
   if YPlusLongDuree>(3*365) then
    begin
     if YSaisie then
752
      E_Saisie6
753
     else
754
755
756
      P_DatesLuesSurFichier;
     YResume:=YResume+B_Date(YDateDebut)+'-'+B_Date(YDateFin);
    end
757
758
759
760
761
762

   {--cas o la priode de rsultats sans lacune est trop courte}
   else
    begin
     YcalculPossibl:=false;
     if YSaisie then
763
      E_Message7;
764
765
766
767
    end;

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

768

769
770
771
772
773
 {--CAS OU ON FAIT LES CALCULS--}
 if YCalculPossibl then

  begin

774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
   {==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);
801
802
   for Yi:=YDateDebut to YDateFin do
    begin
803
804
805
806
807
808
809
     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];
810
    end;
811
812
813
814
   close(YFicDebTemp);
   close(YFicDeb);
   for YrangLac:=1 to YNbLacs do
    YQMmoyen[YRangLac]:=YQMmoyen[YRangLac]/(YDateFin-YDateDebut+1);
815

816
   {==ASSIGNATION DES FICHIERS TYPES DE VOLUME ET MANQUE DE CAPACITE==}
817
818
819
820
821
822
   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);
823
824
825
    end;
   assign(YFicVdef[1],YRepT+YNomFicVdef+'1');
   assign(YFicVdef[2],YRepT+YNomFicVdef+'2');
826

827
828
829
830
   {--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;
831

832
833
834
835
   {--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;
836
837


838
839
840
841
842
   {--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);
843
844

   {--CALCUL DE LA PARTIE FRACTIONNAIRE DES DATES DE VOLUME, pour chaque--}
845
846
847
848
849
850
   {--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.                                     --}
851
852
853
854
855
856
857
858
859
860
   for YRangLac:=1 to YNbLacs do
    begin
     YDecalDateQ[YRanglac]:=YDecalDateQ[YrangLac]+0.5;
     if YdecalDateQ[YrangLac] >=1 then
      begin
       YdateFinLac[YRangLac]:=YdateFinLac[YRanglac]+1;
       YDecalDateQ[YRanglac]:=YDecaldateQ[YRangLac]-1;
      end;
    end;

861
862
863
   {==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    ==}
864
865

   {--criture d'en-tte sur fichier texte de rsultats des chroniques --}
866
867
   {--de volume                                                        --}
   Ici_EcritureEntete_Volumes;
868
869

   {--ouverture des fichiers typs de volume et de manque de capacit pour --}
870
   {--chaque lac, et  premire lecture de volume Vobj  date fractionnaire --}
871
872
873
874
   for YRangLac:=1 to YNbLacs do
    begin
     reset(YFicVobj[YRangLac]);
     reset(YFicVges[YRangLac]);
875
     reset(YFicVop[YRangLac]);
876
     read(YFicVobj[YRanglac],YVobjAnt[YrangLac]);
877
878
     YVobjMoyrel[1,YRangLac]:=0;
     YVobjMoyrel[2,YRangLac]:=0;
879
880
881
882
883
884
885
886
    end;

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

    begin

     {--criture sur fichier texte de date  la station aval--}
887
     write(YFicResult4,B_Date(Yi-1),'   ');
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902

     {--boucle sur les lacs--}
     for YRangLac:=1 to YNbLacs do

      begin

       {--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]);
903
         YVobj[YRangLac]:=YVobjAnt[YrangLac]+(YVobjNouv[YRanglac]-YVobjAnt[Yranglac])*(1-YDecalDateQ[YRangLac]);
904
         seek(YFicVobj[YRangLac],filepos(YFicVobj[YRangLac])-2);
905
         write(YFicVobj[YRanglac],YVobj[YRangLac]);
906
         seek(YficVobj[YRangLac],filepos(YficVobj[Yranglac])+1);
907
         YVobjMoyRel[1,YRangLac]:=YVobjMoyRel[1,YRangLac]+YVobj[YRangLac];
908
909
910
911
912
913
        end

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

917
918
919
920
       {--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];
921

922
923
924
925
926
927
928
929
930
931
932
933
934
       {--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
935
       else
936
937
938
939
940
941
942
943
944
        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;
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
       {--incrmentation du volume  date fractionnaire--}
       YVobjAnt[YRanglac]:=YVobjNouv[YRangLac];

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

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

    end; {fin de boucle sur les temps}

   {--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]);
961
962
     YVobjMoyRel[1,YRangLac]:=YVobjMoyRel[1,YRangLac]/(YVtot[YrangLac]*(YdateFin-YDateDebut));
     YVobjMoyRel[2,YRangLac]:=YVobjMoyRel[2,YRangLac]/(YVtot[YrangLac]*(YdateFin-YDateDebut+1));
963
964
965
966
967
968
969
970
971
972
973
974
    end;
   close(YFicResult4);

   {--remplacement des valeurs  date fractionnaire par les valeurs  24H--}
   {--dans le fichier typ des volumes Vop                               --}
   for YRanglac:=1 to YNbLacs do
    begin
     reset(YficVop[YRangLac]);
     read(YFicVop[YRanglac],YVopAnt);
     for Yi:=YDateDebut+1 to YDateFin do
      begin
       read(YFicVop[YRangLac],YVopNouv);
975
       YVop[YRangLac]:=YVopAnt+(YVopNouv-YVopAnt)*(1-YDecalDateQ[YRangLac]);
976
       seek(YFicVop[YRangLac],filepos(YFicVop[YRangLac])-2);
977
       write(YFicVop[YRanglac],YVop[YRangLac]);
978
979
980
981
982
983
       seek(YficVop[YRangLac],filepos(YficVop[Yranglac])+1);
       YVopAnt:=YVopNouv;
      end;
     close(YFicVop[Yranglac]);
    end;

984
985
986
987
988
989
   {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;
990
991

   {--message cran--}
992
   E_Message9;
993
994
995
996
997
998
999
1000
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013

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

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

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

     {--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;

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

     {--calcul et criture sur fichier texte des rsultats isofrquence   --}
     {--des volumes limites Vobj ( 24 H)  respecter dans le rservoir   --}
1014
     assign(YFicResult4,YRepR+YNomFicResultVobj+YExtension+YExIsofreq);
1015
     if YTypeObjectif=0 then
1016
1017
1018
1019
1020
1021
1022
1023
      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
1024
     else
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
      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 ';
1036
     if YTypeObjectif=0 then
1037
      YChaine:=YChaine+'trop vide)'
1038
     else
1039
1040
      YChaine:=YChaine+'trop plein)';
     C2_ISOFREQ(YdateDebutV,YDateDebutV,YDateFinV,YFicVobj[YRangLac],Ychaine,YNorme[YRangLac],YAbsRel);
1041
1042
1043

     {--calcul et criture sur fichier texte des rsultats isofrquence   --}
     {--des volumes limites Vop ( 24H)  respecter dans le rservoir     --}
1044
     assign(YFicResult4,YRepR+YNomFicResultVop+YExtension+YExIsofreq);
1045
1046
1047
1048
     if YTypeObjectif=0 then
      YChaine:='minimale'
     else
      YChaine:='maximale';
1049
1050
1051
     YChaine:='Valeur '+YChaine+' Vop';
     if YAbsRel=1 then
      YChaine:=YChaine+' (m3)'
1052
     else
1053
1054
1055
1056
1057
1058
      YChaine:=YChaine+'/Vtot';
     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);
1059
1060
1061

    end; {fin de boucle sur les rservoirs}

1062
1063
1064
1065
1066
1067
1068
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
   {--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);

1104
1105
   {--rsultats de dfaillances de dbit invitables absolus (mme avec des--}
   {--capacits de rservoirs sans limite                                  --}
1106
   assign(YFicResult4,YRepR+YNomFicResDef+'0B'+YExIsofreq);
1107
1108
1109
1110
   if YTypeObjectif=0 then
    YChaine:='Manque'
   else
    YChaine:='Exces';
1111
1112
1113
1114
   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);
1115
1116

   {--rsultats de dfaillances de dbit invitables --}
1117
   assign(YFicResult4,YRepR+YNomFicResDef+'1B'+YExIsofreq);
1118
1119
1120
1121
   if YTypeObjectif=0 then
    YChaine:='Manque'
   else
    YChaine:='Exces';
1122
1123
1124
   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);
1125
1126
1127
1128
1129

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

Dorchies David's avatar
Dorchies David committed
1135
  end {fin du cas o on fait les calculs}
1136

Dorchies David's avatar
Dorchies David committed
1137
1138
 else {cas o les calculs n'ont pas pu tre faits}
  if YSaisie then
1139
   E_Message11;
Dorchies David's avatar
Dorchies David committed
1140
1141

 {cration du fichier marquant la fin programme}
1142
1143
1144
1145
1146
 assign(YFitext,YNomFicFin);
 rewrite(YFitext);
 close(YFitext);

END.