SUBROUTINE TAUHF (FRMAX, 1TAUHF.4
*CALL ARGWVAL 
TAUHF.5
*CALL ARGWVCP TAUHF.6
*CALL ARGWVTB TAUHF.7
& icode) TAUHF.8
TAUHF.9
*CALL PARCONS TAUHF.10
*CALL PARWVTB TAUHF.11
TAUHF.12
PARAMETER (JTOT = 50) TAUHF.13
TAUHF.14
*CALL TYPWVCP TAUHF.15
*CALL TYPWVTB TAUHF.16
*CALL TYPWVAL TAUHF.17
TAUHF.18
C ---------------------------------------------------------------------- TAUHF.19
C TAUHF.20
C**** *TAUHF* - COMPUTATION OF HIGH-FREQUENCY STRESS. TAUHF.21
C TAUHF.22
C PETER A.E.M. JANSSEN KNMI OCTOBER 90 TAUHF.23
C TAUHF.24
C* PURPOSE. TAUHF.25
C --------- TAUHF.26
C TAUHF.27
C COMPUTE HIGH-FREQUENCY WAVE STRESS TAUHF.28
C TAUHF.29
C** INTERFACE. TAUHF.30
C ---------- TAUHF.31
C TAUHF.32
C *CALL* *TAUHF (FRMAX)* TAUHF.33
C *FRMAX - LAST MODEL FREQUENCY FR(ML). TAUHF.34
C TAUHF.35
C METHOD. TAUHF.36
C ------- TAUHF.37
C TAUHF.38
C SEE REFERENCE FOR WAVE STRESS CALCULATION. TAUHF.39
C TAUHF.40
C EXTERNALS. TAUHF.41
C ---------- TAUHF.42
C TAUHF.43
C NONE. TAUHF.44
C TAUHF.45
C REFERENCE. TAUHF.46
C ---------- TAUHF.47
C TAUHF.48
C FOR QUASILINEAR EFFECT SEE PETER A.E.M. JANSSEN,1990. TAUHF.49
C TAUHF.50
C ---------------------------------------------------------------------- TAUHF.51
C* 1. PRELIMINARY CALCULATIONS. TAUHF.52
C ------------------------- TAUHF.53
C TAUHF.54
USTARM = 5. TAUHF.55
ALPHAM = 10.*ALPHA TAUHF.56
DELUST = USTARM/FLOAT(IUSTAR) TAUHF.57
DELALP = ALPHAM/FLOAT(IALPHA) TAUHF.58
C TAUHF.59
CONST1 = BETAMAX/XKAPPA**2 TAUHF.60
OMEGAC = ZPI*FRMAX TAUHF.61
C TAUHF.62
DO 1100 L=0,IALPHA TAUHF.63
DO 1200 K=0,IUSTAR TAUHF.64
TAUHFT(K,L) = 0. TAUHF.65
1200 CONTINUE TAUHF.66
1100 CONTINUE TAUHF.67
C TAUHF.68
C* 2. CALCULATE HIGH-FREQUENCY CONTRIBUTION TO STRESS. TAUHF.69
C ------------------------------------------------ TAUHF.70
C TAUHF.71
X0 = 0.05 TAUHF.72
DO 2100 L=0,IALPHA TAUHF.73
DO 2200 K=0,IUSTAR TAUHF.74
UST = MAX(FLOAT(K)*DELUST,0.000001) TAUHF.75
Z0 = UST**2*(ALPHA+FLOAT(L)*DELALP)/G TAUHF.76
OMEGACC = MAX(OMEGAC,X0*G/UST) TAUHF.77
YC = OMEGACC*SQRT(Z0/G) TAUHF.78
DELY = MAX((1.-YC)/FLOAT(JTOT),0.) TAUHF.79
DO 2300 J=1,JTOT TAUHF.80
Y = YC+FLOAT(J-1)*DELY TAUHF.81
OMEGA = Y*SQRT(G/Z0) TAUHF.82
CM = G/OMEGA TAUHF.83
ZX = UST/CM +ZALP TAUHF.84
ZARG = MIN(XKAPPA/ZX,20.) TAUHF.85
ZMU = MIN(G*Z0/CM**2*EXP(ZARG),1.) TAUHF.86
C TAUHF.87
ZLOG = MIN(ALOG(ZMU),0.) TAUHF.88
ZBETA = CONST1*ZMU*ZLOG**4 TAUHF.89
TAUHFT(K,L) = TAUHFT(K,L)+ZBETA/Y*DELY TAUHF.90
2300 CONTINUE TAUHF.91
2200 CONTINUE TAUHF.92
2100 CONTINUE TAUHF.93
TAUHF.94
RETURN TAUHF.95
END TAUHF.96
*ENDIF TAUHF.97