April 5, 2011 Uploaded here under Inundation/bio/ is a sample test run "ww3_swan_2007" This one of the three storms chosen for the extratropical domain. As in the previous test run this is 5 nested grids that start with the large NW Atlantic grid1 (Resolution = 12 min) and zoom in to Scituate Harbour grid5 (resolution = 1.44 sec or ~ 38 m). The main changes from the first test run are: i) improved smoothed bathymetry for grid 5 ii) higher model spectral resolution, now we use 29 frequencies starting at ~0.041 and goes to ~0.594 Hz, (instead of 25 frequencies starting at ~0.041 and goes to ~0.406 Hz) and 36 directions with delta-theta=10 deg. (instead of 24 directions with delta-theta=15 deg.) I use in this test run 1) the latest release of WAVEWATCH III ver. 3.14 for grid1 & grid2, which are run together in ww3-multi grid mode and with two-way nesting between the two grids involved. From this run I output and save the 2D spectra at all of grid3 water boundary points at grid2 resolution (= 2 min). 2) the latest release of SWAN ver. 40.81 for grids 3, 4 & 5. SWAN run over grid3 takes its input boundary data from WW3 grid2, and output boundary data for grid4. Then SWAN run over grid4 takes its input boundary data from SWAN grid3 and output boundary data for grid5. Last SWAN run over grid5 takes its input boundary data from SWAN grid4 and it doesn't have to output any nest information. #------------------------------------------------------------------------------- Below is a description of the grids domain, resolution, mesh size, time step and the wind forcing. WW3 grid1: NW Atlantic - forced by interp. NARR winds -80.00 ==> -50.00 W (long. range = 30.0 deg) +25.00 ==> +50.00 N (lati. range = 25.0 deg) dx = 0.20 deg = 12.0 arc min ==> 151 x 126 mesh WW3 dt = 600 sec (CFL criteria) WW3 grid2: Gulf of Main (~ GOM3) - forced by interp. WRF winds -73.00 ==> - 63.00 W (long. range = 10.0 deg) +38.00 ==> +46.00 N (lati. range = 8.0 deg) dx = 0.05 deg = 3.0 arc min ==> 201 x 161 mesh WW3 dt = 120 sec (CFL criteria) SWAN grid3: Inside Gulf of Main - forced by interp. WRF winds -71.50 ==> -69.50 W (long. range = 2.00 deg) +41.50 ==> +43.00 N (lati. range = 1.50 deg) dx = 0.01 deg = 36.0 arc sec ==> 201 x 151 mesh WW3 dt = 30 sec (CFL criteria) SWAN dt = 360 sec (12 x CFL criteria) SWAN grid4: Scituate Harbour + approaches - forced by interp. WRF winds -70.90 ==> -70.50 W (long. range = 0.40 deg) +42.00 ==> +42.35 N (lati. range = 0.35 deg) dx = 0.002 deg = 7.20 arc sec ==> 201 x 176 mesh WW3 dt = 5 sec (CFL criteria) SWAN dt = 60 sec (12 x CFL criteria) SWAN grid5: Scituate Harbour - forced by interp. WRF winds -70.76 ==> -70.70 W (long. range = 0.06 deg) +42.17 ==> +42.23 N (lati. range = 0.06 deg) dx = 0.0004 deg = 1.44 arc sec ==> 151 x 151 mesh WW3 dt = 1 sec (CFL criteria) SWAN dt = 12 sec (12 x CFL criteria) #------------------------------------------------------------------------------- This test simulation was for ~6-day hindcast for April 2007 storm from 15 April 2007 00:00UTC to 20 April 2007 21:00UTC So, I placed the input/output files under the directory path inundation/bio/ww3_swan_2007/ in two separate directories: input_files/ and output_files/ and kept the packed tar file for each. #------------------------------------------------------------------------------- Under input_files/ directory you'll find all the depth files used for all 5 grids. These are obtained by interp. from ETOPO1 and the unstructured depth files for GOM3 & Scituate Harbour. You'll also find all the wind files used in all 5 grids. For grid1 I use the interp. NARR winds and for all the remaining 4 grids I use the interp. WRF winds. Under output_files directory you'll find 5 subdirectories grd1, grd2, grd3, grd4 & grd5 In WW3 grids (grid1 & grid2) you'll find WW3 type output 1) WW3 gridded fields output in subdirectory outf/ Eight fields were selected all in ASCII format: 1.1) *.hs for significant wave height (m) 1.2) *.dir for mean wave direction (deg) 1.3) *.t for mean wave period (sec) 1.4) *.fp for peak wave frequency (Hz) 1.5) *.dp for peak wave direction (deg) 1.6) *.fpl for wind sea frequency (Hz) 1.7) *.dpl for wind sea direction (deg) 1.8) *.wlv for water level 2) WW3 point output in subdirectory outp/ Three type of point output were done also all in ASCII format: 2.1) 1D energy at the selected location in tab12.ww3 2.2) 2D energy at the selected location in ww3.*.spc 2.3) time series of mean wave parameters at the selected loc in tab22.ww3 3) WW3 nest output is a point output type 2.2 at all water boundary points. The file is very large that I kept it out. I only kept the input cmd file and corresponding output. 2.2) 2D energy at all water boundary points in ww3.*.spc Note: For grd1 & grd2 the selected buoys locations are put in files 'stns_in_grd#' under grd#/ directory #------------------------------------------------------------------------------- In SWAN grids (grid3, grid4 & grid5) you'll find SWAN type output 1) SWAN gridded fields output are available for grd3, grd4 & grd4 Seven fields were selected in ASCII format: 1.1) grd#.hs for significant wave height (m) 1.2) grd#.dir for mean wave direction (deg) 1.3) grd#.tm01 for mean absolute wave period (sec) 1.4) grd#.tm02 for mean absolute zero crossing wave period (sec) 1.5) grd#.pper for peak wave period (sec) 1.6) grd#.pdir for peak wave direction (deg) 1.7) grd#.swl for swell significant wave height (m) Same 7 fields were also outputted in Matlab format: 1.1) grd#_hs.mat for significant wave height (m) 1.2) grd#_dir.mat for mean wave direction (deg) 1.3) grd#_tm01.mat for mean absolute wave period (sec) 1.4) grd#_tm02.mat for mean absolute zero crossing wave period (sec) 1.5) grd#_pper.mat for peak wave period (sec) 1.6) grd#_pdir.mat for peak wave direction (deg) 1.7) grd#_swl.mat for swell significant wave height (m) 2) SWAN point output Three type of point output were done 2.1) 1D energy at the selected location in *.spc1 only at Buoys loc. 2.2) 2D energy at the selected location in *.spc2 only at Buoys loc. 2.3) time series of mean wave parameters at the selected loc in *.tbl at all buoys loc. & watermark stations. 3) SWAN nest output 3.1) grd#.nst for 2D spectra at the boundary of next SWAN nested grid. These are large binary files that I left out. Note: For each grid the selected buoys locations and/or water mark stations are put in file 'stns_in_grd#' under grd#/ directory #------------------------------------------------------------------------------- #===============================================================================