MetOffice / MetOffice/CDDS

bug longitude matrices in files for standard latitude/longitude domains when force_coordinate_rotation = True

Open
#937 0 comments 0 reactions 0 assignees View on GitHub
Dominant language
Python
Stars
8
Forks
7
Avg merge
1d 15h
Merged PRs (30d)
23

Description

When attempting to convert some RCM output that has been run over a non rotated domain over Africa with force_coordinate_rotation = True , the real (rather than rotated) longitude matrices have incorrect values (I think out be out by 180 degrees).

I thought that fixing the issues might be as simple as changing
560 rotated_cs = RotatedGeogCS(90., 0., ellipsoid=cs)
to
560 rotated_cs = RotatedGeogCS(90., 180., ellipsoid=cs)

in mip_convert//mip_convert/new_variable.py . Off line testing showed that the correct latitude longitude matrices get produced, however if I try to run get cdds_convert, it crashes with something about an error in cmor.write.

Note I do not have an urgent need for this to be fixed, just reporting in the hope that it can get picked up at some point.

Contributor guide

No contributing guide indexed for this repository

Research direction

Start in mip_convert/mip_convert/new_variable.py around line 560, where the RotatedGeogCS is created. Reproduce the conversion with non-rotated African RCM output and force_coordinate_rotation = True, then compare the real latitude and longitude matrices. Run cdds_convert and investigate the reported cmor.write failure; done means the matrices are correct and conversion completes successfully.

Written by the indexing model from the issue text.

Assessment

Tech stack
python
Domain
data
Issue type
Bug
Difficulty
3/5
Estimated time
1-2 days
Activity status
Quiet
Clarity
Mostly clear
Newbie friendliness
56/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.