← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/lagr/lagaff.f90

  • Committer: Package Import Robot
  • Author(s): Sylvestre Ledru
  • Date: 2011-11-01 17:43:32 UTC
  • mto: (6.1.7 sid)
  • mto: This revision was merged to the branch mainline in revision 11.
  • Revision ID: package-import@ubuntu.com-20111101174332-tl4vk45no0x3emc3
Tags: upstream-2.1.0
Import upstream version 2.1.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/lagr/lagaff.f90
1
1
!-------------------------------------------------------------------------------
2
2
 
3
 
!     This file is part of the Code_Saturne Kernel, element of the
4
 
!     Code_Saturne CFD tool.
5
 
 
6
 
!     Copyright (C) 1998-2009 EDF S.A., France
7
 
 
8
 
!     contact: saturne-support@edf.fr
9
 
 
10
 
!     The Code_Saturne Kernel is free software; you can redistribute it
11
 
!     and/or modify it under the terms of the GNU General Public License
12
 
!     as published by the Free Software Foundation; either version 2 of
13
 
!     the License, or (at your option) any later version.
14
 
 
15
 
!     The Code_Saturne Kernel is distributed in the hope that it will be
16
 
!     useful, but WITHOUT ANY WARRANTY; without even the implied warranty
17
 
!     of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
18
 
!     GNU General Public License for more details.
19
 
 
20
 
!     You should have received a copy of the GNU General Public License
21
 
!     along with the Code_Saturne Kernel; if not, write to the
22
 
!     Free Software Foundation, Inc.,
23
 
!     51 Franklin St, Fifth Floor,
24
 
!     Boston, MA  02110-1301  USA
 
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.
25
20
 
26
21
!-------------------------------------------------------------------------------
27
22
 
28
23
subroutine lagaff &
29
24
!================
30
25
 
31
 
 ( idbia0 , idbra0 ,                                              &
32
 
   ndim   , ncelet , ncel   , nfac   , nfabor , nfml   , nprfml , &
33
 
   nnod   , lndfac , lndfbr , ncelbr ,                            &
34
 
   nvar   , nscal  , nphas  ,                                     &
 
26
 ( nvar   , nscal  ,                                              &
35
27
   nbpmax , nvp    , nvp1   , nvep   , nivep  ,                   &
36
28
   ntersl , nvlsta , nvisbr ,                                     &
37
 
   nideve , nrdeve , nituse , nrtuse ,                            &
38
 
   ifacel , ifabor , ifmfbr , ifmcel , iprfml ,                   &
39
 
   ipnfac , nodfac , ipnfbr , nodfbr , itepa  ,                   &
40
 
   idevel , ituser , ia     ,                                     &
41
 
   xyzcen , surfac , surfbo , cdgfac , cdgfbo , xyznod , volume , &
 
29
   itepa  ,                                                       &
42
30
   dt     , rtpa   , rtp    , propce , propfa , propfb ,          &
43
31
   coefa  , coefb  ,                                              &
44
 
   ettp   , ettpa  , tepa   , taup   , tlag   , tempct , statis , &
45
 
   w1     , w2     , w3     ,                                     &
46
 
   rdevel , rtuser , ra     )
 
32
   ettp   , ettpa  , tepa   , taup   , tlag   , tempct , statis )
47
33
 
48
34
!===============================================================================
49
35
! FONCTION :
64
50
!__________________.____._____.________________________________________________.
65
51
! name             !type!mode ! role                                           !
66
52
!__________________!____!_____!________________________________________________!
67
 
! idbia0           ! i  ! <-- ! number of first free position in ia            !
68
 
! idbra0           ! i  ! <-- ! number of first free position in ra            !
69
 
! ndim             ! i  ! <-- ! spatial dimension                              !
70
 
! ncelet           ! i  ! <-- ! number of extended (real + ghost) cells        !
71
 
! ncel             ! i  ! <-- ! number of cells                                !
72
 
! nfac             ! i  ! <-- ! number of interior faces                       !
73
 
