search

ecmwf / ecmwf-opendata

A package to download ECMWF open data
Apache License 2.0
160 stars 27 forks source link

Reproject aifs to a regular lat/lon grid #36

Open sailgrib opened 4 months ago

sailgrib commented 4 months ago

I am trying to reproject the reduced-gaussian aifs grib file to a regulat 0.25° lat/lon grid. I have tried 3 approaches, one using nearest neighbor and 2 others that should be better but both failed.

The first one using nearest neighbor : cdo remapnn,target_grid_world_025 input_file output_file That works but is it really appropriate?

The second one with bilenear interpolation: cdo remapbil,target_grid_world_025.txt input_file output_file => cdo remapbil: Bilinear weights from gaussian_reduced (542080) to lonlat (1440x721) grid => cdo remapbil (Abort): Can't do bilinear interpolation if the source grid is not a regular 2D grid!

The third one using first-order conservative remapping: cdo remapycon,target_grid_world_025.txt input_file output_file => cdo remapcon (Abort): Source grid cell corner coordinates missing!

Any idea as to how I should correctly reproject the aifs to a regular lat/lon grid?

igorROIK commented 4 months ago

I've tried several ways to plot ECMWF AIFS data in Python, but without success. Could you help me or provide a script using this data? Thanks.

sailgrib commented 3 months ago

Sorry, I cannot help with this. I have not tried to plot ECMWF AIFS. i am using the data as lat/lon grib files.