...
At present, for different wind input files, different headers are used, which are read in different ways in the sourcecode. Some effeort has been put to in making the files more self-descriptive, readible and to make the input of the different files more consistent. For three wind-input options, a suggestion to standardise the headers is shown below. The suggestion is conform the possible transition NetCDF format in the future. See NetCDF standard names for more information on the used standard NetCDF names. Further below (in section: Integrated standardisation), an An integrated file format is proposed which can be used for all wind input options. Comments are indicated by a # (this is open for discussion).
- The Pseudo ArcInfo wind input file, with keyword in MD-file Filwp
Code Block |
---|
# START OF HEADER
# This file is created by Deltares
n_cols# Additional commments
# General part:
t_ref = 201
n_rows = 19991026.000000 = 173
x0 = -12.000
y0 # = 48.000
value_posReference time
t_ref_unit = centre
#value_posModifiedJulianDate = corner
cellsize = 0.12500 0.083333333
x_unit # Unit used =for reference time degree
#xt_unit = m
t_refmin = 20080312.000000
t_ref_unit = ModifiedJulianDate
t_unit = hrs
#t_unit # =Time unit (sec, ModifiedJulianDate
n_quantitymin, hrs, etc.)
#t_unit = 3
quantityhrs = x_wind y_wind air_pressure
unit
missing_value = m s-1-999.000 m s-1 Pa
missing_value = -999.000
# END OF HEADER
TIME = 0.0
1000.776001 1000.687927 1000.599976 1000.512024 1000.423950 1000.334961 1000.245911 1000.156921 1000.067932
TIME = 120.0
1000.786001 1000.637927 1000.579976 1000.582024 1000.424950 1000.324961 1000.265911 1000.159921 1000.060932
|
- The wind on a curvi-linear grid input file, with keyword in MD-file Fwndgp
Code Block |
---|
# START OF HEADER
# This file is created by Deltares
grid_file = hirlam.grd
first_row = nmax
last_row = 1
first_col = 1
last_col = mmax
t_ref = 20080312.000000
t_ref_unit = ModifiedJulianDate
t_unit = hrs
n_quantity = 1
quantity = air_pressure
unit = Pa
missing_value = -999.0
# END OF HEADER
TIME = 2130
101894.125 101880.625 101867.125 101853.625 101840.625 101827.125 101813.125 101798.625 101784.625 101769.625 101754.125 101739.125 101723.625
TIME = 2273.0
100038.75 100037.25 100035.75 100035.25 100035.25 100035.25 100036.5 100038.0 100039.75 100042.25 100045.0 100048.0 100052.0 100056.0 100060.0
|
- The Spiderweb wind input file, with keyword in MD-file Filweb
Code Block |
---|
# START OF HEADER
# This file is created by Deltares
# Spiders web derived from TRACK file: 19991026_0000.imd
# Wind direction winddir = where it blows from, degrees from true north
n_cols = 400
n_rows = 16
radius = 600000.00
rad_unit = m
x_unit = degree
#x_unit = m
t_ref = 19991026.000000
t_ref_unit = ModifiedJulianDate
t_unit = min
#t_unit = hrs
n_quantity = 3
quantity = wind_speed wind_from_direction p_drop
unit = m s-1 degree Pa
missing_value = -999.000
# END OF HEADER
TIME = 0.0000000E+00
x0 = 115.4
y0 = 18.9
p_drop_eye = -5300
1.388999 1.388999 1.388999 1.388999 1.388999 1.388999 1.388999 1.388999 1.388999 1.388999
155.0000 346.2500 357.5000 8.750000 20.00000 31.25000 42.50000 53.75000 65.00000 76.25000
5300.000 5300.000 5300.000 5300.000 5300.000 5300.000 5300.000 5300.000 5300.000 5300.000
TIME = 1440
x0 = 114.4000
y0 = 18.90000
p_drop_eye = -5300.000
2.435253 2.435253 2.435253 2.435253 2.435253 2.435253 2.435253 2.435253 2.435253 2.435253
65.01978 75.11859 85.25320 95.45370 105.7529 116.1833 126.7762 # Value used for 137.5612undefined or 148.5649 159.8090
5298.544 5298.544 5298.544 5298.544missing data
# For Pseudo ArcInfo or Spiderweb wind:
n_cols = 5298.544 400 5298.544 5298.544 5298.544 5298.544 5298.544
|
Integrated standardisation
An integrated file format, suitable for all wind input, might look as follows:
Code Block |
---|
# START OF HEADER
# This file is created by Deltares
# Additional commments
# General part:
t_ref # Number of columns used for wind datafield
n_rows = 19991026.000000 = 16 # Reference time
t_ref_unit = ModifiedJulianDate # Number of rows used for wind datafield
grid_unit = degree # Unit used for reference time
t_unit = min # Unit of distances on the grid in both x- and y-direction
#grid_unit = m
# Only for #Pseudo Time unit (sec, min, hrs, etc.)
#t_unitArcInfo wind:
x0 = hrs-12.000
missing_value = -999.000 # Xcoordinate of starting point of grid (in units specified in grid_unit)
y0 = 48.000 # Value used for undefined or missing data
# For Pseudo ArcInfo or Spiderweb wind:
n_cols = 400 # Ycoordinate of starting point of grid (in units specified in grid_unit)
value_pos = centre # Number of columns used for wind datafield
n_rows = 16 # Location where the data is specified in a cell (centre or corner)
#value_pos = corner # Number of rows used for wind datafield
x_unit
dx = degree = 0.12500 # Unit of distances on the grid
#x_unit # Grid dimensions: dx and dy (in units =specified in grid_unit)
dy m
# Only for Pseudo ArcInfo= wind:
x0 0.083333333 = -12.000
n_quantity = 3 # Xcoordinate of starting point of grid (in units specified in x_unit)
y0 = 48.000 # Number of quantities prescribed in the file
quantity = x_wind y_wind # Ycoordinate of starting point ofair_pressure grid (in units specified in x_unit)
value_pos # Names of =the quantities
unit centre = m s-1 m s-1 Pa # Location where the data is specified in a# cellUnits (centreof orthe corner)
#value_pos =quantities
# Only for Spiderweb wind:
radius corner = 600000.00
cell_size = 0.12500 0.083333333 # Radius of spiderweb (in units specified in rad_unit)
rad_unit = # Gridm dimensions: dx and dy (in units specified in x_unit)
n_quantity = 3 # Unit of radius
n_quantity = 3 # Number of quantities prescribed in the file
quantity = x_wind y_wind air_pressure # NamesNumber of quantities prescribed in the quantitiesfile
unitquantity1 = wind_speed = m s-1 m s-1 Pa # Name of quantity1
quantity2 # Units of the quantities
#= Only for Spiderweb wind:
radius_from_direction = 600000.00 # Name of quantity2
quantity3 = p_drop # Radius of spiderweb (in units specified in rad_unit)
rad_unit = m # Name of quantity3
unit1 = m s-1 # Unit of radius
n_quantity = 3 # Unit of quantity1
unit2 = degree # Number of quantities prescribed in the file
quantity = wind_speed wind_from_direction p_drop # Unit of quantity2
unit3 # Names of the quantities
unit = Pa = m s-1 degree Pa # UnitsUnit of the quantitiesquantity3
# Only for wind on a separate curvi-linear grid:
grid_file = hirlam.grd # Separate (curvi-linear) grid on which the wind can be specified
first_row = nmax # These four parameters determine the order in which the data
last_row = 1 # is read from the separate wind file. In this way, data which is
first_col = 1 # supplied in a reversed order, can also be read.
last_col = mmax
n_quantity = 1 # Number of quantities prescribed in the file
quantity quantity1 = air_pressure # NamesName of the quantitiesquantity1
unit = Pa # UnitsUnit of the quantitiesquantity1
# END OF HEADER
TIME = 0.0000000E+00 # Time where first dataset is supplied
# Only for Spiderweb wind (every TIME input):
x0x_eye = 115.4 # Xcoordinate of cyclone eye (Spiderweb centre) at specified TIME
y0 y_eye = 18.9 # Ycoordinate of cyclone eye (Spiderweb centre) at specified TIME
p_drop_eye = 5300 # Pressure drop at cyclone eye (Spiderweb centre) at specified TIME
# Actual data field
1.388999 1.388999 1.388999 1.388999 1.388999 1.388999 1.388999 1.388999 1.388999.388999 1.388999
155.0000 346.2500 357.5000 8.750000 20.00000 131.38899925000 142.38899950000
15553.000075000 34665.250000000 35776.500025000
85300.750000000 205300.00000000 315300.25000000 425300.50000000 535300.75000000 655300.00000000 765300.25000000
5300.000 5300.000 5300.000
TIME 5300.000 = 5300.000 1440 5300.000 5300.000 5300.000 5300.000 5300.000
TIME = # Time where 1440next dataset is supplied
# Only for Spiderweb wind (every TIME input): # Time where next dataset is supplied
x0
x_eye = 114.4000 # Xcoordinate of cyclone eye (Spiderweb centre) at specified TIME
y0 y_eye = 18.90000 # Ycoordinate of cyclone eye (Spiderweb centre) at specified TIME
p_drop_eye = 5300.000 # Pressure drop at cyclone eye (Spiderweb centre) at specified TIME
# Actual data field
2.435253 2.435253 2.435253 2.435253 2.435253 2.435253 2.435253 2.435253 2.435253 2.435253
65.01978 75.11859 85.25320 95.45370 105.7529 116.1833 126.7762 137.5612 148.5649 159.8090
5298.544 5298.544 5298.544 5298.544 5298.544 5298.544 5298.544 5298.544 5298.544 5298.544
|