cosake.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\cosake.f
00002 !
00043                      SUBROUTINE COSAKE
00044 !                    *****************
00045 !
00046      &(KARMAN,CMU,C1,C2,SIGMAK,SIGMAE,ESTAR,SCHMIT,KMIN,KMAX,EMIN,EMAX)
00047 !
00048 !***********************************************************************
00049 ! TELEMAC2D   V6P1                                   21/08/2010
00050 !***********************************************************************
00051 !
00052 !
00053 !
00054 !
00055 !
00056 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00057 !| C1             |<--| CONSTANT IN K-EPSILON MODEL
00058 !| C2             |<--| CONSTANT IN K-EPSILON MODEL
00059 !| CMU            |<--| CONSTANT IN K-EPSILON MODEL
00060 !| EMAX           |<--| MAXIMUM EPSILON IF CLIPPING
00061 !| EMIN           |-->| MINIMUM EPSILON IF CLIPPING
00062 !| ESTAR          |-->| CONSTANT IN K-EPSILON MODEL
00063 !| KARMAN         |-->| VON KARMAN CONSTANT
00064 !| KMAX           |<--| MAXIMUM K IF CLIPPING
00065 !| KMIN           |<--| MINIMUM K IF CLIPPING
00066 !| SCHMIT         |-->| SCHMITT NUMBER
00067 !| SIGMAE         |<--| K-EPSILON CONSTANT FOR DIFFUSION OF EPSILON
00068 !| SIGMAK         |<--| K-EPSILON CONSTANT FOR DIFFUSION OF K
00069 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00070 !
00071       IMPLICIT NONE
00072       INTEGER LNG,LU
00073       COMMON/INFO/LNG,LU
00074 !
00075 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00076 !
00077       DOUBLE PRECISION, INTENT(OUT) :: KMIN,KMAX,EMIN,EMAX
00078       DOUBLE PRECISION, INTENT(OUT) :: KARMAN,CMU,C1,C2
00079       DOUBLE PRECISION, INTENT(OUT) :: SIGMAK,SIGMAE,ESTAR,SCHMIT
00080 !
00081 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00082 !
00083 ! KARMAN CONSTANT
00084 !
00085 !     UP TO VERSION 6.0 : 0.41, FROM VERSION 6.1 ON : 0.40
00086       KARMAN = 0.40D0
00087       CMU    = 0.09D0
00088       C1     = 1.44D0
00089       C2     = 1.92D0
00090       SIGMAK = 1.00D0
00091       SIGMAE = 1.30D0
00092       ESTAR  = 0.15D0
00093 !
00094 ! SCHMIDT NUMBER
00095 !
00096       SCHMIT = 0.50D0
00097 !
00098 ! RANGE OF VALUES USED TO CLIP K AND EPSILON
00099 !
00100       KMIN = 1.D-8
00101       EMIN = 1.D-8
00102       KMAX = 1.D10
00103       EMAX = 1.D10
00104 !
00105 !-----------------------------------------------------------------------
00106 !
00107       RETURN
00108       END

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