← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/fuel/furays.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/fuel/furays.f90
1
 
!-------------------------------------------------------------------------------
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
25
 
 
26
 
!-------------------------------------------------------------------------------
27
 
 
28
 
subroutine furays &
29
 
!================
30
 
 
31
 
  ( ivar   , ncelet , ncel   ,                                    &
32
 
    volume , rtpa   , propce , smbrs  , rovsdt )
33
 
 
34
 
 
35
 
!===============================================================================
36
 
!  FONCTION  :
37
 
!  ---------
38
 
 
39
 
! ROUTINE PHYSIQUE PARTICULIERE : COMBUSTION FIOUL LOURD
40
 
!   PRISE EN COMPTE DES TERMES SOURCES RADIATIFS
41
 
!   IMPLICITE ET EXPLICITE DANS L'EQUATION DES PARTICULES
42
 
!   DE LA CLASSE ICLA
43
 
 
44
 
!-------------------------------------------------------------------------------
45
 
! Arguments
46
 
!__________________.____._____.________________________________________________.
47
 
! name             !type!mode ! role                                           !
48
 
!__________________!____!_____!________________________________________________!
49
 
! ivar             ! e  ! <-- ! numero de la variable scalaire                 !
50
 
!                  !    !     !   energie (enthalpie h2) pour le               !
51
 
!                  !    !     !   charbon                                      !
52
 
! ncelet           ! i  ! <-- ! number of extended (real + ghost) cells        !
53
 
! ncel             ! i  ! <-- ! number of cells                                !
54
 
! volume(ncelet    ! tr ! <-- ! volume des cellules                            !
55
 
! rtp, rtpa        ! ra ! <-- ! calculated variables at cell centers           !
56
 
!  (ncelet, *)     !    !     !  (at current and previous time steps)          !
57
 
! propce(ncelet, *)! ra ! <-- ! physical properties at cell centers            !
58
 
! smbrs(ncelet     ! tr ! <-- ! second membre du systeme                       !
59
 
! rovsdt(ncelet    ! tr ! <-- ! diagonale du systeme                           !
60
 
!__________________!____!_____!________________________________________________!
61
 
 
62
 
!     TYPE : E (ENTIER), R (REEL), A (ALPHANUMERIQUE), T (TABLEAU)
63
 
!            L (LOGIQUE)   .. ET TYPES COMPOSES (EX : TR TABLEAU REEL)
64
 
!     MODE : <-- donnee, --> resultat, <-> Donnee modifiee
65
 
!            --- tableau de travail
66
 
!===============================================================================
67
 
 
68
 
implicit none
69
 
 
70
 
!===============================================================================
71
 
! Common blocks
72
 
!===============================================================================
73
 
 
74
 
include "paramx.h"
75
 
include "cstnum.h"
76
 
include "cstphy.h"
77
 
include "entsor.h"
78
 
include "numvar.h"
79
 
include "ppppar.h"
80
 
include "ppthch.h"
81
 
include "ppincl.h"
82
 
include "radiat.h"
83
 
 
84
 
!===============================================================================
85
 
 
86
 
! Arguments
87
 
 
88
 
integer          ivar , ncelet, ncel
89
 
 
90
 
double precision volume(ncelet)
91
 
double precision smbrs(ncelet)
92
 
double precision rovsdt(ncelet)
93
 
double precision rtpa(ncelet,*)
94
 
double precision propce(ncelet,*)
95
 
 
96
 
! Local variables
97
 
 
98
 
integer          iel , numcla , ipcl
99
 
 
100
 
!===============================================================================
101
 
 
102
 
!===============================================================================
103
 
! 1. RECHERCHE DE LA ZONE MEMOIRE (IPH) EN FONCTION DU NUMERO DE PHASE
104
 
!    COURANT IPHAS POUR TROUVER LES BONS TERMES SOURCES
105
 
!===============================================================================
106
 
 
107
 
numcla = ivar-isca(ihlf(1))+1
108
 
ipcl   = 1+numcla
109
 
 
110
 
!===============================================================================
111
 
! 2. PRISE EN COMPTE DES TERMES SOURCES RADIATIFS
112
 
!===============================================================================
113
 
 
114
 
 
115
 
do iel = 1,ncel
116
 
  propce(iel,ipproc(itsri(ipcl))) = max(-propce(iel,ipproc(itsri(ipcl))),zero)
117
 
enddo
118
 
 
119
 
do iel = 1,ncel
120
 
  if ( propce(iel,ipproc(iyfol(numcla))) .gt. epzero ) then
121
 
 
122
 
!--> PARTIE EXPLICITE
123
 
 
124
 
    smbrs(iel)  = smbrs(iel) +  propce(iel,ipproc(itsre(ipcl)))*volume(iel)  &
125
 
                               *rtpa(iel,isca(iyfol(numcla)))
126
 
 
127
 
!--> PARTIE IMPLICITE
128
 
 
129
 
    rovsdt(iel) = rovsdt(iel) + propce(iel,ipproc(itsri(ipcl)))*volume(iel)
130
 
  endif
131
 
 
132
 
enddo
133
 
 
134
 
 
135
 
!----
136
 
! FIN
137
 
!----
138
 
 
139
 
return
140
 
 
141
 
end subroutine