← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/turb/vorimp.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/turb/vorimp.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 vorimp &
 
24
!================
 
25
 
 
26
 ( ient   )
 
27
 
 
28
!===============================================================================
 
29
!  FONCTION  :
 
30
!  ---------
 
31
 
 
32
! IMPRESSION DES PARAMETRES DE LA METHODE DES VORTEX
 
33
! APRES INTERVENTION UTILISATEUR DANS usvort
 
34
!-------------------------------------------------------------------------------
 
35
! Arguments
 
36
!__________________.____._____.________________________________________________.
 
37
! name             !type!mode ! role                                           !
 
38
!__________________!____!_____!________________________________________________!
 
39
! ient             ! e  ! <-- ! numero de l'entree                             !
 
40
! nnent            ! e  ! <-- ! nombre d'entrees utilisant des vortex          !
 
41
! nvort(nentmx)    ! te ! ->  ! nombre de vortex a chaque entree               !
 
42
!__________________!____!_____!________________________________________________!
 
43
 
 
44
!     TYPE : E (ENTIER), R (REEL), A (ALPHANUMERIQUE), T (TABLEAU)
 
45
!            L (LOGIQUE)   .. ET TYPES COMPOSES (EX : TR TABLEAU REEL)
 
46
!     MODE : <-- donnee, --> resultat, <-> Donnee modifiee
 
47
!            --- tableau de travail
 
48
!===============================================================================
 
49
 
 
50
!===============================================================================
 
51
! Module files
 
52
!===============================================================================
 
53
 
 
54
use paramx
 
55
use entsor
 
56
use optcal
 
57
use vorinc
 
58
 
 
59
!===============================================================================
 
60
 
 
61
implicit none
 
62
 
 
63
! Arguments
 
64
 
 
65
integer          ient
 
66
 
 
67
! Local variables
 
68
 
 
69
integer          ipass
 
70
data             ipass /0/
 
71
save             ipass
 
72
 
 
73
!===============================================================================
 
74
! On ne peut pas faire ces impressions dans impini car la geometrie
 
75
!  des entrees n'est pas encore connue.
 
76
 
 
77
!===============================================================================
 
78
! 1. Methode vortex
 
79
!===============================================================================
 
80
 
 
81
ipass = ipass + 1
 
82
 
 
83
write(nfecra,9010) ient, nvort(ient), icvor(ient), icas(ient)
 
84
write(nfecra,9020) dir1(1,ient),dir1(2,ient),dir1(3,ient),        &
 
85
                   dir2(1,ient),dir2(2,ient),dir2(3,ient),        &
 
86
                   dir3(1,ient),dir3(2,ient),dir3(3,ient),        &
 
87
                   cen(1,ient) ,cen(2,ient) ,cen(3,ient)
 
88
if(icas(ient).eq.1) then
 
89
  write(nfecra,9030) iclvor(1,ient),iclvor(2,ient),               &
 
90
                     iclvor(3,ient),iclvor(4,ient)
 
91
endif
 
92
write(nfecra,9040) ymin(ient),ymax(ient),zmin(ient),zmax(ient)
 
93
if(icas(ient).eq.1) then
 
94
  write(nfecra,9050) lly(ient),llz(ient)
 
95
endif
 
96
if(icas(ient).eq.2) then
 
97
  write(nfecra,9060) lld(ient)
 
98
endif
 
99
write(nfecra,9070) itlivo(ient)
 
100
if(itlivo(ient).eq.1) then
 
101
  write(nfecra,9080) tlimvo(ient)
 
102
endif
 
103
write(nfecra,9090) isgmvo(ient)
 
104
if(isgmvo(ient).eq.1) then
 
105
  write(nfecra,9100) xsgmvo(ient)
 
106
endif
 
107
write(nfecra,9110) idepvo(ient)
 
108
if(idepvo(ient).eq.1) then
 
109
  write(nfecra,9120) ud(ient)
 
110
endif
 
111
if(icas(ient).eq.1.or.icas(ient).eq.2.or.icas(ient).eq.3) then
 
112
  write(nfecra,9130) ndat(ient)
 
113
elseif(icas(ient).eq.4) then
 
114
  write(nfecra,9140) udebit(ient),kdebit(ient),edebit(ient)
 
115
endif
 
116
 
 
117
#if defined(_CS_LANG_FR)
 
