← Back to branch summary

~ubuntu-branches/ubuntu/precise/code-saturne/precise

« back to all changes in this revision

Viewing changes to src/comb/cs_fuel_param.f90

  • Committer: Package Import Robot
  • Author(s): Sylvestre Ledru
  • Date: 2011-11-24 00:00:08 UTC
  • mfrom: (6.1.9 sid)
  • Revision ID: package-import@ubuntu.com-20111124000008-2vo99e38267942q5
Tags: 2.1.0-3
Install a missing file

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/comb/cs_fuel_param.f90
 
1
!-------------------------------------------------------------------------------
 
2
 
 
3
! This file is part of Code_Saturne, a general-purpose CFD tool.
 
4
!
 
5
! Copyright (C) 1998-2011 EDF S.A.
 
6
!
 
7
! This program is free software; you can redistribute it and/or modify it under
 
8
! the terms of the GNU General Public License as published by the Free Software
 
9
! Foundation; either version 2 of the License, or (at your option) any later
 
10
! version.
 
11
!
 
12
! This program is distributed in the hope that it will be useful, but WITHOUT
 
13
! ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
 
14
! FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
 
15
! details.
 
16
!
 
17
! You should have received a copy of the GNU General Public License along with
 
18
! this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
 
19
! Street, Fifth Floor, Boston, MA 02110-1301, USA.
 
20
 
 
21
!-------------------------------------------------------------------------------
 
22
 
 
23
subroutine cs_fuel_param
 
24
!=======================
 
25
 
 
26
!===============================================================================
 
27
!  FONCTION  :
 
28
!  ---------
 
29
 
 
30
!       INIT DES OPTIONS DES VARIABLES TRANSPORTEES
 
31
!                     ET DES VARIABLES ALGEBRIQUES
 
32
!               COMBUSTION FUEL
 
33
 
 
34
!-------------------------------------------------------------------------------
 
35
! Arguments
 
36
!__________________.____._____.________________________________________________.
 
37
! name             !type!mode ! role                                           !
 
38
!__________________!____!_____!________________________________________________!
 
39
!__________________!____!_____!________________________________________________!
 
40
 
 
41
!     TYPE : E (ENTIER), R (REEL), A (ALPHANUMERIQUE), T (TABLEAU)
 
42
!            L (LOGIQUE)   .. ET TYPES COMPOSES (EX : TR TABLEAU REEL)
 
43
!     MODE : <-- donnee, --> resultat, <-> Donnee modifiee
 
44
!            --- tableau de travail
 
45
!===============================================================================
 
46
 
 
47
!===============================================================================
 
48
! Module files
 
49
!===============================================================================
 
50
 
 
51
use paramx
 
52
use dimens
 
53
use numvar
 
54
use optcal
 
55
use cstphy
 
56
use entsor
 
57
use cstnum
 
58
use ppppar
 
59
use ppthch
 
60
use coincl
 
61
use cpincl
 
62
use cs_fuel_incl
 
63
use ppincl
 
64
use ppcpfu
 
65
 
 
66
!===============================================================================
 
67
 
 
68
implicit none
 
69
 
 
70
integer          ipp , ii , jj , iok , icla
 
71
integer          isc
 
72
double precision wmolme
 
73
 
 
74
!===============================================================================
 
75
! 0. VERIFICATION ISCALT, ISCSTH
 
76
!===============================================================================
 
77
!     L'utilisateur ne doit pas y avoir touche.
 
78
 
 
79
if(iscalt.ne.-1) then
 
80
  write(nfecra,1000)iscalt
 
81
  call csexit (1)
 
82
  !==========
 
83
endif
 
84
do ii = 1, nscapp
 
85
  if(iscsth(iscapp(ii)).ne.-10) then
 
86
    write(nfecra,1001)ii,iscapp(ii),iscapp(ii),iscsth(iscapp(ii))
 
87
    call csexit (1)
 
88
    !==========
 
89
  endif
 
90
enddo
 
91
 
 
92
!===============================================================================
 
93
! 1. VARIABLES TRANSPORTEES
 
94
!===============================================================================
 
95
 
 
96
! --> Definition des scamin et des scamax des variables transportees
 
