ustar1.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\tomawac\ustar1.f
00002 !
00059                      SUBROUTINE USTAR1
00060 !                    *****************
00061 !
00062      &( USTAR , Z0    , TAUWAV, UV    , VV    , CDRAG , ALPHA , XKAPPA,
00063      &  ZVENT , GRAVIT, NPOIN2)
00064 !
00065 !***********************************************************************
00066 ! TOMAWAC   V6P1                                   29/06/2011
00067 !***********************************************************************
00068 !
00069 !
00070 !
00071 !reference  JANSSEN P.A.E.M (1989) :
00072 !+                     "WIND-INDUCED STRESS AND THE DRAG OF AIR
00073 !+                      FLOW OVER SEA WAVES". JPO, VOL 19, PP 745-754.
00074 !reference JANSSEN P.A.E.M (1991) :
00075 !+                     "QUASI-LINEAR THEORY OF WIND-WAVE GENERATION
00076 !+                      APPLIED TO WAVE FORECASTING". JPO, VOL 21, PP 1631-1642.
00077 !
00078 !
00079 !
00080 !
00081 !
00082 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00083 !| ALPHA          |-->| CHARNOCK CONSTANT
00084 !| CDRAG          |-->| WIND DRAG COEFFICIENT
00085 !| GRAVIT         |-->| GRAVITY ACCELERATION
00086 !| NPOIN2         |-->| NUMBER OF POINTS IN 2D MESH
00087 !| TAUWAV         |-->| STRESS DUE TO WAVES
00088 !| USTAR          |<--| FRICTION VELOCITY
00089 !| UV             |-->| WIND VELOCITY ALONG X AT THE MESH POINTS
00090 !| VV             |-->| WIND VELOCITY ALONG Y AT THE MESH POINTS
00091 !| XKAPPA         |-->| VON KARMAN CONSTANT
00092 !| Z0             |<--| SURFACE ROUGHNESS LENGTH
00093 !| ZVENT          |-->| COTE A LAQUELLE EST MESURE LE VENT (M)
00094 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00095 !
00096       IMPLICIT NONE
00097 !
00098 !.....VARIABLES IN ARGUMENT
00099 !     """"""""""""""""""""
00100       INTEGER  NPOIN2
00101       DOUBLE PRECISION USTAR(NPOIN2) , Z0(NPOIN2) , TAUWAV(NPOIN2)
00102       DOUBLE PRECISION    UV(NPOIN2) , VV(NPOIN2)
00103       DOUBLE PRECISION CDRAG , ALPHA , XKAPPA , ZVENT, GRAVIT
00104 !
00105 !.....LOCAL VARIABLES
00106 !     """""""""""""""""
00107       INTEGER  ITRMIN, ITRMAX, ITR   , IP
00108       DOUBLE PRECISION TAUT  , UVENT , TAUW  , USMIN , SEUIL , X
00109       DOUBLE PRECISION USTEMP
00110 !
00111 !
00112       USMIN =1.D-6
00113       SEUIL =1.D-7
00114       ITRMIN=1
00115       ITRMAX=15
00116 !
00117 !.....MAIN LOOP ON THE NODES OF THE 2D MESH
00118 !     """""""""""""""""""""""""""""""""""""""""""""""""""""
00119       DO IP=1,NPOIN2
00120 !
00121 !.......COMPUTES THE TOTAL STRESS
00122 !       """""""""""""""""""""""""""""""
00123         UVENT=SQRT(UV(IP)**2+VV(IP)**2)
00124         TAUW =TAUWAV(IP)
00125         CALL TAUTOT
00126      &( TAUT  , UVENT , TAUW  , CDRAG , ALPHA , XKAPPA, ZVENT , SEUIL ,
00127      &  GRAVIT, ITR   , ITRMIN, ITRMAX)
00128 !
00129 !.......COMPUTES THE FRICTION VELOCITY
00130 !       """""""""""""""""""""""""""""""""""
00131         USTAR(IP)=SQRT(TAUT)
00132 !
00133 !.......COMPUTES TEH ROUGHNESS LENGTH
00134 !       """"""""""""""""""""""""""""""""""
00135         USTEMP=MAX(USTAR(IP),USMIN)
00136         X     =MIN(TAUWAV(IP)/USTEMP**2,0.999D0)
00137         Z0(IP)=ALPHA*USTEMP**2/(GRAVIT*SQRT(1.D0-X))
00138 !
00139       ENDDO
00140 !
00141       RETURN
00142       END

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