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

Add (OSGB!) HadUK-Grid data to the CEDA-WPS #13

Open
5 of 6 tasks
agstephens opened this issue Jul 30, 2021 · 4 comments
Open
5 of 6 tasks

Add (OSGB!) HadUK-Grid data to the CEDA-WPS #13

agstephens opened this issue Jul 30, 2021 · 4 comments
Assignees

Comments

@agstephens
Copy link
Member

agstephens commented Jul 30, 2021

  • Review the HadUK-Grid datasets - which are the priorities?
    • 1km data (monthly) is the most popular so try that first.
  • Review the requirements for supporting OSGB in:
    • clisops / daops / roocs-utils - we can use lat/lon to specify the bounding box and keep the data in OSGB, thus avoiding regridding
    • phoenix / pywps / owslib - need to test rendering of the domain (UK) in Phoenix
  • Add a new process to flamingo - added for HadUK-Grid 1km mon
  • Test and roll out
@agstephens agstephens self-assigned this Jul 30, 2021
@agstephens
Copy link
Member Author

agstephens commented Jul 30, 2021

NOTE: it is possible to subset on the northings and eastings, as below. However, we are avoiding that route at this stage.

>>> import xarray as xr
>>> df = xr.open_dataset('/badc/ukmo-hadobs/data/insitu/MOHC/HadOBS/HadUK-Grid/v1.0.3.0/5km/tasmin/mon/v20210712/tasmin_hadukgrid_uk_5km_mon_195301-195312.nc')
>>> v = df.sel(projection_x_coordinate=slice(400000, 500000), projection_y_coordinate=slice(400000, 500000))
>>> v.tasmin.shape
(12, 20, 20)

@agstephens
Copy link
Member Author

Test 2: the files include auxiliary coordinate variables for latitude and longitude - so clisops might just work!!!

@agstephens
Copy link
Member Author

It looks like it works. Roll it out and get users to test and comment.

@agstephens
Copy link
Member Author

Why do 180 NC files get read from the FileMapper/DatasetMapper?

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

No branches or pull requests

1 participant