calmet_05_wind_field_opts_params: Set CALMET parameters for wind field options and other...
In rich-iannone/PuffR: A set of tools to help with CALPUFF dispersion modelling.

Description Usage Arguments

This function validates and writes CALMET parameters for wind field options and parameters to the working CALMET.INP file.

calmet_05_wind_field_opts_params(calmet_inp = "calmet_template.txt",
  iwfcod = 1, ifradj = 1, ikine = 0, iobr = 0, islope = 1,
  iextrp = -4, icalm = 0, bias = c(-1, -1, -0.5, -0.3, 0, 1, 1, 1, 1, 1,
  1, 1), rmin2 = -1, iprog = 0, isteppg = 1, igfmet = 0, lvary = TRUE,
  rmax1 = 5, rmax2 = 10, rmax3 = 20, rmin = 0.1, terrad = 5, r1 = 2,
  r2 = 1, rprog = 0, divlim = 5e-06, niter = 50, nsmth = c(2, 4, 4, 4,
  4, 4, 4, 4, 4, 4, 4, 4), nintr2 = 99, critfn = 1, alpha = 0.1,
  fextr2 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), nbar = 0, kbar = 10,
  xbbar = 0, ybbar = 0, xebar = 0, yebar = 0, idiopt1 = 0,
  isurft = 1, idiopt2 = 0, iupt = 1, zupt = 200, idiopt3 = 0,
  iupwnd = -1, zupwnd = c(1, 1000), idiopt4 = 0, idiopt5 = 0,
  llbreze = FALSE, nbox = 0, xg1 = 0, xg2 = 0, yg1 = 0, yg2 = 0,
  xbcst = 0, ybcst = 0, xecst = 0, yecst = 0, nlb = 0, metbxid = 0)

`calmet_inp`	the absolute path and filename for the working CALMET input file.
`iwfcod`	for model, choose betwen: (0) objective analysis only, or (1) diagnostic wind module.
`ifradj`	compute Froude number adjustment effects.
`ikine`	compute kinematic effects.
`iobr`	use O'Brien procedure for adjustment of the vertical velocity.
`islope`	compute slope flow effects.
`iextrp`	extrapolate surface wind observations to upper layers.
`icalm`	extrapolate surface winds even if calm?
`bias`	layer-dependent biases modifying the weights of surface and upper air stations
`rmin2`	minimum distance from nearest upper air station to surface station for which extrapolation of surface winds at surface station will be allowed.
`iprog`	use gridded prognostic wind field model output fields as input to the diagnostic wind field model.
`isteppg`	timestep in hours of the prognostic model input data.
`igfmet`	use coarse CALMET fields as initial guess fields.
`lvary`	use varying radius of influence.
`rmax1`	maximum radius of influence over land in the surface layer.
`rmax2`	maximum radius of influence over land aloft.
`rmax3`	maximum radius of influence over water.
`rmin`	minimum radius of influence used in the wind field interpolation.
`terrad`	radius of influence of terrain features.
`r1`	the relative weighting of the first guess field and observations in the surface layer.
`r2`	the relative weighting of the first guess field and observations in layers aloft.
`rprog`	relative weighting parameter of the prognostic wind field data.
`divlim`	maximum acceptable divergence in the divergence minimization procedure.
`niter`	maximum number of iterations in the divergence min. procedure.
`nsmth`	number of passes in the smoothing procedure.
`nintr2`	maximum number of stations used in each layer for the interpolation of data to a grid point.
`critfn`	the critical Froude number.
`alpha`	an empirical factor controlling the influence of kinematic effects.
`fextr2`	multiplicative scaling factors for extrapolation of surface observations to upper layers.
`nbar`	number of barriers to interpolation of the wind fields.
`kbar`	level up to which barriers apply (must be value in the range of 1 to NZ).
`xbbar`	the x coordinates for the beginning of each barrier.
`ybbar`	the y coordinates for the beginning of each barrier.
`xebar`	the x coordinates for the ending of each barrier.
`yebar`	the y coordinates for the ending of each barrier.
`idiopt1`	method for computation of surface temperatures: (0) compute internally from hourly surface observations, (1) read preprocessed values from a data file (DIAG.DAT).
`isurft`	surface meteorological station to use for the surface temperature (must be value in the range of 1 to nssta).
`idiopt2`	method for computation of domain-averaged temperature lapse rate: (0) compute internally from twice-daily upper air observations, (1) read hourly preprocessed values from a data file (DIAG.DAT).
`iupt`	upper air station to use for the domain-scale lapse rate (must be a value in the range of 1 to NUSTA).
`zupt`	Depth in meters through which the domain-scale lapse rate is computed.
`idiopt3`	method for computation of domain-averaged wind components: (0) compute internally from twice-daily upper air observations, (1) read hourly preprocessed values a data file (DIAG.DAT).
`iupwnd`	upper air station to use for the domain-scale winds (must be a value in the range of -1 to NUSTA).
`zupwnd`	bottom and top of layer through which the domain-scale winds are computed.
`idiopt4`	selection of observed surface wind components for wind field module: (0) read WS and WD values from a surface data file (SURF.DAT), (1) read hourly preprocessed U and V values from a data file (DIAG.DAT).
`idiopt5`	selection of observed upper air wind components for wind field module: (0) read WS and WD values from an upper air data file (UP.DAT), (1) read hourly preprocessed U and V values from a data file (DIAG.DAT).
`llbreze`	use lake breeze module?
`nbox`	number of lake breeze regions.
`xg1`	the x direction grid line 1 defining the region of interest for the lake breeze module.
`xg2`	the x direction grid line 2 defining the region of interest for the lake breeze module.
`yg1`	the y direction grid line 1 defining the region of interest for the lake breeze module.
`yg2`	the y direction grid line 2 defining the region of interest for the lake breeze module.
`xbcst`	the beginning x point (in kilometers) defining the coastline (straight line).
`ybcst`	the beginning y point (in kilometers) defining the coastline (straight line).
`xecst`	the ending x point (in kilometers) defining the coastline (straight line).
`yecst`	the ending y point (in kilometers) defining the coastline (straight line).
`nlb`	the combined number of meteorological and upper air stations in the region.
`metbxid`	station identifiers for the region; include surface stations first, then upper air stations.