--- Welcome to the official ADCIRCWiki site! The site is currently under construction, with limited information. ---
For general information, see the ADCIRC site, ADCIRC FAQ, or the ADCIRC Wikipedia page. For model documentation not yet available on the wiki, see the ADCIRC site. New content is being continuously added to the ADCIRCWiki, and material from the main ADCIRC site will be gradually transitioned over to the wiki.

NWS12

From ADCIRCWiki
Revision as of 17:30, 15 June 2020 by Alexcrosby (talk | contribs) (Created page with "=NWS12= Oceanweather Inc. (OWI) "WIN" and "PRE" file formats as supported in ADCIRC ==fort.22== NWS12 uses the '''fort.22''' file as a control file because the OWI format s...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

NWS12

Oceanweather Inc. (OWI) "WIN" and "PRE" file formats as supported in ADCIRC

fort.22

NWS12 uses the fort.22 file as a control file because the OWI format supports individual files for wind and pressure inputs as well as multiple grid overlays of varying resolutions.

fort.22 layout:

 NWSET
 NWBS
 DWM
 <...>

Where:

  • NWSET is the number of nested grid overlays
  • NWBS is the number of "blank" (constant 0 m/s wind speed and 1013mb pressure) timesteps to insert
  • DWM is a wind multiplier to modify the wind speed magnitudes. Generally this should be set to 1 if using inputs sourced from OWI or if winds can be considered "marine exposure", with an averaging period between 10-min and 3-hrs at a height of 10 meters.

File Formatting

Winds and pressure data formats are similar. The header format is the same, but in the wind file the header is followed by U then V components while in the pressure file the header is followed by just pressures.

The file begins with a header indicating the starting and ending dates and is followed by a grid/date header for each time step and the u and v components of the wind in meters/second or pressures in millibars. Starting/Ending dates are in YYYYMMDDHH format where:

      YYYY    4-character Year
	MM    2-character Month
	DD    2-character Day
	HH    2-character Hour


example WIN:

 OWI WWS Wind Output Ucomp,Vcomp in m/s           Start:1995060600 End:1995060600
 iLat=  67iLong=  67DX= 1.250DY=  .833SWLat=  22.500SWlon= -82.500Dt=199506060000
 -1.16856  -1.06439   -.84875  -1.03460  -1.50047  -2.09462  -2.80243  -3.55863
 -4.24125  -4.84273  -5.59486  -5.37088  -5.30224  -5.12534  -4.89537  -4.67412
 -4.49203  -4.35772  -4.26612  -4.20260  -4.14746  -4.08396  -4.00686  -3.92213
 -3.83615  -3.74765  -3.65182  -3.54998  -3.45299  -3.37660  -3.32959  -3.28037
 -3.10631  -2.67723  -2.08363  -1.53773  -1.12623   -.83526   -.62870   -.47371
 [rest deleted]
  • iLat is the number of parallels
  • iLong is the number of meridions
  • DX is the grid spacing in degrees of longitude
  • DY is the grid spacing in degrees of latitude
  • SWLat is the latitude of the South West corner
  • SWlon is the longitude of the South West corner
  • Dt is the date/time in YYYYMMDDHHmm (same as master header date format but with mm Minutes as well)

The number of grid points is iLat*iLong, the u component of the winds in meters/second is followed by the v component.

Sample code

Sample fortran to read a win file (first time step only):

c     Read in begining/ending dates of win file
10    format (t56,i10,t71,i10)      
      read (20,10) date1,date2
c     Read Grid Specifications/Date
11    format (t6,i4,t16,i4,t23,f6.0,t32,f6.0,t44,f8.0,t58,f8.0,t69,i10,i2)
      read (20,11) iLat, iLong, dx, dy, swlat, swlong, lCYMDH, iMin
c     Read U/V Components of the wind
12    format (8f10.0)
      read (20,12) ((uu(i,j),i=1,ilong),j=1,ilat)
      read (20,12) ((vv(i,j),i=1,ilong),j=1,ilat)

Latitude and Longitude for each point can be calculated as follows:

      do 20  icnt = 1,iLat
        slat(icnt) = SWlat + (icnt - 1) * DY
20    continue
      do 30  jcnt = 1,iLong
        slon(jcnt) = SWlong + (jcnt - 1) * DX
30    continue