Adding diagnostic output variables for HRLDAS-urban simulation

[cenlinhe]

Currently, when users activate the urban physics (>0), there is no urban-related output variables in the default HRLDAS. For example, the T2MB and T2MV are the 2-m temperature for the rural portion of the urban grid instead of the urban 2-m temperature in HRLDAS.

There is a need to add new diagnostic output variables for urban-specific variables and the grid-mean diagnostic variables including rural and urban portions of the grid following the WRF calculations: https://github.com/wrf-model/WRF/blob/release-v4.5.1/phys/module_surface_driver.F#L3369

These new diagnostic output variable calculations should be added in hrldas driver somewhere here: https://github.com/NCAR/hrldas/blob/master/hrldas/IO_code/module_NoahMP_hrldas_driver.F#L1072