*IF DEF,W08_1A GLW1F404.51 C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.15643 C GTS2F400.15644 C Use, duplication or disclosure of this code is subject to the GTS2F400.15645 C restrictions as set forth in the contract. GTS2F400.15646 C GTS2F400.15647 C Meteorological Office GTS2F400.15648 C London Road GTS2F400.15649 C BRACKNELL GTS2F400.15650 C Berkshire UK GTS2F400.15651 C RG12 2SZ GTS2F400.15652 C GTS2F400.15653 C If no contract has been raised with this copy of the code, the use, GTS2F400.15654 C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.15655 C to do so must first be obtained in writing from the Head of Numerical GTS2F400.15656 C Modelling at the above address. GTS2F400.15657 C ******************************COPYRIGHT****************************** GTS2F400.15658 C GTS2F400.15659 SEMEAN.3SUBROUTINE SEMEAN (F3, IJS, IJL, 2SEMEAN.4 *CALL ARGWVAL
SEMEAN.5 *CALL ARGWVFD
SEMEAN.6 *CALL ARGWVMN
SEMEAN.7 & icode) SEMEAN.8 SEMEAN.9 *CALL TYPWVFD
SEMEAN.10 *CALL TYPWVMN
SEMEAN.11 *CALL TYPWVAL
SEMEAN.12 SEMEAN.13 C ---------------------------------------------------------------------- SEMEAN.14 C SEMEAN.15 C**** *SEMEAN* - COMPUTATION OF TOTAL ENERGY AT EACH GRID POINT. SEMEAN.16 C SEMEAN.17 C S.D. HASSELMANN. SEMEAN.18 C OPTIMIZED BY: L. ZAMBRESKY AND H. GUENTHER SEMEAN.19 C SEMEAN.20 C* PURPOSE. SEMEAN.21 C -------- SEMEAN.22 C SEMEAN.23 C TO COMPUTE TOTAL ENERGY AT EACH GRID POINT. SEMEAN.24 C SEMEAN.25 C** INTERFACE. SEMEAN.26 C ---------- SEMEAN.27 C SEMEAN.28 C *CALL* *SEMEAN(F3, IJS, IJL)* SEMEAN.29 C *F3* - SPECTRUM. SEMEAN.30 C *IJS* - INDEX OF FIRST GRIDPOINT SEMEAN.31 C *IJL* - INDEX OF LAST GRIDPOINT SEMEAN.32 C SEMEAN.33 C METHOD. SEMEAN.34 C ------- SEMEAN.35 C SEMEAN.36 C NONE. SEMEAN.37 C SEMEAN.38 C EXTERNALS. SEMEAN.39 C ---------- SEMEAN.40 C SEMEAN.41 C NONE. SEMEAN.42 C SEMEAN.43 C REFERENCE. SEMEAN.44 C ---------- SEMEAN.45 C SEMEAN.46 C NONE. SEMEAN.47 C SEMEAN.48 C ---------------------------------------------------------------------- SEMEAN.49 C SEMEAN.50 DIMENSION F3(0:NIBLO,NANG,NFRE), TEMP(NIBLO) SEMEAN.51 C SEMEAN.52 C ---------------------------------------------------------------------- SEMEAN.53 C SEMEAN.54 C* 1. INITIALISE ENERGY ARRAY. SEMEAN.55 C ------------------------ SEMEAN.56 C SEMEAN.57 1000 CONTINUE SEMEAN.58 DO 1001 IJ=IJS,IJL SEMEAN.59 EMEAN(IJ) = 0. SEMEAN.60 1001 CONTINUE SEMEAN.61 C SEMEAN.62 C ---------------------------------------------------------------------- SEMEAN.63 C SEMEAN.64 C* 2. INTEGRATE OVER FREQUENCIES AND DIRECTION. SEMEAN.65 C ----------------------------------------- SEMEAN.66 C SEMEAN.67 2000 CONTINUE SEMEAN.68 DO 2001 M=1,NFRE SEMEAN.69 DO 2002 IJ=IJS,IJL SEMEAN.70 TEMP(IJ) = 0. SEMEAN.71 2002 CONTINUE SEMEAN.72 DO 2003 K=1,NANG SEMEAN.73 DO 2004 IJ=IJS,IJL SEMEAN.74 TEMP(IJ) = TEMP(IJ)+F3(IJ,K,M) SEMEAN.75 2004 CONTINUE SEMEAN.76 2003 CONTINUE SEMEAN.77 DO 2005 IJ=IJS,IJL SEMEAN.78 EMEAN(IJ) = EMEAN(IJ)+DFIM(M)*TEMP(IJ) SEMEAN.79 2005 CONTINUE SEMEAN.80 2001 CONTINUE SEMEAN.81 C SEMEAN.82 C ---------------------------------------------------------------------- SEMEAN.83 C SEMEAN.84 C* 3. ADD TAIL ENERGY. SEMEAN.85 C ---------------- SEMEAN.86 C SEMEAN.87 3000 CONTINUE SEMEAN.88 DELT25 = FR(NFRE)/4.*DELTH SEMEAN.89 DO 3001 IJ=IJS,IJL SEMEAN.90 EMEAN(IJ) = EMEAN(IJ)+DELT25*TEMP(IJ) SEMEAN.91 3001 CONTINUE SEMEAN.92 SEMEAN.93 RETURN SEMEAN.94 END SEMEAN.95 *ENDIF SEMEAN.96