← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/base/vorin0.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/vorin0.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 vorin0 &
29
 
!================
30
 
 
31
 
 ( nfabor , irepvo )
32
 
 
33
 
!===============================================================================
34
 
!  FONCTION  :
35
 
!  ---------
36
 
 
37
 
! INITIALISATION PAR DEFAUT DES PARAMETRES POUR LA METHODE
38
 
!   DES VORTEX AVANT DE PASSER LA MAIN A L'UTILISATEUR
39
 
 
40
 
!-------------------------------------------------------------------------------
41
 
! Arguments
42
 
!__________________.____._____.________________________________________________.
43
 
! name             !type!mode ! role                                           !
44
 
!__________________!____!_____!________________________________________________!
45
 
! nfabor           ! i  ! <-- ! number of boundary faces                       !
46
 
! irepvo           ! te ! --> ! numero de l'entree associe a chaque            !
47
 
!     (nfabor)     !    !     ! face de bord (=0 si pas de vortex)             !
48
 
!__________________!____!_____!________________________________________________!
49
 
 
50
 
!     TYPE : E (ENTIER), R (REEL), A (ALPHANUMERIQUE), T (TABLEAU)
51
 
!            L (LOGIQUE)   .. ET TYPES COMPOSES (EX : TR TABLEAU REEL)
52
 
!     MODE : <-- donnee, --> resultat, <-> Donnee modifiee
53
 
!            --- tableau de travail
54
 
!===============================================================================
55
 
 
56
 
implicit none
57
 
 
58
 
!===============================================================================
59
 
! Common blocks
60
 
!===============================================================================
61
 
 
62
 
include "vortex.h"
63
 
 
64
 
!===============================================================================
65
 
 
66
 
! Arguments
67
 
 
68
 
integer          nfabor, irepvo(nfabor)
69
 
 
70
 
! Local variables
71
 
 
72
 
integer          ii, jj, ifac
73
 
 
74
 
!===============================================================================
75
 
 
76
 
 
77
 
!===============================================================================
78
 
! 1. PARAMETRES DE LA METHODE
79
 
!===============================================================================
80
 
 
81
 
! ---> Nombre d'entrees
82
 
 
83
 
nnent = -999
84
 
 
85
 
! ---> Nombre de vortex
86
 
 
87
 
do ii = 1, nentmx
88
 
  nvort(ii) = -999
89
 
enddo
90
 
 
91
 
! ---> Defintion du cas
92
 
 
93
 
do ii = 1, nentmx
94
 
  icas(ii)   = - 999
95
 
enddo
96
 
 
97
 
! ---> Tableau de reperage et repere local
98
 
 
99
 
do ifac = 1, nfabor
100
 
  irepvo(ifac) = 0
101
 
enddo
102
 
 
103
 
do ii = 1, nentmx
104
 
  do jj = 1, 3
105
 
    dir1(jj,ii) = 0.d0
106
 
    dir2(jj,ii) = 0.d0
107
 
    cen(jj,ii)  = 0.d0
108
 
  enddo
109
 
enddo
110
 
 
111
 
! ---> Conditions aux limites et grandeurs caracteristiques
112
 
 
113
 
do ii = 1, nentmx
114
 
  do jj = 1, 4
115
 
    iclvor(jj,ii) = -999
116
 
  enddo
117
 
  lly(ii) = -999.d0
118
 
  llz(ii) = -999.d0
119
 
  lld(ii) = -999.d0
120
 
enddo
121
 
 
122
 
! ---> Parametres physiques
123
 
 
124
 
do ii = 1, nentmx
125
 
  itlivo(ii) = -999
126
 
  tlimvo(ii) = -999.d0
127
 
  isgmvo(ii) = -999
128
 
  xsgmvo(ii) = -999.d0
129
 
  idepvo(ii) = -999
130
 
  ud(ii)     = 0.d0
131
 
enddo
132
 
 
133
 
! ---> Donnees utilisateurs
134
 
 
135
 
do ii = 1, nentmx
136
 
  WRITE(FICVOR(II),'(1A6,I2.2)') 'vordat',II
137
 
  udebit(ii) = 0.d0
138
 
  kdebit(ii) = -999.d0
139
 
  edebit(ii) = -999.d0
140
 
enddo
141
 
 
142
 
return
143
 
end subroutine