97
 
 
98
! ---- Variables propres a la phase dispersee
 
99
 
 
100
scamin(ihm) = -grand
 
101
scamax(ihm) = +grand
 
102
 
 
103
do icla=1,nclafu
 
104
  scamin(ing(icla))   = 0.d0
 
105
  scamax(ing(icla))   = +rinfin
 
106
  scamin(iyfol(icla)) = 0.d0
 
107
  scamax(iyfol(icla)) = 4.d-1
 
108
  scamin(ih2(icla))   = h02fol
 
109
  scamax(ih2(icla))   = +grand
 
110
enddo
 
111
 
 
112
! ---- Variables propres a la phase continue
 
113
 
 
114
scamin(ifvap)  = 0.d0
 
115
scamax(ifvap)  = 1.d0
 
116
!
 
117
! Oxydant
 
118
if ( noxyd .ge. 2 ) then
 
119
  scamin(if4m) = 0.d0
 
120
  scamax(if4m) = 1.d0
 
121
endif
 
122
!
 
123
if ( noxyd .eq. 3 ) then
 
124
  scamin(if5m) = 0.d0
 
125
  scamax(if5m) = 1.d0
 
126
endif
 
127
! Combustion heterogene
 
128
scamin(if7m)   = 0.d0
 
129
scamax(if7m)   = 1.d0
 
130
! Variance
 
131
scamin(ifvp2m)  = 0.d0
 
132
scamax(ifvp2m)  = 0.25d0
 
133
! Modele de CO
 
134
if ( ieqco2 .ge. 1 ) then
 
135
  scamin(iyco2) = 0.d0
 
136
  scamax(iyco2) = 1.d0
 
137
endif
 
138
! Modele de Nox
 
139
if ( ieqnox .eq. 1 ) then
 
140
  scamin(iyhcn)  = 0.d0
 
141
  scamax(iyhcn)  = 1.d0
 
142
  scamin(iyno)   = 0.d0
 
143
  scamax(iyno)   = 1.d0
 
144
  scamin(ihox) = -grand
 
145
  scamax(ihox) = +grand
 
146
endif
 
147
 
 
148
! --> Nature des scalaires transportes
 
149
 
 
150
do isc = 1, nscapp
 
151
 
 
152
! ---- Type de scalaire (0 passif, 1 temperature en K
 
153
!                                 -1 temperature en C
 
154
!                                  2 enthalpie)
 
155
!      La distinction -1/1 sert pour le rayonnement
 
156
  iscsth(iscapp(isc)) = 0
 
157
 
 
158
enddo
 
159
 
 
160
! ---- On resout en enthalpie avec un CP constant (Cf. cpvarp)
 
161
 
 
162
iscalt = ihm
 
163
iscsth(ihm)   = 2
 
164
 
 
165
! --> Donnees physiques ou numeriques propres aux scalaires CP
 
166
 
 
167
do isc = 1, nscapp
 
168
 
 
169
  jj = iscapp(isc)
 
170
 
 
171
  if ( iscavr(jj) .le. 0 ) then
 
172
 
 
173
!        En combustion on considere que la viscosite turbulente domine
 
174
!        ON S'INTERDIT DONC LE CALCUL DES FLAMMES LAMINAIRES AVEC Le =/= 1
 
175
 
 
176
    visls0(jj) = viscl0
 
177
 
 
178
  endif
 
179
 
 
180
! ------ Schmidt ou Prandtl turbulent
 
181
 
 
182
  sigmas(jj) = 0.7d0
 
183
 
 
184
! ------ Coeff dissipation des fluctuations
 
185
 
 
186
  rvarfl(jj) = 0.8d0
 
187
 
 
188
  ii = isca(iscapp(isc))
 
189
 
 
190
! ------ Niveau de detail des impressions pour les variables et
 
191
!          donc les scalaires (valeurs 0 ou 1)
 
192
!          Si = -10000 non modifie par l'utilisateur -> niveau 1
 
193
  if(iwarni(ii).eq.-10000) then
 
194
    iwarni(ii) = 1
 
195
  endif
 
196
 
 
197
! ------ Informations relatives a la resolution des scalaires
 
