kepini.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac3d\kepini.f
00002 !
00051                      SUBROUTINE KEPINI
00052 !                    *****************
00053 !
00054      &(AK,EP,U,V,Z,ZF,NPOIN2,NPLAN,DNUVIH,DNUVIV,KARMAN,CMU,KMIN,EMIN)
00055 !
00056 !***********************************************************************
00057 ! TELEMAC3D   V6P1                                   21/08/2010
00058 !***********************************************************************
00059 !
00060 !
00061 !
00062 !
00063 !
00064 !
00065 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00066 !| AK             |<->| TURBULENT ENERGY
00067 !| CMU            |-->| CONSTANT FOR MODELE K-EPSILON MODEL
00068 !| DNUVIH         |-->| COEFFICIENT FOR HORIZONTAL DIFFUSION OF VELOCITIES
00069 !| DNUVIV         |-->| COEFFICIENT FOR VERTICAL DIFFUSION OF VELOCITIES
00070 !| EMIN           |-->| MINIMUM VALUE FOR EPSILON WHEN CLIPPING
00071 !| EP             |<->| TURBULENT DISSIPATION
00072 !| KARMAN         |-->| KARMAN CONSTANT
00073 !| KMIN           |-->| MINIMUM VALUE FOR K WHEN CLIPPING
00074 !| NPLAN          |-->| NUMBER OF PLANES IN THE 3D MESH OF PRISMS
00075 !| NPOIN2         |-->| NUMBER OF POINTS IN 2D
00076 !| U              |-->| COMPONENT OF VELOCITY
00077 !| V              |-->| COMPONENT OF VELOCITY
00078 !| Z              |-->| ELEVATION OF REAL 3D MESH POINTS
00079 !| ZF             |-->| ELEVATION OF BOTTOM
00080 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00081 !
00082       IMPLICIT NONE
00083       INTEGER LNG, LU
00084       COMMON/INFO/ LNG, LU
00085 !
00086 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00087 !
00088       INTEGER, INTENT(IN)            :: NPOIN2,NPLAN
00089       DOUBLE PRECISION, INTENT(INOUT):: AK(NPOIN2,NPLAN)
00090       DOUBLE PRECISION, INTENT(INOUT):: EP(NPOIN2,NPLAN)
00091       DOUBLE PRECISION, INTENT(IN)   :: U(NPOIN2,NPLAN), V(NPOIN2,NPLAN)
00092       DOUBLE PRECISION, INTENT(IN)   :: Z(NPOIN2,NPLAN)
00093       DOUBLE PRECISION, INTENT(IN)   :: ZF(NPOIN2)
00094       DOUBLE PRECISION, INTENT(IN)   :: KARMAN, DNUVIH, DNUVIV
00095       DOUBLE PRECISION, INTENT(IN)   :: CMU
00096       DOUBLE PRECISION, INTENT(IN)   :: KMIN, EMIN
00097 !
00098 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00099 !
00100       INTEGER IPOIN2,IPLAN
00101 !
00102       INTRINSIC LOG, SQRT, MAX
00103 !
00104       DOUBLE PRECISION, PARAMETER :: FICTIFEPS = 2.D0
00105 !
00106 !-----------------------------------------------------------------------
00107 !
00108 !     A THEORY BY VINCENT BOYER MODIFIED BY MARTIN FERRAND
00109 !
00110 !     DO IPOIN2 = 1,NPOIN2
00111 !       DO IPLAN = 1,NPLAN
00112 !
00113 !         ARBITRARY COMPUTATION OF K EXPRESSED AS A PERCENTAGE OF SPEED
00114 !
00115 !         AK(IPOIN2,IPLAN) = 1.D-3*U(IPOIN2,IPLAN)**2
00116 !         AK(IPOIN2,IPLAN) = MAX(AK(IPOIN2,IPLAN),KMIN)
00117 !
00118 !         COMPUTATION OF EPSILON
00119 !
00120 !         EP INITIALISED ACCORDING TO UETOIL**3/KAPPA/Y
00121 !         WHERE UETOIL IS CALCULATED FROM THE VALUE OF K AT THE WALL
00122 !
00123 !         IF(IPLAN.EQ.1) THEN
00124 !           DIST = (Z(IPOIN2,2)-ZF(IPOIN2))/FICTIFEPS
00125 !         ELSE
00126 !           DIST = Z(IPOIN2,IPLAN)-ZF(IPOIN2)
00127 !         ENDIF
00128 !         EP(IPOIN2,IPLAN)=CMU**0.75*SQRT(AK(IPOIN2,1)**3)/KARMAN/DIST
00129 !         EP(IPOIN2,IPLAN)=MAX(EP(IPOIN2,IPLAN),EMIN)
00130 !       ENDDO
00131 !     ENDDO
00132 !
00133 !-----------------------------------------------------------------------
00134 !
00135 !     HERE: NO INITIAL TURBULENCE
00136 !
00137       DO IPOIN2 = 1,NPOIN2
00138         DO IPLAN = 1,NPLAN
00139           AK(IPOIN2,IPLAN) = KMIN
00140           EP(IPOIN2,IPLAN) = EMIN
00141         ENDDO
00142       ENDDO
00143 !
00144 !-----------------------------------------------------------------------
00145 !
00146       RETURN
00147       END

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