SUBROUTINE FEMEAN (F, IJS, IJL, ishallo, 2FEMEAN.4
*CALL ARGWVAL 
FEMEAN.5
*CALL ARGWVFD FEMEAN.6
*CALL ARGWVMN FEMEAN.7
*CALL ARGWVSH FEMEAN.8
& icode) FEMEAN.9
FEMEAN.10
*CALL PARWVSH FEMEAN.11
*CALL PARCONS FEMEAN.12
FEMEAN.13
*CALL TYPWVFD FEMEAN.14
*CALL TYPWVMN FEMEAN.15
*CALL TYPWVSH FEMEAN.16
*CALL TYPWVAL FEMEAN.17
FEMEAN.18
C ---------------------------------------------------------------------- FEMEAN.19
C FEMEAN.20
C**** *FEMEAN* - COMPUTATION OF MEAN FREQUENCY AT EACH GRID POINT. FEMEAN.21
C FEMEAN.22
C S.D. HASSELMANN FEMEAN.23
C MODIFIED : P.JANSSEN (INTEGRATION OF F**-4 TAIL) FEMEAN.24
C OPTIMIZED BY : L. ZAMBRESKY AND H. GUENTHER FEMEAN.25
C FEMEAN.26
C FEMEAN.27
C* PURPOSE. FEMEAN.28
C -------- FEMEAN.29
C FEMEAN.30
C COMPUTE MEAN FREQUENCY AT EACH GRID POINT. FEMEAN.31
C FEMEAN.32
C** INTERFACE. FEMEAN.33
C ---------- FEMEAN.34
C FEMEAN.35
C *CALL* *FEMEAN (F, IJS, IJL)* FEMEAN.36
C *F* - SPECTRUM. FEMEAN.37
C *IJS* - INDEX OF FIRST GRIDPOINT FEMEAN.38
C *IJL* - INDEX OF LAST GRIDPOINT FEMEAN.39
C FEMEAN.40
C METHOD. FEMEAN.41
C ------- FEMEAN.42
C FEMEAN.43
C NONE. FEMEAN.44
C FEMEAN.45
C EXTERNALS. FEMEAN.46
C ---------- FEMEAN.47
C FEMEAN.48
C NONE. FEMEAN.49
C FEMEAN.50
C REFERENCE. FEMEAN.51
C ---------- FEMEAN.52
C FEMEAN.53
C NONE. FEMEAN.54
C FEMEAN.55
C ---------------------------------------------------------------------- FEMEAN.56
C FEMEAN.57
DIMENSION F(0:NIBLO,NANG,NFRE) FEMEAN.58
C FEMEAN.59
C ---------------------------------------------------------------------- FEMEAN.60
C FEMEAN.61
DIMENSION TEMP1(NIBLO), TEMP2(NIBLO) FEMEAN.62
C FEMEAN.63
C ---------------------------------------------------------------------- FEMEAN.64
C FEMEAN.65
C* 1. INITIALISE MEAN FREQUENCY ARRAY AND TAIL FACTOR. FEMEAN.66
C ------------------------------------------------ FEMEAN.67
C FEMEAN.68
1000 CONTINUE FEMEAN.69
DO 1001 IJ=IJS,IJL FEMEAN.70
FMEAN(IJ) = 0. FEMEAN.71
1001 CONTINUE FEMEAN.72
DELT25 = 0.20*DELTH FEMEAN.73
C FEMEAN.74
C ---------------------------------------------------------------------- FEMEAN.75
C FEMEAN.76
C* 2. INTEGRATE OVER FREQUENCIES AND DIRECTIONS. FEMEAN.77
C ------------------------------------------ FEMEAN.78
C FEMEAN.79
2000 CONTINUE FEMEAN.80
FEMEAN.81
IF (ISHALLO.EQ.1) THEN FEMEAN.82
C FEMEAN.83
C* 2.1 DEEP WATER INTEGRATION. FEMEAN.84
C ----------------------- FEMEAN.85
C FEMEAN.86
2100 CONTINUE FEMEAN.87
DO 2101 M=1,NFRE FEMEAN.88
FD = DFIM(M)/FR(M) FEMEAN.89
DO 2102 IJ=IJS,IJL FEMEAN.90
TEMP2(IJ) = 0. FEMEAN.91
2102 CONTINUE FEMEAN.92
DO 2103 K=1,NANG FEMEAN.93
DO 2104 IJ=IJS,IJL FEMEAN.94
TEMP2(IJ) = TEMP2(IJ)+F(IJ,K,M) FEMEAN.95
2104 CONTINUE FEMEAN.96
2103 CONTINUE FEMEAN.97
DO 2105 IJ=IJS,IJL FEMEAN.98
FMEAN(IJ) = FMEAN(IJ)+FD*TEMP2(IJ) FEMEAN.99
2105 CONTINUE FEMEAN.100
2101 CONTINUE FEMEAN.101
CSHALLOW FEMEAN.102
ELSE FEMEAN.103
C FEMEAN.104
C* 2.2 SHALLOW WATER INTEGRATION. FEMEAN.105
C -------------------------- FEMEAN.106
C FEMEAN.107
2200 CONTINUE FEMEAN.108
DO 2201 IJ=IJS,IJL FEMEAN.109
AKMEAN(IJ) = 0. FEMEAN.110
2201 CONTINUE FEMEAN.111
DO 2202 M=1,NFRE FEMEAN.112
FD=DFIM(M)/FR(M) FEMEAN.113
DO 2203 IJ=IJS,IJL FEMEAN.114
TEMP1(IJ) = DFIM(M)/SQRT(TFAK(INDEP(IJ),M)) FEMEAN.115
TEMP2(IJ) = 0. FEMEAN.116
2203 CONTINUE FEMEAN.117
DO 2204 K=1,NANG FEMEAN.118
DO 2205 IJ=IJS,IJL FEMEAN.119
TEMP2(IJ) = TEMP2(IJ)+F(IJ,K,M) FEMEAN.120
2205 CONTINUE FEMEAN.121
2204 CONTINUE FEMEAN.122
DO 2206 IJ=IJS,IJL FEMEAN.123
FMEAN(IJ) = FMEAN(IJ)+FD*TEMP2(IJ) FEMEAN.124
AKMEAN(IJ) = AKMEAN(IJ)+TEMP1(IJ)*TEMP2(IJ) FEMEAN.125
2206 CONTINUE FEMEAN.126
2202 CONTINUE FEMEAN.127
C FEMEAN.128
C ADD TAIL TO MEAN WAVENUMBER AND NORMALIZE WITH TOTAL ENERGY. FEMEAN.129
C FEMEAN.130
DEL2 = DELT25*SQRT(G)/ZPI FEMEAN.131
DO 2207 IJ=IJS,IJL FEMEAN.132
AKMEAN(IJ) = AKMEAN(IJ)+DEL2*TEMP2(IJ) FEMEAN.133
AKMEAN(IJ) = (EMEAN(IJ)/AKMEAN(IJ))**2 FEMEAN.134
2207 CONTINUE FEMEAN.135
ENDIF FEMEAN.136
CSHALLOW FEMEAN.137
C FEMEAN.138
C* 3. ADD TAIL CORRECTION TO MEAN FREQUENCY AND FEMEAN.139
C* NORMALIZE WITH TOTAL ENERGY. FEMEAN.140
C ------------------------------------------ FEMEAN.141
C FEMEAN.142
3000 CONTINUE FEMEAN.143
FEMEAN.144
FEMEAN.145
DO 3001 IJ=IJS,IJL FEMEAN.146
FMEAN(IJ) = FMEAN(IJ)+DELT25*TEMP2(IJ) FEMEAN.147
FMEAN(IJ) = EMEAN(IJ)/FMEAN(IJ) FEMEAN.148
3001 CONTINUE FEMEAN.149
FEMEAN.150
RETURN FEMEAN.151
END FEMEAN.152
*ENDIF FEMEAN.153