ustar2.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\tomawac\ustar2.f
00002 !
00052                      SUBROUTINE USTAR2
00053 !                    *****************
00054 !
00055      &( USTAR , UV    , VV    , NPOIN2)
00056 !
00057 !***********************************************************************
00058 ! TOMAWAC   V6P1                                   29/06/2011
00059 !***********************************************************************
00060 !
00061 !
00062 !reference  WAMDI GROUP (1988) :
00063 !+                     "A THIRD GENERATION OCEAN WAVE PREDICTION MODEL".
00064 !+                      JPO, VOL 18, PP 1775-1810.
00065 !
00066 !
00067 !
00068 !
00069 !
00070 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00071 !| NPOIN2         |-->| NUMBER OF POINTS IN 2D MESH
00072 !| USTAR          |<--| FRICTION VELOCITY
00073 !| UV             |-->| WIND VELOCITY ALONG X AT THE MESH POINTS
00074 !| VV             |-->| WIND VELOCITY ALONG Y AT THE MESH POINTS
00075 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00076 !
00077       IMPLICIT NONE
00078 !
00079 !.....VARIABLES IN ARGUMENT
00080 !     """"""""""""""""""""
00081       INTEGER  NPOIN2
00082       DOUBLE PRECISION USTAR(NPOIN2) , UV(NPOIN2) , VV(NPOIN2)
00083 !
00084 !.....LOCAL VARIABLES
00085 !     """""""""""""""""
00086       INTEGER  IP
00087       DOUBLE PRECISION UVENT , CDRAG
00088 !
00089 !
00090 !.....MAIN LOOP ON THE NODES OF THE 2D MESH
00091 !     """""""""""""""""""""""""""""""""""""""""""""""""""""
00092       DO IP=1,NPOIN2
00093 !
00094 !.......COMPUTES THE WIND SPEED 10 METERS ABOVE WATER
00095 !       """""""""""""""""""""""""""""""""""""""""
00096         UVENT=DSQRT(UV(IP)**2+VV(IP)**2)
00097 !
00098 !.......COMPUTES THE DRAG COEFFICIENT
00099 !       """""""""""""""""""""""""""""""""
00100         CDRAG = 6.5D-5*UVENT + 8.D-4
00101         IF (UVENT.LT.7.5D0) CDRAG = 1.2875D-3
00102 !
00103 !.......COMPUTES THE FRICTION VELOCITY
00104 !       """""""""""""""""""""""""""""""""""
00105         USTAR(IP)=DSQRT(CDRAG)*UVENT
00106 !
00107       ENDDO ! IP
00108 !
00109       RETURN
00110       END

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