← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/base/vert1d.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/base/vert1d.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 vert1d &
29
24
!================
30
25
 
31
 
 ( idbia0 , idbra0 ,                                              &
32
 
   nfabor , nfpt1d , iappel ,                                     &
33
 
   ifpt1d , nppt1d , iclt1d , ia     ,                            &
 
26
 ( nfabor , nfpt1d , iappel ,                                     &
 
27
   ifpt1d , nppt1d , iclt1d ,                                     &
34
28
   rgpt1d , eppt1d ,                                              &
35
 
   xlmbt1 , rcpt1d , dtpt1d , ra     )
 
29
   xlmbt1 , rcpt1d , dtpt1d )
36
30
 
37
31
!===============================================================================
38
32
! FONCTION :
57
51
!__________________.____._____.________________________________________________.
58
52
! name             !type!mode ! role                                           !
59
53
!__________________!____!_____!________________________________________________!
60
 
! idbia0           ! i  ! <-- ! number of first free position in ia            !
61
 
! idbra0           ! i  ! <-- ! number of first free position in ra            !
62
54
! ncelet           ! i  ! <-- ! number of extended (real + ghost) cells        !
63
 
! ncel             ! i  ! <-- ! number of cells                                !
64
 
! nfac             ! i  ! <-- ! number of interior faces                       !
65
55
! nfabor           ! i  ! <-- ! number of boundary faces                       !
66
 
! nfml             ! i  ! <-- ! number of families (group classes)             !
67
 
! nprfml           ! i  ! <-- ! number of properties per family (group class)  !
68
56
! nvar             ! i  ! <-- ! total number of variables                      !
69
57
! nscal            ! i  ! <-- ! total number of scalars                        !
70
 
! nphas            ! i  ! <-- ! number of phases                               !
71
58
! nfpt1d           ! e  ! <-- ! nombre de faces avec module therm 1d           !
72
59
! ifpt1d           ! e  ! <-- ! numero de la face en traitement                !
73
60
!                  !    !     ! thermique en paroi                             !
79
66
! xlmbt1           ! r  ! <-- ! diffusivite thermique                          !
80
67
! rcpt1d           ! r  ! <-- ! rocp                                           !
81
68
! dtpt1d           ! tr ! <-- ! pas de temps                                   !
82
 
! nideve, nrdeve   ! i  ! <-- ! sizes of idevel and rdevel arrays              !
83
 
! nituse, nrtuse   ! i  ! <-- ! sizes of ituser and rtuser arrays              !
84
69
! iappel           ! e  ! <-- ! indique les donnes a renvoyer                  !
85
 
! ifabor(nfabor)   ! ia ! <-- ! boundary faces -> cells connectivity           !
86
 
! ifmfbr(nfabor)   ! ia ! <-- ! boundary face family numbers                   !
87
 
! iprfml           ! ia ! <-- ! property numbers per family                    !
88
 
!  (nfml, nprfml)  !    !     !                                                !
89
 
! idevel(nideve)   ! ia ! <-> ! integer work array for temporary development   !
90
 
! ituser(nituse)   ! ia ! <-> ! user-reserved integer work array               !
91
 
! ia(*)            ! ia ! --- ! main integer work array                        !
92
 
! rdevel(nrdeve)   ! ra ! <-> ! real work array for temporary development      !
93
 
! rtuser(nrtuse)   ! ra ! <-> ! user-reserved real work array                  !
94
 
! ra(*)            ! ra ! --- ! main real work array                           !
95
70
!__________________!____!_____!________________________________________________!
96
71
 
97
72
!     TYPE : E (ENTIER), R (REEL), A (ALPHANUMERIQUE), T (TABLEAU)
98
73
!            L (LOGIQUE)   .. ET TYPES COMPOSES (EX : TR TABLEAU REEL)
99
74
!     MODE : <-- donnee, --> resultat, <-> Donnee modifiee
100
75
!            --- tableau de travail
 
76
!===============================================================================
 
77
 
 
78
!===============================================================================
 
79
! Module files
 
80
!===============================================================================
 
81
 
 
82
use paramx
 
83
use numvar
 
84
use cstnum
 
85
use entsor
 
86
use optcal
101
87
 
102
88
!===============================================================================
103
89
 
104
90
implicit none
105
91
 
106
 
!===============================================================================
107
 
! Common blocks
108
 
!===============================================================================
109
 
 
110
 
include "paramx.h"
111
 
include "numvar.h"
112
 
include "cstnum.h"
113
 
include "entsor.h"
114
 
include "optcal.h"
115
 
 
116
 
!===============================================================================
117
 
 
118
92
! Arguments
119
93
 
120
 
integer          idbia0 , idbra0
121
94
integer          nfabor , nfpt1d
122
95
integer          iappel
123
96
 
124
97
integer          ifpt1d(nfpt1d) , nppt1d(nfpt1d) , iclt1d(nfpt1d)
125
 
integer          ia(*)
126
98
 
127
99
double precision eppt1d(nfpt1d) , rgpt1d(nfpt1d)
128
100
double precision xlmbt1(nfpt1d) , rcpt1d(nfpt1d) , dtpt1d(nfpt1d)
129
 
double precision ra(*)
130
101
 
131
102
! Variables locales
132
 
integer          idebia , idebra
133
103
integer          ii, ifac
134
104
 
135
105
!===============================================================================
136
106
 
137
107
 
138
 
idebia = idbia0
139
 
idebra = idbra0
140
108
 
141
109
 
142
110
if(iappel.eq.1) then