← 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_varpos.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_varpos.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_varpos
 
24
!========================
 
25
!===============================================================================
 
26
!  FONCTION  :
 
27
!  ---------
 
28
 
 
29
!       INIT DES POSITIONS DES VARIABLES TRANSPORTEES POUR
 
30
!                COMBUSTION FUEL
 
31
 
 
32
!-------------------------------------------------------------------------------
 
33
! Arguments
 
34
!__________________.____._____.________________________________________________.
 
35
! name             !type!mode ! role                                           !
 
36
!__________________!____!_____!________________________________________________!
 
37
!__________________!____!_____!________________________________________________!
 
38
 
 
39
!     TYPE : E (ENTIER), R (REEL), A (ALPHANUMERIQUE), T (TABLEAU)
 
40
!            L (LOGIQUE)   .. ET TYPES COMPOSES (EX : TR TABLEAU REEL)
 
41
!     MODE : <-- donnee, --> resultat, <-> Donnee modifiee
 
42
!            --- tableau de travail
 
43
!===============================================================================
 
44
 
 
45
!===============================================================================
 
46
! Module files
 
47
!===============================================================================
 
48
 
 
49
use paramx
 
50
use dimens
 
51
use numvar
 
52
use optcal
 
53
use cstphy
 
54
use entsor
 
55
use cstnum
 
56
use ppppar
 
57
use ppthch
 
58
use coincl
 
59
use cpincl
 
60
use cs_fuel_incl
 
61
use ppincl
 
62
use ppcpfu
 
63
 
 
64
!===============================================================================
 
65
 
 
66
implicit none
 
67
 
 
68
integer        is, isc, icla
 
69
 
 
70
!===============================================================================
 
71
! 1. DEFINITION DES POINTEURS
 
72
!===============================================================================
 
73
 
 
74
! ---> Variables propres a la suspension gaz - particules
 
75
 
 
76
is    = 1
 
77
ihm   = iscapp(is)
 
78
 
 
79
! ---> Variables propres a la phase dispersee
 
80
 
 
81
do icla = 1, nclafu
 
82
  is          = 1+icla
 
83
  ing(icla)   = iscapp(is)
 
84
  is          = 1+1*nclafu+icla
 
85
  iyfol(icla) = iscapp(is)
 
86
  is          = 1+2*nclafu+icla
 
87
  ih2(icla)   = iscapp(is)
 
88
enddo
 
89
 
 
90
! ---> Variables propres a la phase gaz
 
91
 
 
92
is       = is + 1
 
93
ifvap    = iscapp(is)
 
94
! Oxydant 2
 
95
if ( noxyd .ge. 2 ) then
 
96
  is = is+1
 
97
  if4m  = iscapp(is)
 
98
endif
 
99
! Oxydant 3
 
100
if ( noxyd .ge. 3 ) then
 
101
  is = is+1
 
102
  if5m  = iscapp(is)
 
103
endif
 
104
! combustion heterogene
 
105
is       = is + 1
 
106
if7m   = iscapp(is)
 
107
! Variance
 
108
is       = is + 1
 
109
ifvp2m  = iscapp(is)
 
110
! Transport du CO2
 
111
if ( ieqco2 .ge. 1 ) then
 
112
  is = is+1
 
113
  iyco2 = iscapp(is)
 
114
endif
 
115
!  Transport du NOx : HCN, NOx et Hox
 
116
if ( ieqnox .eq. 1 ) then
 
117
  is = is+1
 
118
  iyhcn = iscapp(is)
 
119
  is = is+1
 
120
  iyno  = iscapp(is)
 
121
  is = is+1
 
122
  ihox = iscapp(is)
 
123
endif
 
124
 
 
125
!===============================================================================
 
126
! 2. PROPRIETES PHYSIQUES
 
127
!    A RENSEIGNER OBLIGATOIREMENT (sinon pb dans varpos)
 
128
!    - PROPRES AUX SCALAIRES   : IVISLS, ISCAVR
 
129
!      Rq : pas de variance associee a un scalaire dans notre cas
 
130
!    - PROPRES A LA SUSPENSION : ICP
 
131
!===============================================================================
 
132
 
 
133
do isc = 1, nscapp
 
134
 
 
135
  if ( iscavr(iscapp(isc)) .le. 0 ) then
 
136
 
 
137
! ---- Viscosite dynamique de reference relative au scalaire
 
138
!      ISCAPP(ISC)
 
139
    ivisls(iscapp(isc)) = 0
 
140
 
 
141
  endif
 
142
 
 
143
enddo
 
144
 
 
145
! ---- Bien que l on soit en enthalpie on conserve un CP constant
 
146
 
 
147
icp    = 0
 
148
 
 
149
!----
 
150
! End
 
151
!----
 
152
 
 
153
return
 
154
end subroutine