The TELEMAC-MASCARET system  trunk
coefro_gaia.f
Go to the documentation of this file.
1 ! **********************
2  SUBROUTINE coefro_gaia
3 ! **********************
4 !
5  &(cf,h,ks,npoin,karman)
6 !
7 !***********************************************************************
8 ! GAIA
9 !***********************************************************************
10 !
12 !
13 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
19 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
20 !
21  USE bief
22 !
24  IMPLICIT NONE
25 !
26 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
27 !
28  INTEGER, INTENT(IN):: NPOIN
29  DOUBLE PRECISION,INTENT(IN):: KARMAN
30 !
31  TYPE(bief_obj), INTENT(INOUT) :: CF
32  TYPE(bief_obj),INTENT(IN) :: KS,H
33 !
34 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
35 !
36  INTEGER N
37  DOUBLE PRECISION AUX
38  INTRINSIC max,log
39 !
40 !-----------------------------------------------------------------------
41 !
42 ! CONSTRUCTION OF THE FRICTION COEFFICIENT ACCORDING TO NIKURADSE LAW
43 !
44  DO n=1,npoin
45  IF(ks%R(n).LE.0.d0) THEN
46  WRITE(lu,*) 'FROTTEMENT NON DEFINI DANS COEFRO AU POINT ',n
47  CALL plante(1)
48  stop
49  ENDIF
50  ENDDO
51 !
52 ! AUX=30.D0/EXP(1.D0) =11.036D0
53  DO n=1,npoin
54  aux = max(1.001d0,h%R(n)*11.036d0/ks%R(n))
55  cf%R(n) = 2.d0 / (log( aux)/karman )**2
56  ENDDO
57 !
58 !-----------------------------------------------------------------------
59 !
60  RETURN
61  END
subroutine coefro_gaia(CF, H, KS, NPOIN, KARMAN)
Definition: coefro_gaia.f:7
Definition: bief.f:3