The TELEMAC-MASCARET system  trunk
ustar2.f
Go to the documentation of this file.
1 ! *****************
2  SUBROUTINE ustar2
3 ! *****************
4 !
5  &( ustar, npoin2)
6 !
7 !***********************************************************************
8 ! TOMAWAC V6P1 29/06/2011
9 !***********************************************************************
10 !
11 !brief COMPUTES THE FRICTION VELOCITY U*
12 !+ FOR ALL THE NODES IN THE 2D MESH.
13 !+
14 !+ USES A DRAG COEFFICIENT WHICH VARIES LINEARLY WITH
15 !+ WIND SPEED. THE FORMULATION IS IDENTICAL TO THAT
16 !+ USED IN WAM-CYCLE 3 (WAMDI GROUP, 1988).
17 !
18 !reference WAMDI GROUP (1988) :
19 !+ "A THIRD GENERATION OCEAN WAVE PREDICTION MODEL".
20 !+ JPO, VOL 18, PP 1775-1810.
21 !
22 !history M. BENOIT (EDF/DER/LNH)
23 !+ 26/03/96
24 !+ V1P1
25 !+
26 !
27 !history N.DURAND (HRW), S.E.BOURBAN (HRW)
28 !+ 13/07/2010
29 !+ V6P0
30 !+ Translation of French comments within the FORTRAN sources into
31 !+ English comments
32 !
33 !history N.DURAND (HRW), S.E.BOURBAN (HRW)
34 !+ 21/08/2010
35 !+ V6P0
36 !+ Creation of DOXYGEN tags for automated documentation and
37 !+ cross-referencing of the FORTRAN sources
38 !
39 !history G.MATTAROLO (EDF - LNHE)
40 !+ 29/06/2011
41 !+ V6P1
42 !+ Translation of French names of the variables in argument
43 !
44 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
45 !| NPOIN2 |-->| NUMBER OF POINTS IN 2D MESH
46 !| USTAR |<--| FRICTION VELOCITY
47 !| UV |-->| WIND VELOCITY ALONG X AT THE MESH POINTS
48 !| VV |-->| WIND VELOCITY ALONG Y AT THE MESH POINTS
49 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
50 !
51  USE declarations_tomawac, ONLY : uv, vv
52  USE interface_tomawac, ex_ustar2 => ustar2
53  IMPLICIT NONE
54 !
55 !.....VARIABLES IN ARGUMENT
56 ! """"""""""""""""""""
57  INTEGER, INTENT(IN) :: NPOIN2
58  DOUBLE PRECISION, INTENT(INOUT) :: USTAR(npoin2)
59 !
60 !.....LOCAL VARIABLES
61 ! """""""""""""""""
62  INTEGER IP
63  DOUBLE PRECISION UVENT , CDRAG
64 !
65 !
66 !.....MAIN LOOP ON THE NODES OF THE 2D MESH
67 ! """""""""""""""""""""""""""""""""""""""""""""""""""""
68  DO ip=1,npoin2
69 !
70 !.......COMPUTES THE WIND SPEED 10 METERS ABOVE WATER
71 ! """""""""""""""""""""""""""""""""""""""""
72  uvent=sqrt(uv(ip)**2+vv(ip)**2)
73 !
74 !.......COMPUTES THE DRAG COEFFICIENT
75 ! """""""""""""""""""""""""""""""""
76  cdrag = 6.5d-5*uvent + 8.d-4
77  IF (uvent.LT.7.5d0) cdrag = 1.2875d-3
78 !
79 !.......COMPUTES THE FRICTION VELOCITY
80 ! """""""""""""""""""""""""""""""""""
81  ustar(ip)=sqrt(cdrag)*uvent
82 !
83  ENDDO ! IP
84 !
85  RETURN
86  END
double precision, dimension(:), pointer uv
subroutine ustar2(USTAR, NPOIN2)
Definition: ustar2.f:7
double precision, dimension(:), pointer vv