198
 
 
199
!         - Facteur multiplicatif du pas de temps
 
200
  cdtvar(ii) = 1.d0
 
201
 
 
202
!         - Schema convectif % schema 2ieme ordre
 
203
!           = 0 : upwind
 
204
!           = 1 : second ordre
 
205
  blencv(ii) = 0.d0
 
206
 
 
207
!         - Type de schema convetif second ordre (utile si BLENCV > 0)
 
208
!           = 0 : Second Order Linear Upwind
 
209
!           = 1 : Centre
 
210
  ischcv(ii) = 1
 
211
 
 
212
!         - Test de pente pour basculer d'un schema centre vers l'upwind
 
213
!           = 0 : utlisation automatique du test de pente
 
214
!           = 1 : calcul sans test de pente
 
215
  isstpc(ii) = 0
 
216
 
 
217
!         - Reconstruction des flux de convetion et de diffusion aux faces
 
218
!           = 0 : pas de reconstruction
 
219
  ircflu(ii) = 0
 
220
 
 
221
enddo
 
222
 
 
223
 
 
224
!---> Variable courante : nom, sortie chrono, suivi listing, sortie histo
 
225
 
 
226
!     Comme pour les autres variables,
 
227
!       si l'on n'affecte pas les tableaux suivants,
 
228
!       les valeurs par defaut seront utilisees
 
229
 
 
230
!     NOMVAR( ) = nom de la variable
 
231
!     ICHRVR( ) = sortie chono (oui 1/non 0)
 
232
!     ILISVR( ) = suivi listing (oui 1/non 0)
 
233
!     IHISVR( ) = sortie historique (nombre de sondes et numeros)
 
234
!     si IHISVR(.,1)  = -1 sortie sur toutes les sondes
 
235
 
 
236
!     NB : Les 8 premiers caracteres du noms seront repris dans le
 
237
!          listing 'developpeur'
 
238
 
 
239
 
 
240
! ---- Variables propres a la suspension gaz - particules
 
241
 
 
242
ipp = ipprtp(isca(ihm))
 
243
NOMVAR(IPP)  = 'Enthalpy'
 
244
ichrvr(ipp)  = 1
 
245
ilisvr(ipp)  = 1
 
246
ihisvr(ipp,1)= -1
 
247
 
 
248
! ---- Variables propres a la phase dispersee
 
249
 
 
250
do icla = 1, nclafu
 
251
  ipp = ipprtp(isca(iyfol(icla)))
 
252
  WRITE(NOMVAR(IPP),'(A8,I2.2)')'YFOL_FOL' ,ICLA
 
253
  ichrvr(ipp)  = 1
 
254
  ilisvr(ipp)  = 1
 
255
  ihisvr(ipp,1)= -1
 
256
  ipp = ipprtp(isca(ih2(icla)))
 
257
  WRITE(NOMVAR(IPP),'(A7,I2.2)')'HLF_FOL' ,ICLA
 
258
  ichrvr(ipp)  = 1
 
259
  ilisvr(ipp)  = 1
 
260
  ihisvr(ipp,1)= -1
 
261
  ipp = ipprtp(isca(ing(icla)))
 
262
  WRITE(NOMVAR(IPP),'(A6,I2.2)')'NG_FOL' ,ICLA
 
263
  ichrvr(ipp)  = 1
 
264
  ilisvr(ipp)  = 1
 
265
  ihisvr(ipp,1)= -1
 
266
enddo
 
267
 
 
268
! ---- Variables propres a la phase gaz
 
269
 
 
270
ipp = ipprtp(isca(ifvap))
 
271
NOMVAR(IPP)  = 'Fr_VAP'
 
272
ichrvr(ipp)  = 1
 
273
ilisvr(ipp)  = 1
 
274
ihisvr(ipp,1)= -1
 
275
ipp = ipprtp(isca(if7m))
 
276
NOMVAR(IPP)  = 'Fr_HET'
 
277
ichrvr(ipp)  = 1
 
278
ilisvr(ipp)  = 1
 
279
ihisvr(ipp,1)= -1
 
280
ipp = ipprtp(isca(ifvp2m))
 