118
 
 
119
 9010 format(                                                           &
 
120
'-------------                                                ',/,&
 
121
'  -- Entree : ',I10                                            /,&
 
122
'-------------                                                ',/,&
 
123
'       NVORT  = ',4X,I10,    ' (Nombre de vortex            )',/,&
 
124
'       ICVOR  = ',4X,I10,   ' (Nombre de faces a l''entree  )',/,&
 
125
'       ICAS   = ',4X,I10,    ' (1 : conduite rectangulaire   ',/,&
 
126
'                ',14X,       '  2 : conduite circulaire      ',/,&
 
127
'                ',14X,       '  3 : sans CL mais avec fichier',/,&
 
128
'                ',14X,       '  4 : sans CL ni fichier      )',/)
 
129
 9020 format(                                                           &
 
130
' --- Directions principales du repere local                  ',/,&
 
131
'                  ---- X ----    ---- Y ----    ---- Z ----  ',/,&
 
132
'       DIR1  = ', E14.5,' ',E14.5,' ',E14.5                    /,&
 
133
'       DIR2  = ', E14.5,' ',E14.5,' ',E14.5                    /,&
 
134
'       DIR3  = ', E14.5,' ',E14.5,' ',E14.5                    /,&
 
135
'                                                             ',/,&
 
136
' --- Coordonnees du centre de l''entree                      ',/,&
 
137
'       CEN   = ', E14.5,' ',E14.5,' ',E14.5,                   /)
 
138
 9030 format(                                                           &
 
139
' --- Conditions aux limites dans le repere local             ',/,&
 
140
' Plan y = -LLY/2 ',4X,I10,    ' (1 : paroi                   ',/,&
 
141
' Plan z =  LLZ/2 ',4X,I10,    '  2 : symetrie                ',/,&
 
142
' Plan y =  LLY/2 ',4X,I10,    '  3 : periodicite            )',/,&
 
143
' Plan z = -LLZ/2 ',4X,I10,    '                              ',/)
 
144
 9040 format(                                                           &
 
145
' --- Dimensions de l''entree dans le repere local            ',/,&
 
146
'                  ---- min ----    ---- max ----             ',/,&
 
147
'       Y       = ',E14.5,'  ',E14.5,'                        ',/,&
 
148
'       Z       = ',E14.5,'  ',E14.5,'                        ',/)
 
149
 9050 format(                                                           &
 
150
'       LLY     = ',E14.5,    ' (longueur de la conduite dans ',/,&
 
151
'       LLZ     = ',E14.5,    '  les directions DIR1 et DIR2) ',/)
 
152
 9060 format(                                                           &
 
153
'       LLD     = ',E14.5,    ' (diametre de la conduite    ) ',/)
 
154
 9070 format(                                                           &
 
155
' --- Duree de vie des vortex                                 ',/,&
 
156
'       ITLIVO  = ',4X,I10,   ' (1 : constante                ',/,&
 
157
'                 ',14X,      '  2 : en k^(3/2).U/epsilon   ) ',/)
 
158
 9080 format(                                                           &
 
159
'       TLIMVO  = ',E14.5,    ' (1 : duree de vie imposee   ) ',/)
 
160
 9090 format(                                                           &
 
161
' --- Taille des vortex                                       ',/,&
 
162
'       ISGMVO  = ',4X,I10,   ' (1 : taille constante         ',/,&
 
163
'                 ',14X,      '  2 : en k^(3/2)/epsilon       ',/,&
 
164
'                 ',14X,      '  2 : en max[nu.k/eps,200.Lk]) ',/)
 
165
 9100 format(                                                           &
 
166
'       XSGMVO  = ',E14.5,    ' (1 : taille imposee         ) ',/)
 
167
 9110 format(                                                           &
 
168
' --- Marche en temps                                         ',/,&
 
169
'       IDEPVO  = ',4X,I10,   ' (1 : deplacement aleatoire    ',/,&
 
170
'                 ',14X,      '  2 : convection des vortex    ',/)
 
171
 9120 format(                                                           &
 
172
'       UD      = ',E14.5,    ' (1 : vit. de deplacement max) ',/)
 
173
 9130 format(                                                           &
 
174
' --- Fichier de donnees                                      ',/,&
 
175
'       NDAT    = ',4X,I10,   ' (Nombre de lignes du fichier )',/)
 
176
 9140 format(                                                           &
 
177
' --- Donnees a l''entree                                     ',/,&
 
178
'       UDEBIT  = ',E14.5,    ' (vitesse debitante imposee)   ',/,&
 
179
'       KDEBIT  = ',E14.5,    ' (energie cinetique imposee)   ',/,&
 
180
'       EDEBIT  = ',E14.5,    ' (dissipation imposee)         ',/)
 
181
 
 
182
#else
 
