← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/cplv/cpini1.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/cplv/cpini1.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
 
49
44
!            --- tableau de travail
50
45
!===============================================================================
51
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 ppincl
 
63
use ihmpre
 
64
use ppcpfu
 
65
 
 
66
!===============================================================================
 
67
 
52
68
implicit none
53
69
 
54
 
!===============================================================================
55
 
! Common blocks
56
 
!===============================================================================
57
 
 
58
 
include "paramx.h"
59
 
include "dimens.h"
60
 
include "numvar.h"
61
 
include "optcal.h"
62
 
include "cstphy.h"
63
 
include "entsor.h"
64
 
include "cstnum.h"
65
 
include "ppppar.h"
66
 
include "ppthch.h"
67
 
include "coincl.h"
68
 
include "cpincl.h"
69
 
include "ppincl.h"
70
 
include "ihmpre.h"
71
 
include "ppcpfu.h"
72
 
 
73
 
!===============================================================================
74
 
 
75
70
integer          ipp , icla , ii , jj , iok
76
 
integer          icha , isc , is, iphas
 
71
integer          icha , isc , is
77
72
double precision wmolme
78
73
 
79
74
!===============================================================================
82
77
!===============================================================================
83
78
!     L'utilisateur ne doit pas y avoir touche.
84
79
 
85
 
do iphas = 1, nphas
86
 
  if(iscalt(iphas).ne.-1) then
87
 
    write(nfecra,1000)iphas,iscalt(iphas)
88
 
    call csexit (1)
89
 
    !==========
90
 
  endif
91
 
enddo
 
80
if(iscalt.ne.-1) then
 
81
  write(nfecra,1000)iscalt
 
82
  call csexit (1)
 
83
  !==========
 
84
endif
92
85
do ii = 1, nscapp
93
86
  if(iscsth(iscapp(ii)).ne.-10) then
94
87
    write(nfecra,1001)ii,iscapp(ii),iscapp(ii),iscsth(iscapp(ii))
188
181
 
189
182
! ---- On resout en enthalpie avec un CP constant (Cf. cpvarp)
190
183
 
191
 
iphas = iphsca(ihm)
192
 
iscalt(iphas) = ihm
 
184
iscalt = ihm
193
185
iscsth(ihm) = 2
194
186
 
195
187
! --> Donnees physiques ou numeriques propres aux scalaires CP
203
195
!        En combustion on considere que la viscosite turbulente domine
204
196
!        ON S'INTERDIT DONC LE CALCUL DES FLAMMES LAMINAIRES AVEC Le =/= 1
205
197
 
206
 
    visls0(jj) = viscl0(iphsca(jj))
 
198
    visls0(jj) = viscl0
207
199
 
208
200
  endif
209
201
 
556
548
          +wmole(ih2o)*oxyh2o(1)+wmole(ico2)*oxyco2(1))           &
557
549
        /(oxyo2(1)+oxyn2(1)+oxyh2o(1)+oxyco2(1))
558
550
 
559
 
ro0(iphas) = p0(iphas)*wmolme / (rr*t0(iphas))
 
551
ro0 = p0*wmolme / (rr*t0)
560
552
 
561
553
! ---- Initialisation pour la masse volumique du coke
562
554
 
580
572
diftl0 =-grand
581
573
 
582
574
! ---> Masse volumique variable et viscosite constante (pour les suites)
583
 
irovar(iphas) = 1
584
 
ivivar(iphas) = 0
 
575
irovar = 1
 
576
ivivar = 0
585
577
 
586
578
!===============================================================================
587
579
! 4. ON REDONNE LA MAIN A L'UTLISATEUR
630
622
'@                                                            ',/,&
631
623
'@  L''utilisateur ne doit pas la renseigner dans usini1, or  ',/,&
632
624
'@    elle a ete affectee comme suit :                        ',/,&
633
 
'@    ISCALT(',I10   ,') = ',I10                               ,/,&
 
625
'@    ISCALT = ',I10                                           ,/,&
634
626
'@                                                            ',/,&
635
627
'@  Le calcul ne sera pas execute.                            ',/,&
636
628
'@                                                            ',/,&