! nfabor           ! i  ! <-- ! number of boundary faces                       !
74
 
! nfml             ! i  ! <-- ! number of families (group classes)             !
75
 
! nprfml           ! i  ! <-- ! number of properties per family (group class)  !
76
 
! nnod             ! i  ! <-- ! number of vertices                             !
77
 
! lndfac           ! i  ! <-- ! size of nodfac indexed array                   !
78
 
! lndfbr           ! i  ! <-- ! size of nodfbr indexed array                   !
79
 
! ncelbr           ! i  ! <-- ! number of cells with faces on boundary         !
80
53
! nvar             ! i  ! <-- ! total number of variables                      !
81
54
! nscal            ! i  ! <-- ! total number of scalars                        !
82
 
! nphas            ! i  ! <-- ! number of phases                               !
83
55
! nbpmax           ! e  ! <-- ! nombre max de particulies autorise             !
84
56
! nvp              ! e  ! <-- ! nombre de variables particulaires              !
85
57
! nvp1             ! e  ! <-- ! nvp sans position, vfluide, vpart              !
88
60
! ntersl           ! e  ! <-- ! nbr termes sources de couplage retour          !
89
61
! nvlsta           ! e  ! <-- ! nombre de var statistiques lagrangien          !
90
62
! nvisbr           ! e  ! <-- ! nombre de statistiques aux frontieres          !
91
 
! nideve, nrdeve   ! i  ! <-- ! sizes of idevel and rdevel arrays              !
92
 
! nituse, nrtuse   ! i  ! <-- ! sizes of ituser and rtuser arrays              !
93
 
! ifacel(2, nfac)  ! ia ! <-- ! interior faces -> cells connectivity           !
94
 
! ifabor(nfabor)   ! ia ! <-- ! boundary faces -> cells connectivity           !
95
 
! ifmfbr(nfabor)   ! ia ! <-- ! boundary face family numbers                   !
96
 
! ifmcel(ncelet)   ! ia ! <-- ! cell family numbers                            !
97
 
! iprfml           ! ia ! <-- ! property numbers per family                    !
98
 
!  (nfml, nprfml)  !    !     !                                                !
99
 
! ipnfac           ! te ! <-- ! position du premier noeud de chaque            !
100
 
!   (nfac+1)       !    !     !  face interne dans nodfac (optionnel)          !
101
 
! nodfac           ! te ! <-- ! connectivite faces internes/noeuds             !
102
 
!   (lndfac)       !    !     !  (optionnel)                                   !
103
 
! ipnfbr           ! te ! <-- ! position du premier noeud de chaque            !
104
 
!  (nfabor+1)      !    !     !  face de bord dans nodfbr (optionnel)          !
105
 
! nodfbr           ! te ! <-- ! connectivite faces de bord/noeuds              !
106
 
