1
!-------------------------------------------------------------------------------
3
! This file is part of the Code_Saturne Kernel, element of the
4
! Code_Saturne CFD tool.
6
! Copyright (C) 1998-2009 EDF S.A., France
8
! contact: saturne-support@edf.fr
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.
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.
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
26
!-------------------------------------------------------------------------------
29
!===============================================================================
31
! METHODE DES VORTEX POUR CONDITIONS AUX LIMITES D'ENTREE EN L.E.S.
36
integer nentmx, ndatmx
37
parameter (nentmx = 10)
38
parameter (ndatmx = 10000)
40
! NENTMX : NOMBRE D'ENTREE MAX
41
! NDATMX : NOMBRE DE POINTS MAX POUR LE FICHIER DES DONNEES
46
integer icvor(nentmx) , icvor2(nentmx) , &
49
common /idimvo/ icvor , icvor2 , icvmax , nvomax
51
! ICVOR : NOMBRE DE FACES (GLOBAL) UTILISANT DES VORTEX
53
! ICVOR2 : COMPTEUR DU NOMBRE LOCAL DE FACES UTILISANT DES VORTEX
54
! ICVMAX : NOMBRE MAX DE FACES UTILISANT DES VORTEX (SUR TOUTES ENTREES
56
! NVOMAX : NOMBRE MAX DE VORTEX UTILISE (TOUTES ENTREES CONFONDUES)
62
integer iirepv , iifagl , iivrce , &
64
iyzcel , iuvort , ivvort , iwvort , &
65
iyzvor , iyzvoa , isignv , ixsigm , &
66
ixgamm , ixtmp , ixtmpl , &
67
iw1x , iw1y , iw1z , iw1v , &
68
iw2x , iw2y , iw2z , iw2v
70
common /iivort/ iirepv , iifagl , iivrce , &
72
iyzcel , iuvort , ivvort , iwvort , &
73
iyzvor , iyzvoa , isignv , ixsigm , &
74
ixgamm , ixtmp , ixtmpl , &
75
iw1x , iw1y , iw1z , iw1v , &
76
iw2x , iw2y , iw2z , iw2v
78
! IIREPV : DEBUT DU TABLEAU ASSOCIANT AUX FACES DE BORD
79
! LE NUMERO D'UNE ENTREE
80
! IIFAGL : DEBUT DU TABLEAU DE CONNECTIVITE
81
! IIVRCE : DEBUT DU TABLEAU REPERANT LA CELLULE LA PLUS VOISINE
83
! IXYZV : DEBUT DU TABLEAUX CONTENANT LES COORDONNEES DE
84
! TOUTES LES FACES D'ENTREE
85
! IVISV : DEBUT DU TABLEAU CONTENANT LA VISCOSITE SUR
86
! TOUTES LES FACES D'ENTREE
87
! IYZCEL : DEBUT DU TABLEAU CONTENANT LES COORDONNEES DES
88
! FACES D'ENTREE DANS LE REPERE LOCAL
89
! IUVORT,...: DEBUTS DES TABLEAUX CONTENANT LES COMPOSANTES DE VITESSE
90
! IYZVOR : DEBUT DU TABLEAU CONTENANT LA POSITION DES VORTEX
91
! DANS LE REPERE LOCAL
92
! IYZVOA : DEBUT DU TABLEAU CONTENANT LA POSITION DES VORTEX
93
! DANS LE REPERE LOCAL AU PAS DE TEMPS PRECEDENT
94
! ISIGNV : DEBUT DU TABLEAU CONTENANT LE SENS DE ROTATION DES
96
! IXSIGM : DEBUT DU TABLEAU CONTENANT LA TAILLE DES VORTEX
97
! IXGAMM : DEBUT DU TABLEAU CONTENANT L'INTENSITE DES VORTEX
98
! IXTMP : DEBUT DU TABLEAU CONTENANT LE TEMPS CUMULE
99
! IXTMPL : DEBUT DU TABLEAU CONTENANT LE TEMPS DE VIE DES VORTEX
100
! IW1X,.. : DEBUT DES TABLEAUX DE TRAVAILS SERVANT A COMMUNIQUER
101
! LES DONNEES AUX ENTREES A TOUS LES PROCESSEURS
102
! (PLUS UTILISE APRES VORPRE)
108
integer nnent , nvort(nentmx) , &
110
icas(nentmx) , itlivo(nentmx) , &
111
isgmvo(nentmx) , idepvo(nentmx) , &
112
iclvor(4,nentmx), ndat(nentmx)
114
common /ioptvo/ nnent , nvort , initvo , &
115
icas , itlivo , isgmvo , idepvo , &
118
! NNENT : NOMBRE D ENTREES UTILISEES
119
! NVORT : NOMBRE DE VORTEX
120
! INITVO : INDICATEUR DE REINITIALISATION
121
! ICAS : TYPE DE GEOMETRIE POUR L'ENTREE
122
! ITLIVO : TYPE DE MODELE POUR LA DUREE DE VIE
123
! ISGMVO : TYPE DE MODELE POUR LA TAILLE DES VORTEX
124
! IDEPVO : TYPE DE MODELE POUR LA MARCHE EN TEMPS
125
! ICLVOR : TYPE DE CONDITION AUX LIMITES
126
! NDAT : NOMBRE DE LIGNES DU FICHIER DE DONNEES
132
double precision tlimvo(nentmx), xsgmvo(nentmx), ud(nentmx), &
133
xdat(ndatmx,nentmx), &
134
ydat(ndatmx,nentmx), zdat(ndatmx,nentmx), &
135
udat(ndatmx,nentmx), &
136
vdat(ndatmx,nentmx), wdat(ndatmx,nentmx), &
137
dudat(ndatmx,nentmx), &
138
kdat(ndatmx,nentmx), epsdat(ndatmx,nentmx), &
139
udebit(nentmx), kdebit(nentmx), edebit(nentmx), &
140
dir1(3,nentmx), dir2(3,nentmx), dir3(3,nentmx), &
141
cen(3,nentmx) , surf(3,nentmx), &
142
ymax(nentmx) , ymin(nentmx), &
143
zmax(nentmx) , zmin(nentmx), &
144
xsurfv(nentmx), llz(nentmx), &
145
lly(nentmx) , lld(nentmx)
147
common /roptvo/ tlimvo , xsgmvo , ud , &
148
xdat , ydat , zdat , &
149
udat , vdat , wdat , &
151
epsdat , udebit , kdebit , edebit , &
152
dir1 , dir2 , dir3 , cen , surf , &
153
ymax , ymin , zmax , zmin , &
154
xsurfv , llz , lly , lld
157
! TLIMVO : TEMPS DE VIE MAX DES VORTEX IMPOSE PAR L'UTILISATEUR
158
! XSGMVO : DIAMETRE DES VORTEX IMPOSE PAR L'UTILISATEUR
159
! UD : VITESSE DE DEPLACEMENT (MAX) IMPOSEE PAR L'UTILISATEUR
160
! XDAT, ... : COORDONNEES DES POINTS OU SONT CONNUES LES DONNEES
161
! UDAT : VITESSE MOYENNE PRINCIPALE (FICHIER DE DONNEES)
162
! VDAT,WDAT : VITESSE MOYENNE TRANSVERSE (FICHIER DE DONNEES)
163
! DUDAT : DERIVE NORMALE DE LA VITESSE PRINCIPALE (FICHIER D'ENTREE)
164
! KDAT : EC MOYENNE (FICHIER D'ENTREE)
165
! EPSDAT : DISSIPATION (FICHIER D'ENTREE)
166
! UDEBIT : VITESSE MOYENNE IMPOSEE PAR L'UTILISATEUR EN ENTREE
167
! KDEBIT : EC IMPOSEE PAR L'UTILISATEUR EN ENTREE
168
! EDEBIT : DISSIPATION IMPOSEE PAR L'UTILISATEUR EN ENTREE
169
! DIR1,... : VECTEURS DEFINISSANT LE REPERE LOCAL DANS LE PLAN D'ENTREE
170
! CEN : COORDONNEES DU CENTRE DE L'ENTREE
171
! SURF : VECTEUR SURFACE DU PLAN D'ENTREE (SUPPOSEE PLANE)
172
! XMAX,... : DIMENSIONS MAX DE L'ENTREE DANS LE REPERE LOCAL
173
! LLZ,LLY,LLD : DIMENSIONS DE L'ENTREE DANS LE CALCUL
175
character*50 ficvor(nentmx)
176
common /coptvo/ ficvor
178
! FICVOR : NOM DU FICHIER DE DONNEE