Update sfc_climo_gen to output fractional vegetation and soil type

[GeorgeGayno-NOAA]

Currently, this program outputs a dominate category at each model point. Update the program to output the percentage of each category within a model point.

[GeorgeGayno-NOAA]

Tried some timing and memory tests on Hera using 4adcb64. I used a small (301x200) CONUS ESG grid and the VIIRS vegetation data.

• Test 1 - use the vegetation_type.viirs.igbp.conus.30s.nc as input. Using one node and 24 tasks, it ran in about 60 seconds.
• Test 2 - use the vegetation_type.viirs.igbp.nh.30s.nc as input. I had to request the two 'bigmem' nodes, 12 tasks per node. It ran in about six minutes.

So, the fractional option can take a lot of memory. We are pushing the limit using the NH version of the file. Using the global version may not be feasible on some machines.

[GeorgeGayno-NOAA]

Tried some timing and memory tests on Hera using 4adcb64. I used a small (301x200) CONUS ESG grid and the VIIRS vegetation data.

• Test 1 - use the vegetation_type.viirs.igbp.conus.30s.nc as input. Using one node and 24 tasks, it ran in about 60 seconds.
• Test 2 - use the vegetation_type.viirs.igbp.nh.30s.nc as input. I had to request the two 'bigmem' nodes, 12 tasks per node. It ran in about six minutes.

So, the fractional option can take a lot of memory. We are pushing the limit using the NH version of the file. Using the global version may not be feasible on some machines.

I repeated this test on WCOSS2 using the NH file. It was able to run on one node/24 tasks, but I had to request 500 GB of memory: #PBS -l select=1:ncpus=24:mem=500GB. It took about 6 minutes to run.

[GeorgeGayno-NOAA]

When processing a fractional grid (points can be a mix of land and non-land), how will the dominate category be defined? Suppose a point is:

• 40% land
• 60% non-land
• 35% forest
• 5% grassland

Is the dominate category forest or water?

[GeorgeGayno-NOAA]

In order to read the water flag from the file, the file attributes and missing values need to be changed. This may be done using netCDF tools. On WCOSS2, the vegetation type files were changed as follows:

module load gsl
module load nco

ncatted -h -a water_category,vegetation_type,c,i,17 file.nc
ncatted -h -a missing_value,vegetation_type,o,i,17 file.nc
ncap2 -O -s 'where(vegetation_type < 0) vegetation_type=17' file.nc file.nc