w3src2md.F90 File Reference

Contains MODULE W3SRC2MD. More...

Go to the source code of this file.

Modules
module	w3src2md
	Tolman and Chalikov (1996) input and dissipation source terms.

Functions/Subroutines
subroutine	w3src2md::w3spr2 (A, CG, WN, DEPTH, FPI, U, USTAR, EMEAN, FMEAN, WNMEAN, AMAX, ALFA, FP)
	Calculate mean wave parameters for the use in the source term routines (Tolman and Chalikov). More...
subroutine	w3src2md::w3sin2 (A, CG, K, U, UDIR, CD, Z0, FPI, S, D)
	Calculate input source term. More...
subroutine	w3src2md::w3sds2 (A, CG, K, FPI, USTAR, ALFA, S, D)
	Calculate whitecapping source term and diagonal term of derivative. More...
subroutine	w3src2md::inptab
	Generate an interpolation table for the air-sea interaction parameter of Chalikov and Belevich (1993). More...
real function	w3beta (OMA, CL, NDST)
	Calculate wind-wave interaction parameter beta. More...

Detailed Description

Contains MODULE W3SRC2MD.

Author

H. L. Tolman

Date

29-May-2009

Definition in file w3src2md.F90.

Function/Subroutine Documentation

w3beta()

real function inptab::w3beta	(	real, intent(in)	OMA,
		real, intent(in)	CL,
		integer, intent(in)	NDST
	)

Calculate wind-wave interaction parameter beta.

Parameters

OMA	Non-dimensional apparent frequency.
CL	Drag coefficient at height l.
NDST

Returns

W3BETA Wind-wave interaction parameter multiplied by density ratio.

Author

H. L. Tolman

D. Chalikov

Date

21-Feb-2004

Definition at line 1082 of file w3src2md.F90.

      !/
      !/                  +-----------------------------------+
      !/                  | WAVEWATCH III           NOAA/NCEP |
      !/                  |           H. L. Tolman            |
      !/                  |            D.Chalikov             |
      !/                  |                        FORTRAN 90 |
      !/                  | Last update :         21-Feb-2004 |
      !/                  +-----------------------------------+
      !/
      !/    06-Dec-1996 : Final version 1.18 / FORTRAN 77 version.
      !/    06-Dec-1999 : Upgrade to FORTRAN 90               ( version 2.00 )
      !/    21-Feb-2004 : Multiple model version.             ( version 3.06 )
      !/
      !  1. Purpose :
      !
      !     Calculate wind-wave interaction parameter beta.
      !
      !  2. Method :
      !
      !     Chalikov and Belevich (1992), see also manual.
      !
      !  3. Parameters :
      !
      !     Parameter list
      !     ----------------------------------------------------------------
      !       W3BETA  Real  O   Wind-wave interaction parameter multiplied
      !                         by density ratio.
      !       OMA     Real  I   Non-dimensional apparent frequency.
      !
      !                         OMA = OMEGA | U | cos(theta-theta ) / g
      !                                        l                 w
      !
      !       CL      Real  I   Drag coefficient at height l
      !     ----------------------------------------------------------------
      !
      !  4. Subroutines used :
      !
      !  5. Called by :
      !
      !  6. Error messages :
      !
      !  7. Remarks :
      !
      !  8. Structure :
      !
      !     See source code.
      !
      !  9. Switches :
      !
      !     !/S   Enable subroutine tracing.
      !     !/T0  Enable test output.
      !
      ! 10. Source code :
      !
      !/ ------------------------------------------------------------------- /
      IMPLICIT NONE
      !/
      !/ ------------------------------------------------------------------- /
      !/ Parameter list
      !/
      INTEGER, INTENT(IN)     :: NDST
      REAL, INTENT(IN)        :: OMA, CL
      !/
      !/ ------------------------------------------------------------------- /
      !/ Local parameters
      !/
#ifdef W3_S
      INTEGER, SAVE           :: IENT = 0
#endif
      REAL                    :: OM1, OM2, A0, A1, A2, A3, A4, A5,    &
           A6, A7, A8, A9, A10
      !/
      !/ ------------------------------------------------------------------- /
      !/
#ifdef W3_S
      CALL strace (ient, 'W3BETA')
#endif
      !
#ifdef W3_T0
      WRITE (ndst,9000) oma, cl
#endif
      !
      ! calculate Omegas
      !
      om1    =  1.075 +  75.*cl
      om2    =  1.2   + 300.*cl
      !
      ! calculate factors a
      !
      a1     =  0.25  + 395.*cl
      a2     =  0.35  + 150.*cl
      a4     =  0.3   + 300.*cl
      a9     =  0.35  + 240.*cl
      a10    = -0.06  + 470.*cl
      !
      a5     =  a4 * om1
      a0     =  0.25 * a5**2 / a4
      a3     = (a0-a2-a1) / (a0+a4+a5)
      a6     =  a0 * (1.-a3)
      a7     = (a9*(om2-1)**2+a10) / (om2-om1)
      a8     =  a7 * om1
      !
#ifdef W3_T0
      WRITE (ndst,9001) om1, om2
      WRITE (ndst,9002) a0, a1, a2, a3, a4, a5, a6, a7, a8, a9, a10
#endif
      !
      ! calculate beta * 1.e4
      !
      IF  ( oma .LT. -1. ) THEN
        w3beta = -a1 * oma**2 - a2
      ELSE IF (oma .LT. om1/2.) THEN
        w3beta =  a3 * oma * ( a4 * oma - a5 ) - a6
      ELSE IF (oma .LT. om1) THEN
        w3beta =       oma * ( a4 * oma - a5 )
      ELSE IF (oma .LT. om2) THEN
        w3beta = a7 * oma - a8
      ELSE
        w3beta = a9 * (oma-1.)**2 + a10
      END IF
      !
      ! beta * dwat / dair
      !
      w3beta = w3beta * 1.e-4
#ifdef W3_T0
      WRITE (ndst,9003) w3beta
#endif
      !
      RETURN
      !
      ! Formats
      !
#ifdef W3_T0
9000  FORMAT ( ' TEST W3BETA : INPUT : ',2e10.3)
9001  FORMAT ( ' TEST W3BETA : OM1-2 : ',2e10.3)
9002  FORMAT ( ' TEST W3BETA : A0-10 : ',5e10.3/                   &
           '             ',6e10.3)
9003  FORMAT ( ' TEST W3BETA : BETA  : ',e10.3)
#endif
      !/
      !/ End of W3BETA ----------------------------------------------------- /
      !/

References w3servmd::strace().

Referenced by w3src2md::inptab().