281
NOMVAR(IPP)  = 'Var_CB'
 
282
ichrvr(ipp)  = 1
 
283
ilisvr(ipp)  = 1
 
284
ihisvr(ipp,1)= -1
 
285
 
 
286
if ( ieqco2 .ge. 1 ) then
 
287
  ipp = ipprtp(isca(iyco2))
 
288
  NOMVAR(IPP)  = 'FR_CO2'
 
289
  ichrvr(ipp)  = 1
 
290
  ilisvr(ipp)  = 1
 
291
  ihisvr(ipp,1)= -1
 
292
endif
 
293
 
 
294
if ( ieqnox .eq. 1 ) then
 
295
  ipp = ipprtp(isca(iyhcn))
 
296
  NOMVAR(IPP)  = 'FR_HCN'
 
297
  ichrvr(ipp)  = 1
 
298
  ilisvr(ipp)  = 1
 
299
  ihisvr(ipp,1)= -1
 
300
  ipp = ipprtp(isca(iyno))
 
301
  NOMVAR(IPP)  = 'FR_NO'
 
302
  ichrvr(ipp)  = 1
 
303
  ilisvr(ipp)  = 1
 
304
  ihisvr(ipp,1)= -1
 
305
  ipp = ipprtp(isca(ihox))
 
306
  NOMVAR(IPP)  = 'Enth_Ox'
 
307
  ichrvr(ipp)  = 1
 
308
  ilisvr(ipp)  = 1
 
309
  ihisvr(ipp,1)= -1
 
310
endif
 
311
 
 
312
!===============================================================================
 
313
! 2. VARIABLES ALGEBRIQUES OU D'ETAT
 
314
!===============================================================================
 
315
 
 
316
 
 
317
! ---> Variables algebriques propres a la suspension gaz - particules
 
318
 
 
319
ipp = ipppro(ipproc(immel))
 
320
NOMVAR(IPP)   = 'XM'
 
321
ichrvr(ipp)   = 0
 
322
ilisvr(ipp)   = 0
 
323
ihisvr(ipp,1) = -1
 
324
 
 
325
! ---> Variables algebriques propres a la phase dispersee
 
326
 
 
327
do icla = 1, nclafu
 
328
  ipp = ipppro(ipproc(itemp2(icla)))
 
329
  write(nomvar(ipp),'(a7,i2.2)')'Tem_FOL' ,icla
 
330
  ichrvr(ipp)   = 1
 
331
  ilisvr(ipp)   = 1
 
332
  ihisvr(ipp,1) = -1
 
333
  ipp = ipppro(ipproc(irom2(icla)))
 
334
  write(nomvar(ipp),'(a7,i2.2)')'Rho_FOL' ,icla
 
335
  ichrvr(ipp)   = 1
 
336
  ilisvr(ipp)   = 1
 
337
  ihisvr(ipp,1) = -1
 
338
  ipp = ipppro(ipproc(idiam2(icla)))
 
339
  write(nomvar(ipp),'(a6,i2.2)')'Dia_gt' ,icla
 
340
  ichrvr(ipp)   = 1
 
341
  ilisvr(ipp)   = 1
 
342
  ihisvr(ipp,1) = -1
 
343
  ipp = ipppro(ipproc(ih1hlf(icla)))
 
344
  write(nomvar(ipp),'(a6,i2.2)')'H1-Hlf' ,icla
 
345
  ichrvr(ipp)   = 1
 
346
  ilisvr(ipp)   = 1
 
347
  ihisvr(ipp,1) = -1
 
348
  ipp = ipppro(ipproc(igmeva(icla)))
 
349
  write(nomvar(ipp),'(a6,i2.2)')'Ga_EVA' ,icla
 
350
  ichrvr(ipp)   = 1
 
351
  ilisvr(ipp)   = 1
 
352
  ihisvr(ipp,1) = -1
 
353
  ipp = ipppro(ipproc(igmhtf(icla)))
 
354
  write(nomvar(ipp),'(a6,i2.2)')'Ga_HET' ,icla
 
355
  ichrvr(ipp)   = 1
 
356
  ilisvr(ipp)   = 1
 
