Interpolation algorithm for ERA5 data to ICON grid changed to bilinear #440
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.
Description
When using the Tegen aerosol climatology in ICON, a horizontally distributed reference pressure field is required for the vertical distribution of the aerosol optical depths (AOD). With the option
itype_vegetation_cycle >= 2, the reference pressure is computed with the ERA5 orography (resolution: ~30 km). For high-resolution ICON runs, the resolution of this orography is more similar to the very coarse one from the Tegen data.So far, the ERA5 orography was remapped to the ICON grid with a first-order conservative scheme. For an ICON simulation at 1km, this leads to a very "blocky" orography, and the artefacts are finally also visible in surface downward direct shortwave radiation (ASWDIR_S) through the computation of reference pressure and vertically resolved AODs:
I discussed this issue with Daniel Rieger and Günther Zängl (DWD). Günther suggested to change the interpolation algorithm for ERA5 orography. I implemented this by changing the interpolation algorithm for ERA5 orography (and other ERA5 variables) to bilinear interpolation. This resolves the block-artefacts in ASWDIR_S:
Workflow for merging PRs
Please read these instructions carefully and follow the steps below before requesting a review by the maintainers. This way we can ensure a smoother review process and your changes will be merged sooner.
Additionally, if this is the first PR you open in EXTPAR make sure to read the "Information for EXTPAR Developers" section in the documentation.
Checklist
docs/directory.If all the points above are satisfied you can ask for a review by selecting
stelliomas a reviewer.For any questions please ping
@stelliomon this PR.Testing and debugging
The most important test for PRs is the one labeled "EXTPAR Testsuite on Jenkins". This checks that the results of all testcases (described by the namelists in
test/testsuite/data) did not change compared to the references.You can launch the testsuite by writing
launch jenkinsas a comment in the PR that you want to test. Once completed, the result of the testsuite will be shown on the PR (failure or success).If you need more details on the testsuite results (e.g., if you are trying to debug an error or you are simply unsure why the tests fail) you can launch the testsuite with
launch jenkins(debug). This will run the tests as usual, but, once completed, you will be given a URL (via a comment on the PR) to access the logfiles and namelists of all tests that were run.