Page 1 Source Listing 2014-09-16 17:02 multiwavespec.f90 multiwavespec.f90(757): remark #8290: Recommended relationship between field width 'W' and the number of fractional digits 'D' in... 75X,'n : Number of fields with Hs > ',f4.2, & ------------------------------------------------^ multiwavespec.f90(759): remark #8290: Recommended relationship between field width 'W' and the number of fractional digits 'D' in... 75X,'x : Number of fields with Hs > ',f4.2, & ------------------------------------------------^ multiwavespec.f90(2053): remark #8290: Recommended relationship between field width 'W' and the number of fractional digits 'D' i... 2931 FORMAT (1X,F6.4,6E11.3) ------------------^ multiwavespec.f90(2052): remark #8290: Recommended relationship between field width 'W' and the number of fractional digits 'D' i... 2930 FORMAT (1X,F6.4,2X,6E11.3) ------------------^ Page 2 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1 !/ ------------------------------------------------------------------- / 2 PROGRAM WAVESPEC 3 !/ 4 !/ +-----------------------------------+ 5 !/ | WAVEWATCH III NOAA/NCEP | 6 !/ | H. L. Tolman | 7 !/ | J.H. Alves | 8 !/ | A. Chawla | 9 !/ | F. Ardhuin | 10 !/ | E. Rogers | 11 !/ | T. Campbell | 12 !/ | FORTRAN 90 | 13 !/ | Last update : 26-Dec-2012 | 14 !/ +-----------------------------------+ 15 !/ 16 !/ 14-Jan-1999 : Final FORTRAN 77 ( version 1.18 ) 17 !/ 21-Jan-2000 : Upgrade to FORTRAN 90 ( version 2.00 ) 18 !/ 14-Feb-2000 : Exact nonlinear interactions ( version 2.01 ) 19 !/ 09-Jan-2001 : U* bug fix in tabular output ( version 2.05 ) 20 !/ 25-Jan-2001 : Flat grid version ( version 2.06 ) 21 !/ 02-Feb-2001 : Xnl version 3.0 ( version 2.07 ) 22 !/ 11-Jun-2001 : Clean up ( version 2.11 ) 23 !/ 11-Oct-2001 : Clean up, X*, Y* in tables ( version 2.14 ) 24 !/ 13-Nov-2002 : Add stress vector ( version 3.00 ) 25 !/ 27-Nov-2002 : First version of VDIA and MDIA ( version 3.01 ) 26 !/ 24-Dec-2004 : Multiple grid version. ( version 3.06 ) 27 !/ 17-Apr-2006 : Filter for directional spread. ( version 3.09 ) 28 !/ 23-Jun-2006 : Linear input added. ( version 3.09 ) 29 !/ 28-Jun-2006 : Adding file name preamble. ( version 3.09 ) 30 !/ 03-Jul-2006 : Separate flux modules. ( version 3.09 ) 31 !/ 28-Oct-2006 : Add partitioning option. ( version 3.10 ) 32 !/ 24-Mar-2007 : Add pars for entire spectrum. ( version 3.11 ) 33 !/ 25-Apr-2007 : Battjes-Janssen Sdb added. ( version 3.11 ) 34 !/ (J. H. Alves) 35 !/ 08-Aug-2007 : Creation of buoy log file added ( version 3.12 ) 36 !/ (switch O14 -- A. Chawla) 37 !/ 09-Oct-2007 : WAM 4+ Sin and Sds added. ( version 3.13 ) 38 !/ (F. Ardhuin) 39 !/ 09-Oct-2007 : Experimental Sbs (BS1) added. ( version 3.13 ) 40 !/ (F. Ardhuin) 41 !/ 09-Apr-2008 : Adding an additional output for ( version 3.12 ) 42 !/ WMO standard (A. Chawla) 43 !/ 29-Apr-2008 : Adjust format partition output. ( version 3.14 ) 44 !/ 29-May-2009 : Preparing distribution version. ( version 3.14 ) 45 !/ 30-Oct-2009 : Implement run-time grid selection. ( version 3.14 ) 46 !/ (W. E. Rogers & T. J. Campbell, NRL) 47 !/ 04-Mar-2010 : Added partitions bulletin output. ( version 3.14 ) 48 !/ (J. H. Alves) 49 !/ 20-Apr-2010 : Fix initialization of USTAR. ( version 3.14.1 ) 50 !/ 16-Jul-2012 : Move GMD (SNL3) and nonlinear filter (SNLS) 51 !/ from 3.15 (HLT). ( version 4.08 ) 52 !/ 23-Aug-2012 : Adding movable bed friction BT4 ( version 4.08 ) 53 !/ 26-Dec-2012 : Modified obsolete declarations. ( version 4.11 ) 54 !/ 55 !/ Copyright 2009-2012 National Weather Service (NWS), 56 !/ National Oceanic and Atmospheric Administration. All rights 57 !/ reserved. WAVEWATCH III is a trademark of the NWS. Page 3 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 58 !/ No unauthorized use without permission. 59 !/ 60 ! 1. Purpose : 61 ! 62 ! Post-processing of point output. 63 ! 64 ! 2. Method : 65 ! 66 ! Data is read from the grid output file out_pnt.ww3 (raw data) 67 ! and from the file multiwavespec.inp ( NDSI, output requests ). 68 ! Model definition and raw data files are read using WAVEWATCH III 69 ! subroutines. 70 ! 71 ! Output types ITYPE : Sub-type OTYPE : 72 ! -------------------- ----------------- 73 ! 0 : Check file. 74 ! 1 : Spectra. 75 ! 1 : Print plots. 76 ! 2 : Table of 1-D spectra 77 ! 3 : Transfer file 78 ! 2 : Table of mean wave parameters 79 ! 1 : Depth, current, wind 80 ! 2 : Mean wave pars. 81 ! 3 : Nondimensional pars. (U*) 82 ! 4 : Nondimensional pars. (U10) 83 ! 5 : Validation table 84 ! 6 : WMO standard output 85 ! 3 : Source terms 86 ! 1 : Print plots. 87 ! 2 : Table of 1-D S(f). 88 ! 3 : Table of 1-D time scales. 89 ! 4 : Transfer file. 90 ! 91 ! 4 : Partitioning and bulletins 92 ! 1 : Spectral partitions table 93 ! 2 : Bulletins ASCII format 94 ! 3 : Bulletins CSV format 95 ! 4 : Bulletins CSV & ASCII format 96 ! 3. Parameters : 97 ! 98 ! 4. Subroutines used : 99 ! 100 ! Name Type Module Description 101 ! ---------------------------------------------------------------- 102 ! W3NMOD Subr. W3GDATMD Set number of model. 103 ! W3SETG Subr. Id. Point to selected model. 104 ! W3NDAT Subr. W3WDATMD Set number of model for wave data. 105 ! W3SETW Subr. Id. Point to selected model for wave data. 106 ! W3NAUX Subr. W3ADATMD Set number of model for aux data. 107 ! W3SETA Subr. Id. Point to selected model for aux data. 108 ! W3NOUT Subr. W3ODATMD Set number of model for output. 109 ! W3SETO Subr. Id. Point to selected model for output. 110 ! ITRACE Subr. W3SERVMD Subroutine tracing initialization. 111 ! STRACE Subr. Id. Subroutine tracing. 112 ! NEXTLN Subr. Id. Get next line from input filw 113 ! EXTCDE Subr. Id. Abort program as graceful as possible. 114 ! STME21 Subr. W3TIMEMD Convert time to string. Page 4 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 115 ! TICK21 Subr. Id. Advance time. 116 ! DSEC21 Func. Id. Difference between times. 117 ! W3IOGR Subr. W3IOGRMD Reading/writing model definition file. 118 ! W3IOPO Subr. W3IOPOMD Reading/writing raw point output file. 119 ! W3EXPO Subr. Internal Execute point output. 120 ! W3BULL Subr. W3BULLMD Generate buletins from spectral part. 121 ! ---------------------------------------------------------------- 122 ! 123 ! 5. Called by : 124 ! 125 ! None, stand-alone program. 126 ! 127 ! 6. Error messages : 128 ! 129 ! Checks on input, checks in W3IOxx. 130 ! 131 ! 7. Remarks : 132 ! 133 ! - Tables written to file 'tabNN.ww3', where NN is the 134 ! unit umber (NDSTAB). 135 ! - Transfder file written to ww3.yymmddhh.spc with multiple 136 ! spectra and times in file. yymmddhh relates to first 137 ! output (NDSTAB). 138 ! - !/IC1 !/IC2 !/IC3 are not included in dissipation term 139 ! FIXME: ICE is a dummy variable at the moment 140 ! Include ice parameters in point output file out_pnt.ww3 141 ! Ice coupling to SIN, SDS and SIC similar to w3srcemd.ftn 142 ! 143 ! 8. Structure : 144 ! 145 ! See source code. 146 ! 147 ! 9. Switches : 148 ! 149 ! !/S Enable subroutine tracing. 150 ! 151 ! !/NCO NCEP NCO modifications for operational implementation. 152 ! 153 ! !/O14 Buoy log file generation. 154 ! 155 ! 10. Source code : 156 ! 157 !/ ------------------------------------------------------------------- / 158 USE CONSTANTS 159 !/ 160 ! USE W3GDATMD, ONLY: W3NMOD, W3SETG 161 USE W3WDATMD, ONLY: W3SETW, W3NDAT 162 USE W3ADATMD, ONLY: W3SETA, W3NAUX 163 USE W3ODATMD, ONLY: W3SETO, W3NOUT 164 USE W3IOGRMD, ONLY: W3IOGR 165 USE W3IOPOMD, ONLY: W3IOPO 166 USE W3SERVMD, ONLY : ITRACE, NEXTLN, EXTCDE 167 USE W3TIMEMD, ONLY: STME21, TICK21, DSEC21 168 !/ 169 USE W3GDATMD 170 USE W3WDATMD, ONLY: TIME 171 USE W3ODATMD, ONLY: NDSE, NDST, NDSO, NOPTS, PTLOC, PTNME, & Page 5 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 172 DPO, WAO, WDO, ASO, CAO, CDO, SPCO, FNMPRE 173 USE W3BULLMD, ONLY: NPTAB, NFLD, NPMAX, BHSMIN, BHSDROP, IYY, & 174 HST, TPT, DMT, ASCBLINE, CSVBLINE 175 USE W3BULLMD, ONLY: CASCBLINE 176 USE W3ODATMD, ONLY: GRDID 177 ! 178 IMPLICIT NONE 179 !/ 180 !/ ------------------------------------------------------------------- / 181 !/ Local parameters 182 !/ 183 INTEGER :: NDSI, NDSM, NDSOP, NDSTRC, NTRACE, & 184 IERR, I, TOUT(2), NOUT, TDUM(2), & 185 NREQ, IPOINT, ITYPE, OTYPE, NDSTAB, & 186 IOTEST, IK, ITH, IOUT, J, DIMXP, & 187 NDSBUL, NDSCSV, ICSV, IJ 188 INTEGER :: NDSCBUL 189 INTEGER :: ISCALE = 0 190 INTEGER :: TIMEV(2) 191 INTEGER :: NDBO 192 REAL :: DTREQ, SCALE1, SCALE2, DTEST 193 REAL :: M2KM 194 REAL, ALLOCATABLE :: XPART(:,:) 195 LOGICAL :: FLFORM, FLSRCE(6) 196 LOGICAL, ALLOCATABLE :: FLREQ(:) 197 CHARACTER :: COMSTR*1, IDTIME*23, IDDDAY*11, & 198 TABNME*9, TFNAME*16 199 CHARACTER(LEN=25) :: IDSRCE(6) 200 CHARACTER :: HSTR*6, HTYPE*3 201 !/ 202 !/ ------------------------------------------------------------------- / 203 !/ 204 DATA IDSRCE / 'Spectrum ' , & 205 'Wind-wave interactions ' , & 206 'Nonlinear interactions ' , & 207 'Dissipation ' , & 208 'Wave-bottom interactions ' , & 209 'Sum of selected sources ' / 210 FLSRCE = .FALSE. 211 ! 212 ! CALL W3TAGB('WAVESPEC',1998,0007,0050,'NP21 ') 213 ! 214 !--- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 215 ! 1. IO set-up. 216 ! 217 CALL W3NMOD ( 1, 6, 6 ) 218 CALL W3SETG ( 1, 6, 6 ) 219 CALL W3NDAT ( 6, 6 ) 220 CALL W3SETW ( 1, 6, 6 ) 221 CALL W3NAUX ( 6, 6 ) 222 CALL W3SETA ( 1, 6, 6 ) 223 CALL W3NOUT ( 6, 6 ) 224 CALL W3SETO ( 1, 6, 6 ) 225 ! 226 NDSI = 10 227 NDSM = 20 228 NDSOP = 20 Page 6 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 229 NDSBUL = 0 230 NDSCBUL = 0 231 ! 232 NDSTRC = 6 233 NTRACE = 10 234 CALL ITRACE ( NDSTRC, NTRACE ) 235 236 ! 237 ! Redo according to NCO 238 ! 239 NDSI = 11 240 NDSO = 6 241 NDSE = NDSO 242 NDST = NDSO 243 NDSM = 12 244 NDSOP = 13 245 NDBO = 14 246 NDSTRC = NDSO 247 ! 248 WRITE (NDSO,900) 249 ! 250 J = LEN_TRIM(FNMPRE) 251 OPEN (NDSI,FILE=FNMPRE(:J)//'multiwavespec.inp',STATUS='OLD', & 252 ERR=800,IOSTAT=IERR) 253 READ (NDSI,'(A)',END=801,ERR=802) COMSTR 254 IF (COMSTR.EQ.' ') COMSTR = '$' 255 WRITE (NDSO,901) COMSTR 256 ! 257 !--- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 258 ! 2. Read model definition file. 259 ! 260 CALL W3IOGR ( 'READ', NDSM ) 261 WRITE (NDSO,920) GNAME 262 ! 263 IF ( FLAGLL ) THEN 264 M2KM = 1. 265 ELSE 266 M2KM = 1.E-3 267 END IF 268 ! 269 DIMXP = ((NK+1)/2) * ((NTH-1)/2) 270 ALLOCATE ( XPART(6,0:DIMXP) ) 271 ! 272 !--- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 273 ! 3. Read general data and first fields from file 274 ! 275 CALL W3IOPO ( 'READ', NDSOP, IOTEST ) 276 ! 277 WRITE (NDSO,930) 278 DO I=1, NOPTS 279 IF ( FLAGLL ) THEN 280 WRITE (NDSO,931) PTNME(I), M2KM*PTLOC(1,I), M2KM*PTLOC(2,I) 281 ELSE 282 WRITE (NDSO,932) PTNME(I), M2KM*PTLOC(1,I), M2KM*PTLOC(2,I) 283 END IF 284 END DO 285 ! Page 7 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 286 !--- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 287 ! 4. Read requests from input file. 288 ! Output times 289 ! 290 CALL NEXTLN ( COMSTR , NDSI , NDSE ) 291 READ (NDSI,*,END=801,ERR=802) TOUT, DTREQ, NOUT 292 DTREQ = MAX ( 0. , DTREQ ) 293 IF ( DTREQ.EQ.0 ) NOUT = 1 294 NOUT = MAX ( 1 , NOUT ) 295 ! 296 CALL STME21 ( TOUT , IDTIME ) 297 WRITE (NDSO,940) IDTIME 298 ! 299 TDUM = 0 300 CALL TICK21 ( TDUM , DTREQ ) 301 CALL STME21 ( TDUM , IDTIME ) 302 IF ( DTREQ .GE. 86400. ) THEN 303 WRITE (IDDDAY,'(I10,1X)') INT(DTREQ/86400.) 304 ELSE 305 IDDDAY = ' ' 306 END IF 307 IDTIME(1:11) = IDDDAY 308 IDTIME(21:23) = ' ' 309 WRITE (NDSO,941) IDTIME, NOUT 310 ! 311 ! ... Output points 312 ! 313 ALLOCATE ( FLREQ(NOPTS) ) 314 FLREQ = .FALSE. 315 NREQ = 0 316 ! 317 410 CONTINUE 318 CALL NEXTLN ( COMSTR , NDSI , NDSE ) 319 READ (NDSI,*,END=801,ERR=802) IPOINT 320 IF ( IPOINT .GT. 0 ) THEN 321 IF ( IPOINT .LE. NOPTS ) THEN 322 IF ( .NOT. FLREQ(IPOINT) ) NREQ = NREQ + 1 323 FLREQ(IPOINT) = .TRUE. 324 END IF 325 GOTO 410 326 END IF 327 ! 328 ! ... Output type 329 ! 330 CALL NEXTLN ( COMSTR , NDSI , NDSE ) 331 READ (NDSI,*,END=801,ERR=802) ITYPE 332 ! 333 ! ... ITYPE = 0 334 ! 335 IF ( ITYPE .EQ. 0 ) THEN 336 ! 337 WRITE (NDSO,942) ITYPE, 'Generating buoy log file' 338 OPEN (NDBO,FILE=FNMPRE(:J)//'buoy_log.ww3', & 339 STATUS='NEW',ERR=805,IOSTAT=IERR) 340 DO I = 1,NOPTS 341 WRITE(NDBO,945) I, PTNME(I), PTLOC(1,I), & 342 PTLOC(2,I), GRDID(I) Page 8 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 343 END DO 344 CLOSE(NDBO) 345 ! 346 WRITE (NDSO,942) ITYPE, 'Checking contents of file' 347 DO 348 CALL STME21 ( TIME , IDTIME ) 349 WRITE (NDSO,948) IDTIME 350 CALL W3IOPO ( 'READ', NDSOP, IOTEST ) 351 IF ( IOTEST .EQ. -1 ) THEN 352 WRITE (NDSO,949) 353 GOTO 888 354 END IF 355 END DO 356 ! 357 ! ... ITYPE = 1 358 ! 359 ELSE IF (ITYPE .EQ. 1) THEN 360 WRITE (NDSO,942) ITYPE, '1-D and/or 2-D spectra' 361 CALL NEXTLN ( COMSTR , NDSI , NDSE ) 362 READ (NDSI,*,END=801,ERR=802) OTYPE, SCALE1, SCALE2, & 363 NDSTAB, FLFORM 364 NDSTAB = 51 365 IF (OTYPE .EQ. 1) THEN 366 WRITE (NDSO,943) 'print plots' 367 IF ( SCALE1 .LT. 0. ) THEN 368 WRITE (NDSO,1940) '1-D' 369 ELSE IF ( SCALE1 .EQ. 0. ) THEN 370 WRITE (NDSO,1941) '1-D' 371 ELSE 372 WRITE (NDSO,1942) '1-D', SCALE1 373 END IF 374 IF ( SCALE2 .LT. 0. ) THEN 375 WRITE (NDSO,1940) '2-D' 376 ELSE IF ( SCALE2 .EQ. 0. ) THEN 377 WRITE (NDSO,1941) '2-D' 378 ELSE 379 WRITE (NDSO,1942) '2-D', SCALE2 380 END IF 381 ELSE IF ( OTYPE .EQ. 2 ) THEN 382 WRITE (NDSO,943) 'Table of 1-D spectral data' 383 TABNME = 'tab--.ww3' 384 IF ( NDSTAB.LE.0 .OR. NDSTAB.GT.99 ) NDSTAB = 51 385 WRITE ( TABNME(4:5) , '(I2.2)' ) NDSTAB 386 J = LEN_TRIM(FNMPRE) 387 OPEN (NDSTAB,FILE=FNMPRE(:J)//TABNME,ERR=803,IOSTAT=IERR) 388 WRITE (NDSO,1947) TABNME 389 ELSE IF ( OTYPE .EQ. 3 ) THEN 390 TFNAME = 'ww3.--------.spc' 391 WRITE (TFNAME(5:12),'(I6.6,I2.2)') & 392 MOD(TOUT(1),1000000), TOUT(2)/10000 393 WRITE (NDSO,943) 'Transfer file' 394 IF ( FLFORM ) THEN 395 WRITE (NDSO,1943) TFNAME, 'UNFORMATTED' 396 J = LEN_TRIM(FNMPRE) 397 OPEN (NDSTAB,FILE=FNMPRE(:J)//TFNAME,ERR=804, & 398 IOSTAT=IERR,FORM='UNFORMATTED') 399 WRITE (NDSTAB) 'WAVEWATCH III SPECTRA', & Page 9 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 400 NK, NTH, NREQ, GNAME 401 WRITE (NDSTAB) (SIG(IK)*TPIINV,IK=1,NK) 402 ! 403 ! conversion of directions from trignonmetric to nautical (still uses directions TO ) 404 ! 405 WRITE (NDSTAB) (MOD(2.5*PI-TH(ITH),TPI),ITH=1,NTH) 406 407 ELSE 408 WRITE (NDSO,1943) TFNAME, 'FORMATTED' 409 J = LEN_TRIM(FNMPRE) 410 OPEN (NDSTAB,FILE=FNMPRE(:J)//TFNAME,ERR=804, & 411 IOSTAT=IERR,FORM='FORMATTED') 412 WRITE (NDSTAB,1944) 'WAVEWATCH III SPECTRA', & 413 NK, NTH, NREQ, GNAME 414 WRITE (NDSTAB,1945) (SIG(IK)*TPIINV,IK=1,NK) 415 WRITE (NDSTAB,1946) & 416 (MOD(2.5*PI-TH(ITH),TPI),ITH=1,NTH) 417 END IF 418 ELSE 419 WRITE (NDSE,1011) OTYPE 420 CALL EXTCDE ( 10 ) 421 END IF 422 ! 423 ! ... ITYPE = 2 424 ! 425 ELSE IF (ITYPE .EQ. 2) THEN 426 WRITE (NDSO,942) ITYPE, 'Table of mean wave parameters' 427 CALL NEXTLN ( COMSTR , NDSI , NDSE ) 428 READ (NDSI,*,END=801,ERR=802) OTYPE, NDSTAB 429 NDSTAB = 51 430 TABNME = 'tab--.ww3' 431 IF ( NDSTAB.LE.0 .OR. NDSTAB.GT.99 ) NDSTAB = 51 432 WRITE ( TABNME(4:5) , '(I2.2)' ) NDSTAB 433 J = LEN_TRIM(FNMPRE) 434 OPEN (NDSTAB,FILE=FNMPRE(:J)//TABNME,ERR=803,IOSTAT=IERR) 435 IF ( OTYPE .EQ. 1 ) THEN 436 WRITE (NDSO,2940) 'Depth, current and wind', TABNME 437 ELSE IF ( OTYPE .EQ. 2 ) THEN 438 WRITE (NDSO,2940) 'Mean wave parameters', TABNME 439 ELSE IF ( OTYPE .EQ. 3 ) THEN 440 WRITE (NDSO,2940) 'Nondimensional parameters (U*)', & 441 TABNME 442 ELSE IF ( OTYPE .EQ. 4 ) THEN 443 WRITE (NDSO,2940) 'Nondimensional parameters (U10)', & 444 TABNME 445 ELSE IF ( OTYPE .EQ. 5 ) THEN 446 WRITE (NDSO,2940) 'Validation parameters', TABNME 447 ELSE IF ( OTYPE .EQ. 6 ) THEN 448 WRITE (NDSO,2940) 'WMO standard mean parameters', TABNME 449 ELSE 450 WRITE (NDSE,1011) OTYPE 451 CALL EXTCDE ( 20 ) 452 END IF 453 ! 454 ! ... ITYPE = 3 455 ! 456 ELSE IF (ITYPE .EQ. 3) THEN Page 10 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 457 WRITE (NDSO,942) ITYPE, 'Source terms' 458 CALL NEXTLN ( COMSTR , NDSI , NDSE ) 459 READ (NDSI,*,END=801,ERR=802) OTYPE, SCALE1, SCALE2, & 460 NDSTAB, FLSRCE, ISCALE, FLFORM 461 NDSTAB = 51 462 ISCALE = MAX ( 0 , MIN ( 5 , ISCALE ) ) 463 IF ( OTYPE .EQ. 1 ) THEN 464 WRITE (NDSO,943) 'Print plots' 465 ELSE IF ( OTYPE .EQ. 2 ) THEN 466 IF ( ISCALE .LE. 2) THEN 467 WRITE (NDSO,943) 'Tables as a function of freq.' 468 ELSE 469 WRITE (NDSO,943) 'Tables as a function of f/fp.' 470 END IF 471 IF ( MOD(ISCALE,3) .EQ. 1 ) THEN 472 WRITE (NDSO,944) '(nondimensional based on U10)' 473 ELSE IF ( MOD(ISCALE,3) .EQ. 2) THEN 474 WRITE (NDSO,944) '(nondimensional based on U*)' 475 END IF 476 ELSE IF ( OTYPE .EQ. 3 ) THEN 477 IF ( ISCALE .LE. 2) THEN 478 WRITE (NDSO,943) 'Time scales as a function of freq.' 479 ELSE 480 WRITE (NDSO,943) 'Time scales as a function of f/fp.' 481 END IF 482 IF ( ISCALE .EQ. 1 ) THEN 483 WRITE (NDSO,944) '(nondimensional based on U10)' 484 ELSE IF ( ISCALE .EQ. 2) THEN 485 WRITE (NDSO,944) '(nondimensional based on U*)' 486 END IF 487 ELSE IF ( OTYPE .EQ. 4 ) THEN 488 TFNAME = 'ww3.--------.src' 489 WRITE (TFNAME(5:12),'(I6.6,I2.2)') & 490 MOD(TOUT(1),1000000), TOUT(2)/10000 491 WRITE (NDSO,943) 'Transfer file' 492 IF ( FLFORM ) THEN 493 WRITE (NDSO,3943) TFNAME, 'UNFORMATTED' 494 J = LEN_TRIM(FNMPRE) 495 OPEN (NDSTAB,FILE=FNMPRE(:J)//TFNAME,ERR=804, & 496 IOSTAT=IERR,FORM='UNFORMATTED') 497 WRITE (NDSTAB) 'WAVEWATCH III SOURCES', & 498 NK, NTH, NREQ, FLSRCE 499 WRITE (NDSTAB) (SIG(IK)*TPIINV,IK=1,NK) 500 WRITE (NDSTAB) (MOD(2.5*PI-TH(ITH),TPI),ITH=1,NTH) 501 502 ELSE 503 WRITE (NDSO,3943) TFNAME, 'FORMATTED' 504 J = LEN_TRIM(FNMPRE) 505 OPEN (NDSTAB,FILE=FNMPRE(:J)//TFNAME,ERR=804, & 506 IOSTAT=IERR,FORM='FORMATTED') 507 WRITE (NDSTAB,3944) 'WAVEWATCH III SOURCES', & 508 NK, NTH, NREQ, FLSRCE 509 WRITE (NDSTAB,3945) (SIG(IK)*TPIINV,IK=1,NK) 510 WRITE (NDSTAB,3946) & 511 (MOD(2.5*PI-TH(ITH),TPI),ITH=1,NTH) 512 END IF 513 ELSE Page 11 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 514 WRITE (NDSE,1011) OTYPE 515 CALL EXTCDE ( 30 ) 516 END IF 517 ! 518 DO I=1, 6 519 IF ( FLSRCE(I) ) WRITE (NDSO,3940) IDSRCE(I) 520 END DO 521 WRITE (NDSO,*) ' ' 522 ! 523 IF ( OTYPE .EQ. 1 ) THEN 524 IF ( SCALE1 .LT. 0. ) THEN 525 WRITE (NDSO,1940) '1-D' 526 ELSE IF ( SCALE1 .EQ. 0. ) THEN 527 WRITE (NDSO,1941) '1-D' 528 ELSE 529 WRITE (NDSO,1942) '1-D', SCALE1 530 END IF 531 IF ( SCALE2 .LT. 0. ) THEN 532 WRITE (NDSO,1940) '2-D' 533 ELSE IF ( SCALE2 .EQ. 0. ) THEN 534 WRITE (NDSO,1941) '2-D' 535 ELSE 536 WRITE (NDSO,1942) '2-D', SCALE2 537 END IF 538 END IF 539 ! 540 IF ( OTYPE.EQ.2 .OR. OTYPE.EQ.3 ) THEN 541 TABNME = 'tab--.ww3' 542 IF ( NDSTAB.LE.0 .OR. NDSTAB.GT.99 ) NDSTAB = 51 543 WRITE ( TABNME(4:5) , '(I2.2)' ) NDSTAB 544 J = LEN_TRIM(FNMPRE) 545 OPEN (NDSTAB,FILE=FNMPRE(:J)//TABNME,ERR=803,IOSTAT=IERR) 546 WRITE (NDSO,3941) TABNME 547 END IF 548 ! 549 ! ... ITYPE = 4 550 ! 551 ELSE IF (ITYPE .EQ. 4) THEN 552 WRITE (NDSO,942) ITYPE, 'Spectral partitions or bulletins' 553 CALL NEXTLN ( COMSTR , NDSI , NDSE ) 554 READ (NDSI,*,END=801,ERR=802) OTYPE, NDSTAB, TIMEV, HTYPE 555 NDSTAB = 51 556 IF ( OTYPE .EQ. 1 ) THEN 557 WRITE (NDSO,943) 'Partitioning of spectra' 558 TABNME = 'tab--.ww3' 559 IF ( NDSTAB.LE.0 .OR. NDSTAB.GT.99 ) NDSTAB = 51 560 WRITE ( TABNME(4:5) , '(I2.2)' ) NDSTAB 561 J = LEN_TRIM(FNMPRE) 562 OPEN (NDSTAB,FILE=FNMPRE(:J)//TABNME,ERR=803,IOSTAT=IERR) 563 WRITE (NDSO,1947) TABNME 564 565 ELSEIF ( OTYPE .GE. 2 ) THEN 566 IF (OTYPE .EQ. 2 .OR. OTYPE .EQ. 4 ) THEN 567 WRITE (NDSO,943) 'Bulletins, ASCII format' 568 J = LEN_TRIM(FNMPRE) 569 DO IJ = 1,NOPTS 570 IF ( COUNT(FLREQ .EQV. .TRUE. ) .GT. 1 ) THEN Page 12 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 571 ! ... This version only allows single point output for bulletins 572 WRITE (NDSE,1012) OTYPE 573 CALL EXTCDE ( 45 ) 574 ENDIF 575 IF (FLREQ(IJ)) THEN 576 NDSBUL = NDSTAB + (IJ - 1) 577 OPEN(NDSBUL,FILE=TRIM(PTNME(IJ))//'.bull',ERR=803,IOSTAT=IERR) 578 WRITE (NDSO,1947) TRIM(PTNME(IJ))//'.bull' 579 NDSCBUL = NDSTAB + (IJ - 1) + NOPTS 580 OPEN(NDSCBUL,FILE=TRIM(PTNME(IJ))//'.cbull',ERR=803,IOSTAT=IERR) 581 WRITE (NDSO,1947) TRIM(PTNME(IJ))//'.cbull' 582 ENDIF 583 ENDDO 584 ENDIF 585 IF ( OTYPE .EQ. 3 .OR. OTYPE .EQ. 4 ) THEN 586 WRITE (NDSO,943) 'Bulletins, CSV format' 587 J = LEN_TRIM(FNMPRE) 588 DO IJ = 1,NOPTS 589 IF (FLREQ(IJ)) THEN 590 ICSV = 0 591 IF ( NDSBUL .GT. 0 ) ICSV = NDSBUL 592 IF ( NDSCBUL .GT. 0 ) ICSV = NDSCBUL 593 NDSCSV = NDSTAB + (IJ - 1) + ICSV 594 OPEN(NDSCSV,FILE=TRIM(PTNME(IJ))//'.csv',ERR=803,IOSTAT=IERR) 595 WRITE (NDSO,1947) TRIM(PTNME(IJ))//'.csv' 596 ENDIF 597 ENDDO 598 ENDIF 599 ELSE 600 WRITE (NDSE,1011) OTYPE 601 CALL EXTCDE ( 50 ) 602 END IF 603 ! 604 ! ... ITYPE = ILLEGAL 605 ! 606 ELSE 607 WRITE (NDSE,1010) ITYPE 608 CALL EXTCDE ( 1 ) 609 ! 610 END IF 611 ! 612 ! ... Output of output points 613 ! 614 WRITE (NDSO,950) NREQ 615 DO I=1, NOPTS 616 IF (FLREQ(I)) THEN 617 IF ( FLAGLL ) THEN 618 WRITE (NDSO,951) PTNME(I), M2KM*PTLOC(1,I), & 619 M2KM*PTLOC(2,I) 620 ELSE 621 WRITE (NDSO,953) PTNME(I), M2KM*PTLOC(1,I), & 622 M2KM*PTLOC(2,I) 623 END IF 624 END IF 625 END DO 626 ! 627 IF ( ITYPE.EQ.3 .AND. OTYPE.EQ.4 ) WRITE (NDSO,952) Page 13 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 628 ! 629 !--- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 630 ! 5. Time management. 631 ! 632 IOUT = 0 633 ! 634 635 DO 636 DTEST = DSEC21 ( TIME , TOUT ) 637 IF ( DTEST .GT. 0. ) THEN 638 CALL W3IOPO ( 'READ', NDSOP, IOTEST ) 639 IF ( IOTEST .EQ. -1 ) THEN 640 WRITE (NDSO,949) 641 EXIT 642 END IF 643 CYCLE 644 END IF 645 IF ( DTEST .LT. 0. ) THEN 646 CALL TICK21 ( TOUT , DTREQ ) 647 CYCLE 648 END IF 649 ! 650 IOUT = IOUT + 1 651 CALL STME21 ( TOUT , IDTIME ) 652 IF ( ( ITYPE.EQ.1 .AND. OTYPE.EQ.1 ) .OR. & 653 ( ITYPE.EQ.3 .AND. OTYPE.EQ.1 ) & 654 ) WRITE (NDSO,960) IDTIME 655 CALL W3EXPO 656 CALL TICK21 ( TOUT , DTREQ ) 657 IF ( IOUT .GE. NOUT ) EXIT 658 END DO 659 ! 660 ! ... ITYPE=4 & OTYPES=[2,4] requires adding lines at bottom of 661 ! bulletin output for compatibility with version 2.22 662 ! 663 IF (ITYPE .EQ. 4 .AND. ( OTYPE .EQ. 2 .OR. OTYPE .EQ. 4 ) ) THEN 664 DO IJ = 1,NOPTS 665 IF (FLREQ(IJ)) THEN 666 NDSBUL = NDSTAB + (IJ - 1) 667 WRITE(NDSBUL,971) 668 WRITE(NDSBUL,974) BHSDROP, BHSMIN 669 NDSCBUL = NDSTAB + (IJ - 1) + NOPTS 670 WRITE(NDSCBUL,961) 671 WRITE(NDSCBUL,962) 672 ENDIF 673 ENDDO 674 ENDIF 675 ! 676 GOTO 888 677 ! 678 ! Escape locations read errors : 679 ! 680 800 CONTINUE 681 WRITE (NDSE,1000) IERR 682 CALL EXTCDE ( 40 ) 683 ! 684 801 CONTINUE Page 14 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 685 WRITE (NDSE,1001) 686 CALL EXTCDE ( 41 ) 687 ! 688 802 CONTINUE 689 WRITE (NDSE,1002) IERR 690 CALL EXTCDE ( 42 ) 691 ! 692 803 CONTINUE 693 WRITE (NDSE,1003) IERR 694 CALL EXTCDE ( 43 ) 695 ! 696 804 CONTINUE 697 WRITE (NDSE,1004) IERR 698 CALL EXTCDE ( 44 ) 699 ! 700 805 CONTINUE 701 WRITE (NDSE,1005) IERR 702 CALL EXTCDE ( 45 ) 703 ! 704 888 CONTINUE 705 ! 706 WRITE (NDSO,999) 707 ! 708 ! CALL W3TAGE('WAVESPEC') 709 ! 710 ! Formats 711 ! 712 900 FORMAT (/15X,' *** WAVEWATCH III Point output post.*** '/ & 713 15X,'==============================================='/) 714 901 FORMAT ( ' Comment character is ''',A,''''/) 715 ! 716 920 FORMAT ( ' Grid name : ',A/) 717 ! 718 930 FORMAT ( ' Points in file : '/ & 719 ' ------------------------------------') 720 931 FORMAT ( ' ',A,2F10.2) 721 932 FORMAT ( ' ',A,2(F8.1,'E3')) 722 ! 723 940 FORMAT (/' Output time data : '/ & 724 ' --------------------------------------------------'/ & 725 ' First time : ',A) 726 941 FORMAT ( ' Interval : ',A/ & 727 ' Number of requests : ',I4) 728 942 FORMAT (/' Output type ',I2,' :'/ & 729 ' --------------------------------------------------'/ & 730 ' ',A/) 731 943 FORMAT ( ' Subtype : ',A) 732 944 FORMAT ( ' ',A) 733 945 FORMAT ( ' ',I5,3X,A,2F10.2,3X,A) 734 948 FORMAT ( ' Data for ',A) 735 949 FORMAT (/' End of file reached '/) 736 ! 737 950 FORMAT (/' Requested output for',I3,' points : '/ & 738 ' --------------------------------------------------') 739 951 FORMAT ( ' ',A,2F10.2) 740 953 FORMAT ( ' ',A,2(F8.1,'E3')) 741 952 FORMAT (/' Output times :'/ & Page 15 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 742 ' --------------------------------------------------') 743 961 FORMAT ('----------------------------------------', & 744 '---------------------------') 745 962 FORMAT ( 'DD = Day of Month'/ & 746 'HH = Hour of Day'/ & 747 'HS = Total Significant Wave Height (feet)'/ & 748 'SS = Significant Wave Height of separate system (feet)'/ & 749 'PP = Peak Period of separate system (whole seconds)'/ & 750 'DDD = Mean Direction of separate system (degrees,"from")') 751 971 FORMAT (' +-------+-----------+-----------------+', & 752 '-----------------+-----------------+----', & 753 '-------------+-----------------+--------', & 754 '---------+')! 755 974 FORMAT ( & 756 75X,'Hst : Total sigificant wave height.'/ & 757 75X,'n : Number of fields with Hs > ',f4.2, & ................................................1 (1) Recommended relationship between field width 'W' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+3'. 758 ' in 2-D spectrum.'/ & 759 75X,'x : Number of fields with Hs > ',f4.2, & ................................................1 (1) Recommended relationship between field width 'W' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+3'. 760 ' not in table.'/ & 761 75X,'Hs : Significant wave height of separate wave field.'/ & 762 75X,'Tp : Peak period of separate wave field.'/ & 763 75X,'dir : Mean direction of separate wave field.'/ & 764 75X,'* : Wave generation due to local wind probable.') 765 766 1940 FORMAT ( ' ',A,' print plots not requested.') 767 1941 FORMAT ( ' ',A,' print plots normalized.') 768 1942 FORMAT ( ' Scale factor ',A,' spectrum : ',E10.3) 769 1943 FORMAT ( ' File name : ',A,' (',A,')') 770 1944 FORMAT ('''',A,'''',1X,3I6,1X,'''',A,'''') 771 1945 FORMAT (8E10.3) 772 1946 FORMAT (7E11.3) 773 1947 FORMAT ( ' File name : ',A) 774 ! 775 2940 FORMAT ( ' Table output : ',A/ & 776 ' File name : ',A) 777 ! 778 3940 FORMAT ( ' ',A) 779 3941 FORMAT ( ' File name : ',A) 780 3943 FORMAT ( ' File name : ',A,' (',A,')') 781 3944 FORMAT ('''',A,'''',1X,3I6,6L2) 782 3945 FORMAT (8E10.3) 783 3946 FORMAT (7E11.3) 784 ! 785 960 FORMAT (//' Output for ',A/ & 786 ' --------------------------------------------------') 787 ! 788 999 FORMAT (/' End of program '/ & 789 ' ========================================='/ & 790 ' WAVEWATCH III Point output '/) 791 ! 792 1000 FORMAT (/' *** WAVEWATCH III ERROR IN WAVESPEC : '/ & Page 16 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 793 ' ERROR IN OPENING INPUT FILE'/ & 794 ' IOSTAT =',I5/) 795 ! 796 1001 FORMAT (/' *** WAVEWATCH III ERROR IN WAVESPEC : '/ & 797 ' PREMATURE END OF INPUT FILE'/) 798 ! 799 1002 FORMAT (/' *** WAVEWATCH III ERROR IN WAVESPEC : '/ & 800 ' ERROR IN READING FROM INPUT FILE'/ & 801 ' IOSTAT =',I5/) 802 ! 803 1003 FORMAT (/' *** WAVEWATCH III ERROR IN WAVESPEC : '/ & 804 ' ERROR IN OPENING TABLE FILE'/ & 805 ' IOSTAT =',I5/) 806 ! 807 1004 FORMAT (/' *** WAVEWATCH III ERROR IN WAVESPEC : '/ & 808 ' ERROR IN OPENING IDL FILE'/ & 809 ' IOSTAT =',I5/) 810 ! 811 1005 FORMAT (/' *** WAVEWATCH III ERROR IN WAVESPEC : '/ & 812 ' ERROR IN OPENING BUOY LOG FILE'/ & 813 ' IOSTAT =',I5/) 814 ! 815 1010 FORMAT (/' *** WAVEWATCH III ERROR IN WAVESPEC : '/ & 816 ' ILLEGAL TYPE, ITYPE =',I4/) 817 ! 818 1011 FORMAT (/' *** WAVEWATCH III ERROR IN WAVESPEC : '/ & 819 ' ILLEGAL TYPE, OTYPE =',I4/) 820 ! 821 1012 FORMAT (/' *** WAVEWATCH III ERROR IN WAVESPEC : '/ & 822 ' MULTIPLE OUTPUT POINTS DEFINED, ITYPE =',I4,/ & 823 ' ONLY SINGLE POINT ALLOWED IN THIS VERSION'/) 824 ! 825 !/ 826 !/ Internal subroutine W3EXPO ---------------------------------------- / 827 !/ 828 CONTAINS 829 !/ ------------------------------------------------------------------- / 830 SUBROUTINE W3EXPO 831 !/ 832 !/ +-----------------------------------+ 833 !/ | WAVEWATCH III NOAA/NCEP | 834 !/ | H. L. Tolman | 835 !/ | J.H. Alves | 836 !/ | F. Ardhuin | 837 !/ | A. Chawla | 838 !/ | FORTRAN 90 | 839 !/ | Last update : 26-Dec-2012 | 840 !/ +-----------------------------------+ 841 !/ 842 !/ 08-Jun-1999 : Final FORTRAN 77 ( version 1.18 ) 843 !/ 21-Jan-2000 : Upgrade to FORTRAN 90 ( version 2.00 ) 844 !/ Massive changes to logistics 845 !/ 09-Jan-2001 : U* bug fix in tabular output ( version 2.05 ) 846 !/ 25-Jan-2001 : Flat grid version ( version 2.06 ) 847 !/ 02-Feb-2001 : Xnl version 3.0 ( version 2.07 ) 848 !/ 11-Jun-2001 : Clean up ( version 2.11 ) 849 !/ 11-Oct-2001 : Clean up, X*, Y* in tables ( version 2.14 ) Page 17 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 850 !/ 24-Dec-2004 : Multiple grid version. ( version 3.06 ) 851 !/ 17-Apr-2006 : Filter for directional spread. ( version 3.09 ) 852 !/ 23-Jun-2006 : Linear input added. ( version 3.09 ) 853 !/ 03-Jul-2006 : Separate flux modules. ( version 3.09 ) 854 !/ 28-Oct-2006 : Add partitioning option. ( version 3.10 ) 855 !/ 24-Mar-2007 : Add pars for entire spectrum. ( version 3.11 ) 856 !/ 25-Apr-2007 : Battjes-Janssen Sdb added. ( version 3.11 ) 857 !/ (J. H. Alves) 858 !/ 09-Oct-2007 : WAM 4+ Sin and Sds added. ( version 3.13 ) 859 !/ (F. Ardhuin) 860 !/ 09-Oct-2007 : Experimental Sbs (BS1) added. ( version 3.13 ) 861 !/ (F. Ardhuin) 862 !/ 09-Apr-2008 : Adding an additional output for ( version 3.12 ) 863 !/ WMO standard (A. Chawla) 864 !/ 29-Apr-2008 : Adjust format partition output. ( version 3.14 ) 865 !/ 01-Jul-2011 : Adding BT4 ( version 4.01 ) 866 !/ 16-Jul-2012 : Move GMD (SNL3) and nonlinear filter (SNLS) 867 !/ from 3.15 (HLT). ( version 4.08 ) 868 !/ 26-Dec-2012 : Modified obsolete declarations. ( version 4.11 ) 869 !/ 870 ! 1. Purpose : 871 ! 872 ! Perform actual point output. 873 ! 874 ! 3. Parameters : 875 ! 876 ! 4. Subroutines used : 877 ! 878 ! Name Type Module Description 879 ! ---------------------------------------------------------------- 880 ! W3SPRn Subr. W3SRCnMD Mean wave parameters for use in 881 ! source terms. 882 ! W3FLXn Subr. W3FLXnMD Flux/stress computation. 883 ! W3SLNn Subr. W3SLNnMD Linear input. 884 ! W3SINn Subr. W3SRCnMD Input source term. 885 ! W3SDSn Subr. W3SRCnMD Whitecapping source term 886 ! W3SNLn Subr. W3SNLnMD Nonlinear interactions. 887 ! W3SBTn Subr. W3SBTnMD Bottom friction source term. 888 ! W3SDBn Subr. W3SBTnMD Depth induced breaking source term. 889 ! W3STRn Subr. W3STRnMD Triad interaction source term. 890 ! W3SBSn Subr. W3SBSnMD Bottom scattering source term. 891 ! W3SXXn Subr. W3SXXnMD Unclassified source term. 892 ! W3PART Sunr. W3PARTMD Spectral partitioning routine. 893 ! STRACE Subr. W3SERVMD Subroutine tracing. 894 ! STME21 Subr. W3TIMEMD Convert time to string. 895 ! PRT1DS Subr. W3ARRYMD Print plot of 1-D spectrum. 896 ! PRT1DM Subr. Id. Print plot of several 1-D spectra. 897 ! PRT2DS Subr. Id. Print plot of 2-D spectrum. 898 ! WAVNU1 Subr. W3DISPMD Solve dispersion relation. 899 ! ---------------------------------------------------------------- 900 ! 901 ! 5. Called by : 902 ! 903 ! Main program in which it is contained, 904 ! 905 ! 6. Error messages : 906 ! Page 18 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 907 ! None. 908 ! 909 ! 7. Remarks : 910 ! 911 ! - Spectra are relative frequency energy spectra. 912 ! - Note that arrays CX and CY of the main program now contain 913 ! the absolute current speed and direction respectively. 914 ! 915 ! 8. Structure : 916 ! 917 ! See source code. 918 ! 919 ! 9. Switches : 920 ! 921 ! !/S Enable subroutine tracing. 922 ! !/T Enable test output. 923 ! 924 ! !/FLXx Flux/stress computation. 925 ! !/LNx Linear input package 926 ! !/STx Source term package 927 ! !/NLx Nonlinear interaction package 928 ! !/BTx Bottom friction package 929 ! !/ICx S_ice source term package 930 ! !/DBx Depth-induced breaking package 931 ! !/TRx Triad interaction package 932 ! !/BSx Bottom scattering package 933 ! !/XXx Arbitrary adittional source term package 934 ! 935 ! !/STAB2 Stability correction for !/ST2 936 ! 937 ! 10. Source code : 938 ! 939 !/ ------------------------------------------------------------------- / 940 USE W3SLN1MD 941 USE W3SRC4MD 942 USE W3SNL1MD 943 USE W3SBT1MD 944 USE W3SDB1MD 945 USE W3PARTMD, ONLY: W3PART 946 USE W3DISPMD, ONLY: WAVNU1 947 !/ 948 USE W3ARRYMD, ONLY: PRT1DS, PRT2DS, PRT1DM 949 USE W3DISPMD, ONLY: NAR1D, DFAC, N1MAX, ECG1, EWN1, DSIE 950 USE W3BULLMD, ONLY: W3BULL 951 !/ 952 !/ ------------------------------------------------------------------- / 953 !/ Local parameters 954 !/ 955 INTEGER :: J, I1, I2, ISP, IKM, IKL, IKH, ITH, & 956 IK, IH, IM, IS, IYR, IMTH, IDY, ITT, & 957 I, NPART, IP, IX, IY, ISEA 958 INTEGER, SAVE :: IPASS = 0 959 REAL :: DEPTH, SQRTH, CDIR, SIX, R1, R2, & 960 UDIR, UDIRR, UABS, XL, XH, XL2, XH2, & 961 ET, EWN, ETR, ETX, ETY, EBND, EBX, & 962 EBY, HSIG, WLEN, TMEAN, THMEAN, & 963 THSPRD, EMAX, EL, EH, DENOM, FP, THP,& Page 19 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 964 SPP, CD, USTAR, FACTOR, UNORM, ESTAR,& 965 FPSTAR, FACF, FACE, FACS, HMAT, WNA, & 966 XYZ, AGE1, AFR, AGE2, FACT, XSTAR, & 967 YSTAR, FHIGH, ZWND, Z0, USTD, EMEAN, & 968 FMEAN, WNMEAN, UDIRCA, X, Y, CHARN, M2KM 969 REAL :: AMAX, FMEANS, FMEANWS, TAUWX, TAUWY, & 970 TAUWNX, TAUWNY, FMEAN1, WHITECAP(1:4) 971 REAL :: ICE 972 REAL, SAVE :: HSMIN = 0.05 973 REAL :: WN(NK), CG(NK) 974 REAL :: E(NK,NTH), E1(NK), APM(NK), & 975 THBND(NK), SPBND(NK), A(NTH,NK), & 976 WN2(NTH,NK) 977 REAL :: DIA(NTH,NK), SWN(NK,NTH), SNL(NK,NTH),& 978 SDS(NK,NTH), SBT(NK,NTH), XLN(NTH,NK),& 979 XIN(NTH,NK), XNL(NTH,NK), XTR(NTH,NK),& 980 XDS(NTH,NK), XDB(NTH,NK), XBT(NTH,NK),& 981 XBS(NTH,NK), XXX(NTH,NK), STT(NK,NTH),& 982 XWL(NTH,NK) 983 REAL :: SIN1(NK), SNL1(NK), SDS1(NK), & 984 SBT1(NK), STT1(NK), E1ALL(NK,5) 985 LOGICAL :: LLWS(NSPEC) 986 CHARACTER :: DTME21*23 987 CHARACTER(LEN=4) VAR1(5) 988 CHARACTER(LEN=1) IDLAT, IDLON 989 ! 990 DATA VAR1 / 'Sin ' , 'Snl ', 'Sds ' , 'Sbt ' , 'Stot' / 991 !/ 992 !/ ------------------------------------------------------------------- / 993 !/ 994 ! 1. Initialisations 995 ! 996 IF ( FLAGLL ) THEN 997 M2KM = 1. 998 ELSE 999 M2KM = 1.E-3 1000 END IF 1001 ! 1002 XL = 1./XFR - 1. 1003 XH = XFR - 1. 1004 XL2 = XL**2 1005 XH2 = XH**2 1006 IPASS = IPASS + 1 1007 ! 1008 IF ( ITYPE .EQ. 3 ) THEN 1009 XLN = 0. 1010 XIN = 0. 1011 XNL = 0. 1012 XTR = 0. 1013 XDS = 0. 1014 XDB = 0. 1015 XBT = 0. 1016 XBS = 0. 1017 XWL = 0. 1018 XXX = 0. 1019 END IF 1020 ! Page 20 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1021 ! Output of time 1022 ! 1023 IF ( ( ITYPE.EQ.1 .AND. OTYPE.EQ.3 ) .OR. & 1024 ( ITYPE.EQ.3 .AND. OTYPE.EQ.4 ) ) THEN 1025 IF ( FLFORM ) THEN 1026 WRITE (NDSTAB) TIME 1027 ELSE 1028 WRITE (NDSTAB,900) TIME 1029 END IF 1030 END IF 1031 ! 1032 IF (ITYPE.EQ.2) THEN 1033 IF ( NREQ.EQ.1 .AND. IPASS.EQ.1 ) THEN 1034 IF ( OTYPE.EQ.1 ) THEN 1035 WRITE (NDSTAB,1901) 1036 ELSE IF ( OTYPE.EQ.2 ) THEN 1037 WRITE (NDSTAB,1902) 1038 ELSE IF ( OTYPE.EQ.3 ) THEN 1039 WRITE (NDSTAB,1903) 1040 ELSE IF ( OTYPE.EQ.4 ) THEN 1041 WRITE (NDSTAB,1904) 1042 ELSE IF ( OTYPE.EQ.5 ) THEN 1043 WRITE (NDSTAB,1905) 1044 ELSE IF ( OTYPE.EQ.6 ) THEN 1045 WRITE (NDSTAB,1906) 1046 END IF 1047 END IF 1048 IF ( NREQ.NE.1 ) THEN 1049 CALL STME21 ( TIME , DTME21 ) 1050 IF ( IPASS .NE. 1 ) WRITE (NDSTAB,1910) 1051 IF ( OTYPE.EQ.1 ) THEN 1052 IF ( FLAGLL ) THEN 1053 WRITE (NDSTAB,1911) DTME21 1054 ELSE 1055 WRITE (NDSTAB,1711) DTME21 1056 END IF 1057 ELSE IF ( OTYPE.EQ.2 ) THEN 1058 IF ( FLAGLL ) THEN 1059 WRITE (NDSTAB,1912) DTME21 1060 ELSE 1061 WRITE (NDSTAB,1712) DTME21 1062 END IF 1063 ELSE IF ( OTYPE.EQ.3 ) THEN 1064 WRITE (NDSTAB,1913) DTME21 1065 ELSE IF ( OTYPE.EQ.4 ) THEN 1066 WRITE (NDSTAB,1914) DTME21 1067 ELSE IF ( OTYPE.EQ.5 ) THEN 1068 IF ( FLAGLL ) THEN 1069 WRITE (NDSTAB,1915) DTME21 1070 ELSE 1071 WRITE (NDSTAB,1715) DTME21 1072 END IF 1073 ELSE IF ( OTYPE.EQ.6 ) THEN 1074 IF ( FLAGLL ) THEN 1075 WRITE (NDSTAB,1916) DTME21 1076 ELSE 1077 WRITE (NDSTAB,1716) DTME21 Page 21 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1078 END IF 1079 END IF 1080 END IF 1081 END IF 1082 ! 1083 IF (ITYPE.EQ.3) THEN 1084 IF ( OTYPE .EQ. 4 ) THEN 1085 CALL STME21 ( TIME , DTME21 ) 1086 WRITE (NDSO,905) DTME21 1087 END IF 1088 END IF 1089 ! 1090 !--- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 1091 ! Loop over output points. 1092 ! 1093 DO J=1, NOPTS 1094 IF ( FLREQ(J) ) THEN 1095 ! 1096 ! 2. Calculate grid parameters using and inlined version of WAVNU1. 1097 ! 1098 DEPTH = MAX ( DMIN, DPO(J) ) 1099 SQRTH = SQRT ( DEPTH ) 1100 UDIR = MOD ( 270. - WDO(J)*RADE , 360. ) 1101 UDIRCA = WDO(J)*RADE 1102 UDIRR = WDO(J) 1103 UABS = MAX ( 0.001 , WAO(J) ) 1104 CDIR = MOD ( 270. - CDO(J)*RADE , 360. ) 1105 ! 1106 DO IK=1, NK 1107 SIX = SIG(IK) * SQRTH 1108 I1 = INT(SIX/DSIE) 1109 IF (I1.LE.N1MAX) THEN 1110 I2 = I1 + 1 1111 R1 = SIX/DSIE - REAL(I1) 1112 R2 = 1. - R1 1113 WN(IK) = ( R2*EWN1(I1) + R1*EWN1(I2) ) / DEPTH 1114 CG(IK) = ( R2*ECG1(I1) + R1*ECG1(I2) ) * SQRTH 1115 ELSE 1116 WN(IK) = SIG(IK)*SIG(IK)/GRAV 1117 CG(IK) = 0.5 * GRAV / SIG(IK) 1118 END IF 1119 ! 1120 END DO 1121 ! 1122 ! 3. Prepare spectra etc. 1123 ! 3.a Mean wave parameters. 1124 ! 1125 ET = 0. 1126 EWN = 0. 1127 ETR = 0. 1128 ETX = 0. 1129 ETY = 0. 1130 DO IK=1, NK 1131 EBND = 0. 1132 EBX = 0. 1133 EBY = 0. 1134 DO ITH=1, NTH Page 22 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1135 ISP = ITH + (IK-1)*NTH 1136 E(IK,ITH) = SPCO(ISP,J) 1137 EBND = EBND + SPCO(ISP,J) 1138 EBX = EBX + SPCO(ISP,J)*ECOS(ITH) 1139 EBY = EBY + SPCO(ISP,J)*ESIN(ITH) 1140 END DO 1141 E1(IK) = EBND * DTH 1142 APM(IK)= E1(IK) / ( TPI * GRAV**2 / SIG(IK)**5 ) 1143 IF ( E1(IK) .GT. 1.E-5) THEN 1144 THBND(IK) = MOD(630.- RADE*ATAN2(EBY,EBX),360.) 1145 SPBND(IK) = RADE * SQRT ( MAX ( 0. , 2.*( 1. - & 1146 SQRT( MAX(0.,(EBX**2+EBY**2)/EBND**2) ) ) ) ) 1147 ELSE 1148 THBND(IK) = -999.9 1149 SPBND(IK) = -999.9 1150 END IF 1151 EBND = E1(IK) * DSII(IK) * TPIINV 1152 ET = ET + EBND 1153 EWN = EWN + EBND / WN(IK) 1154 ETR = ETR + EBND / SIG(IK) 1155 ETX = ETX + EBX * DSII(IK) 1156 ETY = ETY + EBY * DSII(IK) 1157 END DO 1158 ! 1159 ! tail factors for radian action etc ...! 1160 ! 1161 EBND = E1(NK) * TPIINV / ( SIG(NK) * DTH ) 1162 ET = ET + FTE *EBND 1163 EWN = EWN + FTWL*EBND 1164 ETR = ETR + FTTR*EBND 1165 ETX = DTH*ETX*TPIINV + FTE*EBX*TPIINV/SIG(NK) 1166 ETY = DTH*ETY*TPIINV + FTE*EBY*TPIINV/SIG(NK) 1167 ! 1168 HSIG = 4. * SQRT ( ET ) 1169 IF ( HSIG .GT. HSMIN ) THEN 1170 WLEN = EWN / ET * TPI 1171 TMEAN = ETR / ET * TPI 1172 THMEAN = MOD ( 630. - RADE*ATAN2(ETY,ETX) , 360. ) 1173 THSPRD = RADE * SQRT ( MAX ( 0. , 2.*( 1. - SQRT( & 1174 MAX(0.,(ETX**2+ETY**2)/ET**2) ) ) ) ) 1175 IF ( THSPRD .LT. 0.01*RADE*DTH ) THSPRD = 0. 1176 ELSE 1177 WLEN = 0. 1178 TMEAN = 0. 1179 THMEAN = 0. 1180 THSPRD = 0. 1181 DO IK=1, NK 1182 E1(IK) = 0. 1183 DO ITH=1, NTH 1184 E(IK,ITH) = 0. 1185 END DO 1186 END DO 1187 END IF 1188 ! 1189 ! peak frequency 1190 ! 1191 EMAX = E1(NK) Page 23 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1192 IKM = NK 1193 ! 1194 DO IK=NK-1, 1, -1 1195 IF ( E1(IK) .GT. EMAX ) THEN 1196 EMAX = E1(IK) 1197 IKM = IK 1198 END IF 1199 END DO 1200 ! 1201 IKL = MAX ( 1 , IKM-1 ) 1202 IKH = MIN ( NK , IKM+1 ) 1203 EL = E1(IKL) - E1(IKM) 1204 EH = E1(IKH) - E1(IKM) 1205 DENOM = XL*EH - XH*EL 1206 ! 1207 IF ( HSIG .GE. HSMIN ) THEN 1208 FP = SIG(IKM) * ( 1. + 0.5 * ( XL2*EH - XH2*EL ) & 1209 / SIGN ( MAX(ABS(DENOM),1.E-15) , DENOM ) ) 1210 THP = THBND(IKM) 1211 SPP = SPBND(IKM) 1212 IF ( SPP .LT. 0.01*RADE*DTH ) SPP = 0. 1213 ELSE 1214 FP = 0. 1215 THP = 0. 1216 SPP = 0. 1217 END IF 1218 ! 1219 ! spectral partitioning 1220 ! 1221 IF ( ITYPE.EQ.4 ) CALL W3PART & 1222 ( E, UABS, UDIRCA, DEPTH, WN, NPART, XPART, DIMXP ) 1223 ! 1224 ! nondimensional parameters 1225 ! 1226 IF ( ( ITYPE.EQ.2 .AND. (OTYPE.EQ.3.OR.OTYPE.EQ.4) ) .OR. & 1227 ( ITYPE.EQ.1 .AND. (OTYPE.EQ.2) ) ) THEN 1228 ! 1229 DO IK=1, NK 1230 FACTOR = TPIINV * CG(IK) / SIG(IK) 1231 DO ITH=1, NTH 1232 ISP = ITH + (IK-1)*NTH 1233 A(ITH,IK) = FACTOR * SPCO(ISP,J) 1234 WN2(ITH,IK) = WN(IK) 1235 END DO 1236 END DO 1237 ! 1238 LLWS(:) = .TRUE. 1239 ZWND = ZZWND 1240 TAUWX = 0. 1241 TAUWY = 0. 1242 USTAR = 1. 1243 ! 1244 CALL W3SPR4 (A, CG, WN, EMEAN, FMEAN, FMEAN1, & 1245 WNMEAN, AMAX, UABS, UDIRR, USTAR, USTD,& 1246 TAUWX, TAUWY, CD, Z0, CHARN, LLWS, FMEANWS ) 1247 ! 1248 DO ITT=1, 3 Page 24 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1249 IX=1 1250 IY=1 1251 CALL W3SIN4 ( A, CG, WN2, UABS, USTAR, DAIR/DWAT,& 1252 ASO(J), UDIRR, Z0, CD, TAUWX, TAUWY,& 1253 TAUWNX, TAUWNY, ICE, XIN, DIA, LLWS, IX, IY ) 1254 CALL W3SPR4 (A, CG, WN, EMEAN, FMEAN, FMEAN1, & 1255 WNMEAN, AMAX, UABS, UDIRR, USTAR, USTD,& 1256 TAUWX, TAUWY, CD, Z0, CHARN, LLWS, FMEANWS ) 1257 END DO 1258 ! 1259 ! Add alternative flux calculations here as part of !/ST2 option .... 1260 ! Also add before actual source term calculation !!! 1261 ! 1262 IF ( WAO(J) .LT. 0.01 ) THEN 1263 UNORM = 0. 1264 ESTAR = 0. 1265 FPSTAR = 0. 1266 ELSE 1267 IF ( OTYPE.EQ.3 ) THEN 1268 UNORM = USTAR 1269 ELSE 1270 UNORM = WAO(J) 1271 END IF 1272 ESTAR = ET * GRAV**2 / UNORM**4 1273 FPSTAR = FP * TPIINV * UNORM / GRAV 1274 XSTAR = PTLOC(1,J) * GRAV / UNORM**2 1275 YSTAR = PTLOC(2,J) * GRAV / UNORM**2 1276 IF ( FLAGLL ) THEN 1277 XSTAR = XSTAR * DERA * RADIUS & 1278 * COS(PTLOC(2,J)*DERA) 1279 YSTAR = YSTAR * DERA * RADIUS 1280 END IF 1281 END IF 1282 ! 1283 END IF 1284 ! 1285 ! 3.4 source terms 1286 ! 1287 IF ( ITYPE.EQ.3 ) THEN 1288 ! 1289 DO IK=1, NK 1290 FACTOR = TPIINV * CG(IK) / SIG(IK) 1291 DO ITH=1, NTH 1292 A(ITH,IK) = FACTOR * SPCO(ITH+(IK-1)*NTH,J) 1293 WN2(ITH,IK) = WN(IK) 1294 END DO 1295 END DO 1296 ! 1297 ZWND = ZZWND 1298 USTAR = 0. 1299 USTD = 0. 1300 TAUWX = 0. 1301 TAUWY = 0. 1302 ! 1303 CALL W3SPR4 (A, CG, WN, EMEAN, FMEAN, FMEAN1, & 1304 WNMEAN, AMAX, UABS, UDIRR, USTAR, USTD,& 1305 TAUWX, TAUWY, CD, Z0, CHARN, LLWS, FMEANWS ) Page 25 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1306 ! 1307 DO ITT=1, 3 1308 CALL W3SPR4 (A, CG, WN, EMEAN, FMEAN, FMEAN1, & 1309 WNMEAN, AMAX, UABS, UDIRR, USTAR, USTD,& 1310 TAUWX, TAUWY, CD, Z0, CHARN, LLWS, FMEANWS ) 1311 CALL W3SIN4 ( A, CG, WN2, UABS, USTAR, DAIR/DWAT,& 1312 ASO(J), UDIRR, Z0, CD,TAUWX, TAUWY, & 1313 TAUWNX, TAUWNY, ICE, XIN, DIA, LLWS, IX, IY ) 1314 END DO 1315 ! 1316 IF ( FLSRCE(2) ) THEN 1317 CALL W3SLN1 (WN, FHIGH, USTAR, UDIRR, XLN ) 1318 ! 1319 CALL W3SIN4 ( A, CG, WN2, UABS, USTAR, & 1320 DAIR/DWAT, ASO(J), UDIRR, & 1321 Z0, CD, TAUWX, TAUWY,TAUWNX, TAUWNY, & 1322 ICE, XIN, DIA, LLWS, IX, IY ) 1323 END IF 1324 IF ( FLSRCE(3) ) THEN 1325 CALL W3SNL1 ( A, CG, WNMEAN*DEPTH, XNL, DIA ) 1326 ! 1327 END IF 1328 IF ( FLSRCE(4) ) THEN 1329 CALL W3SDS4 ( A, WN, CG, USTAR, USTD, DEPTH, XDS, & 1330 DIA, IX, IY, WHITECAP ) 1331 ! 1332 CALL W3SDB1 ( A, WN2, DEPTH, EMEAN, FMEAN, & 1333 WNMEAN, XDB, DIA ) 1334 END IF 1335 IF ( FLSRCE(5) ) THEN 1336 CALL W3SBT1 ( A, CG, WN, DEPTH, XBT, DIA ) 1337 1338 1339 1340 1341 1342 ! 1343 END IF 1344 ! 1345 IF ( ISCALE.EQ.0 .OR. ISCALE.EQ.3 ) THEN 1346 FACF = TPIINV 1347 FACE = 1. 1348 FACS = 1. 1349 ELSE IF ( ISCALE.EQ.1 .OR. ISCALE.EQ.4 ) THEN 1350 FACF = TPIINV * UABS / GRAV 1351 FACE = GRAV**3 / UABS**5 1352 FACS = GRAV**2 / UABS**4 1353 ELSE IF ( ISCALE.EQ.2 .OR. ISCALE.EQ.5 ) THEN 1354 FACF = TPIINV * USTAR / GRAV 1355 FACE = GRAV**3 / USTAR**5 1356 FACS = GRAV**2 / USTAR**4 1357 END IF 1358 ! 1359 DO IK=1, NK 1360 FACTOR = TPI / CG(IK) * SIG(IK) 1361 E1 (IK) = 0. 1362 SIN1(IK) = 0. Page 26 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1363 SNL1(IK) = 0. 1364 SDS1(IK) = 0. 1365 SBT1(IK) = 0. 1366 STT1(IK) = 0. 1367 DO ITH=1, NTH 1368 ISP = ITH + (IK-1)*NTH 1369 E (IK,ITH) = SPCO(ISP,J) 1370 SWN(IK,ITH) = ( XLN(ITH,IK) + XIN(ITH,IK) ) * FACTOR 1371 SNL(IK,ITH) = ( XNL(ITH,IK) + XTR(ITH,IK) ) * FACTOR 1372 SDS(IK,ITH) = ( XDS(ITH,IK) + XDB(ITH,IK) ) * FACTOR 1373 SBT(IK,ITH) = ( XBT(ITH,IK) * XBS(ITH,IK) ) * FACTOR 1374 STT(IK,ITH) = SWN(IK,ITH) + SNL(IK,ITH) + SDS(IK,ITH) + & 1375 SBT(IK,ITH) + XXX(ITH,IK) * FACTOR 1376 E1 (IK) = E1 (IK) + E(IK,ITH) 1377 SIN1(IK) = SIN1(IK) + SWN(IK,ITH) 1378 SNL1(IK) = SNL1(IK) + SNL(IK,ITH) 1379 SDS1(IK) = SDS1(IK) + SDS(IK,ITH) 1380 SBT1(IK) = SBT1(IK) + SBT(IK,ITH) 1381 END DO 1382 E1 (IK) = E1(IK) * DTH * FACE 1383 SIN1(IK) = SIN1(IK) * DTH * FACS 1384 SNL1(IK) = SNL1(IK) * DTH * FACS 1385 SDS1(IK) = SDS1(IK) * DTH * FACS 1386 SBT1(IK) = SBT1(IK) * DTH * FACS 1387 END DO 1388 ! 1389 STT1 = SIN1 + SNL1 + SDS1 + SBT1 1390 E1ALL(:,1) = SIN1 1391 E1ALL(:,2) = SNL1 1392 E1ALL(:,3) = SDS1 1393 E1ALL(:,4) = SBT1 1394 E1ALL(:,5) = STT1 1395 ! 1396 END IF 1397 ! 1398 ! 4.a Perform output type 1 ( print plots / tables / file ) 1399 ! 1400 IF ( ITYPE .EQ. 1 ) THEN 1401 ! 1402 IF ( OTYPE .EQ. 1 ) THEN 1403 ! 1404 IF ( SCALE1 .GE. 0. ) & 1405 CALL PRT1DS (NDSO, NK, E1, SIG(1:NK), 'RAD/S',& 1406 17, SCALE1, 'E(f)', 'm^2s', PTNME(J) ) 1407 IF ( SCALE2 .GE. 0. ) & 1408 CALL PRT2DS (NDSO, NK, NK, NTH, E, SIG(1:NK), & 1409 'RAD/S', 1., SCALE2, 0.0001, 'E(f,th)', & 1410 'm^2s', PTNME(J) ) 1411 WRITE (NDSO,910) DPO(J), UABS 1412 IF ( WAO(J) .GT. 0. ) WRITE (NDSO,911) UDIR 1413 WRITE (NDSO,912) ASO(J), CAO(J) 1414 IF ( CAO(J) .GT. 0. ) WRITE (NDSO,913) CDIR 1415 WRITE (NDSO,914) HSIG, WLEN, TMEAN, THMEAN, THSPRD 1416 ! 1417 ELSE IF ( OTYPE .EQ. 2 ) THEN 1418 ! 1419 CALL STME21 ( TIME , DTME21 ) Page 27 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1420 IF ( FLAGLL ) THEN 1421 WRITE (NDSTAB,920) DTME21, PTNME(J), & 1422 M2KM*PTLOC(1,J), M2KM*PTLOC(2,J), & 1423 DPO(J), USTAR, WAO(J), UDIR 1424 ELSE 1425 WRITE (NDSTAB,720) DTME21, PTNME(J), & 1426 M2KM*PTLOC(1,J), M2KM*PTLOC(2,J), & 1427 DPO(J), USTAR, WAO(J), UDIR 1428 END IF 1429 IF ( FP .EQ. 0. ) FP = SIG(NK) 1430 DO IK=1, NK 1431 WRITE (NDSTAB,921) TPIINV*SIG(IK), SIG(IK)/FP, & 1432 E1(IK), THBND(IK), SPBND(IK), APM(IK) 1433 END DO 1434 IF ( FP .EQ. SIG(NK) ) FP = 0. 1435 WRITE (NDSTAB,922) 1436 ! 1437 ELSE IF ( OTYPE .EQ. 3 ) THEN 1438 ! 1439 IF ( FLFORM ) THEN 1440 WRITE (NDSTAB) PTNME(J), PTLOC(2,J), & 1441 PTLOC(1,J), DPO(J), WAO(J), & 1442 UDIR, CAO(J), CDIR 1443 WRITE (NDSTAB) ((E(IK,ITH),IK=1,NK),ITH=1,NTH) 1444 ELSE 1445 WRITE (NDSTAB,901) PTNME(J), M2KM*PTLOC(2,J), & 1446 M2KM*PTLOC(1,J), DPO(J), & 1447 WAO(J), UDIR, CAO(J), CDIR 1448 WRITE (NDSTAB,902) & 1449 ((E(IK,ITH),IK=1,NK),ITH=1,NTH) 1450 END IF 1451 ! 1452 END IF 1453 ! 1454 ! 4.b Perform output type 2 ( tables ) 1455 ! 1456 ELSE IF ( ITYPE .EQ. 2 ) THEN 1457 ! 1458 IF ( NREQ .EQ. 1 ) THEN 1459 ! 1460 IYR = TIME(1) / 10000 1461 IMTH = MOD ( TIME(1) , 10000 ) / 100 1462 IDY = MOD ( TIME(1) , 100 ) 1463 IH = TIME(2) / 10000 1464 IM = MOD ( TIME(2) , 10000 ) / 100 1465 IS = MOD ( TIME(2) , 100 ) 1466 IF ( OTYPE .EQ. 1 ) THEN 1467 WRITE (NDSTAB,1921) TIME(1), IH, IM, IS, & 1468 DPO(J), CAO(J), CDIR, WAO(J), UDIR 1469 ELSE IF ( OTYPE .EQ. 2 ) THEN 1470 WRITE (NDSTAB,1922) TIME(1), IH, IM, IS, & 1471 HSIG, WLEN, TMEAN, THMEAN, THSPRD, & 1472 FP*TPIINV, THP, SPP 1473 ELSE IF ( OTYPE.EQ.3 ) THEN 1474 WRITE (NDSTAB,1923) TIME(1), IH, IM, IS, & 1475 UNORM, ESTAR, FPSTAR, CD*1000., APM(NK)*100. 1476 ELSE IF ( OTYPE.EQ.4 ) THEN Page 28 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1477 WRITE (NDSTAB,1924) TIME(1), IH, IM, IS, & 1478 UNORM, ESTAR, FPSTAR, CD*1000., APM(NK)*100. 1479 ELSE IF ( OTYPE.EQ.5 ) THEN 1480 HMAT = MIN ( 100. , 3.33*GRAV*HSIG/UABS**2 ) 1481 IF ( HSIG .GE. HSMIN ) THEN 1482 CALL WAVNU1 ( FP, DPO(J), WNA, XYZ ) 1483 AGE1 = MIN ( 100. , FP / WNA / UABS ) 1484 AFR = TPI / TMEAN 1485 CALL WAVNU1 ( AFR, DPO(J), WNA, XYZ ) 1486 AGE2 = MIN ( 100. , AFR / WNA / UABS ) 1487 ELSE 1488 AGE1 = -9.99 1489 AGE2 = -9.99 1490 END IF 1491 WRITE (NDSTAB,1925) TIME(1), IH, IM, IS, & 1492 WAO(J), UDIR, HSIG, HMAT, AGE1, AGE2, & 1493 ASO(J) 1494 ELSE IF ( OTYPE.EQ.6 ) THEN 1495 IF ( HSIG .GE. HSMIN ) THEN 1496 WRITE (NDSTAB,1926) IYR, IMTH, IDY, IH, & 1497 WAO(J), NINT(UDIR), HSIG, TPI / FP 1498 ELSE 1499 WRITE (NDSTAB,1926) IYR, IMTH, IDY, IH, & 1500 WAO(J), NINT(UDIR), HSIG, 0.0 1501 END IF 1502 END IF 1503 ! 1504 ELSE 1505 ! 1506 IF ( OTYPE .EQ. 1 ) THEN 1507 IF ( FLAGLL ) THEN 1508 WRITE (NDSTAB,1931) M2KM*PTLOC(1,J), & 1509 M2KM*PTLOC(2,J), DPO(J), CAO(J), & 1510 CDIR, WAO(J), UDIR 1511 ELSE 1512 WRITE (NDSTAB,1731) M2KM*PTLOC(1,J), & 1513 M2KM*PTLOC(2,J), DPO(J), CAO(J), & 1514 CDIR, WAO(J), UDIR 1515 END IF 1516 ELSE IF ( OTYPE .EQ. 2 ) THEN 1517 IF ( FLAGLL ) THEN 1518 WRITE (NDSTAB,1932) M2KM*PTLOC(1,J), & 1519 M2KM*PTLOC(2,J), HSIG, WLEN, & 1520 TMEAN, THMEAN, THSPRD, FP*TPIINV, & 1521 THP, SPP 1522 ELSE 1523 WRITE (NDSTAB,1732) M2KM*PTLOC(1,J), & 1524 M2KM*PTLOC(2,J), HSIG, WLEN, & 1525 TMEAN, THMEAN, THSPRD, FP*TPIINV, & 1526 THP, SPP 1527 END IF 1528 ELSE IF ( OTYPE .EQ. 3 ) THEN 1529 WRITE (NDSTAB,1933) 1.E-4*XSTAR, & 1530 1.E-4*YSTAR, UNORM, ESTAR, FPSTAR, & 1531 CD*1000., APM(NK)*100. 1532 ELSE IF ( OTYPE .EQ. 4 ) THEN 1533 WRITE (NDSTAB,1934) XSTAR, YSTAR, UNORM, & Page 29 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1534 ESTAR, FPSTAR, CD*1000., APM(NK)*100. 1535 ELSE IF ( OTYPE .EQ. 5 ) THEN 1536 HMAT = MIN ( 100. , 3.33*GRAV*HSIG/UABS**2 ) 1537 CALL WAVNU1 ( FP, DPO(J), WNA, XYZ ) 1538 AGE1 = MIN ( 100. , FP / WNA / UABS ) 1539 AFR = TPI / TMEAN 1540 CALL WAVNU1 ( AFR, DPO(J), WNA, XYZ ) 1541 AGE2 = MIN ( 100. , AFR / WNA / UABS ) 1542 IF ( FLAGLL ) THEN 1543 WRITE (NDSTAB,1935) M2KM*PTLOC(1,J), & 1544 M2KM*PTLOC(2,J), WAO(J), UDIR, & 1545 HSIG, HMAT, AGE1, AGE2, ASO(J) 1546 ELSE 1547 WRITE (NDSTAB,1735) M2KM*PTLOC(1,J), & 1548 M2KM*PTLOC(2,J), WAO(J), UDIR, & 1549 HSIG, HMAT, AGE1, AGE2, ASO(J) 1550 END IF 1551 ELSE IF ( OTYPE .EQ. 6 ) THEN 1552 IF ( HSIG .GE. HSMIN ) THEN 1553 IF ( FLAGLL ) THEN 1554 WRITE (NDSTAB,1936) M2KM*PTLOC(1,J), & 1555 M2KM*PTLOC(2,J), WAO(J), NINT(UDIR),& 1556 HSIG, TPI / FP 1557 ELSE 1558 WRITE (NDSTAB,1736) M2KM*PTLOC(1,J), & 1559 M2KM*PTLOC(2,J), WAO(J), NINT(UDIR),& 1560 HSIG, TPI / FP 1561 END IF 1562 ELSE 1563 IF ( FLAGLL ) THEN 1564 WRITE (NDSTAB,1936) M2KM*PTLOC(1,J), & 1565 M2KM*PTLOC(2,J), WAO(J), NINT(UDIR),& 1566 HSIG, 0.0 1567 ELSE 1568 WRITE (NDSTAB,1736) M2KM*PTLOC(1,J), & 1569 M2KM*PTLOC(2,J), WAO(J), NINT(UDIR),& 1570 HSIG, 0.0 1571 END IF 1572 END IF 1573 END IF 1574 ! 1575 END IF 1576 ! 1577 ! 4.c Perform output type 3 ( source terms ) 1578 ! 1579 ELSE IF ( ITYPE .EQ. 3 ) THEN 1580 ! 1581 IF ( OTYPE .EQ. 1 ) THEN 1582 ! 1583 IF ( SCALE1 .GE. 0. ) THEN 1584 IF ( FLSRCE(1) ) & 1585 CALL PRT1DS (NDSO, NK, E1, SIG(1:NK), & 1586 'RAD/S', 17, 0., 'E(f)', 'm^2s', & 1587 PTNME(J) ) 1588 IF (FLSRCE(2) .OR. FLSRCE(3) .OR. & 1589 FLSRCE(4) .OR. FLSRCE(5) .OR. & 1590 FLSRCE(6) ) & Page 30 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1591 CALL PRT1DM (NDSO, NK, 5, E1ALL, SIG(1:NK),& 1592 'RAD/S', 17, SCALE1, VAR1, 'M2', & 1593 PTNME(J) ) 1594 END IF 1595 IF ( SCALE2 .GE. 0. ) THEN 1596 IF ( FLSRCE(1) ) & 1597 CALL PRT2DS (NDSO, NK, NK, NTH, E, & 1598 SIG(1:NK), 'RAD/S', 1., 0., 0.0001, & 1599 'E(f,th)', 'm^2s', PTNME(J) ) 1600 IF ( FLSRCE(2) ) & 1601 CALL PRT2DS (NDSO, NK, NK, NTH, SWN, & 1602 SIG(1:NK), 'RAD/S', 1., SCALE2, 0.0001,& 1603 'Sin(f,th)', 'm^2', PTNME(J) ) 1604 IF ( FLSRCE(3) ) & 1605 CALL PRT2DS (NDSO, NK, NK, NTH, SNL, & 1606 SIG(1:NK), 'RAD/S', 1., SCALE2, 0.0001,& 1607 'Snl(f,th)', 'm^2', PTNME(J) ) 1608 IF ( FLSRCE(4) ) & 1609 CALL PRT2DS (NDSO, NK, NK, NTH, SDS, & 1610 SIG(1:NK), 'RAD/S', 1., SCALE2, 0.0001,& 1611 'Sds(f,th)', 'm^2', PTNME(J) ) 1612 IF ( FLSRCE(5) ) & 1613 CALL PRT2DS (NDSO, NK, NK, NTH, SBT, & 1614 SIG(1:NK), 'RAD/S', 1., SCALE2, 0.0001,& 1615 'Sbt(f,th)', 'm^2', PTNME(J) ) 1616 IF ( FLSRCE(6) ) & 1617 CALL PRT2DS (NDSO, NK, NK, NTH, STT, & 1618 SIG(1:NK), 'RAD/S', 1., SCALE2, 0.0001,& 1619 'Stot(f,th)', 'm^2', PTNME(J) ) 1620 END IF 1621 ! 1622 ELSE IF ( OTYPE .EQ. 2 ) THEN 1623 ! 1624 CALL STME21 ( TIME , DTME21 ) 1625 IF ( FLAGLL ) THEN 1626 WRITE (NDSTAB,2920) DTME21, PTNME(J), & 1627 M2KM*PTLOC(1,J), M2KM*PTLOC(2,J), & 1628 DPO(J), USTAR, WAO(J) 1629 ELSE 1630 WRITE (NDSTAB,2720) DTME21, PTNME(J), & 1631 M2KM*PTLOC(1,J), M2KM*PTLOC(2,J), & 1632 DPO(J), USTAR, WAO(J) 1633 END IF 1634 IF ( ISCALE.EQ.0 ) THEN 1635 WRITE (NDSTAB,2921) 1636 ELSE IF ( ISCALE.EQ.1 .OR. ISCALE.EQ.2 ) THEN 1637 WRITE (NDSTAB,2922) 1638 ELSE IF ( ISCALE.EQ.3 ) THEN 1639 WRITE (NDSTAB,2923) 1640 ELSE IF ( ISCALE.EQ.4 .OR. ISCALE.EQ.5 ) THEN 1641 WRITE (NDSTAB,2924) 1642 END IF 1643 IF ( ISCALE.GE.3 ) FACF = 1. / FP 1644 DO IK=1, NK 1645 WRITE (NDSTAB,2930) FACF*SIG(IK), E1(IK), & 1646 SIN1(IK), SNL1(IK), SDS1(IK), SBT1(IK), STT1(IK) 1647 END DO Page 31 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1648 WRITE (NDSTAB,2940) 1649 ! 1650 ELSE IF ( OTYPE .EQ. 3 ) THEN 1651 ! 1652 CALL STME21 ( TIME , DTME21 ) 1653 IF ( FLAGLL ) THEN 1654 WRITE (NDSTAB,2920) DTME21, PTNME(J), & 1655 M2KM*PTLOC(1,J), M2KM*PTLOC(2,J), & 1656 DPO(J), USTAR, WAO(J) 1657 ELSE 1658 WRITE (NDSTAB,2720) DTME21, PTNME(J), & 1659 M2KM*PTLOC(1,J), M2KM*PTLOC(2,J), & 1660 DPO(J), USTAR, WAO(J) 1661 END IF 1662 IF ( ISCALE.EQ.0 ) THEN 1663 WRITE (NDSTAB,2925) 1664 ELSE IF ( ISCALE.EQ.1 .OR. ISCALE.EQ.2 ) THEN 1665 WRITE (NDSTAB,2926) 1666 ELSE IF ( ISCALE.EQ.3 ) THEN 1667 WRITE (NDSTAB,2927) 1668 ELSE IF ( ISCALE.EQ.4 .OR. ISCALE.EQ.5 ) THEN 1669 WRITE (NDSTAB,2928) 1670 END IF 1671 ! 1672 IF ( ISCALE.GE.3 ) FACF = 1. / FP 1673 DO IK=1, NK 1674 FACT = 1. / MAX ( 1.E-10 , E1(IK) ) 1675 IF ( E1(IK) .GT. 1.E-10 ) THEN 1676 WRITE (NDSTAB,2931) FACF*SIG(IK), E1(IK), & 1677 FACT*SIN1(IK), FACT*SNL1(IK), & 1678 FACT*SDS1(IK), FACT*SBT1(IK), FACT*STT1(IK) 1679 ELSE 1680 WRITE (NDSTAB,2931) FACF*SIG(IK), E1(IK) 1681 END IF 1682 END DO 1683 WRITE (NDSTAB,2940) 1684 ! 1685 ELSE IF ( OTYPE .EQ. 4 ) THEN 1686 ! 1687 IF ( FLFORM ) THEN 1688 WRITE (NDSTAB) PTNME(J), PTLOC(2,J), & 1689 PTLOC(1,J), DPO(J), WAO(J), & 1690 UDIR, CAO(J), CDIR 1691 IF ( FLSRCE(1) ) WRITE (NDSTAB) & 1692 ((E(IK,ITH),IK=1,NK),ITH=1,NTH) 1693 IF ( FLSRCE(2) ) WRITE (NDSTAB) & 1694 ((SWN(IK,ITH),IK=1,NK),ITH=1,NTH) 1695 IF ( FLSRCE(3) ) WRITE (NDSTAB) & 1696 ((SNL(IK,ITH),IK=1,NK),ITH=1,NTH) 1697 IF ( FLSRCE(4) ) WRITE (NDSTAB) & 1698 ((SDS(IK,ITH),IK=1,NK),ITH=1,NTH) 1699 IF ( FLSRCE(5) ) WRITE (NDSTAB) & 1700 ((SBT(IK,ITH),IK=1,NK),ITH=1,NTH) 1701 IF ( FLSRCE(6) ) WRITE (NDSTAB) & 1702 ((STT(IK,ITH),IK=1,NK),ITH=1,NTH) 1703 ELSE 1704 IF ( FLAGLL ) THEN Page 32 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1705 WRITE (NDSTAB,901) PTNME(J), & 1706 M2KM*PTLOC(2,J), M2KM*PTLOC(1,J), & 1707 DPO(J), WAO(J), UDIR, CAO(J), CDIR 1708 ELSE 1709 WRITE (NDSTAB,701) PTNME(J), & 1710 M2KM*PTLOC(2,J), M2KM*PTLOC(1,J), & 1711 DPO(J), WAO(J), UDIR, CAO(J), CDIR 1712 END IF 1713 IF ( FLSRCE(1) ) WRITE (NDSTAB,902) & 1714 ((E(IK,ITH),IK=1,NK),ITH=1,NTH) 1715 IF ( FLSRCE(2) ) WRITE (NDSTAB,902) & 1716 ((SWN(IK,ITH),IK=1,NK),ITH=1,NTH) 1717 IF ( FLSRCE(3) ) WRITE (NDSTAB,902) & 1718 ((SNL(IK,ITH),IK=1,NK),ITH=1,NTH) 1719 IF ( FLSRCE(4) ) WRITE (NDSTAB,902) & 1720 ((SDS(IK,ITH),IK=1,NK),ITH=1,NTH) 1721 IF ( FLSRCE(5) ) WRITE (NDSTAB,902) & 1722 ((SBT(IK,ITH),IK=1,NK),ITH=1,NTH) 1723 IF ( FLSRCE(6) ) WRITE (NDSTAB,902) & 1724 ((STT(IK,ITH),IK=1,NK),ITH=1,NTH) 1725 END IF 1726 ! 1727 END IF 1728 ! 1729 ! 4.d Perform output type 4 ( Spectral partitions and bulletins ) 1730 ! 1731 ELSE IF ( ITYPE .EQ. 4 ) THEN 1732 ! 1733 IF ( OTYPE .EQ. 1 ) THEN 1734 ! 1735 IF ( FLAGLL ) THEN 1736 IF ( PTLOC(1,J) .LT. 0. ) & 1737 PTLOC(1,J) = PTLOC(1,J) + 360. 1738 WRITE (NDSTAB,940) TIME, M2KM*PTLOC(2,J), & 1739 M2KM*PTLOC(1,J), PTNME(J), NPART, DEPTH, & 1740 WAO(J), UDIR, CAO(J), CDIR 1741 ELSE 1742 WRITE (NDSTAB,943) TIME, M2KM*PTLOC(1,J), & 1743 M2KM*PTLOC(2,J), PTNME(J), NPART, DEPTH, & 1744 WAO(J), UDIR, CAO(J), CDIR 1745 END IF 1746 ! WRITE (NDSTAB,941) 1747 DO I=0, NPART 1748 WRITE (NDSTAB,942) I, XPART(:,I) 1749 END DO 1750 ! 1751 ELSEIF ( OTYPE .GE. 2 ) THEN 1752 CALL W3BULL (NPART, XPART, DIMXP, UABS, UDIR, & 1753 J, IOUT, TIMEV ) 1754 ! 1755 IF ( FLAGLL ) THEN 1756 X = M2KM * PTLOC(1,J) 1757 Y = M2KM * PTLOC(2,J) 1758 1759 X = MOD ( X+720. , 360. ) 1760 IF ( X .LE. 180. ) THEN 1761 IDLON = 'E' Page 33 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1762 ELSE 1763 X = 360. - X 1764 IDLON = 'W' 1765 ENDIF 1766 !IF ( ABS(Y) .LE. 0.0049 ) THEN 1767 !IDLAT = '-' 1768 IF ( Y .GE. 0. ) THEN 1769 IDLAT = 'N' 1770 ELSE 1771 IDLAT = 'S' 1772 Y = -Y 1773 ENDIF 1774 ELSE 1775 IDLAT = ' ' 1776 IDLON = ' ' 1777 ENDIF 1778 IF ( OTYPE .EQ. 2 .OR. OTYPE .EQ. 4 ) THEN 1779 NDSBUL=NDSTAB + (J - 1) 1780 NDSCBUL=NDSTAB + (J - 1) + NOPTS 1781 IF (IOUT .EQ. 1) THEN 1782 WRITE(HSTR,'(I2,1X,A)') TIMEV(2)/10000, & 1783 HTYPE 1784 WRITE (NDSBUL,970) PTNME(J), Y, IDLAT, X, & 1785 IDLON, GNAME, TIMEV(1), & 1786 HSTR 1787 WRITE (NDSBUL,971) 1788 WRITE (NDSBUL,972) 1789 WRITE (NDSBUL,971) 1790 WRITE (NDSCBUL,960) PTNME(J), Y, IDLAT, & 1791 X, IDLON, GNAME, TIMEV(1), HSTR 1792 WRITE (NDSCBUL,961) 1793 ENDIF 1794 1795 WRITE (NDSBUL,973) ASCBLINE 1796 WRITE (NDSCBUL,963) CASCBLINE 1797 ENDIF 1798 IF ( OTYPE .EQ. 3 .OR. OTYPE .EQ. 4 ) THEN 1799 ICSV = 0 1800 IF ( NDSBUL .GT. 0 ) ICSV = NDSBUL 1801 IF ( NDSCBUL .GT. 0 ) ICSV = NDSCBUL 1802 NDSCSV = NDSTAB + (J - 1) + ICSV 1803 WRITE (NDSCSV,'(A664)') CSVBLINE 1804 ENDIF 1805 END IF 1806 ! 1807 END IF 1808 ! ... End of fields loop 1809 ! 1810 END IF 1811 END DO 1812 ! 1813 RETURN 1814 ! 1815 ! Formats 1816 ! 1817 900 FORMAT (I8.8,I7.6) 1818 901 FORMAT ('''',A10,'''',2F7.2,F10.1,2(F7.2,F6.1)) Page 34 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1819 701 FORMAT ('''',A10,'''',2(F8.1,'E3'),F10.1,2(F7.2,F6.1)) 1820 902 FORMAT (7E11.3) 1821 905 FORMAT (9X,A) 1822 910 FORMAT (/15X,' Water depth :',F7.1,' (m)'/ & 1823 15X,' Wind speed :',F8.2,' (m/s)') 1824 911 FORMAT ( 15X,' Wind direction :',F7.1,' (degr)') 1825 912 FORMAT ( 15X,' Air-sea temp. dif.:',F7.1,' (degr)'/ & 1826 15X,' Current speed :',F8.2,' (m/s)') 1827 913 FORMAT ( 15X,' Current direction :',F7.1,' (degr)') 1828 914 FORMAT ( 15X,' Wave height :',F8.2,' (m)'/ & 1829 15X,' Mean wave length :',F6.0,' (m)'/ & 1830 15X,' Mean wave period :',F7.1,' (s)'/ & 1831 15X,' Mean wave direct. :',F7.1,' (degr)'/ & 1832 15X,' Direct. spread :',F7.1,' (degr)'/) 1833 920 FORMAT (' Time : ',A/ & 1834 ' Location : ',A,' (',2F8.2,' )'/ & 1835 ' depth : ',F7.1,' m'/ & 1836 ' U* : ',F9.3,' m/s'/ & 1837 ' U10 : ',F7.1,' m/s'/ & 1838 ' Dir U10 : ',F7.1,' degr'// & 1839 ' f f/fp F(f) theta spr alpha '/ & 1840 ' (Hz) (-) (m2s) (deg) (deg) (-) '/ & 1841 ' --------------------------------------------------') 1842 720 FORMAT (' Time : ',A/ & 1843 ' Location : ',A,' (',2(F8.1,'E3'),' )'/ & 1844 ' depth : ',F7.1,' m'/ & 1845 ' U* : ',F9.3,' m/s'/ & 1846 ' U10 : ',F7.1,' m/s'/ & 1847 ' Dir U10 : ',F7.1,' degr'// & 1848 ' f f/fp F(f) theta spr alpha '/ & 1849 ' (Hz) (-) (m2s) (deg) (deg) (-) '/ & 1850 ' --------------------------------------------------') 1851 921 FORMAT (1x,F8.5,F7.3,E11.3,2F8.1,F8.4) 1852 922 FORMAT (' '/' ') 1853 ! 1854 940 FORMAT (1X,I8.8,1X,I6.6,2F8.3,2X,'''',A10,'''', & 1855 1X,I2,F7.1,F5.1,f6.1,F5.2,F6.1) 1856 943 FORMAT (1X,I8.8,1X,I6.6,2(F8.1,'E3'),2X,'''',A10,'''', & 1857 1X,I2,F7.1,F5.1,f6.1,F5.2,F6.1) 1858 941 FORMAT (' hs tp lp theta sp wf') 1859 942 FORMAT (I3,3F8.2,2F9.2,F7.2) 1860 ! 1861 960 FORMAT ( 'Location : ',A,' (',F5.2,A,1X,F6.2,A,')'/ & 1862 'Model : ',A/ & 1863 'Cycle : ',I8,1X,A// & 1864 'DDHH HS SS PP DDD SS PP DDD SS PP DDD', & 1865 ' SS PP DDD SS PP DDD SS PP DDD') 1866 961 FORMAT ('----------------------------------------', & 1867 '---------------------------') 1868 963 FORMAT (A) 1869 ! 1870 970 FORMAT ( ' Location : ',A,' (',F5.2,A,1X,F6.2,A,')'/ & 1871 ' Model : ',A/ & 1872 ' Cycle : ',I8,1X,A) 1873 971 FORMAT (' +-------+-----------+-----------------+', & 1874 '-----------------+-----------------+----', & 1875 '-------------+-----------------+--------', & Page 35 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1876 '---------+') 1877 972 FORMAT (' | day & | Hst n x | Hs Tp dir |', & 1878 ' Hs Tp dir |', & 1879 ' Hs Tp dir |', & 1880 ' Hs Tp dir |', & 1881 ' Hs Tp dir |', & 1882 ' Hs Tp dir |'/ & 1883 ' | hour | (m) - - | (m) (s) (d) |', & 1884 ' (m) (s) (d) |', & 1885 ' (m) (s) (d) |', & 1886 ' (m) (s) (d) |', & 1887 ' (m) (s) (d) |', & 1888 ' (m) (s) (d) |') 1889 973 FORMAT (1X,A) 1890 ! 1891 1901 FORMAT ( & 1892 ' Date Time d Uc Dir. U10 Dir. '/ & 1893 ' h m s (m) (m/s) (d.N) (m/s) (d.N) '/ & 1894 ' ---------------------------------------------------------') 1895 1902 FORMAT ( & 1896 ' Date Time Hs L Tr Dir. Spr. ', & 1897 ' fp p_dir p_spr'/ & 1898 ' h m s (m) (m) (s) (d.N) (deg)', & 1899 ' (Hz) (d.N) (deg)'/ & 1900 ' -------------------------------------------------------', & 1901 '-----------------------') 1902 1903 FORMAT ( & 1903 ' Date Time U* E* fp* Cd alpha'/& 1904 ' h m s (m/s) (-) (-) *1000 *100'/ & 1905 ' --------------------------------------------------------------') 1906 1904 FORMAT ( & 1907 ' Date Time U10 E* fp* Cd alpha'/& 1908 ' h m s (m/s) (-) (-) *1000 *100'/ & 1909 ' --------------------------------------------------------------') 1910 1905 FORMAT ( & 1911 ' Date Time U10 Dir. Hs H* cp/U ', & 1912 ' cm/U Dt'/ & 1913 ' (m/s) (d.N) (m) (-) (-) ', & 1914 ' (-) (deg)'/ & 1915 ' --------------------------------------------------', & 1916 '---------------------') 1917 1906 FORMAT ( & 1918 ' Time U10 Dir. Hs Tp '/ & 1919 ' yr mth dy h (m/s) (d.N) (m) (s) '/ & 1920 ' ----------------------------------') 1921 1910 FORMAT ( ' '/' ' ) 1922 1911 FORMAT (' Time : ',A// & 1923 ' Long. Lat. d Uc Dir. U10 Dir. '/ & 1924 ' (m) (m/s) (d.N) (m/s) (d.N) '/ & 1925 ' --------------------------------------------------------') 1926 1912 FORMAT (' Time : ',A// & 1927 ' Long. Lat. Hs L Tr Dir. Spr. ', & 1928 ' fp p_dir p_spr'/ & 1929 ' (m) (m) (s) (d.N) (deg)', & 1930 ' (Hz) (d.N) (deg)'/ & 1931 ' ------------------------------------------------------', & 1932 '-----------------------') Page 36 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1933 1711 FORMAT (' Time : ',A// & 1934 ' X Y d Uc Dir. U10 Dir. '/ & 1935 ' (m) (m) (m) (m/s) (d.N) (m/s) (d.N) '/ & 1936 ' ----------------------------------------------------------') 1937 1712 FORMAT (' Time : ',A// & 1938 ' X Y Hs L Tr Dir. Spr. ', & 1939 ' fp p_dir p_spr'/ & 1940 ' (m) (m)) (m) (m) (s) (d.N) (deg)', & 1941 ' (Hz) (d.N) (deg)'/ & 1942 ' ------------------------------------------------------', & 1943 '-------------------------') 1944 1913 FORMAT (' Time : ',A// & 1945 ' X* Y* U* E* fp* Cd alpha'/& 1946 ' (-) (-) (m/s) (-) (-) *1000 *100'/ & 1947 ' --------------------------------------------------------------') 1948 1914 FORMAT (' Time : ',A// & 1949 ' X* Y* U10 E* fp* Cd alpha'/ & 1950 ' (-) (-) (m/s) (-) (-) *1000 *100 '/ & 1951 ' --------------------------------------------------------------') 1952 1915 FORMAT (' Time : ',A// & 1953 ' Long. Lat. U10 Dir. Hs H* cp/U ', & 1954 ' cm/U Dt'/ & 1955 ' (m/s) (d.N) (m) (-) (-) ', & 1956 ' (-) (deg)'/ & 1957 ' -------------------------------------------------', & 1958 '---------------------') 1959 1715 FORMAT (' Time : ',A// & 1960 ' X Y U10 Dir. Hs H* cp/U ', & 1961 ' cm/U Dt'/ & 1962 ' (m) (m) (m/s) (d.N) (m) (-) (-) ', & 1963 ' (-) (deg)'/ & 1964 ' ---------------------------------------------------', & 1965 '---------------------') 1966 1916 FORMAT (' Time : ',A// & 1967 ' Long. Lat. U10 Dir. Hs Tp '/ & 1968 ' (m/s) (d.N) (m) (s) '/ & 1969 '-----------------------------------------------') 1970 1716 FORMAT (' Time : ',A// & 1971 ' X Y U10 Dir. Hs Tp '/ & 1972 ' (m) (m) (m/s) (d.N) (m) (s) '/ & 1973 '---------------------------------------------') 1974 1921 FORMAT ( 2X,I8.8,I3,2(1X,I2.2),F10.1,F6.2,F7.1,F6.2,F7.1) 1975 1922 FORMAT ( 2X,I8.8,I3,2(1X,I2.2),F9.3,F7.1,F7.2,F7.1,F7.2, & 1976 F8.4,F7.1,F7.2) 1977 1923 FORMAT ( 2X,I8.8,I3,2(1X,I2.2),F8.4,2E11.3,2F7.3) 1978 1924 FORMAT ( 2X,I8.8,I3,2(1X,I2.2),F7.1,2E11.3,2F7.3) 1979 1925 FORMAT ( 2X,I8.8,I3,2(1X,I2.2),F7.2,F7.1,2F7.2,2F8.2,F7.1) 1980 1926 FORMAT ( 2X,I4,3(1X,I2),F5.1,1X,I3,2F5.1) 1981 1931 FORMAT ( 2X,2F8.3,F10.1,F6.2,F7.1,F6.2,F7.1) 1982 1932 FORMAT ( 2X,2F8.3,F9.3,F7.1,F7.2,F7.1,F7.2, & 1983 F8.4,F7.1,F7.2) 1984 1731 FORMAT ( 2X,2(F7.1,'E3'),F10.1,F6.2,F7.1,F6.2,F7.1) 1985 1732 FORMAT ( 2X,2(F7.1,'E3'),F9.3,F7.1,F7.2,F7.1,F7.2, & 1986 F8.4,F7.1,F7.2) 1987 1933 FORMAT ( 2X,2(F7.1,'E4'),F8.4,2E11.3,2F7.3) 1988 1934 FORMAT ( 2X,2F9.1,F7.1,2E11.3,2F7.3) 1989 1935 FORMAT ( 2X,2F8.3,F7.2,F7.1,2F7.2,2F8.2,F7.1) Page 37 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 1990 1735 FORMAT ( 2X,2(F7.1,'E3'),F7.2,F7.1,2F7.2,2F8.2,F7.1) 1991 1936 FORMAT ( 2X,2F8.3,F5.1,1X,I3,2F5.1) 1992 1736 FORMAT ( 2X,2(F7.1,'E3'),F5.1,1X,I3,2F5.1) 1993 ! 1994 2920 FORMAT (' Time : ',A/ & 1995 ' Location : ',A,' (',2F8.2,' )'/ & 1996 ' depth : ',F7.1,' m'/ & 1997 ' U* : ',F9.3,' m/s'/ & 1998 ' U10 : ',F7.1,' m/s'/) 1999 2720 FORMAT (' Time : ',A/ & 2000 ' Location : ',A,' (',2(F8.1,'E3'),' )'/ & 2001 ' depth : ',F7.1,' m'/ & 2002 ' U* : ',F9.3,' m/s'/ & 2003 ' U10 : ',F7.1,' m/s'/) 2004 2921 FORMAT (' f E ', & 2005 ' Sin Snl Sds Sbt Stot'/ & 2006 ' (Hz) (m2s) ', & 2007 ' (m2) (m2) (m2) (m2) (m2)'/ & 2008 ' -------------------------------------', & 2009 '--------------------------------------') 2010 2922 FORMAT (' f* E* ', & 2011 ' Sin* Snl* Sds* Sbt* Stot*'/ & 2012 ' (-) (-) ', & 2013 ' (-) (-) (-) (-) (-)'/ & 2014 ' -------------------------------------', & 2015 '--------------------------------------') 2016 2923 FORMAT (' f/fp E ', & 2017 ' Sin Snl Sds Sbt Stot'/ & 2018 ' (-) (m2s) ', & 2019 ' (m2) (m2) (m2) (m2) (m2)'/ & 2020 ' -------------------------------------', & 2021 '--------------------------------------') 2022 2924 FORMAT (' f/fp E* ', & 2023 ' Sin* Snl* Sds* Sbt* Stot*'/ & 2024 ' (-) (-) ', & 2025 ' (-) (-) (-) (-) (-)'/ & 2026 ' -------------------------------------', & 2027 '--------------------------------------') 2028 2925 FORMAT (' f E ', & 2029 ' Tini Tnli Tdsi Tbti Ttoti'/ & 2030 ' (Hz) (m2s) ', & 2031 ' (1/s) (1/s) (1/s) (1/s) (1/s)'/ & 2032 ' -----------------------------------', & 2033 '--------------------------------------') 2034 2926 FORMAT (' f* E* ', & 2035 ' Tini* Tnli* Tdsi* Tbti* Ttoti*'/ & 2036 ' (-) (-) ', & 2037 ' (-) (-) (-) (-) (-)'/ & 2038 ' -----------------------------------', & 2039 '--------------------------------------') 2040 2927 FORMAT (' f/fp E ', & 2041 ' Tini Tnli Tdsi Tbti Ttoti'/ & 2042 ' (-) (m2s) ', & 2043 ' (1/s) (1/s) (1/s) (1/s) (1/s)'/ & 2044 ' -----------------------------------', & 2045 '--------------------------------------') 2046 2928 FORMAT (' f/fp E* ', & Page 38 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 2047 ' Tini* Tnli* Tdsi* Tbti* Ttoti*'/ & 2048 ' (-) (-) ', & 2049 ' (-) (-) (-) (-) (-)'/ & 2050 ' -----------------------------------', & 2051 '--------------------------------------') 2052 2930 FORMAT (1X,F6.4,2X,6E11.3) ..................1 (1) Recommended relationship between field width 'W' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+3'. 2053 2931 FORMAT (1X,F6.4,6E11.3) ..................1 (1) Recommended relationship between field width 'W' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+3'. 2054 2940 FORMAT ( ' '/' ' ) 2055 ! 2056 !/ 2057 !/ End of W3EXPO ----------------------------------------------------- / 2058 !/ 2059 END SUBROUTINE W3EXPO ENTRY POINTS Name wavespec_IP_w3expo_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 1711 Label 1933 1055 1712 Label 1937 1061 1715 Label 1959 1071 1716 Label 1970 1077 1731 Label 1984 1512 1732 Label 1985 1523 1735 Label 1990 1547 1736 Label 1992 1558,1568 1901 Label 1891 1035 1902 Label 1895 1037 1903 Label 1902 1039 1904 Label 1906 1041 1905 Label 1910 1043 1906 Label 1917 1045 1910 Label 1921 1050 1911 Label 1922 1053 1912 Label 1926 1059 1913 Label 1944 1064 1914 Label 1948 1066 1915 Label 1952 1069 1916 Label 1966 1075 1921 Label 1974 1467 1922 Label 1975 1470 1923 Label 1977 1474 1924 Label 1978 1477 Page 39 Source Listing W3EXPO 2014-09-16 17:02 Symbol Table multiwavespec.f90 Name Object Declared Type Bytes Dimen Elements Attributes References 1925 Label 1979 1491 1926 Label 1980 1496,1499 1931 Label 1981 1508 1932 Label 1982 1518 1933 Label 1987 1529 1934 Label 1988 1533 1935 Label 1989 1543 1936 Label 1991 1554,1564 2720 Label 1999 1630,1658 2920 Label 1994 1626,1654 2921 Label 2004 1635 2922 Label 2010 1637 2923 Label 2016 1639 2924 Label 2022 1641 2925 Label 2028 1663 2926 Label 2034 1665 2927 Label 2040 1667 2928 Label 2046 1669 2930 Label 2052 1645 2931 Label 2053 1676,1680 2940 Label 2054 1648,1683 701 Label 1819 1709 720 Label 1842 1425 900 Label 1817 1028 901 Label 1818 1445,1705 902 Label 1820 1448,1713,1715,1717,1719,1721,1723 905 Label 1821 1086 910 Label 1822 1411 911 Label 1824 1412 912 Label 1825 1413 913 Label 1827 1414 914 Label 1828 1415 920 Label 1833 1421 921 Label 1851 1431 922 Label 1852 1435 940 Label 1854 1738 941 Label 1858 942 Label 1859 1748 943 Label 1856 1742 960 Label 1861 1790 961 Label 1866 1792 963 Label 1868 1796 970 Label 1870 1784 971 Label 1873 1787,1789 972 Label 1877 1788 973 Label 1889 1795 A Local 975 R(4) 4 2 0 1233,1244,1251,1254,1292,1303,1308 ,1311,1319,1325,1329,1332,1336 ABS Func 1209 scalar 1209 AFR Local 966 R(4) 4 scalar 1484,1485,1486,1539,1540,1541 AGE1 Local 966 R(4) 4 scalar 1483,1488,1492,1538,1545,1549 AGE2 Local 966 R(4) 4 scalar 1486,1489,1492,1541,1545,1549 AMAX Local 969 R(4) 4 scalar 1245,1255,1304,1309 APM Local 974 R(4) 4 1 0 1142,1432,1475,1478,1531,1534 ASCBLINE Local 1795 CHAR 129 scalar 174,1795 Page 40 Source Listing W3EXPO 2014-09-16 17:02 Symbol Table multiwavespec.f90 Name Object Declared Type Bytes Dimen Elements Attributes References ASO Local 1252 R(4) 4 1 1 PTR 172,1252,1312,1320,1413,1493,1545, 1549 ATAN2 Func 1144 scalar 1144,1172 CAO Local 1413 R(4) 4 1 1 PTR 172,1413,1414,1442,1447,1468,1509, 1513,1690,1707,1711,1740,1744 CASCBLINE Local 1796 CHAR 67 scalar 175,1796 CD Local 964 R(4) 4 scalar 1246,1252,1256,1305,1310,1312,1321 ,1475,1478,1531,1534 CDIR Local 959 R(4) 4 scalar 1104,1414,1442,1447,1468,1510,1514 ,1690,1707,1711,1740,1744 CDO Local 1104 R(4) 4 1 1 PTR 172,1104 CG Local 973 R(4) 4 1 0 1114,1117,1230,1244,1251,1254,1290 ,1303,1308,1311,1319,1325,1329,133 6,1360 CHARN Local 968 R(4) 4 scalar 1246,1256,1305,1310 COS Func 1278 scalar 1278 CSVBLINE Local 1803 CHAR 664 scalar 174,1803 DAIR Param 1251 R(4) 4 scalar 1251,1311,1320 DENOM Local 963 R(4) 4 scalar 1205,1209 DEPTH Local 959 R(4) 4 scalar 1098,1099,1113,1222,1325,1329,1332 ,1336,1739,1743 DERA Param 1277 R(4) 4 scalar 1277,1278,1279 DFAC Param 949 R(4) 4 scalar 949 DIA Local 977 R(4) 4 2 0 1253,1313,1322,1325,1330,1333,1336 DIMXP Local 1222 I(4) 4 scalar 269,270,1222,1752 DMIN Local 1098 R(4) 4 scalar PTR 1098 DPO Local 1098 R(4) 4 1 1 PTR 172,1098,1411,1423,1427,1441,1446, 1468,1482,1485,1509,1513,1537,1540 ,1628,1632,1656,1660,1689,1707,171 1 DSIE Local 949 R(4) 4 scalar 949,1108,1111 DSII Local 1151 R(4) 4 1 1 PTR 1151,1155,1156 DTH Local 1141 R(4) 4 scalar PTR 1141,1161,1165,1166,1175,1212,1382 ,1383,1384,1385,1386 DTME21 Local 986 CHAR 23 scalar 1049,1053,1055,1059,1061,1064,1066 ,1069,1071,1075,1077,1085,1086,141 9,1421,1425,1624,1626,1630,1652,16 54,1658 DWAT Param 1251 R(4) 4 scalar 1251,1311,1320 E Local 974 R(4) 4 2 0 1136,1184,1222,1369,1376,1408,1443 ,1449,1597,1692,1714 E1 Local 974 R(4) 4 1 0 1141,1142,1143,1151,1161,1182,1191 ,1195,1196,1203,1204,1361,1376,138 2,1405,1432,1585,1645,1674,1675,16 76,1680 E1ALL Local 984 R(4) 4 2 0 1390,1391,1392,1393,1394,1591 EBND Local 961 R(4) 4 scalar 1131,1137,1141,1146,1151,1152,1153 ,1154,1161,1162,1163,1164 EBX Local 961 R(4) 4 scalar 1132,1138,1144,1146,1155,1165 EBY Local 962 R(4) 4 scalar 1133,1139,1144,1146,1156,1166 ECG1 Local 949 R(4) 4 1 122 949,1114 ECOS Local 1138 R(4) 4 1 1 PTR 1138 EH Local 963 R(4) 4 scalar 1204,1205,1208 EL Local 963 R(4) 4 scalar 1203,1205,1208 EMAX Local 963 R(4) 4 scalar 1191,1195,1196 Page 41 Source Listing W3EXPO 2014-09-16 17:02 Symbol Table multiwavespec.f90 Name Object Declared Type Bytes Dimen Elements Attributes References EMEAN Local 967 R(4) 4 scalar 1244,1254,1303,1308,1332 ESIN Local 1139 R(4) 4 1 1 PTR 1139 ESTAR Local 964 R(4) 4 scalar 1264,1272,1475,1478,1530,1534 ET Local 961 R(4) 4 scalar 1125,1152,1162,1168,1170,1171,1174 ,1272 ETR Local 961 R(4) 4 scalar 1127,1154,1164,1171 ETX Local 961 R(4) 4 scalar 1128,1155,1165,1172,1174 ETY Local 961 R(4) 4 scalar 1129,1156,1166,1172,1174 EWN Local 961 R(4) 4 scalar 1126,1153,1163,1170 EWN1 Local 949 R(4) 4 1 122 949,1113 FACE Local 965 R(4) 4 scalar 1347,1351,1355,1382 FACF Local 965 R(4) 4 scalar 1346,1350,1354,1643,1645,1672,1676 ,1680 FACS Local 965 R(4) 4 scalar 1348,1352,1356,1383,1384,1385,1386 FACT Local 966 R(4) 4 scalar 1674,1677,1678 FACTOR Local 964 R(4) 4 scalar 1230,1233,1290,1292,1360,1370,1371 ,1372,1373,1375 FHIGH Local 967 R(4) 4 scalar 1317 FLAGLL Local 996 L(4) 4 scalar 263,279,617,996,1052,1058,1068,107 4,1276,1420,1507,1517,1542,1553,15 63,1625,1653,1704,1735,1755 FLFORM Local 1025 L(4) 4 scalar 363,394,460,492,1025,1439,1687 FLREQ Local 1094 L(4) 4 1 1 ALC 313,314,322,323,570,575,589,616,66 5,1094 FLSRCE Local 1316 L(4) 4 1 6 210,460,498,508,519,1316,1324,1328 ,1335,1584,1588,1589,1590,1596,160 0,1604,1608,1612,1616,1691,1693,16 95,1697,1699,1701,1713,1715,1717,1 719,1721,1723 FMEAN Local 968 R(4) 4 scalar 1244,1254,1303,1308,1332 FMEAN1 Local 970 R(4) 4 scalar 1244,1254,1303,1308 FMEANS Local 969 R(4) 4 scalar FMEANWS Local 969 R(4) 4 scalar 1246,1256,1305,1310 FP Local 963 R(4) 4 scalar 1208,1214,1273,1429,1431,1434,1472 ,1482,1483,1497,1520,1525,1537,153 8,1556,1560,1643,1672 FPSTAR Local 965 R(4) 4 scalar 1265,1273,1475,1478,1530,1534 FTE Local 1162 R(4) 4 scalar PTR 1162,1165,1166 FTTR Local 1164 R(4) 4 scalar PTR 1164 FTWL Local 1163 R(4) 4 scalar PTR 1163 GNAME Local 1785 CHAR 30 scalar PTR 261,400,413,1785,1791 GRAV Param 1116 R(4) 4 scalar 1116,1117,1142,1272,1273,1274,1275 ,1350,1351,1352,1354,1355,1356,148 0,1536 HMAT Local 965 R(4) 4 scalar 1480,1492,1536,1545,1549 HSIG Local 962 R(4) 4 scalar 1168,1169,1207,1415,1471,1480,1481 ,1492,1495,1497,1500,1519,1524,153 6,1545,1549,1552,1556,1560,1566,15 70 HSMIN Local 972 R(4) 4 scalar 972,1169,1207,1481,1495,1552 HSTR Local 1782 CHAR 6 scalar 1782,1786,1791 HTYPE Local 1783 CHAR 3 scalar 554,1783 I Local 957 I(4) 4 scalar 1747,1748 I1 Local 955 I(4) 4 scalar 1108,1109,1110,1111,1113,1114 I2 Local 955 I(4) 4 scalar 1110,1113,1114 Page 42 Source Listing W3EXPO 2014-09-16 17:02 Symbol Table multiwavespec.f90 Name Object Declared Type Bytes Dimen Elements Attributes References ICE Local 971 R(4) 4 scalar 1253,1313,1322 ICSV Local 1799 I(4) 4 scalar 590,591,592,593,1799,1800,1801,180 2 IDLAT Local 988 CHAR 1 scalar 1769,1771,1775,1784,1790 IDLON Local 988 CHAR 1 scalar 1761,1764,1776,1785,1791 IDY Local 956 I(4) 4 scalar 1462,1496,1499 IH Local 956 I(4) 4 scalar 1463,1467,1470,1474,1477,1491,1496 ,1499 IK Local 956 I(4) 4 scalar 1106,1107,1113,1114,1116,1117,1130 ,1135,1136,1141,1142,1143,1144,114 5,1148,1149,1151,1153,1154,1155,11 56,1181,1182,1184,1194,1195,1196,1 197,1229,1230,1232,1233,1234,1289, 1290,1292,1293,1359,1360,1361,1362 ,1363,1364,1365,1366,1368,1369,137 0,1371,1372,1373,1374,1375,1376,13 77,1378,1379,1380,1382,1383,1384,1 385,1386,1430,1431,1432,1443,1449, 1644,1645,1646,1673,1674,1675,1676 ,1677,1678,1680,1692,1694,1696,169 8,1700,1702,1714,1716,1718,1720,17 22,1724 IKH Local 955 I(4) 4 scalar 1202,1204 IKL Local 955 I(4) 4 scalar 1201,1203 IKM Local 955 I(4) 4 scalar 1192,1197,1201,1202,1203,1204,1208 ,1210,1211 IM Local 956 I(4) 4 scalar 1464,1467,1470,1474,1477,1491 IMTH Local 956 I(4) 4 scalar 1461,1496,1499 INT Func 1108 scalar 303,1108 IOUT Local 1753 I(4) 4 scalar 632,650,657,1753,1781 IP Local 957 I(4) 4 scalar IPASS Local 958 I(4) 4 scalar 958,1006,1033,1050 IS Local 956 I(4) 4 scalar 1465,1467,1470,1474,1477,1491 ISCALE Local 1345 I(4) 4 scalar 189,460,462,466,471,473,477,482,48 4,1345,1349,1353,1634,1636,1638,16 40,1643,1662,1664,1666,1668,1672 ISEA Local 957 I(4) 4 scalar ISP Local 955 I(4) 4 scalar 1135,1136,1137,1138,1139,1232,1233 ,1368,1369 ITH Local 955 I(4) 4 scalar 1134,1135,1136,1138,1139,1183,1184 ,1231,1232,1233,1234,1291,1292,129 3,1367,1368,1369,1370,1371,1372,13 73,1374,1375,1376,1377,1378,1379,1 380,1443,1449,1692,1694,1696,1698, 1700,1702,1714,1716,1718,1720,1722 ,1724 ITT Local 956 I(4) 4 scalar 1248,1307 ITYPE Local 1008 I(4) 4 scalar 331,335,337,346,359,360,425,426,45 6,457,551,552,607,627,652,653,663, 1008,1023,1024,1032,1083,1221,1226 ,1227,1287,1400,1456,1579,1731 IX Local 957 I(4) 4 scalar 1249,1253,1313,1322,1330 IY Local 957 I(4) 4 scalar 1250,1253,1313,1322,1330 IYR Local 956 I(4) 4 scalar 1460,1496,1499 J Local 955 I(4) 4 scalar 1093,1094,1098,1100,1101,1102,1103 Page 43 Source Listing W3EXPO 2014-09-16 17:02 Symbol Table multiwavespec.f90 Name Object Declared Type Bytes Dimen Elements Attributes References ,1104,1136,1137,1138,1139,1233,125 2,1262,1270,1274,1275,1278,1292,13 12,1320,1369,1406,1410,1411,1412,1 413,1414,1421,1422,1423,1425,1426, 1427,1440,1441,1442,1445,1446,1447 ,1468,1482,1485,1492,1493,1497,150 0,1508,1509,1510,1512,1513,1514,15 18,1519,1523,1524,1537,1540,1543,1 544,1545,1547,1548,1549,1554,1555, 1558,1559,1564,1565,1568,1569,1587 ,1593,1599,1603,1607,1611,1615,161 9,1626,1627,1628,1630,1631,1632,16 54,1655,1656,1658,1659,1660,1688,1 689,1690,1705,1706,1707,1709,1710, 1711,1736,1737,1738,1739,1740,1742 ,1743,1744,1753,1756,1757,1779,178 0,1784,1790,1802 LLWS Local 985 L(4) 4 1 0 1238,1246,1253,1256,1305,1310,1313 ,1322 M2KM Local 968 R(4) 4 scalar 997,999,1422,1426,1445,1446,1508,1 509,1512,1513,1518,1519,1523,1524, 1543,1544,1547,1548,1554,1555,1558 ,1559,1564,1565,1568,1569,1627,163 1,1655,1659,1706,1710,1738,1739,17 42,1743,1756,1757 MAX Func 1098 scalar 292,294,462,1098,1103,1145,1146,11 73,1174,1201,1209,1674 MIN Func 1202 scalar 462,1202,1480,1483,1486,1536,1538, 1541 MOD Func 1100 scalar 392,405,416,471,473,490,500,511,11 00,1104,1144,1172,1461,1462,1464,1 465,1759 N1MAX Local 949 I(4) 4 scalar 949,1109 NAR1D Param 949 I(4) 4 scalar 949 NDSBUL Local 1779 I(4) 4 scalar 229,576,577,591,666,667,668,1779,1 784,1787,1788,1789,1795,1800 NDSCBUL Local 1780 I(4) 4 scalar 230,579,580,592,669,670,671,1780,1 790,1792,1796,1801 NDSCSV Local 1802 I(4) 4 scalar 593,594,1802,1803 NDSO Local 1086 I(4) 4 scalar PTR 171,240,241,242,246,248,255,261,27 7,280,282,297,309,337,346,349,352, 360,366,368,370,372,375,377,379,38 2,388,393,395,408,426,436,438,440, 443,446,448,457,464,467,469,472,47 4,478,480,483,485,491,493,503,519, 521,525,527,529,532,534,536,546,55 2,557,563,567,578,581,586,595,614, 618,621,627,640,654,706,1086,1405, 1408,1411,1412,1413,1414,1415,1585 ,1591,1597,1601,1605,1609,1613,161 7 NDSTAB Local 1026 I(4) 4 scalar 363,364,384,385,387,397,399,401,40 5,410,412,414,415,428,429,431,432, 434,460,461,495,497,499,500,505,50 7,509,510,542,543,545,554,555,559, Page 44 Source Listing W3EXPO 2014-09-16 17:02 Symbol Table multiwavespec.f90 Name Object Declared Type Bytes Dimen Elements Attributes References 560,562,576,579,593,666,669,1026,1 028,1035,1037,1039,1041,1043,1045, 1050,1053,1055,1059,1061,1064,1066 ,1069,1071,1075,1077,1421,1425,143 1,1435,1440,1443,1445,1448,1467,14 70,1474,1477,1491,1496,1499,1508,1 512,1518,1523,1529,1533,1543,1547, 1554,1558,1564,1568,1626,1630,1635 ,1637,1639,1641,1645,1648,1654,165 8,1663,1665,1667,1669,1676,1680,16 83,1688,1691,1693,1695,1697,1699,1 701,1705,1709,1713,1715,1717,1719, 1721,1723,1738,1742,1748,1779,1780 ,1802 NINT Func 1497 scalar 1497,1500,1555,1559,1565,1569 NK Local 973 I(4) 4 scalar PTR 269,400,401,413,414,498,499,508,50 9,973,974,975,976,977,978,979,980, 981,982,983,984,1106,1130,1161,116 5,1166,1181,1191,1192,1194,1202,12 29,1289,1359,1405,1408,1429,1430,1 434,1443,1449,1475,1478,1531,1534, 1585,1591,1597,1598,1601,1602,1605 ,1606,1609,1610,1613,1614,1617,161 8,1644,1673,1692,1694,1696,1698,17 00,1702,1714,1716,1718,1720,1722,1 724 NOPTS Local 1093 I(4) 4 scalar PTR 171,278,313,321,340,569,579,588,61 5,664,669,1093,1780 NPART Local 957 I(4) 4 scalar 1222,1739,1743,1747,1752 NREQ Local 1033 I(4) 4 scalar 315,322,400,413,498,508,614,1033,1 048,1458 NSPEC Local 985 I(4) 4 scalar PTR 985 NTH Local 974 I(4) 4 scalar PTR 269,400,405,413,416,498,500,508,51 1,974,975,976,977,978,979,980,981, 982,1134,1135,1183,1231,1232,1291, 1292,1367,1368,1408,1443,1449,1597 ,1601,1605,1609,1613,1617,1692,169 4,1696,1698,1700,1702,1714,1716,17 18,1720,1722,1724 OTYPE Local 1023 I(4) 4 scalar 362,365,381,389,419,428,435,437,43 9,442,445,447,450,459,463,465,476, 487,514,523,540,554,556,565,566,57 2,585,600,627,652,653,663,1023,102 4,1034,1036,1038,1040,1042,1044,10 51,1057,1063,1065,1067,1073,1084,1 226,1227,1267,1402,1417,1437,1466, 1469,1473,1476,1479,1494,1506,1516 ,1528,1532,1535,1551,1581,1622,165 0,1685,1733,1751,1778,1798 PRT1DM Subr 948 948,1591 PRT1DS Subr 948 948,1405,1585 PRT2DS Subr 948 948,1408,1597,1601,1605,1609,1613, 1617 PTLOC Local 1274 R(4) 4 2 1 PTR 171,280,282,341,342,618,619,621,62 2,1274,1275,1278,1422,1426,1440,14 Page 45 Source Listing W3EXPO 2014-09-16 17:02 Symbol Table multiwavespec.f90 Name Object Declared Type Bytes Dimen Elements Attributes References 41,1445,1446,1508,1509,1512,1513,1 518,1519,1523,1524,1543,1544,1547, 1548,1554,1555,1558,1559,1564,1565 ,1568,1569,1627,1631,1655,1659,168 8,1689,1706,1710,1736,1737,1738,17 39,1742,1743,1756,1757 PTNME Local 1406 CHAR 10 1 1 PTR 171,280,282,341,577,578,580,581,59 4,595,618,621,1406,1410,1421,1425, 1440,1445,1587,1593,1599,1603,1607 ,1611,1615,1619,1626,1630,1654,165 8,1688,1705,1709,1739,1743,1784,17 90 R1 Local 959 R(4) 4 scalar 1111,1112,1113,1114 R2 Local 959 R(4) 4 scalar 1112,1113,1114 RADE Param 1100 R(4) 4 scalar 1100,1101,1104,1144,1145,1172,1173 ,1175,1212 RADIUS Param 1277 R(4) 4 scalar 1277,1279 REAL Func 1111 scalar 1111 SBT Local 978 R(4) 4 2 0 1373,1375,1380,1613,1700,1722 SBT1 Local 984 R(4) 4 1 0 1365,1380,1386,1389,1393,1646,1678 SCALE1 Local 1404 R(4) 4 scalar 362,367,369,372,459,524,526,529,14 04,1406,1583,1592 SCALE2 Local 1407 R(4) 4 scalar 362,374,376,379,459,531,533,536,14 07,1409,1595,1602,1606,1610,1614,1 618 SDS Local 978 R(4) 4 2 0 1372,1374,1379,1609,1698,1720 SDS1 Local 983 R(4) 4 1 0 1364,1379,1385,1389,1392,1646,1678 SIG Local 1107 R(4) 4 1 1 PTR 401,414,499,509,1107,1116,1117,114 2,1154,1161,1165,1166,1208,1230,12 90,1360,1405,1408,1429,1431,1434,1 585,1591,1598,1602,1606,1610,1614, 1618,1645,1676,1680 SIGN Func 1209 scalar 1209 SIN1 Local 983 R(4) 4 1 0 1362,1377,1383,1389,1390,1646,1677 SIX Local 959 R(4) 4 scalar 1107,1108,1111 SNL Local 977 R(4) 4 2 0 1371,1374,1378,1605,1696,1718 SNL1 Local 983 R(4) 4 1 0 1363,1378,1384,1389,1391,1646,1677 SPBND Local 975 R(4) 4 1 0 1145,1149,1211,1432 SPCO Local 1136 R(4) 4 2 1 PTR 172,1136,1137,1138,1139,1233,1292, 1369 SPP Local 964 R(4) 4 scalar 1211,1212,1216,1472,1521,1526 SQRT Func 1099 scalar 1099,1145,1146,1168,1173 SQRTH Local 959 R(4) 4 scalar 1099,1107,1114 STME21 Subr 1049 167,296,301,348,651,1049,1085,1419 ,1624,1652 STT Local 981 R(4) 4 2 0 1374,1617,1702,1724 STT1 Local 984 R(4) 4 1 0 1366,1389,1394,1646,1678 SWN Local 977 R(4) 4 2 0 1370,1374,1377,1601,1694,1716 TAUWNX Local 970 R(4) 4 scalar 1253,1313,1321 TAUWNY Local 970 R(4) 4 scalar 1253,1313,1321 TAUWX Local 969 R(4) 4 scalar 1240,1246,1252,1256,1300,1305,1310 ,1312,1321 TAUWY Local 969 R(4) 4 scalar 1241,1246,1252,1256,1301,1305,1310 ,1312,1321 THBND Local 975 R(4) 4 1 0 1144,1148,1210,1432 Page 46 Source Listing W3EXPO 2014-09-16 17:02 Symbol Table multiwavespec.f90 Name Object Declared Type Bytes Dimen Elements Attributes References THMEAN Local 962 R(4) 4 scalar 1172,1179,1415,1471,1520,1525 THP Local 963 R(4) 4 scalar 1210,1215,1472,1521,1526 THSPRD Local 963 R(4) 4 scalar 1173,1175,1180,1415,1471,1520,1525 TIME Local 1026 I(4) 4 1 1 PTR 170,348,636,1026,1028,1049,1085,14 19,1460,1461,1462,1463,1464,1465,1 467,1470,1474,1477,1491,1624,1652, 1738,1742 TIMEV Local 1753 I(4) 4 1 2 554,1753,1782,1785,1791 TMEAN Local 962 R(4) 4 scalar 1171,1178,1415,1471,1484,1520,1525 ,1539 TPI Param 1142 R(4) 4 scalar 405,416,500,511,1142,1170,1171,136 0,1484,1497,1539,1556,1560 TPIINV Param 1151 R(4) 4 scalar 401,414,499,509,1151,1161,1165,116 6,1230,1273,1290,1346,1350,1354,14 31,1472,1520,1525 UABS Local 960 R(4) 4 scalar 1103,1222,1245,1251,1255,1304,1309 ,1311,1319,1350,1351,1352,1411,148 0,1483,1486,1536,1538,1541,1752 UDIR Local 960 R(4) 4 scalar 1100,1412,1423,1427,1442,1447,1468 ,1492,1497,1500,1510,1514,1544,154 8,1555,1559,1565,1569,1690,1707,17 11,1740,1744,1752 UDIRCA Local 968 R(4) 4 scalar 1101,1222 UDIRR Local 960 R(4) 4 scalar 1102,1245,1252,1255,1304,1309,1312 ,1317,1320 UNORM Local 964 R(4) 4 scalar 1263,1268,1270,1272,1273,1274,1275 ,1475,1478,1530,1533 USTAR Local 964 R(4) 4 scalar 1242,1245,1251,1255,1268,1298,1304 ,1309,1311,1317,1319,1329,1354,135 5,1356,1423,1427,1628,1632,1656,16 60 USTD Local 967 R(4) 4 scalar 1245,1255,1299,1304,1309,1329 VAR1 Local 987 CHAR 4 1 5 990,1592 W3ARRYMD Module 948 948 W3BULL Subr 950 950,1752 W3BULLMD Module 950 950 W3DISPMD Module 946 946,949 W3EXPO Subr 830 655 W3PART Subr 945 945,1221 W3PARTMD Module 945 945 W3SBT1 Subr 1336 1336 W3SBT1MD Module 943 943 W3SDB1 Subr 1332 1332 W3SDB1MD Module 944 944 W3SDS4 Subr 1329 1329 W3SIN4 Subr 1251 1251,1311,1319 W3SLN1 Subr 1317 1317 W3SLN1MD Module 940 940 W3SNL1 Subr 1325 1325 W3SNL1MD Module 942 942 W3SPR4 Subr 1244 1244,1254,1303,1308 W3SRC4MD Module 941 941 WAO Local 1103 R(4) 4 1 1 PTR 172,1103,1262,1270,1412,1423,1427, 1441,1447,1468,1492,1497,1500,1510 ,1514,1544,1548,1555,1559,1565,156 Page 47 Source Listing W3EXPO 2014-09-16 17:02 Symbol Table multiwavespec.f90 Name Object Declared Type Bytes Dimen Elements Attributes References 9,1628,1632,1656,1660,1689,1707,17 11,1740,1744 WAVNU1 Subr 946 946,1482,1485,1537,1540 WDO Local 1100 R(4) 4 1 1 PTR 172,1100,1101,1102 WHITECAP Local 970 R(4) 4 1 4 1330 WLEN Local 962 R(4) 4 scalar 1170,1177,1415,1471,1519,1524 WN Local 973 R(4) 4 1 0 1113,1116,1153,1222,1234,1244,1254 ,1293,1303,1308,1317,1329,1336 WN2 Local 976 R(4) 4 2 0 1234,1251,1293,1311,1319,1332 WNA Local 965 R(4) 4 scalar 1482,1483,1485,1486,1537,1538,1540 ,1541 WNMEAN Local 968 R(4) 4 scalar 1245,1255,1304,1309,1325,1333 X Local 968 R(4) 4 scalar 1756,1759,1760,1763,1784,1791 XBS Local 981 R(4) 4 2 0 1016,1373 XBT Local 980 R(4) 4 2 0 1015,1336,1373 XDB Local 980 R(4) 4 2 0 1014,1333,1372 XDS Local 980 R(4) 4 2 0 1013,1329,1372 XFR Local 1002 R(4) 4 scalar PTR 1002,1003 XH Local 960 R(4) 4 scalar 1003,1005,1205 XH2 Local 960 R(4) 4 scalar 1005,1208 XIN Local 979 R(4) 4 2 0 1010,1253,1313,1322,1370 XL Local 960 R(4) 4 scalar 1002,1004,1205 XL2 Local 960 R(4) 4 scalar 1004,1208 XLN Local 978 R(4) 4 2 0 1009,1317,1370 XNL Local 979 R(4) 4 2 0 1011,1325,1371 XPART Local 1222 R(4) 4 2 1 ALC 270,1222,1748,1752 XSTAR Local 966 R(4) 4 scalar 1274,1277,1529,1533 XTR Local 979 R(4) 4 2 0 1012,1371 XWL Local 982 R(4) 4 2 0 1017 XXX Local 981 R(4) 4 2 0 1018,1375 XYZ Local 966 R(4) 4 scalar 1482,1485,1537,1540 Y Local 968 R(4) 4 scalar 1757,1768,1772,1784,1790 YSTAR Local 967 R(4) 4 scalar 1275,1279,1530,1533 Z0 Local 967 R(4) 4 scalar 1246,1252,1256,1305,1310,1312,1321 ZWND Local 967 R(4) 4 scalar 1239,1297 ZZWND Local 1239 R(4) 4 scalar PTR 1239,1297 Page 48 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 2060 !/ 2061 !/ End of WAVESPEC ----------------------------------------------------- / 2062 !/ 2063 END PROGRAM WAVESPEC ENTRY POINTS Name MAIN__ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 1000 Label 792 681 1001 Label 796 685 1002 Label 799 689 1003 Label 803 693 1004 Label 807 697 1005 Label 811 701 1010 Label 815 607 1011 Label 818 419,450,514,600 1012 Label 821 572 1940 Label 766 368,375,525,532 1941 Label 767 370,377,527,534 1942 Label 768 372,379,529,536 1943 Label 769 395,408 1944 Label 770 412 1945 Label 771 414 1946 Label 772 415 1947 Label 773 388,563,578,581,595 2940 Label 775 436,438,440,443,446,448 3940 Label 778 519 3941 Label 779 546 3943 Label 780 493,503 3944 Label 781 507 3945 Label 782 509 3946 Label 783 510 410 Label 317 325 800 Label 680 252 801 Label 684 253,291,319,331,362,428,459,554 802 Label 688 253,291,319,331,362,428,459,554 803 Label 692 387,434,545,562,577,580,594 804 Label 696 397,410,495,505 805 Label 700 339 888 Label 704 353,676 900 Label 712 248 901 Label 714 255 920 Label 716 261 930 Label 718 277 931 Label 720 280 932 Label 721 282 940 Label 723 297 941 Label 726 309 Page 49 Source Listing W3EXPO 2014-09-16 17:02 Symbol Table multiwavespec.f90 Name Object Declared Type Bytes Dimen Elements Attributes References 942 Label 728 337,346,360,426,457,552 943 Label 731 366,382,393,464,467,469,478,480,49 1,557,567,586 944 Label 732 472,474,483,485 945 Label 733 341 948 Label 734 349 949 Label 735 352,640 950 Label 737 614 951 Label 739 618 952 Label 741 627 953 Label 740 621 960 Label 785 654 961 Label 743 670 962 Label 745 671 971 Label 751 667 974 Label 755 668 999 Label 788 706 BHSDROP Param 173 R(4) 4 scalar 173,668 BHSMIN Param 173 R(4) 4 scalar 173,668 COMSTR Local 197 CHAR 1 scalar 253,254,255,290,318,330,361,427,45 8,553 CONSTANTS Module 158 158 COUNT Func 570 scalar 570 DMT Local 174 R(4) 4 2 12 174 DSEC21 Func 167 R(4) 4 scalar 167,636 DTEST Local 192 R(4) 4 scalar 636,637,645 DTREQ Local 192 R(4) 4 scalar 291,292,293,300,302,303,646,656 EXTCDE Subr 166 166,420,451,515,573,601,608,682,68 6,690,694,698,702 FNMPRE Local 172 CHAR 80 scalar 172,250,251,338,386,387,396,397,40 9,410,433,434,494,495,504,505,544, 545,561,562,568,587 GRDID Local 176 CHAR 10 1 1 PTR 176,342 HST Local 174 R(4) 4 2 12 174 I Local 184 I(4) 4 scalar 278,280,282,340,341,342,518,519,61 5,616,618,619,621,622 IDDDAY Local 197 CHAR 11 scalar 303,305,307 IDSRCE Local 199 CHAR 25 1 6 204,519 IDTIME Local 197 CHAR 23 scalar 296,297,301,307,308,309,348,349,65 1,654 IERR Local 184 I(4) 4 scalar 252,339,387,398,411,434,496,506,54 5,562,577,580,594,681,689,693,697, 701 IJ Local 187 I(4) 4 scalar 569,575,576,577,578,579,580,581,58 8,589,593,594,595,664,665,666,669 IK Local 186 I(4) 4 scalar 401,414,499,509 IOTEST Local 186 I(4) 4 scalar 275,350,351,638,639 IPOINT Local 185 I(4) 4 scalar 319,320,321,322,323 ITH Local 186 I(4) 4 scalar 405,416,500,511 ITRACE Subr 166 166,234 IYY Local 173 L(4) 4 1 80 173 J Local 186 I(4) 4 scalar 250,251,338,386,387,396,397,409,41 0,433,434,494,495,504,505,544,545, 561,562,568,587 LEN_TRIM Func 250 scalar 250,386,396,409,433,494,504,544,56 Page 50 Source Listing W3EXPO 2014-09-16 17:02 Symbol Table multiwavespec.f90 Name Object Declared Type Bytes Dimen Elements Attributes References 1,568,587 M2KM Local 193 R(4) 4 scalar 264,266,280,282,618,619,621,622 NDBO Local 191 I(4) 4 scalar 245,338,341,344 NDSE Local 171 I(4) 4 scalar PTR 171,241,290,318,330,361,419,427,45 0,458,514,553,572,600,607,681,685, 689,693,697,701 NDSI Local 183 I(4) 4 scalar 226,239,251,253,290,291,318,319,33 0,331,361,362,427,428,458,459,553, 554 NDSM Local 183 I(4) 4 scalar 227,243,260 NDSOP Local 183 I(4) 4 scalar 228,244,275,350,638 NDST Local 171 I(4) 4 scalar PTR 171,242 NDSTRC Local 183 I(4) 4 scalar 232,234,246 NEXTLN Subr 166 166,290,318,330,361,427,458,553 NFLD Param 173 I(4) 4 scalar 173 NOUT Local 184 I(4) 4 scalar 291,293,294,309,657 NPMAX Param 173 I(4) 4 scalar 173 NPTAB Param 173 I(4) 4 scalar 173 NTRACE Local 183 I(4) 4 scalar 233,234 PI Param 405 R(4) 4 scalar 405,416,500,511 TABNME Local 198 CHAR 9 scalar 383,385,387,388,430,432,434,436,43 8,441,444,446,448,541,543,545,546, 558,560,562,563 TDUM Local 184 I(4) 4 1 2 299,300,301 TFNAME Local 198 CHAR 16 scalar 390,391,395,397,408,410,488,489,49 3,495,503,505 TH Local 405 R(4) 4 1 1 PTR 405,416,500,511 TICK21 Subr 167 167,300,646,656 TOUT Local 184 I(4) 4 1 2 291,296,392,490,636,646,651,656 TPT Local 174 R(4) 4 2 12 174 TRIM Func 577 scalar 577,578,580,581,594,595 W3ADATMD Module 162 162 W3BULLMD Module 173 173,175 W3GDATMD Module 169 169 W3IOGR Subr 164 164,260 W3IOGRMD Module 164 164 W3IOPO Subr 165 165,275,350,638 W3IOPOMD Module 165 165 W3NAUX Subr 162 162,221 W3NDAT Subr 161 161,219 W3NMOD Subr 217 217 W3NOUT Subr 163 163,223 W3ODATMD Module 163 163,171,176 W3SERVMD Module 166 166 W3SETA Subr 162 162,222 W3SETG Subr 218 218 W3SETO Subr 163 163,224 W3SETW Subr 161 161,220 W3TIMEMD Module 167 167 W3WDATMD Module 161 161,170 WAVESPEC Prog 2 Page 51 Source Listing W3EXPO 2014-09-16 17:02 Subprograms/Common Blocks multiwavespec.f90 SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References W3EXPO Subr 830 655 WAVESPEC Prog 2 COMPILER OPTIONS BEING USED -align nocommons -align nodcommons -align noqcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume byterecl -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_unit_star -assume old_ldout_format -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume nostd_mod_proc_name -assume norealloc_lhs -assume underscore -assume no2underscores -auto no -auto_scalar no -bintext -ccdefault default -check noargs -check noarg_temp_created -check nobounds -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check power -check noshape -check nounderflow -check nouninitialized -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1210 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __pentium4 -D __pentium4__ -D __tune_pentium4__ -D __SSE2__ -D __SSE3__ -D __SSSE3__ -D __SSE4_1__ -D __SSE4_2__ -D __SSE__ -D __MMX__ -D __AVX__ -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals no -fpconstant Page 52 Source Listing W3EXPO 2014-09-16 17:02 multiwavespec.f90 -fpe3 -fprm nearest no -ftz -fp_model noprecise -fp_model fast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -heap_arrays 0 no -threadprivate_compat -free -g0 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude -O2 no -pad_source -real_size 32 no -recursive -reentrancy none no -sharable_localsaves -vec=simd -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w noargument_checking -w nodeclarations -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage -includepath : /gpfs/gp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/,.f,./.f,/usrx/local/intel/composerxe/mkl/include/.f, /usrx/local/intel/composerxe/tbb/include/.f,/gpfs/gp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/.f, /gpfs/gp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/.f,/usr/local/include/.f,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/.f, /usr/include/.f,/usr/include/.f -list filename : multiwavespec.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100