MPAS-Dev / MPAS-Model

Repository for MPAS models and shared framework releases.
235 stars 312 forks source link

Enable supersampling of 30" terrestrial fields: ter, ivgtyp, isltyp, greenfrac #1120

Closed mgduda closed 10 months ago

mgduda commented 10 months ago

This PR enables the supersampling of the 30" terrain elevation, land use category, soil category, and green vegetation fraction datasets through the introduction of a new namelist option, config_30s_supersample_factor, in the &data_sources namelist group in the namelist.init_atmosphere file.

The config_30s_supersample_factor option defaults to 1 (i.e., to no supersampling), which appears to be reasonable for MPAS meshes with >= 3 km cells. Experimentally, a supersampling factor of 2 was found to be sufficient for a 1-km quasi-uniform regional mesh centered over Colorado.

The init_atm_map_static_data routine, which is used by the interp_terrain, interp_landuse, and interp_soilcat routines in the init_atmosphere core, already supported supersampling, and all that was required for these fields was to pass the namelist supersampling value down the subroutine call stack. The code to remap the MODIS FPAR 30-arc-second greenness fraction data did not support supersampling, and so was modified to use the same approach as for terrain, land use, and soil category.