!   (lndfbr  )     !    !     !  (optionnel)                                   !
107
63
! itepa            ! te ! <-- ! info particulaires (entiers)                   !
108
64
! (nbpmax,nivep    !    !     !   (cellule de la particule,...)                !
109
 
! idevel(nideve)   ! ia ! <-> ! integer work array for temporary development   !
110
 
! ituser(nituse)   ! ia ! <-> ! user-reserved integer work array               !
111
 
! ia(*)            ! ia ! --- ! main integer work array                        !
112
 
! xyzcen           ! ra ! <-- ! cell centers                                   !
113
 
!  (ndim, ncelet)  !    !     !                                                !
114
 
! surfac           ! ra ! <-- ! interior faces surface vectors                 !
115
 
!  (ndim, nfac)    !    !     !                                                !
116
 
! surfbo           ! ra ! <-- ! boundary faces surface vectors                 !
117
 
!  (ndim, nfabor)  !    !     !                                                !
118
 
! cdgfac           ! ra ! <-- ! interior faces centers of gravity              !
119
 
!  (ndim, nfac)    !    !     !                                                !
120
 
! cdgfbo           ! ra ! <-- ! boundary faces centers of gravity              !
121
 
!  (ndim, nfabor)  !    !     !                                                !
122
 
! xyznod           ! tr ! <-- ! coordonnes des noeuds                          !
123
 
! (ndim,nnod)      !    !     !                                                !
124
 
! volume(ncelet    ! tr ! <-- ! volume d'un des ncelet elements                !
125
65
! dt(ncelet)       ! ra ! <-- ! time step (per cell)                           !
126
66
! rtp, rtpa        ! ra ! <-- ! calculated variables at cell centers           !
127
67
!  (ncelet, *)     !    !     !  (at current and previous time steps)          !
142
82
!  (nbpmax,2)      !    !     !                                                !
143
83
! statis(ncelet    ! tr ! <-- ! cumul des statistiques volumiques              !
144
84
!    nvlsta)       !    !     !                                                !
145
 
! w1..w3(ncelet    ! tr ! --- ! tableaux de travail                            !
146
 
! rdevel(nrdeve)   ! ra ! <-> ! real work array for temporary development      !
147
 
! rtuser(nrtuse)   ! ra ! <-> ! user-reserved real work array                  !
148
 
! ra(*)            ! ra ! --- ! main real work array                           !
149
85
!__________________!____!_____!________________________________________________!
150
86
 
151
87
!     TYPE : E (ENTIER), R (REEL), A (ALPHANUMERIQUE), T (TABLEAU)
152
88
!            L (LOGIQUE)   .. ET TYPES COMPOSES (EX : TR TABLEAU REEL)
153
89
!     MODE : <-- donnee, --> resultat, <-> Donnee modifiee
154
90
!            --- tableau de travail
 
91
!===============================================================================
 
92
 
 
93
!===============================================================================
 
94
! Module files
 
95
!===============================================================================
 
96
 
 
97
use paramx
 
98
use numvar
 
99
use cstnum
 
100
use optcal
 
101
use pointe
 
102
use entsor
 
103
use parall
 
104
use lagpar
 
105
use lagran
 
106
use cstphy
 
107
use mesh
155
108
 
156
109
!===============================================================================
157
110
 
158
111
implicit none
159
112
 
160
 
!===============================================================================
161
 
! Common blocks
162
 
!===============================================================================
163
 
 
164
 
include "paramx.h"
165
 
include "numvar.h"
166
 
include "cstnum.h"
167
 
include "optcal.h"
168
 
include "pointe.h"
169
 
include "entsor.h"
170
 
include "parall.h"
171
 
include "lagpar.h"
172
 
include "lagran.h"
173
 
include "cstphy.h"
174
 
 
175
 
!===============================================================================
176
 
 
177
113
! Arguments
178
114
 
179
 
integer          idbia0 , idbra0
180
 
integer          ndim   , ncelet , ncel   , nfac   , nfabor
181
 
integer          nfml   , nprfml
182
 
integer          nnod   , lndfac , lndfbr , ncelbr
183
 
integer          nvar   , nscal  , nphas
 
115
integer          nvar   , nscal
184
116
integer          nbpmax , nvp    , nvp1   , nvep  , nivep
185
117
integer          ntersl , nvlsta , nvisbr
186
 
integer          nideve , nrdeve , nituse , nrtuse
187
 
integer          ifacel(2,nfac) , ifabor(nfabor)
188
 
integer          ifmfbr(nfabor) , ifmcel(ncelet)
189
 
integer          iprfml(nfml,nprfml)
190
 
integer          ipnfac(nfac+1) , nodfac(lndfac)
191
 
integer          ipnfbr(nfabor+1) , nodfbr(lndfbr)
 
118
 
192
119
integer          itepa(nbpmax,nivep)
193
 
integer          idevel(nideve), ituser(nituse)
194
 
integer          ia(*)
195
120
 
196
 
double precision xyzcen(ndim,ncelet)
197
 
double precision surfac(ndim,nfac) , surfbo(ndim,nfabor)
198
 
double precision cdgfac(ndim,nfac) , cdgfbo(ndim,nfabor)
199
 
double precision xyznod(ndim,nnod) , volume(ncelet)
200
121
double precision dt(ncelet) , rtp(ncelet,*) , rtpa(ncelet,*)
201
122
double precision propce(ncelet,*)
202
123
double precision propfa(nfac,*) , propfb(nfabor,*)
205
126
double precision tepa(nbpmax,nvep)
206
127
double precision taup(nbpmax) , tlag(nbpmax,3) , tempct(nbpmax,2)
207
128
double precision statis(ncelet,nvlsta)
208
 
double precision w1(ncelet), w2(ncelet), w3(ncelet)
209
 
double precision rdevel(nrdeve) , rtuser(nrtuse)
210
 
double precision ra(*)
211
129
 
212
130
! Local variables
213
131
 
214
 
integer          idebia , idebra
215
132
 
216
 
integer          iphas
217
133
double precision dnbpr
218
134
 
219
135
! NOMBRE DE PASSAGES DANS LA ROUTINE
228
144
! 0.  GESTION MEMOIRE
229
145
!===============================================================================
230
146
 
231
 
idebia = idbia0
232
 
idebra = idbra0
233
147
 
234
148
!===============================================================================
235
149
! 1. INITIALISATION
236
150
!===============================================================================
237
151
 
238
 
iphas = ilphas
239
152
ipass = ipass + 1
240
153
 
241
154
!===============================================================================
279
192
         nbpart        , dnbpar        ,                          &
280
193
         nbpnew        , dnbpnw        ,                          &
281
194
         nbpout-nbperr , dnbpou-dnbper ,                          &
 
195
         nbpdep        , dnbdep        ,                          &
282
196
         nbperr        , dnbper        ,                          &
283
197
         dnbpr         ,                                          &
284
198
         npcsup        , dnpcsu        ,                          &
293
207
         nbpart        , dnbpar        ,                          &
294
208
         nbpnew        , dnbpnw        ,                          &
295
209
         nbpout-nbperr , dnbpou-dnbper ,                          &
 
210
         nbpdep        , dnbdep        ,                          &
296
211
         nbperr        , dnbper        ,                          &
297
212
         dnbpr         ,                                          &
298
213
         npcsup        , dnpcsu        ,                          &
306
221
         nbpart        , dnbpar        ,                          &
307
222
         nbpnew        , dnbpnw        ,                          &
308
223
         nbpout-nbperr , dnbpou-dnbper ,                          &
 
224
         nbpdep        , dnbdep        ,                          &
309
225
         nbperr        , dnbper        ,                          &
310
226
         dnbpr         ,                                          &
311
227
         npencr        , dnpenc
316
232
         nbpart        , dnbpar        ,                          &
317
233
         nbpnew        , dnbpnw        ,                          &
318
234
         nbpout-nbperr , dnbpou-dnbper ,                          &
 
235
         nbpdep        , dnbdep        ,                          &
319
236
         nbperr        , dnbper        ,                          &
320
237
         dnbpr
321
238
 
331
248
!--------
332
249
 
333
250
 1000 format('# ** INFORMATIONS SUR LE CALCUL LAGRANGIEN     ',/, &
334
 
       '#    -------------------------------------     ',/, &
335
 
       '#                                         ',/,      &
336
 
       '# colonne  1 : numero de pas de temps      ',/,     &
337
 
       '# colonne  2 : temps physique              ',/,     &
338
 
       '# colonne  3 : nbre inst. de part.    ',/,          &
339
 
       '# colonne  4 : nbre inst. de part. (avec poids)',/, &
340
 
       '# colonne  5 : nbre inst. de part. injectees',/,    &
341
 
       '# colonne  6 : nbre inst. de part. injectees',            &
342
 
       ' (avec poids)',/,                                   &
343
 
       '# colonne  7 : nbre inst. de part. sorties ou deposees',/,&
344
 
       '# colonne  8 : nbre inst. de part. sorties ou deposees',  &
345
 
       ' (avec poids)',/,                                   &
346
 
       '# colonne  9 : nbre inst. de part. perdues (reperage) ',/,&
347
 
       '# colonne 10 : nbre inst. de part. perdues',              &
348
 
        ' (reperage, avec poids)',/,                        &
349
 
       '# colonne 11 : % de part. perdues',/,               &
350
 
       '# colonne 12 : nbre inst. de part. qui ont subi le',      &
351
 
       ' clonage',/,                                        &
352
 
       '# colonne 13 : nbre inst. de part. qui ont subi le',      &
353
 
       ' clonage (avec poids)',/,                           &
354
 
       '# colonne 14 : nbre inst. de nouvel. part. par clonage',/,&
355
 
       '# colonne 15 : nbre inst. de nouvel. part. par clonage',  &
356
 
       ' (avec poids)',/,                                   &
357
 
       '# colonne 16 : nbre inst. de nouvel. part. eliminees par',&
358
 
       ' roulette russe ',/,                                &
359
 
       '# colonne 17 : nbre inst. de nouvel. part. eliminees par',&
360
 
       ' roulette russe (avec poids)',/,                    &
361
 
       '# colonne 18 : nbre inst. de part encrassees',            &
362
 
       ' (Charbon)) '/,                                     &
363
 
       '# colonne 19 : nbre inst. de part encrassees',            &
364
 
       ' (Charbon, avec poids))',/,                         &
 
251
       '#    -------------------------------------     '      ,/, &
 
252
       '#                                              '      ,/, &
 
253
       '# colonne  1 : numero de pas de temps          '      ,/, &
 
254
       '# colonne  2 : temps physique                  '      ,/, &
 
255
       '# colonne  3 : nbre inst. de part.             '      ,/, &
 
256
       '# colonne  4 : nbre inst. de part. (avec poids)'      ,/, &
 
257
       '# colonne  5 : nbre inst. de part. injectees   '      ,/, &
 
258
       '# colonne  6 : nbre inst. de part. injectees   '      ,   &
 
259
       ' (avec poids)'                                        ,/, &
 
260
       '# colonne  7 : nbre inst. de part. sorties, ou deposees et supprimees',/,&
 
261
       '# colonne  8 : nbre inst. de part. sorties, ou deposees et supprimees',  &
 
262
       ' (avec poids)'                                        ,/, &
 
263
       '# colonne  9 : nbre inst. de part. deposees '         ,/, &
 
264
       '# colonne 10 : nbre inst. de part. deposees '         ,   &
 
265
       ' (avec poids)'                                        ,/, &
 
266
       '# colonne 11 : nbre inst. de part. perdues (reperage)',/, &
 
267
       '# colonne 12 : nbre inst. de part. perdues'           ,   &
 
268
       ' (reperage, avec poids)'                              ,/, &
 
269
       '# colonne 13 : % de part. perdues'                    ,/, &
 
270
       '# colonne 14 : nbre inst. de part. qui ont subi le'   ,   &
 
271
       ' clonage'                                             ,/, &
 
272
       '# colonne 15 : nbre inst. de part. qui ont subi le'   ,   &
 
273
       ' clonage (avec poids)'                                ,/, &
 
274
       '# colonne 16 : nbre inst. de nouvel. part. par clonage',/,&
 
275
       '# colonne 17 : nbre inst. de nouvel. part. par clonage',  &
 
276
       ' (avec poids)'                                        ,/, &
 
277
       '# colonne 18 : nbre inst. de nouvel. part. eliminees par',&
 
278
       ' roulette russe '                                     ,/, &
 
279
       '# colonne 19 : nbre inst. de nouvel. part. eliminees par',&
 
280
       ' roulette russe (avec poids)'                         ,/, &
 
281
       '# colonne 20 : nbre inst. de part encrassees'         ,   &
 
282
       ' (Charbon) '                                          ,/, &
 
283
       '# colonne 21 : nbre inst. de part encrassees'         ,   &
 
284
       ' (Charbon, avec poids)'                               ,/, &
365
285
       '# ')
366
286
 
367
287
 1001 format('# ** INFORMATIONS SUR LE CALCUL LAGRANGIEN     ',/, &
368
 
       '#    -------------------------------------     ',/, &
369
 
       '#                                         ',/,      &
370
 
       '# colonne  1 : numero de pas de temps      ',/,     &
371
 
       '# colonne  2 : temps physique              ',/,     &
372
 
       '# colonne  3 : nbre inst. de part.    ',/,          &
373
 
       '# colonne  4 : nbre inst. de part. (avec poids)',/, &
374
 
       '# colonne  5 : nbre inst. de part. injectees',/,    &
375
 
       '# colonne  6 : nbre inst. de part. injectees',            &
376
 
       ' (avec poids)',/,                                   &
377
 
       '# colonne  7 : nbre inst. de part. sorties ou deposees',/,&
378
 
       '# colonne  8 : nbre inst. de part. sorties ou deposees',  &
379
 
       ' (avec poids)',/,                                   &
380
 
       '# colonne  9 : nbre inst. de part. perdues (reperage) ',/,&
381
 
       '# colonne 10 : nbre inst. de part. perdues',              &
382
 
        ' (reperage, avec poids)',/,                        &
383
 
       '# colonne 11 : % de part. perdues',/,               &
384
 
       '# colonne 12 : nbre inst. de part. qui ont subi le',      &
385
 
       ' clonage',/,                                        &
386
 
       '# colonne 13 : nbre inst. de part. qui ont subi le',      &
387
 
       ' clonage (avec poids)',/,                           &
388
 
       '# colonne 14 : nbre inst. de nouvel. part. par clonage',/,&
389
 
       '# colonne 15 : nbre inst. de nouvel. part. par clonage',  &
390
 
       ' (avec poids)',/,                                   &
391
 
       '# colonne 16 : nbre inst. de nouvel. part. eliminees par',&
392
 
       ' roulette russe ',/,                                &
393
 
       '# colonne 17 : nbre inst. de nouvel. part. eliminees par',&
394
 
       ' roulette russe (avec poids)',/,                    &
 
288
       '#    -------------------------------------     '      ,/, &
 
289
       '#                                              '      ,/, &
 
290
       '# colonne  1 : numero de pas de temps          '      ,/, &
 
291
       '# colonne  2 : temps physique                  '      ,/, &
 
292
       '# colonne  3 : nbre inst. de part.             '      ,/, &
 
293
       '# colonne  4 : nbre inst. de part. (avec poids)'      ,/, &
 
294
       '# colonne  5 : nbre inst. de part. injectees   '      ,/, &
 
295
       '# colonne  6 : nbre inst. de part. injectees   '        , &
 
296
       ' (avec poids)'                                        ,/, &
 
297
       '# colonne  7 : nbre inst. de part. sorties, ou deposees et supprimees',/,&
 
298
       '# colonne  8 : nbre inst. de part. sorties, ou deposees et supprimees',  &
 
299
       ' (avec poids)'                                        ,/, &
 
300
       '# colonne  9 : nbre inst. de part. deposees    '      ,/, &
 
301
       '# colonne 10 : nbre inst. de part. deposees    '        , &
 
302
       ' (avec poids)'                                        ,/, &
 
303
       '# colonne 11 : nbre inst. de part. perdues (reperage) ',/,&
 
304
       '# colonne 12 : nbre inst. de part. perdues'             , &
 
305
       ' (reperage, avec poids)'                              ,/, &
 
306
       '# colonne 13 : % de part. perdues'                    ,/, &
 
307
       '# colonne 14 : nbre inst. de part. qui ont subi le'     , &
 
308
       ' clonage'                                             ,/, &
 
309
       '# colonne 15 : nbre inst. de part. qui ont subi le'     , &
 
310
       ' clonage (avec poids)'                                ,/, &
 
311
       '# colonne 16 : nbre inst. de nouvel. part. par clonage',/,&
 
312
       '# colonne 17 : nbre inst. de nouvel. part. par clonage'  ,&
 
313
       ' (avec poids)'                                        ,/, &
 
314
       '# colonne 18 : nbre inst. de nouvel. part. eliminees par',&
 
315
       ' roulette russe '                                     ,/, &
 
316
       '# colonne 19 : nbre inst. de nouvel. part. eliminees par',&
 
317
       ' roulette russe (avec poids) '                        ,/, &
395
318
       '# ')
396
319
 
397
320
 1002 format('# ** INFORMATIONS SUR LE CALCUL LAGRANGIEN     ',/, &
398
 
       '#    -------------------------------------     ',/, &
399
 
       '#                                         ',/,      &
400
 
       '# colonne  1 : numero de pas de temps      ',/,     &
401
 
       '# colonne  2 : temps physique              ',/,     &
402
 
       '# colonne  3 : nbre inst. de part.    ',/,          &
403
 
       '# colonne  4 : nbre inst. de part. (avec poids)',/, &
404
 
       '# colonne  5 : nbre inst. de part. injectees',/,    &
405
 
       '# colonne  6 : nbre inst. de part. injectees',            &
406
 
       ' (avec poids)',/,                                   &
407
 
       '# colonne  7 : nbre inst. de part. sorties ou deposees',/,&
408
 
       '# colonne  8 : nbre inst. de part. sorties ou deposees',  &
409
 
       ' (avec poids)',/,                                   &
410
 
       '# colonne  9 : nbre inst. de part. perdues (reperage)',/, &
411
 
       '# colonne 10 : nbre inst. de part. perdues',              &
412
 
        ' (reperage, avec poids)',/,                        &
413
 
       '# colonne 11 : % de part. perdues ',/,              &
414
 
       '# colonne 12 : nbre inst. de part. encrassees',           &
415
 
       ' (Charbon)) '/,                                     &
416
 
       '# colonne 13 : nbre inst. de part. encrassees',           &
417
 
       ' (Charbon, avec poids))',/,                         &
 
321
       '#    -------------------------------------     '      ,/, &
 
322
       '#                                              '      ,/, &
 
323
       '# colonne  1 : numero de pas de temps          '      ,/, &
 
324
       '# colonne  2 : temps physique                  '      ,/, &
 
325
       '# colonne  3 : nbre inst. de part.             '      ,/, &
 
326
       '# colonne  4 : nbre inst. de part. (avec poids)'      ,/, &
 
327
       '# colonne  5 : nbre inst. de part. injectees   '      ,/, &
 
328
       '# colonne  6 : nbre inst. de part. injectees   '        , &
 
329
       ' (avec poids)'                                        ,/, &
 
330
       '# colonne  7 : nbre inst. de part. sorties, ou deposees et supprimees',/,&
 
331
       '# colonne  8 : nbre inst. de part. sorties, ou deposees et supprimees',  &
 
332
       ' (avec poids)'                                        ,/, &
 
333
       '# colonne  9 : nbre inst. de part. deposees'          ,/, &
 
334
       '# colonne 10 : nbre inst. de part. deposees'            , &
 
335
       ' (avec poids)'                                        ,/, &
 
336
       '# colonne 11 : nbre inst. de part. perdues (reperage)',/, &
 
337
       '# colonne 12 : nbre inst. de part. perdues'             , &
 
338
       ' (reperage, avec poids)'                              ,/, &
 
339
       '# colonne 13 : % de part. perdues'                    ,/, &
 
340
       '# colonne 14 : nbre inst. de part. encrassees'          , &
 
341
       ' (Charbon) '                                          ,/, &
 
342
       '# colonne 15 : nbre inst. de part. encrassees'          , &
 
343
       ' (Charbon, avec poids)'                               ,/, &
418
344
       '# ')
419
345
 
420
346
 1003 format('# ** INFORMATIONS SUR LE CALCUL LAGRANGIEN     ',/, &
421
 
       '#    -------------------------------------     ',/, &
422
 
       '#                                         ',/,      &
423
 
       '# colonne  1 : numero de pas de temps      ',/,     &
424
 
       '# colonne  2 : temps physique              ',/,     &
425
 
       '# colonne  3 : nbre inst. de part.    ',/,          &
426
 
       '# colonne  4 : nbre inst. de part. (avec poids)',/, &
427
 
       '# colonne  5 : nbre inst. de part. injectees',/,    &
428
 
       '# colonne  6 : nbre inst. de part. injectees',            &
429
 
       ' (avec poids)',/,                                   &
430
 
       '# colonne  7 : nbre inst. de part. sorties ou deposees',/,&
431
 
       '# colonne  8 : nbre inst. de part. sorties ou deposees',  &
432
 
       ' (avec poids)',/,                                   &
433
 
       '# colonne  9 : nbre inst. de part. perdues (reperage)',/, &
434
 
       '# colonne 10 : nbre inst. de part. perdues',              &
435
 
        ' (reperage, avec poids)',/,                        &
436
 
       '# colonne 11 : % de part. perdues ',/,              &
 
347
       '#    -------------------------------------     '      ,/, &
 
348
       '#                                              '      ,/, &
 
349
       '# colonne  1 : numero de pas de temps          '      ,/, &
 
350
       '# colonne  2 : temps physique                  '      ,/, &
 
351
       '# colonne  3 : nbre inst. de part.             '      ,/, &
 
352
       '# colonne  4 : nbre inst. de part. (avec poids)'      ,/, &
 
353
       '# colonne  5 : nbre inst. de part. injectees   '      ,/, &
 
354
       '# colonne  6 : nbre inst. de part. injectees   '        , &
 
355
       ' (avec poids)'                                        ,/, &
 
356
       '# colonne  7 : nbre inst. de part. sorties, ou deposees et supprimees',/,&
 
357
       '# colonne  8 : nbre inst. de part. sorties, ou deposees et supprimees',  &
 
358
       ' (avec poids)'                                        ,/, &
 
359
       '# colonne  9 : nbre inst. de part. deposees'          ,/, &
 
360
       '# colonne 10 : nbre inst. de part. deposees'            , &
 
361
       ' (avec poids)'                                        ,/, &
 
362
       '# colonne 11 : nbre inst. de part. perdues (reperage) ',/,&
 
363
       '# colonne 12 : nbre inst. de part. perdues'             , &
 
364
        ' (reperage, avec poids)'                             ,/, &
 
365
       '# colonne 13 : % de part. perdues'                    ,/, &
437
366
       '# ')
438
367
 
439
 
 2000 format(1x,i8,2x,e10.4,2x,4(i8,2x,e10.4),2x,e10.4,4(i8,2x,e10.4))
440
 
 2001 format(1x,i8,2x,e10.4,2x,4(i8,2x,e10.4),2x,e10.4,3(i8,2x,e10.4))
441
 
 2002 format(1x,i8,2x,e10.4,2x,4(i8,2x,e10.4),2x,e10.4,1(i8,2x,e10.4))
442
 
 2003 format(1x,i8,2x,e10.4,2x,4(i8,2x,e10.4),2x,e10.4)
 
368
 2000 format(1x,i8,2x,e10.4,2x,5(i8,2x,e10.4),2x,e10.4,4(i8,2x,e10.4))
 
369
 2001 format(1x,i8,2x,e10.4,2x,5(i8,2x,e10.4),2x,e10.4,3(i8,2x,e10.4))
 
370
 2002 format(1x,i8,2x,e10.4,2x,5(i8,2x,e10.4),2x,e10.4,1(i8,2x,e10.4))
 
371
 2003 format(1x,i8,2x,e10.4,2x,5(i8,2x,e10.4),2x,e10.4)
443
372
 
444
373
!====
445
374
! FIN