cstkep.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac3d\cstkep.f
00002 !
00095                      SUBROUTINE CSTKEP
00096 !                    *****************
00097 !
00098 !***********************************************************************
00099 ! TELEMAC3D   V7P0                                   21/08/2010
00100 !***********************************************************************
00101 !
00102 !
00103 !
00104 !
00105 !
00106 !
00107 !
00108 !
00109 !
00110 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00111 !| ALPHA          |<->| K-OMEGA CONSTANT
00112 !| BETA           |<->| K-OMEGA CONSTANT
00113 !| BETAS          |<->| K-OMEGA CONSTANT
00114 !| C1             |<->| K-EPSILON CONSTANT
00115 !| C2             |<->| K-EPSILON CONSTANT
00116 !| CMU            |<->| K-EPSILON CONSTANT
00117 !| EMAX           |<->| EPSILON MAXIMUM
00118 !| EMIN           |<->| EPSILON MINIMUM
00119 !| ITURBV         |-->| TURBULENCE MODEL (3:K-EPSILON 7:K-OMEGA)
00120 !| KARMAN         |<->| VON KARMAN CONSTANT
00121 !| KMAX           |<->| K MAXIMUM
00122 !| KMIN           |<->| K MINIMUM
00123 !| OMSTAR         |<->| K-OMEGA CONSTANT
00124 !| PRANDTL        |<->| PRANDTL NUMBER
00125 !| SCHMIT         |<->| SCHMIT NUMBER
00126 !| SIGMAE         |<->| K-EPSILON OR K-OMEGA CONSTANT
00127 !| SIGMAK         |<->| K-EPSILON OR K-OMEGA CONSTANT
00128 !| VIRT           |<->| VIRTUAL ORIGIN FOR EPSILON
00129 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00130 !
00131       USE DECLARATIONS_TELEMAC3D, ONLY : CMU,C1,C2,SIGMAK,SIGMAE,
00132      &                                   VIRT,SCHMIT,KMIN,KMAX,
00133      &                                   EMIN,EMAX,ALPHA,
00134      &                                   BETA,BETAS,OMSTAR,ITURBV,
00135      &                                   CLIPK,CLIPE,WSIK,YAP,
00136      &                                   PERNORM2,PERPROD,RIMIN,RIMAX,
00137      &                                   OPTPROD,LIMKF,LIMEF,LIMKS,LIMES
00138       IMPLICIT NONE
00139       INTEGER LNG,LU
00140       COMMON/INFO/LNG,LU
00141 !
00142 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00143 !
00144 !     INTEGER,          INTENT(IN )   :: ITURBV
00145 !     DOUBLE PRECISION, INTENT(INOUT) :: KMIN,KMAX,EMIN,EMAX
00146 !     DOUBLE PRECISION, INTENT(INOUT) :: KARMAN,CMU,C1,C2,SIGMAK,SIGMAE
00147 !     DOUBLE PRECISION, INTENT(INOUT) :: VIRT,PRANDTL,SCHMIT
00148 !     DOUBLE PRECISION, INTENT(INOUT) :: ALPHA,BETA,BETAS,OMSTAR
00149 !
00150 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00151 !
00152 !
00153 !
00154 !=======================================================================
00155 !
00156 !     ALL MODELS
00157 !
00158 !=======================================================================
00159 !
00160 !     VON KARMAN CONSTANT
00161 !
00162 !     UP TO VERSION 6.0, 0.41  FROM NOW ON : 0.40
00163 !     FROM 7.0 USE KEYWORDS
00164 !      KARMAN = 0.40D0
00165 !
00166 !     SCHMIDT NUMBER (not used)
00167 !
00168 !      SCHMIT = 1.D0
00169 !
00170 !     PRANDTL NUMBER (BETWEEN 0.8 AND 0.9 FOR TEMPERATURE)
00171 !
00172 !###>TBE - Prandtl number should be 1 by default for sediment
00173 !cc      PRANDTL = 0.71D0
00174 !      PRANDTL = 1.D0
00175 !###<TBE
00176 !
00177 !     K-EPSILON OR K-OMEGA MODEL
00178 !
00179       IF(ITURBV.EQ.3) THEN
00180         SIGMAK = 1.D0
00181         SIGMAE = 1.3D0
00182       ELSEIF(ITURBV.EQ.7) THEN
00183         SIGMAK = 2.D0
00184         SIGMAE = 2.D0
00185       ENDIF
00186 !
00187 !=======================================================================
00188 !
00189 !     K-EPSILON MODEL
00190 !
00191 !=======================================================================
00192 !
00193       CMU    = 0.09D0
00194       C1     = 1.44D0
00195       C2     = 1.92D0
00196 !
00197 !-----------------------------------------------------------------------
00198 !
00199 !     BOUNDARY CONDITIONS AT BOTTOM AND FREE SURFACE
00200 !
00201 !-----------------------------------------------------------------------
00202 !
00203 !     K : K
00204 !     E : EPSILON
00205 !     F : BOTTOM
00206 !     S : FREE SURFACE
00207 !     1 : NEUMANN
00208 !     2 : DIRICHLET (VALUES COMPUTED IN KEPCL3)
00209 !
00210       LIMKF = 1
00211       LIMEF = 1
00212       LIMKS = 1
00213       LIMES = 1
00214 !
00215 !-----------------------------------------------------------------------
00216 !
00217 !     PARAMETERS USED IN SUBROUTINE SOUKEP
00218 !
00219 !-----------------------------------------------------------------------
00220 !
00221 !     LIMITATION OF K AND EPSILON WITH PHYSICAL CRITERIA
00222 !
00223       CLIPK    = .TRUE.
00224       CLIPE    = .TRUE.
00225       PERNORM2 = 0.5D0
00226       PERPROD  = 0.1D0
00227 !
00228 !     MIN AND MAX OF RICHARDSON NUMBER
00229 !
00230       RIMIN=0.D0
00231       RIMAX=100.D0
00232 !
00233 !     OPTION FOR PRODUCTION
00234 !     1: LINEAR (NOT STANDARD)
00235 !     2: QUADRATIC (STANDARD)
00236 !
00237       OPTPROD=2
00238 !
00239 !     WIND STRESS IN K
00240 !
00241       WSIK=.TRUE.
00242 !
00243 !     YAP CORRECTION
00244 !
00245       YAP=.FALSE.
00246 !
00247 !=======================================================================
00248 !
00249 !     K-OMEGA MODEL
00250 !
00251 !=======================================================================
00252 !
00253       ALPHA  = 5.D0/9.D0
00254       BETA   = 3.D0/40.D0
00255       BETAS  = 0.09D0
00256 !
00257 !     TO COMPUTE THE FREE SURFACE VALUE OF OMEGA
00258 !
00259       OMSTAR  = 100.D0
00260 !
00261 !     VIRTUAL ORIGIN FOR EPSILON
00262 !
00263       VIRT = 0.07D0
00264 !
00265 !     MINIMA AND MAXIMA FOR CLIPPING
00266 !
00267       IF(ITURBV.EQ.3) THEN
00268         KMIN = 1.D-10
00269         EMIN = 1.D-16
00270         KMAX = 1.D4
00271         EMAX = 1.D10
00272       ELSEIF(ITURBV.EQ.7) THEN
00273         KMIN = 1.D-8
00274         EMIN = 1.D-3
00275         KMAX = 1.D-1
00276         EMAX = 1.D4
00277       ENDIF
00278 !
00279 !-----------------------------------------------------------------------
00280 !
00281       RETURN
00282       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0