ugrid: Define an urban grid

In sebschub/dcepucp: Generate an urban canopy parameters file for DCEP

Description

Function to define an urban grid including COSMO-CLM grid parameters. Some consistency checks are done.

Usage

ugrid(pollat, pollon, polgam = 0, dlat, dlon, startlat_tot, startlon_tot,
  ie_tot, je_tot, ke_uhl, hhl_uhl, n_uclass = 1, n_udir = 2,
  angle_udir = c(0, 90))

Arguments

pollat	Geographical latitude of the rotated north pole (in degrees, north >0); for a non-rotated lat-lon grid set pollat = 90.
pollon	Geographical longitude of the rotated north pole (in degrees, east >0); for a non-rotated lat-lon grid set pollon = -180.
polgam	Angle between the north poles of two rotated grids (in degrees, east > 0)
dlat	'Meridional' (rotated lat-direction) grid spacing (in degrees).
dlon	'Zonal' (rotated lon-direction) grid spacing (in degrees).
startlat_tot	Latitude of the lower left scalar grid point of the total domain (in degrees, north >0, rotated coordinates).
startlon_tot	Longitude of the lower left scalar grid point of the total domain (in degrees, east >0, rotated coordinates).
ie_tot	Number of gridpoints of the total domain in 'west-east' direction of the rotated coordinates.
je_tot	Number of gridpoints of the total domain in 'south-north' direction of the rotated coordinates.
ke_uhl	Number of main urban height levels corresponding to the number of wall elements.
hhl_uhl	Vector of the height of urban half levels. The ke_uhl+1 values are the heights of the roof levels.
n_uclass	Number of urban classes within a grid cell.
n_udir	Number considered street orientations for each urban class.
angle_udir	Vector of street orientations relative to the south-north axis.

Value

An object of class ugrid.

Note

In upar2nc, dlon and startlon_tot as well as dlat and startlat_tot are used to define the dimensions "rlon" and "rlat", respectively; ke_uhl and hhl_uhl for "uheight1"; n_udir and angle_udir for "udir"; n_uclass for "uclass".

Examples

grid <- ugrid(pollat = 40., pollon = -70.6,
              startlat_tot = -0.882, startlon_tot = -0.882,
              ie_tot = 195, je_tot = 195,
              dlat = 0.009, dlon = 0.009,
              ke_uhl = 8, hhl_uhl = c(0, 5, 10, 15, 20, 25, 30, 35,40),
              n_uclass = 1,
              n_udir = 2, angle_udir = c(0., 90.))

sebschub/dcepucp documentation built on May 20, 2019, 2:23 p.m.