Retire ExtState%WLI (land-water-ice index) #170
Merged
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This PR is 1 of 2 for super-projects. See instructions here:
geoschem/GCHP#264
geoschem/GCClassic#26
This pull request retires use of external meteorology field land-water-ice (LWI) index and the ExtState array it is stored in (
ExtState%WLI
). We are removing it as a dependency because it will not be in the GEOS-IT dataset and determining whether a grid cell is over ocean, ice, or land is possible through other input meteorology fields, specifically GMAO meteorology constantsFRLAND
,FRLANDIC
,FRLAKE
, andFROCEAN
, and time-varying fieldsFRSEAICE
andFRSNO
.The
ExtState%WLI
field is currently used in HEMCO to retrieve land type via subroutineHCO_LANDTYPE
inhco_geotools_mod.F90
. That subroutine is called in the following extension modules, which will be impacted by this update:Note that the LWI input field used for WLI in HEMCO is currently computed during GEOS-Chem input pre-processing from the other GMAO surface type met-fields. Changing to using these fields at run-time will result in differences in the model. In addition, the new way for determining if mostly land, ice, or ocean is to evaluate which surface type fraction is the greatest. This is a slightly different algorithm than what was used to compute LWI and this will also cause small differences.
There is an accompanying pull request in GEOS-Chem that should be merged at the same time:
geoschem/geos-chem#1435.