*IF DEF,W06_1A FLDOUT.2
C *****************************COPYRIGHT****************************** FLDOUT.3
C (c) CROWN COPYRIGHT 1996, METEOROLOGICAL OFFICE, All Rights Reserved. FLDOUT.4
C FLDOUT.5
C Use, duplication or disclosure of this code is subject to the FLDOUT.6
C restrictions as set forth in the contract. FLDOUT.7
C FLDOUT.8
C Meteorological Office FLDOUT.9
C London Road FLDOUT.10
C BRACKNELL FLDOUT.11
C Berkshire UK FLDOUT.12
C RG12 2SZ FLDOUT.13
C FLDOUT.14
C If no contract has been raised with this copy of the code, the use, FLDOUT.15
C duplication or disclosure of it is strictly prohibited. Permission FLDOUT.16
C to do so must first be obtained in writing from the Head of Numerical FLDOUT.17
C Modelling at the above address. FLDOUT.18
C ******************************COPYRIGHT****************************** FLDOUT.19
C*LLL FLDOUT.20
C 12. SUBROUTINE FLDOUT FLDOUT.21
C FLDOUT.22
! FLDOUT.23
! Description: FLDOUT.24
! FLDOUT.25
! Method: FLDOUT.26
! FLDOUT.27
! FLDOUT.28
! FLDOUT.29
! Current Code Owner: Martin Holt FLDOUT.30
! FLDOUT.31
! History: FLDOUT.32
! Version Date Comment FLDOUT.33
! ------- ---- ------- FLDOUT.34
! UM4.1 June 1996 Original code. M Holt FLDOUT.35
! FLDOUT.36
! Code Description: FLDOUT.37
! Language: FORTRAN 77 + common extensions. FLDOUT.38
! FLDOUT.39
!- End of header FLDOUT.40
! FLDOUT.41
C DOCUMENTATION FLDOUT.42
C FLDOUT.43
C SEE WAVE MODEL DOCUMENTATION PAPER. FLDOUT.44
C FLDOUT.45
C DESCRIPTION FLDOUT.46
C FLDOUT.47
C THIS ROUTINE OUTPUTS INTEGRATED WAVE PARAMETERS CALCULATED FLDOUT.48
C BY CALLING WAVEH. Arrays are prepared for use by STASH FLDOUT.49
C FLDOUT.50
C ********************************************************************* FLDOUT.51
FLDOUT.52
SUBROUTINE FLDOUT(ia1, 1,3FLDOUT.53
*CALL ARGWVAL 
FLDOUT.54
*CALL ARGWVFD FLDOUT.55
*CALL ARGWVGD FLDOUT.56
+ l_wvtra,len1,energy,uwind,vwind,mdata,ndata,nfldmax FLDOUT.57
+ ,icode) FLDOUT.58
FLDOUT.59
FLDOUT.60
FLDOUT.61
*CALL TYPWVAL FLDOUT.62
*CALL TYPWVFD FLDOUT.63
*CALL TYPWVGD FLDOUT.64
FLDOUT.65
C FLDOUT.66
LOGICAL IA1(ngx,ngy) FLDOUT.67
REAL WH(ngx*ngy,nfldmax) FLDOUT.68
real windsp(ndata),windir(ndata) FLDOUT.69
REAL uwind(mdata),vwind(mdata) ! input U10 wind components FLDOUT.70
FLDOUT.71
LOGICAL l_wvtra ! IN T if wavetrains required FLDOUT.72
FLDOUT.73
INTEGER ICODE FLDOUT.74
C FLDOUT.75
*CALL C_MDI FLDOUT.76
C FLDOUT.77
C ---------------------------------------------------------------------- FLDOUT.78
C VALUES STORED IN WH FIELD NUMBER PP FLDOUT.79
C 1 TOTAL WAVE HEIGHT 94 387 FLDOUT.80
C 2 MEAN DIRECTION - ? FLDOUT.81
C 3 PRINCIPAL DIRECTION 96 394 FLDOUT.82
C 4 ZERO UP CROSSING PERIOD 95 393 FLDOUT.83
C 5 MEAN PERIOD - ? FLDOUT.84
C 6 PEAK PERIOD 116 392 FLDOUT.85
C 7 TOTAL WINDSEA WAVE HEIGHT 97 385 FLDOUT.86
C 8 MEAN DIRECTION - ? FLDOUT.87
C 9 PRINCIPAL DIRECTION 99 389 FLDOUT.88
C 10 ZERO UP CROSSING PERIOD 98 388 FLDOUT.89
C 11 MEAN PERIOD - ? FLDOUT.90
C 12 PEAK PERIOD - ? FLDOUT.91
C 13 TOTAL SWELL WAVE HEIGHT 100 386 FLDOUT.92
C 14 MEAN DIRECTION - ? FLDOUT.93
C 15 PRINCIPAL DIRECTION 102 391 FLDOUT.94
C 16 ZERO UP CROSSING PERIOD 101 390 FLDOUT.95
C 17 MEAN PERIOD - ? FLDOUT.96
C 18 PEAK PERIOD - ? FLDOUT.97
FLDOUT.98
C ?? WAVE STRESS - ABSOLUTE VALUE 211 366 FLDOUT.99
C ?? ROUGHNESS LENGTH WAVE DEPENDENT 212 367 FLDOUT.100
C FLDOUT.101
c 19 swh ) first FLDOUT.102
c 20 swh ) second wave train FLDOUT.103
c 21 swh ) third FLDOUT.104
c 22 swh ) fourth FLDOUT.105
c FLDOUT.106
c 23 period ) first FLDOUT.107
c 24 period ) second wave train FLDOUT.108
c 25 period ) third FLDOUT.109
c 26 period ) fourth FLDOUT.110
c FLDOUT.111
c 27 dir ) first FLDOUT.112
c 28 dir ) second wave train FLDOUT.113
c 29 dir ) third FLDOUT.114
c 30 dir ) fourth FLDOUT.115
c FLDOUT.116
c 31 number of wave trains present FLDOUT.117
C ---------------------------------------------------------------------- FLDOUT.118
C FLDOUT.119
C still need somewhere to put windsp windir and other data arrays into FLDOUT.120
C WH for passing to stash eg stress and zo FLDOUT.121
FLDOUT.122
C FLDOUT.123
C need to have all diagnostic arrays on the full grid for STASH FLDOUT.124
C FLDOUT.125
FLDOUT.126
C calculate wind speed and direction from components FLDOUT.127
FLDOUT.128
do i=1,ndata FLDOUT.129
FLDOUT.130
windsp(i) = sqrt(uwind(i)*uwind(i) + vwind(i)*vwind(i) ) FLDOUT.131
windir(i) = atan2(uwind(i),vwind(i)) FLDOUT.132
FLDOUT.133
enddo FLDOUT.134
FLDOUT.135
C FLDOUT.136
C *** CALCULATE WAVE INTEGRATED FIELDS FLDOUT.137
C FLDOUT.138
ngrid=ngx*ngy FLDOUT.139
FLDOUT.140
call waveh(ia1, FLDOUT.141
*CALL ARGWVAL FLDOUT.142
*CALL ARGWVFD FLDOUT.143
& l_wvtra,ndata,ngrid,nfldmax,len1,energy,windsp,windir,wh,icode) FLDOUT.144
FLDOUT.145
RETURN FLDOUT.146
END FLDOUT.147
*ENDIF FLDOUT.148