Explicit urban vars

[Peter9192]

Add spatially explicit radiation parameters to SLUCM

TYPE: new feature

KEYWORDS: slucm, 2d fields, radiation, albedo, emissivity

SOURCE: Peter Kalverla (Netherlands eScience Center)

DESCRIPTION OF CHANGES:
Problem:
Currently, urban albedo and emissivity can only be provided per urban class. We are investigating methods to obtain spatially explicit radiation fields and want to assess its impact in WRF.

Solution:
We added 6 parameters to the model code that are propagated into the land surface schemes and eventually the urban scheme. See attached for details: urban_flow.pdf

LIST OF MODIFIED FILES: list of changed files (use git diff --name-status master to get formatted list)

M Registry/Registry.EM_COMMON
M dyn_em/module_first_rk_step_part1.F
M phys/module_sf_clm.F
M phys/module_sf_noahdrv.F
M phys/module_sf_urban.F
M phys/module_surface_driver.F
M phys/noahmp (submodule)

TESTS CONDUCTED:

1. Do mods fix problem? How can that be demonstrated, and was that test conducted?
2. Are the Jenkins tests all passing?

RELEASE NOTE: Include a stand-alone message suitable for the inclusion in the minor and annual releases. A publication citation is appropriate.

[Peter9192]

To do:

• Don't abuse merge as ternary? --> now just passing values directly, without presence check
• Check -999. treated correctly and no more reliance on presence (since it's always present now)
• Check whether 2D input fields ingested directly yield same results as when derived from LCZ param table in urban scheme (use nearest neighbor)?
• Document namelist options
• Document WPS procedure
• Make PR to noah subrepo ?

[Peter9192]

Currently, if new vars present in met_em, will use that (with -999. as fillvalue). If not present in met_em, real.exe generates new fields anyway, but fills them with 0.

[Peter9192]

Added two namelist options in the physics section:

urb2d_additional_vars      = 'ALBB','ALBR','ALBG',   ! specify one or more from 'ALBB','ALBR','ALBG','EPSB','EPSR','EPSG'
num_urb2d_additional_vars  = 3                       ! specify the number of variable provided above (e.g. if set to 2, ALBG is skipped.)

Alternatively you could do:

urb2d_additional_vars     = 'ALBB','ALBR','ALBG','EPSB','EPSR','EPSG'
num_urb2d_additional_vars = 1  ! only ALBB will be considered