PyAPS icon indicating copy to clipboard operation
PyAPS copied to clipboard

Convert the ECMWF model level products to pressure level products

Open xhuang-jpl opened this issue 2 years ago • 3 comments

@scottstanie @mirzaees @vbrancat

To make the pyAPS accept the ECMWF model level products, this PR is to convert the model level ECMWF products in NetCDF format to the pressure level products. Both the theory and partial codes for the conversion can be found in ECMWF website. This PR is tested comparing with the ECMWF pressure level products. Some tests results are shown as,

1. Temperature

image

2. Specific Humidity

image

3. Geopotential

image

xhuang-jpl avatar Jan 30 '24 22:01 xhuang-jpl

@scottstanie @mirzaaes can you please review this PR? We would like to have a quick turn around on this and I believe this PR should be also beneficial for your projects. Thank you.

vbrancat avatar Feb 14 '24 15:02 vbrancat

I am a little concerned about how to deal with those NaNs that are out range of the pressure, interestingly, the downloaded pressure level products do not have this issue.

I am thinking two options, Option 1: Interpolating them, in this way, we might introduce more delays. Option 2: Skip the NaNs when calculate the delay, in which we might need more code change of the pyAPS to deal with NaNs.

xhuang-jpl avatar Feb 22 '24 01:02 xhuang-jpl

Since there is already GRIB test data in the repo, it should be quick to add a test of the main conversion function

scottstanie avatar Feb 22 '24 13:02 scottstanie