357
  ihisvr(ipp,1) = -1
 
358
enddo
 
359
 
 
360
! ---> Variables algebriques propres a la phase continue
 
361
 
 
362
ipp = ipppro(ipproc(itemp1))
 
363
NOMVAR(IPP)   = 'Temp_GAZ'
 
364
ichrvr(ipp)   = 1
 
365
ilisvr(ipp)   = 1
 
366
ihisvr(ipp,1) = -1
 
367
ipp = ipppro(ipproc(irom1))
 
368
NOMVAR(IPP)   = 'ROM_GAZ'
 
369
ichrvr(ipp)   = 1
 
370
ilisvr(ipp)   = 1
 
371
ihisvr(ipp,1) = -1
 
372
ipp = ipppro(ipproc(iym1(1)))
 
373
NOMVAR(IPP)   = 'YM_FO0'
 
374
ichrvr(ipp)   =  0
 
375
ilisvr(ipp)   =  0
 
376
ihisvr(ipp,1) = -1
 
377
ipp = ipppro(ipproc(iym1(2)))
 
378
NOMVAR(IPP)   = 'YM_FOV'
 
379
ichrvr(ipp)   = 1
 
380
ilisvr(ipp)   = 1
 
381
ihisvr(ipp,1) = -1
 
382
ipp = ipppro(ipproc(iym1(3)))
 
383
NOMVAR(IPP)   = 'YM_CO'
 
384
ichrvr(ipp)   = 1
 
385
ilisvr(ipp)   = 1
 
386
ihisvr(ipp,1) = -1
 
387
ipp = ipppro(ipproc(iym1(4)))
 
388
NOMVAR(IPP)   = 'YM_H2S'
 
389
ichrvr(ipp)   = 1
 
390
ilisvr(ipp)   = 1
 
391
ihisvr(ipp,1) = -1
 
392
ipp = ipppro(ipproc(iym1(5)))
 
393
NOMVAR(IPP)   = 'YM_H2'
 
394
ichrvr(ipp)   = 1
 
395
ilisvr(ipp)   = 1
 
396
ihisvr(ipp,1) = -1
 
397
ipp = ipppro(ipproc(iym1(6)))
 
398
NOMVAR(IPP)   = 'YM_HCN'
 
399
ichrvr(ipp)   = 1
 
400
ilisvr(ipp)   = 1
 
401
ihisvr(ipp,1) = -1
 
402
ipp = ipppro(ipproc(iym1(7)))
 
403
NOMVAR(IPP)   = 'YM_NH3'
 
404
ichrvr(ipp)   = 1
 
405
ilisvr(ipp)   = 1
 
406
ihisvr(ipp,1) = -1
 
407
ipp = ipppro(ipproc(iym1(8)))
 
408
NOMVAR(IPP)   = 'YM_O2'
 
409
ichrvr(ipp)   = 1
 
410
ilisvr(ipp)   = 1
 
411
ihisvr(ipp,1) = -1
 
412
ipp = ipppro(ipproc(iym1(9)))
 
413
NOMVAR(IPP)   = 'YM_CO2'
 
414
ichrvr(ipp)   = 1
 
415
ilisvr(ipp)   = 1
 
416
ihisvr(ipp,1) = -1
 
417
ipp = ipppro(ipproc(iym1(10)))
 
418
NOMVAR(IPP)   = 'YM_H2O'
 
419
ichrvr(ipp)   = 1
 
420
ilisvr(ipp)   = 1
 
421
ihisvr(ipp,1) = -1
 
422
ipp = ipppro(ipproc(iym1(11)))
 
423
NOMVAR(IPP)   = 'YM_SO2'
 
424
ichrvr(ipp)   = 1
 
425
ilisvr(ipp)   = 1
 
426
ihisvr(ipp,1) = -1
 
427
ipp = ipppro(ipproc(iym1(12)))
 
428
NOMVAR(IPP)   = 'YM_N2'
 
429
ichrvr(ipp)   = 1
 
430
ilisvr(ipp)   = 1
 
431
ihisvr(ipp,1) = -1
 
432
!
 
433
if ( ieqnox .eq. 1 ) then
 
