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

Dataset problem: CORDEX datasets in a Lambert Conformal coordinate system #1879

Open
sloosvel opened this issue Dec 20, 2022 · 2 comments · May be fixed by #1956
Open

Dataset problem: CORDEX datasets in a Lambert Conformal coordinate system #1879

sloosvel opened this issue Dec 20, 2022 · 2 comments · May be fixed by #1956
Assignees

Comments

@sloosvel
Copy link
Contributor

Describe the dataset issue
PR #1765 fixed datasets that use rotated pole coordinates, but not Lambert conformal ones. I see two possibilities to fix the data:

1- Regrid the lambert conformal datasets to the rotated pole coordinates

2- Take the standardised rlat and rlot coordinates, convert their coordinate system to LambertConformal and interpolate the points to match the shape of the data. For the geographical coordinates, just interpolate to match the shape of the data.

I am not sure if either approach is very precise, so more ideas are welcome.

FYI @pepcos

Assign the right people
If you are already a member of the ESMValTool GitHub project, please assign Valeriu Predoi (valeriupredoi) and
Klaus Zimmermann (zklaus) to the issue. They will then check the issue raised and propagate it further to the
data model developers.

Please attach

  • The recipe that you are trying to run, you can find a copy in the run directory in the output directory
  • The main_log_debug.txt file, this can also be found in the run directory in the output directory
  • Attach a text file with the output of ncdump -h /path/to/file.nc

Cheers 🍺

@nhsavage
Copy link

see also #1156 and iris ticket #4812 (I am not clear if the only issue is the lack of area weights or if there are other problems)

@sloosvel
Copy link
Contributor Author

After discussing with @zklaus and @Neah-Ko, and considering the cordex domain descriptions, we will store in a json file the values of the domain boundaries as given in the specifications:

image

We can start by checking how does the boundary domain compare to domain in each dataset, by building a mesh grid and plotting the grid points. It is likely that the model domains will expand beyond the standardised domains, as this is done to avoid issues in the boundaries when running the model

With that information, then fixing the grid would be a matter of transforming the lambert and geographical coordinate systems back and forth, whether that is to fix coordinate points or generate the bounds. If relevant, and depending on how the mentioned plots look like, we could consider cropping the cube to fit the standardised domain.

@Neah-Ko Neah-Ko linked a pull request Mar 3, 2023 that will close this issue
9 tasks
@Neah-Ko Neah-Ko linked a pull request Mar 3, 2023 that will close this issue
9 tasks
@zklaus zklaus removed their assignment Feb 7, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging a pull request may close this issue.

4 participants