Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Bugfix: mean land elevation + v13 bcs #1049

Open
wants to merge 67 commits into
base: develop
Choose a base branch
from

Conversation

biljanaorescanin
Copy link
Contributor

@biljanaorescanin biljanaorescanin commented Jan 31, 2025

This PR does 2 things:

1. Update Target Mean Elevation for Land Tiles

make_bcs records the elevation of each land tile in the catchment.def file, based on the topography data underlying the Pfafstetter catchments. To ensure that the mean elevation of the land tiles matches that computed from the topography used for the atmosphere (NCAR GMTED TOPO 30arcsec), a scaling factor is applied to the land tile elevations before recording them in the catchment.def file.

Previously, make_bcs used a hardcoded value of 614.649 m for the target mean land tile elevation computed from NCAR GMTED TOPO 30arcsec. This value is the straight average across the tiles, whereas the mean adjustment uses the tile-area weighted mean elevation, which resulted in an error in the mean elevation scaling in all existing bcs versions.

This PR corrects the hardcoded value for the target mean land elevation to 656.83 m, which reflects the cosine-lat-weighted mean elevation from the NCAR GMTED TOPO 30arcsec dataset.

Impact: Ensures that land elevation scaling is based on the correct target mean elevation. Results in differences in the land tile elevations in catchment.def when reproducing existing bcs versions, hence the PR is non-zero-diff for make_bcs. Since the land tile elevations from catchment.def are not used by the GCM or the LDAS, this PR is trivially zero-diff for the GCM and the LDAS.

2. Release of v13 BCS set for GCM v12

Defines new v13 bcs dataset for use with GCM v12. The v13 bcs data are identical to v12 data except for the aforementioned land tile elevations and the associated topography data for the atmosphere. The latter are added to the bcs manually at this time. Work is in progress to include the processing of the topography for the atmosphere into make_bcs. Adding a new set of bcs only impacts make_bcs and is trivially zero-diff for the GCM and the LDAS.

Contingencies:
Before merging this PR, both #909 and #1028 have to be merged.
Base of PR should eventually be develop. (Currently, it is feature/wjiang/nc4_tilefile of #1028.)

@wmputman @rdkoster @weiyuan-jiang

weiyuan-jiang and others added 30 commits February 23, 2024 10:51
…Sgcm_GridComp into feature/wjiang/nc4_tilefile
…ram.F90, mkEASETilesParam.F90, rmTinyCatchParaMod.F90)
…butes() [formerly catchment_def()] (rmTinyCatchParaMod.F90)
@biljanaorescanin biljanaorescanin added Not 0-diff for BCS 0 diff The changes in this pull request have verified to be zero-diff with the target branch. labels Jan 31, 2025
@gmao-rreichle gmao-rreichle changed the base branch from feature/wjiang/nc4_tilefile to develop January 31, 2025 17:09
@gmao-rreichle gmao-rreichle changed the base branch from develop to feature/wjiang/nc4_tilefile January 31, 2025 17:09
@biljanaorescanin biljanaorescanin changed the base branch from feature/wjiang/nc4_tilefile to develop February 3, 2025 18:06
@biljanaorescanin biljanaorescanin marked this pull request as ready for review February 3, 2025 18:06
@biljanaorescanin biljanaorescanin requested review from a team as code owners February 3, 2025 18:06
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
0 diff The changes in this pull request have verified to be zero-diff with the target branch. Not 0-diff for BCS
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants