--- 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: Difference between revisions

From ADCIRCWiki
Jump to navigation Jump to search
(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...")
 
(add wtiminc and fort.22 and fort.221-4 explainations)
Line 12: Line 12:
   NWBS
   NWBS
   DWM
   DWM
   <...>
   <...to be expanded in new verions>


Where:  
Where:  
Line 18: Line 18:
* '''NWBS''' is the number of "blank" (constant 0 m/s wind speed and 1013mb pressure) timesteps to insert
* '''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.
* '''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.
NWS12 WIN/PRE text files containing wind and pressure data should be named fort.22[1-4], as follows in order of increasing precedence in application to the ADCIRC mesh:
* ''fort.221'': Largest scale pressure (mb) fields
* ''fort.222'': Largest scale wind vector (m/s) fields
* ''fort.223'': Finer scale pressure (mb) fields
* ''fort.224'': Finer scale wind vector (m/s) fields
* ...to be expanded in new versions
==WTIMINC in fort.15==
Although the WIN/PRE files contain datetime and timestepping information for the input data, ''WTIMINC'' must be set to the timestep of the input files. This limitation locks all nested grid inputs to the same input timesteps, and requires that the input timesteps are constant for the duration of the input files.


== File Formatting==
== File Formatting==

Revision as of 13:40, 15 June 2020

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
 <...to be expanded in new verions>

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.

NWS12 WIN/PRE text files containing wind and pressure data should be named fort.22[1-4], as follows in order of increasing precedence in application to the ADCIRC mesh:

  • fort.221: Largest scale pressure (mb) fields
  • fort.222: Largest scale wind vector (m/s) fields
  • fort.223: Finer scale pressure (mb) fields
  • fort.224: Finer scale wind vector (m/s) fields
  • ...to be expanded in new versions

WTIMINC in fort.15

Although the WIN/PRE files contain datetime and timestepping information for the input data, WTIMINC must be set to the timestep of the input files. This limitation locks all nested grid inputs to the same input timesteps, and requires that the input timesteps are constant for the duration of the input files.

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