Learn R Programming

PuffR (version 0.1)

calmet_05_wind_field_opts_params: Set CALMET parameters for wind field options and other options

Description

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

Usage

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)

Arguments

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.