434
  ipp = ipppro(ipproc(ighcn1))
 
435
  NOMVAR(IPP)   = 'EXP1'
 
436
  ichrvr(ipp)   = 1
 
437
  ilisvr(ipp)   = 1
 
438
  ihisvr(ipp,1) = -1
 
439
  ipp = ipppro(ipproc(ighcn2))
 
440
  NOMVAR(IPP)   = 'EXP2'
 
441
  ichrvr(ipp)   = 1
 
442
  ilisvr(ipp)   = 1
 
443
  ihisvr(ipp,1) = -1
 
444
  ipp = ipppro(ipproc(ignoth))
 
445
  NOMVAR(IPP)   = 'EXP3'
 
446
  ichrvr(ipp)   = 1
 
447
  ilisvr(ipp)   = 1
 
448
  ihisvr(ipp,1) = -1
 
449
endif
 
450
!
 
451
ipp = ipppro(ipproc(ibcarbone))
 
452
NOMVAR(IPP)   = 'Bilan_C'
 
453
ichrvr(ipp)   = 1
 
454
ilisvr(ipp)   = 1
 
455
ihisvr(ipp,1) = -1
 
456
ipp = ipppro(ipproc(iboxygen))
 
457
NOMVAR(IPP)   = 'Bilan_O'
 
458
ichrvr(ipp)   = 1
 
459
ilisvr(ipp)   = 1
 
460
ihisvr(ipp,1) = -1
 
461
ipp = ipppro(ipproc(ibhydrogen))
 
462
NOMVAR(IPP)   = 'Bilan_H'
 
463
ichrvr(ipp)   = 1
 
464
ilisvr(ipp)   = 1
 
465
ihisvr(ipp,1) = -1
 
466
!
 
467
!===============================================================================
 
468
! 3. INFORMATIONS COMPLEMENTAIRES
 
469
!===============================================================================
 
470
 
 
471
! ---> Initialisation
 
472
 
 
473
! ---- Calcul de RO0 a partir de T0 et P0
 
