← Back to branch summary

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

« back to all changes in this revision

Viewing changes to users/rayt/usray1.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:
users/rayt/usray1.f90
2
2
 
3
3
!VERS
4
4
 
5
 
 
6
 
!     This file is part of the Code_Saturne Kernel, element of the
7
 
!     Code_Saturne CFD tool.
8
 
 
9
 
!     Copyright (C) 1998-2009 EDF S.A., France
10
 
 
11
 
!     contact: saturne-support@edf.fr
12
 
 
13
 
!     The Code_Saturne Kernel is free software; you can redistribute it
14
 
!     and/or modify it under the terms of the GNU General Public License
15
 
!     as published by the Free Software Foundation; either version 2 of
16
 
!     the License, or (at your option) any later version.
17
 
 
18
 
!     The Code_Saturne Kernel is distributed in the hope that it will be
19
 
!     useful, but WITHOUT ANY WARRANTY; without even the implied warranty
20
 
!     of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
21
 
!     GNU General Public License for more details.
22
 
 
23
 
!     You should have received a copy of the GNU General Public License
24
 
!     along with the Code_Saturne Kernel; if not, write to the
25
 
!     Free Software Foundation, Inc.,
26
 
!     51 Franklin St, Fifth Floor,
27
 
!     Boston, MA  02110-1301  USA
 
5
! This file is part of Code_Saturne, a general-purpose CFD tool.
 
6
!
 
7
! Copyright (C) 1998-2011 EDF S.A.
 
8
!
 
9
! This program is free software; you can redistribute it and/or modify it under
 
10
! the terms of the GNU General Public License as published by the Free Software
 
11
! Foundation; either version 2 of the License, or (at your option) any later
 
12
! version.
 
13
!
 
14
! This program is distributed in the hope that it will be useful, but WITHOUT
 
15
! ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
 
16
! FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
 
17
! details.
 
18
!
 
19
! You should have received a copy of the GNU General Public License along with
 
20
! this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
 
21
! Street, Fifth Floor, Boston, MA 02110-1301, USA.
28
22
 
29
23
!-------------------------------------------------------------------------------
30
24
 
94
88
!     mode: <-- input, --> output, <-> modifies data, --- work array
95
89
!===============================================================================
96
90
 
 
91
!===============================================================================
 
92
! Module files
 
93
!===============================================================================
 
94
 
 
95
use paramx
 
96
use dimens
 
97
use numvar
 
98
use entsor
 
99
use optcal
 
100
use cstphy
 
101
use parall
 
102
use period
 
103
use ppppar
 
104
use radiat
 
105
 
 
106
!===============================================================================
 
107
 
97
108
implicit none
98
109
 
99
 
!===============================================================================
100
 
! Common blocks
101
 
!===============================================================================
102
 
 
103
 
include "paramx.h"
104
 
include "dimens.h"
105
 
include "numvar.h"
106
 
include "entsor.h"
107
 
include "optcal.h"
108
 
include "cstphy.h"
109
 
include "parall.h"
110
 
include "period.h"
111
 
include "ppppar.h"
112
 
include "radiat.h"
113
 
 
114
 
!===============================================================================
115
 
 
116
 
integer          iphas, ipp
 
110
integer          irphas, ipp
117
111
character*2      num
118
112
 
119
113
integer       ipass
246
240
  ihisvr(ipp,1) = -1
247
241
 
248
242
 
249
 
  do iphas = 1,nphasc
 
243
  do irphas = 1,nrphas
250
244
 
251
 
    write(num,'(I1)') iphas
 
245
    write(num,'(I1)') irphas
252
246
 
253
247
    !--> TERME SOURCE RADIATIF (ANALYTIQUE/CONSERVATIF/SEMI-ANALYTIQUE)
254
248
 
255
 
    ipp = ipppro(ipproc(itsre(iphas)))
 
249
    ipp = ipppro(ipproc(itsre(irphas)))
256
250
    nomvar(ipp)   = 'Srad'//num
257
251
    ichrvr(ipp)   = 0
258
252
    ilisvr(ipp)   = 0
260
254
 
261
255
    !--> PART DE L'ABSORPTION DANS LE TERME SOURCE RADIATIF
262
256
 
263
 
    ipp = ipppro(ipproc(iabs(iphas)))
 
257
    ipp = ipppro(ipproc(iabs(irphas)))
264
258
    nomvar(ipp)   = 'Absorp'//num
265
259
    ichrvr(ipp)   = 0
266
260
    ilisvr(ipp)   = 0
268
262
 
269
263
    !--> PART DE L'EMISSION DANS LE TERME SOURCE RADIATIF
270
264
 
271
 
    ipp = ipppro(ipproc(iemi(iphas)))
 
265
    ipp = ipppro(ipproc(iemi(irphas)))
272
266
    nomvar(ipp)   = 'Emiss'//num
273
267
    ichrvr(ipp)   = 0
274
268
    ilisvr(ipp)   = 0
276
270
 
277
271
    !--> COEFFICIENT D'ABSORPTION DU MILIEU SEMI-TRANSPARENT
278
272
 
279
 
    ipp = ipppro(ipproc(icak(iphas)))
 
273
    ipp = ipppro(ipproc(icak(irphas)))
280
274
    nomvar(IPP)   = 'CoefAb_'//num
281
275
    ichrvr(ipp)   = 0
282
276
    ilisvr(ipp)   = 0
324
318
 
325
319
  !--> FLUX CONVECTIF AUX FACES FRONTIERES DE PAROIS
326
320
 
327
 
  nbrvaf(ifconp) = 'Flux_convectif'
 
321
  nbrvaf(ifconp) = 'Convective_flux'
328
322
  irayvf(ifconp) = 0
329
323
 
330
324
  !--> COEFFICIENT D'ECHANGE CONVECTIF AUX FACES FRONTIERES DE PAROIS