Project

General

Profile

Regrid/interpolate CERRA dataset to regular lat/lon

Added by Alis A 16 days ago

Good morning, I have a netcdf file containing a CERRA dataset in Lambert Conformal Conic projection (details about the grid here: https://confluence.ecmwf.int/display/CKB/Copernicus+European+Regional+ReAnalysis+%28CERRA%29%3A+product+user+guide#heading-515TheCERRAgriddescription).
I need to transform it into a regular lat/lon grid, in order to be able to visualise it with an external tool (https://adaguc.knmi.nl/).
I am going in circles unfortunately. I read it should be possible with cdo, can you tell me how to?
The output of ncdump -h is

dimensions:
    values = 1142761 ;
    time = UNLIMITED ; // (11 currently)
variables:
    int64 time(time) ;
        time:units = "seconds since 2020-01-01 12:00:00" ;
        time:long_name = "time" ;
        time:calendar = "gregorian" ;
    float latitude(values) ;
        latitude:units = "degrees_north" ;
        latitude:long_name = "latitude" ;
    float longitude(values) ;
        longitude:units = "degrees_east" ;
        longitude:long_name = "longitude" ;
    float \10u(time, values) ;
    float \10v(time, values) ;
    float \2d(time, values) ;
    float \2t(time, values) ;
    float cos_julian_day(time, values) ;
    float cos_latitude(time, values) ;
    float cos_local_time(time, values) ;
    float cos_longitude(time, values) ;
    float insolation(time, values) ;
    float lsm(time, values) ;
    float msl(time, values) ;
    float q_100(time, values) ;
    float q_1000(time, values) ;
    float q_150(time, values) ;
    float q_200(time, values) ;
    float q_250(time, values) ;
    float q_300(time, values) ;
    float q_400(time, values) ;
    float q_50(time, values) ;
    float q_500(time, values) ;
    float q_700(time, values) ;
    float q_850(time, values) ;
    float q_925(time, values) ;
    float sin_julian_day(time, values) ;
    float sin_latitude(time, values) ;
    float sin_local_time(time, values) ;
    float sin_longitude(time, values) ;
    float skt(time, values) ;
    float sp(time, values) ;
    float t_100(time, values) ;
    float t_1000(time, values) ;
    float t_150(time, values) ;
    float t_200(time, values) ;
    float t_250(time, values) ;
    float t_300(time, values) ;
    float t_400(time, values) ;
    float t_50(time, values) ;
    float t_500(time, values) ;
    float t_700(time, values) ;
    float t_850(time, values) ;
    float t_925(time, values) ;
    float u_100(time, values) ;
    float u_1000(time, values) ;
    float u_150(time, values) ;
    float u_200(time, values) ;
    float u_250(time, values) ;
    float u_300(time, values) ;
    float u_400(time, values) ;
    float u_50(time, values) ;
    float u_500(time, values) ;
    float u_700(time, values) ;
    float u_850(time, values) ;
    float u_925(time, values) ;
    float v_100(time, values) ;
    float v_1000(time, values) ;
    float v_150(time, values) ;
    float v_200(time, values) ;
    float v_250(time, values) ;
    float v_300(time, values) ;
    float v_400(time, values) ;
    float v_50(time, values) ;
    float v_500(time, values) ;
    float v_700(time, values) ;
    float v_850(time, values) ;
    float v_925(time, values) ;
    float z_100(time, values) ;
    float z_1000(time, values) ;
    float z_150(time, values) ;
    float z_200(time, values) ;
    float z_250(time, values) ;
    float z_300(time, values) ;
    float z_400(time, values) ;
    float z_50(time, values) ;
    float z_500(time, values) ;
    float z_700(time, values) ;
    float z_850(time, values) ;
    float z_925(time, values) ;
    float z(time, values) ;
    float tp(time, values) ;
}

The file does not seem to contain metadata, which might be why some generic regridding I tried fails?
Commands such as
cdo -remapbil,global_0.1 myfile.nc myfile_regridded.nc
Give errors
cdo remapbil (Abort): Unsupported generic coordinates (Variable: latitude)!

I am new to this, so any help is greatly appreciated!