The TELEMAC-MASCARET system  trunk
tobw_gaia.f
Go to the documentation of this file.
1 ! ********************
2  SUBROUTINE tobw_gaia
3 ! ********************
4 !
5  &(tobw ,cf, fw, uw,tw,hn,npoin,xmve)
6 !
7 !***********************************************************************
8 ! GAIA
9 !***********************************************************************
10 !
13 !
14 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
23 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
24 !
26  IMPLICIT NONE
27 !
28 !
29 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
30 !
31  INTEGER, INTENT(IN) :: NPOIN
32 !
33  DOUBLE PRECISION, INTENT(IN) :: CF(npoin)
34  DOUBLE PRECISION, INTENT(IN) :: UW(npoin),TW(npoin),HN(npoin)
35  DOUBLE PRECISION, INTENT(IN) :: XMVE
36  DOUBLE PRECISION, INTENT(INOUT) :: TOBW(npoin),FW(npoin)
37 !
38 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
39 !
40  INTEGER I
41  DOUBLE PRECISION KS,AUX
42  DOUBLE PRECISION PI,AW
43  DOUBLE PRECISION, PARAMETER :: KARMAN = 0.4d0
44 !
45 !-----------------------------------------------------------------------
46 !
47  pi = 4.d0 * atan( 1.d0 )
48 !
49 !-----------------------------------------------------------------------
50 !
51  DO i=1,npoin
52 ! KS : NIKURADSE COEFFICIENT (TOTAL FRICTION)
53  aux=1.d0+karman*sqrt(2.d0/max(cf(i),1.d-10))
54  ks=30.d0*max(hn(i),1.d-8)*exp(-aux)
55  aw= uw(i)*tw(i) / (2.d0*pi)
56  IF(aw/ks.GT.1.59d0) THEN
57  fw(i)=exp( -6.d0 + 5.2d0 * (aw/ks)**(-0.19d0) )
58  ELSE
59  fw(i)=0.3d0
60  ENDIF
61  tobw(i)=0.5d0 * xmve * fw(i) * uw(i)*uw(i)
62  ENDDO
63 !
64 !-----------------------------------------------------------------------
65 !
66  RETURN
67  END SUBROUTINE tobw_gaia
subroutine tobw_gaia(TOBW, CF, FW, UW, TW, HN, NPOIN, XMVE)
Definition: tobw_gaia.f:7