tautot.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\tomawac\tautot.f
00002 !
00056                      SUBROUTINE TAUTOT
00057 !                    *****************
00058 !
00059      &( TAUT  , UVENT , TAUW  , CDRAG , ALPHA , XKAPPA, ZVENT , SEUIL ,
00060      &  GRAVIT, ITR   , ITRMIN, ITRMAX)
00061 !
00062 !***********************************************************************
00063 ! TOMAWAC   V6P1                                   28/06/2011
00064 !***********************************************************************
00065 !
00066 !
00067 !
00068 !reference  JANSSEN P.A.E.M (1989) :
00069 !+                     "WIND-INDUCED STRESS AND THE DRAG OF AIR FLOW
00070 !+                      OVER SEA WAVES". JPO, VOL 19, PP 745-754.
00071 !reference JANSSEN P.A.E.M (1991) :
00072 !+                     "QUASI-LINEAR THEORY OF WIND-WAVE GENERATION
00073 !+                      APPLIED TO WAVE FORECASTING". JPO, VOL 21, PP 1631-1642.
00074 !
00075 !
00076 !
00077 !
00078 !
00079 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00080 !| ALPHA          |-->| CHARNOCK CONSTANT
00081 !| CDRAG          |-->| WIND DRAG COEFFICIENT
00082 !| GRAVIT         |-->| GRAVITY ACCELERATION
00083 !| ITR            |<--| NOMBRE D'ITERATIONS EFFECTUES
00084 !| ITRMAX         |-->| MAXIMUM NUMBER OF ITERATIONS
00085 !| ITRMIN         |-->| MINIMUM NUMBER OF ITERATIONS
00086 !| SEUIL          |-->| CONVERGENCY THRESHOLD - NEWTON METHOD
00087 !| TAUT           |<--| TOTAL STRESS
00088 !| TAUW           |-->| STRESS DUE TO THE WAVES
00089 !| UVENT          |-->| WIND SPEED AT HEIGHT ZVENT
00090 !| XKAPPA         |-->| VON KARMAN CONSTANT
00091 !| ZVENT          |-->| WIND MEASUREMENT LEVEL
00092 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00093 !
00094       IMPLICIT NONE
00095 !
00096 !.....VARIABLES IN ARGUMENT
00097 !     """"""""""""""""""""
00098       INTEGER  ITRMIN, ITRMAX, ITR
00099       DOUBLE PRECISION TAUT  , UVENT , TAUW  , ALPHA , XKAPPA , ZVENT
00100       DOUBLE PRECISION SEUIL , CDRAG , GRAVIT
00101 !
00102 !.....LOCAL VARIABLES
00103 !     """""""""""""""""
00104       DOUBLE PRECISION TAUMIN, XNUAIR, AUX   , USTO  , TAUO   , TAUN
00105       DOUBLE PRECISION USTN  , X     , ZE    , DIFF  , FOLD   , DFOLD
00106 !
00107 !
00108       ITR   =0
00109       TAUMIN=1.D-5
00110       XNUAIR=1.D-5
00111 !
00112 !.....INITIAL VALUES
00113 !     """"""""""""""""""
00114       USTO  =UVENT*SQRT(CDRAG)
00115       TAUO  =MAX(USTO**2,TAUW+TAUMIN)
00116 !
00117   190 CONTINUE
00118       ITR   = ITR+1
00119 !
00120 !.....ITERATION BY THE METHOD OF NEWTON
00121 !     """""""""""""""""""""""""""""""""""
00122       USTO  = SQRT(TAUO)
00123       X     = TAUW/TAUO
00124       ZE    = MAX(0.1D0*XNUAIR/USTO,ALPHA*TAUO/(GRAVIT*SQRT(1.D0-X)))
00125       AUX   = DLOG(ZVENT/ZE)
00126       FOLD  = USTO-XKAPPA*UVENT/AUX
00127       DFOLD = 1.D0-2.D0*XKAPPA*UVENT*(1.D0-1.5D0*X)/AUX**2/USTO/(1.D0-X)
00128       USTN  = USTO-FOLD/DFOLD
00129       TAUN  = MAX(USTN**2,TAUW+TAUMIN)
00130 !
00131 !.....CONVERGENCE CRITERIA
00132 !     """"""""""""""""""""""""
00133       DIFF=ABS(TAUN-TAUO)/TAUO
00134       TAUO=TAUN
00135       IF (ITR.LT.ITRMIN) GOTO 190
00136       IF ((DIFF.GT.SEUIL).AND.(ITR.LT.ITRMAX)) GOTO 190
00137 !
00138 !.....APPLIES THE SOLUTION
00139 !     """""""""""""""""""""""""""""""""""
00140       TAUT=TAUN
00141 !
00142       RETURN
00143       END

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