183
 
 
184
 9010 format(                                                           &
 
185
'-----------                                                  ',/,&
 
186
'  -- Inlet: ',I10                                              /,&
 
187
'-----------                                                  ',/,&
 
188
'       NVORT  = ',4X,I10,    ' (Number of vortices          )',/,&
 
189
'       ICVOR  = ',4X,I10,    ' (Number of faces at the inlet)',/,&
 
190
'       ICAS   = ',4X,I10,    ' (1 : rectangular duct         ',/,&
 
191
'                ',14X,       '  2 : pipe                     ',/,&
 
192
'                ',14X,       '  3 : wo BC but with file      ',/,&
 
193
'                ',14X,       '  4 : wo BC neither file      )',/)
 
194
 9020 format(                                                           &
 
195
' --- Main directions for the local frame                     ',/,&
 
196
'                  ---- X ----    ---- Y ----    ---- Z ----  ',/,&
 
197
'       DIR1  = ', E14.5,' ',E14.5,' ',E14.5                    /,&
 
198
'       DIR2  = ', E14.5,' ',E14.5,' ',E14.5                    /,&
 
199
'       DIR3  = ', E14.5,' ',E14.5,' ',E14.5                    /,&
 
200
'                                                             ',/,&
 
201
' --- Inlet center coordinates                                ',/,&
 
202
'       CEN   = ', E14.5,' ',E14.5,' ',E14.5,                   /)
 
203
 9030 format(                                                           &
 
204
' --- Boundary conditions in the local frame                  ',/,&
 
205
' Y plane = -LLY/2 ',4X,I10,    ' (1 : wall                   ',/,&
 
206
' Z plane =  LLZ/2 ',4X,I10,    '  2 : symmetry               ',/,&
 
207
' Y plane =  LLY/2 ',4X,I10,    '  3 : periodicity           )',/,&
 
208
' Z plane = -LLZ/2 ',4X,I10,    '                             ',/)
 
209
 9040 format(                                                           &
 
210
' --- Inlet dimensions in the local framae                    ',/,&
 
211
'                  ---- min ----    ---- max ----             ',/,&
 
212
'       Y       = ',E14.5,'  ',E14.5,'                        ',/,&
 
213
'       Z       = ',E14.5,'  ',E14.5,'                        ',/)
 
214
 9050 format(                                                           &
 
215
'       LLY     = ',E14.5,    ' (duct length in the directions',/,&
 
216
'       LLZ     = ',E14.5,    '  DIR1 and DIR2              ) ',/)
 
217
 9060 format(                                                           &
 
218
'       LLD     = ',E14.5,    ' (pipe diameter              ) ',/)
 
219
 9070 format(                                                           &
 
220
' --- Vortices life time                                      ',/,&
 
221
'       ITLIVO  = ',4X,I10,   ' (1 : constant                 ',/,&
 
222
'                 ',14X,      '  2 : in k^(3/2).U/epsilon   ) ',/)
 
223
 9080 format(                                                           &
 
224
'       TLIMVO  = ',E14.5,    ' (1 : given life time        ) ',/)
 
225
 9090 format(                                                           &
 
226
' --- Vortices size                                           ',/,&
 
227
'       ISGMVO  = ',4X,I10,   ' (1 : constant                 ',/,&
 
228
'                 ',14X,      '  2 : in k^(3/2)/epsilon       ',/,&
 
229
'                 ',14X,      '  2 : in max[nu.k/eps,200.Lk]) ',/)
 
230
 9100 format(                                                           &
 
231
'       XSGMVO  = ',E14.5,    ' (1 : given size             ) ',/)
 
232
 9110 format(                                                           &
 
233
' --- Time marching                                           ',/,&
 
234
'       IDEPVO  = ',4X,I10,   ' (1 : random displacement      ',/,&
 
235
'                 ',14X,      '  2 : vortices convection      ',/)
 
236
 9120 format(                                                           &
 
237
'       UD      = ',E14.5,    ' (1 : max displacement vel.  ) ',/)
 
238
 9130 format(                                                           &
 
239
' --- Data file                                               ',/,&
 
240
'       NDAT    = ',4X,I10,   ' (Number of lines in the file )',/)
 
241
 9140 format(                                                           &
 
242
' --- Inlet date                                              ',/,&
 
243
'       UDEBIT  = ',E14.5,    ' (given velocity      )        ',/,&
 
244
'       KDEBIT  = ',E14.5,    ' (given kinetic energy)        ',/,&
 
245
'       EDEBIT  = ',E14.5,    ' (given dissipation   )        ',/)
 
246
 
 
247
#endif
 
248
 
 
249
!===============================================================================
 
250
! 2. FIN
 
251
!===============================================================================
 
252
 
 
253
return
 
254
end subroutine