474
!        (loi des gaz parfaits applliquee a l'air)
 
475
!    On initialise RO0 avec l'oxydant 1 qui est sense etre
 
476
!    l'oxydant majoritaire
 
477
 
 
478
wmolme = ( wmole(io2) *oxyo2(1) +wmole(in2) *oxyn2(1)             &
 
479
          +wmole(ih2o)*oxyh2o(1)+wmole(ico2)*oxyco2(1))           &
 
480
        /(oxyo2(1)+oxyn2(1)+oxyh2o(1)+oxyco2(1))
 
481
ro0 = p0*wmolme / (rr*t0)
 
482
 
 
483
! ---- Initialisation pour la masse volumique du coke
 
484
!    Maintenant c'est fait dans FULECD
 
485
!     RHOKF = RHO0FL
 
486
 
 
487
! ---> Coefficient de relaxation de la masse volumique
 
488
!      RHO(n+1) = SRROM * RHO(n) + (1-SRROM) * RHO(n+1)
 
489
srrom = 0.90d0
 
490
 
 
491
! ---> Viscosite laminaire associee au scalaire enthalpie
 
492
!       DIFTL0 (diffusivite dynamique en kg/(m s))
 
493
!     C'est cette valeur par defaut qui est TOUJOURS utilisee dans les
 
494
!       calculs charbon (un peu etonnant de ne pas prendre en
 
495
!       compte les variations de ce parametre physique si on
 
496
!       recherche des informations sur les flux thermiques aux parois)
 
497
 
 
498
diftl0      = 4.25d-5
 
499
visls0(ihm) = diftl0
 
500
 
 
501
! ---> Masse volumique variable et viscosite constante (pour les suites)
 
502
irovar = 1
 
503
ivivar = 0
 
504
 
 
505
!===============================================================================
 
506
! 4. ON REDONNE LA MAIN A L'UTLISATEUR
 
507
!===============================================================================
 
508
 
 
509
call user_fuel_ini1
 
510
!==================
 
511
 
 
512
!===============================================================================
 
513
! 5. VERIFICATION DES DONNERS FOURNIES PAR L'UTLISATEUR
 
514
!===============================================================================
 
515
 
 
516
iok = 0
 
517
call cs_fuel_verify(iok)
 
518
!==================
 
519
 
 
520
if(iok.gt.0) then
 
521
  write(nfecra,9999)iok
 
522
  call csexit (1)
 
523
  !==========
 
524
else
 
525
  write(nfecra,9998)
 
526
endif
 
527
 
 
528
 
 
529
 1000 format(                                                     &
 
530
'@                                                            ',/,&
 
531
'@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,&
 
532
'@                                                            ',/,&
 
533
'@ @@ ATTENTION : ARRET A L''ENTREE DES DONNEES               ',/,&
 
534
'@    =========                                               ',/,&
 
535
'@    PHYSIQUE PARTICULIERE (FUEL) DEMANDEE                   ',/,&
 
536
'@                                                            ',/,&
 
537
'@  La valeur de ISCALT est renseignee automatiquement.       ',/,&
 
538
'@                                                            ',/,&
 
539
'@  L''utilisateur ne doit pas la renseigner dans usini1, or  ',/,&
 
540
'@    elle a ete affectee comme suit :                        ',/,&
 
541
'@    ISCALT = ',I10                                           ,/,&
 
542
'@                                                            ',/,&
 
543
'@  Le calcul ne sera pas execute.                            ',/,&
 
544
'@                                                            ',/,&
 
545
'@  Verifier usini1.                                          ',/,&
 
546
'@                                                            ',/,&
 
547
'@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,&
 
548
'@                                                            ',/)
 
549
 1001 format(                                                     &
 
550
'@                                                            ',/,&
 
551
'@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,&
 
552
'@                                                            ',/,&
 
553
'@ @@ ATTENTION : ARRET A L''ENTREE DES DONNEES               ',/,&
 
554
'@    =========                                               ',/,&
 
555
'@    PHYSIQUE PARTICULIERE (FUEL) DEMANDEE                   ',/,&
 
556
'@                                                            ',/,&
 
557
'@  Les valeurs de ISCSTH sont renseignees automatiquement.   ',/,&
 
558
'@                                                            ',/,&
 
559
'@  L''utilisateur ne doit pas les renseigner dans usini1, or ',/,&
 
560
'@    pour le scalaire ',I10   ,' correspondant au scalaire   ',/,&
 
561
'@    physique particuliere ',I10   ,' on a                   ',/,&
 
562
'@    ISCSTH(',I10   ,') = ',I10                               ,/,&
 
563
'@                                                            ',/,&
 
564
'@  Le calcul ne sera pas execute.                            ',/,&
 
565
'@                                                            ',/,&
 
566
'@  Verifier usini1.                                          ',/,&
 
567
'@                                                            ',/,&
 
568
'@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,&
 
569
'@                                                            ',/)
 
570
 9998 format(                                                     &
 
571
'                                                             ',/,&
 
572
' Pas d erreur detectee lors de la verification des donnees   ',/,&
 
573
'                                            (user_fuel_ini1).',/)
 
574
 9999 format(                                                     &
 
575
'@                                                            ',/,&
 
576
'@                                                            ',/,&
 
577
'@                                                            ',/,&
 
578
'@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,&
 
579
'@                                                            ',/,&
 
580
'@ @@ ATTENTION : ARRET A L''ENTREE DES DONNEES               ',/,&
 
581
'@    =========                                               ',/,&
 
582
'@    LES PARAMETRES DE CALCUL SONT INCOHERENTS OU INCOMPLETS ',/,&
 
583
'@                                                            ',/,&
 
584
'@  Le calcul ne sera pas execute (',I10,' erreurs).          ',/,&
 
585
'@                                                            ',/,&
 
586
'@  Se reporter aux impressions precedentes pour plus de      ',/,&
 
587
'@    renseignements.                                         ',/,&
 
588
'@  Verifier user_fuel_ini1.                                  ',/,&
 
589
'@                                                            ',/,&
 
590
'@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,&
 
591
'@                                                            ',/)
 
592
 
 
593
!----
 
594
! End
 
595
!----
 
596
 
 
597
return
 
598
end subroutine