akepin.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\akepin.f
00002 !
00080                      SUBROUTINE AKEPIN
00081 !                    *****************
00082 !
00083      &(AK,EP,U,V,H,NPOIN,KFROT,CMU,C2,ESTAR,SCHMIT,KMIN,EMIN,CF)
00084 !
00085 !***********************************************************************
00086 ! TELEMAC2D   V6P1                                   21/08/2010
00087 !***********************************************************************
00088 !
00089 !
00090 !
00091 !
00092 !
00093 !
00094 !
00095 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00096 !| AK             |<--| TURBULENT KINETIC ENERGY
00097 !| C2             |-->| CONSTANT IN K-EPSILON MODEL
00098 !| CF             |-->| ADIMENSIONAL FRICTION COEFFICIENT
00099 !| CMU            |-->| CONSTANT IN K-EPSILON MODEL
00100 !| EMIN           |-->| MINIMUM EPSILON IF CLIPPING
00101 !| EP             |<--| TURBULENT DISSIPATION
00102 !| ESTAR          |-->| CONSTANT IN K-EPSILON MODEL
00103 !| H              |-->| WATER DEPTH
00104 !| KFROT          |-->| KEY-WORD: "LAW OF BOTTOM FRICTION"
00105 !|                |   | 1:CHEZY 2:LINEAIRE 3:STRICKLER
00106 !| KMIN           |-->| MINIMUM K IF CLIPPING
00107 !| NPOIN          |-->| NUMBER OF POINTS
00108 !| SCHMIT         |-->| SCHMITT NUMBER
00109 !| U              |-->| X-COMPONENT OF VELOCITY
00110 !| V              |-->| Y-COMPONENT OF VELOCITY
00111 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00112 !
00113       IMPLICIT NONE
00114       INTEGER LNG,LU
00115       COMMON/INFO/LNG,LU
00116 !
00117 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00118 !
00119       INTEGER, INTENT(IN) :: NPOIN,KFROT
00120       DOUBLE PRECISION, INTENT(INOUT) :: AK(NPOIN),EP(NPOIN)
00121       DOUBLE PRECISION, INTENT(IN) :: KMIN,EMIN,CMU,C2,ESTAR,SCHMIT
00122       DOUBLE PRECISION, INTENT(IN) :: U(NPOIN),V(NPOIN),H(NPOIN),CF(*)
00123 !
00124 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00125 !
00126       INTEGER K
00127 !
00128       DOUBLE PRECISION TIERS,HAUT,USTAR,CEPS
00129 !
00130       INTRINSIC SQRT,MAX
00131 !
00132 !-----------------------------------------------------------------------
00133 !
00134       TIERS = 1.D0/3.D0
00135 !
00136 !     INITIALISATION OF K AND EPSILON
00137 !
00138 !     *******************
00139       IF(KFROT.EQ.0) THEN
00140 !     *******************
00141 !
00142         IF(LNG.EQ.1) WRITE(LU,100)
00143         IF(LNG.EQ.2) WRITE(LU,101)
00144 100     FORMAT(1X,'AKEPIN N''EST PAS PREVU SANS FROTTEMENT SUR LE FOND')
00145 101     FORMAT(1X,'AKEPIN IS NOT PROVIDED WITHOUT BOTTOM FRICTION')
00146         CALL PLANTE(1)
00147         STOP
00148 !
00149 !     ****
00150       ELSE
00151 !     ****
00152 !
00153         DO K=1,NPOIN
00154            HAUT  = MAX(H(K),1.D-4)
00155            USTAR = SQRT( 0.5D0 * CF(K) * ( U(K)**2 + V(K)**2 ) )
00156            CEPS  = C2*SQRT(CMU)/SQRT(ESTAR*SCHMIT)/(0.5D0*CF(K))**0.75D0
00157            AK(K) = C2*USTAR**2/(0.5D0*CF(K)*CEPS)
00158            EP(K) = MAX( USTAR**3/(HAUT*SQRT(0.5D0*CF(K))) , EMIN )
00159         ENDDO
00160 !
00161 !     *****
00162       ENDIF
00163 !     *****
00164 !
00165 !-----------------------------------------------------------------------
00166 !
00167       RETURN
00168       END

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