CDO

Hi all,

I am trying to remapbil SODA_Y1985_to_2015_DJF_climatol_salt.nc to my ROMS grid atleq3_1985to2015_salt_DJF_media_climatol.nc, to make the differences between them need to leave at the same resolution, but the ROMS grid is different from the SODA grid. As some users before, I encounter the problem of the unsupported generic grid.

ncdump -h atleq3_1985to2015_salt_DJF_media_climatol.nc (ROMS)
dimensions:
    xi_rho = 488 ;
    eta_rho = 292 ;
    s_rho = 30 ;
    ocean_time = UNLIMITED ; // (1 currently)
    bnds = 2 ;
variables:
    double s_rho(s_rho) ;
        s_rho:long_name = "S-coordinate at RHO-points" ;
        s_rho:positive = "up" ;
        s_rho:axis = "Z" ;
        s_rho:field = "s_rho, scalar" ;
    double ocean_time(ocean_time) ;
        ocean_time:standard_name = "time" ;
        ocean_time:long_name = "averaged time since initialization" ;
        ocean_time:bounds = "ocean_time_bnds" ;
        ocean_time:units = "seconds since 1900-1-1 00:00:00" ;
        ocean_time:calendar = "standard" ;
        ocean_time:axis = "T" ;
    double ocean_time_bnds(ocean_time, bnds) ;
    float salt(ocean_time, s_rho, eta_rho, xi_rho) ;
        salt:long_name = "time-averaged salinity" ;
        salt:_FillValue = 1.e+37f ;
        salt:missing_value = 1.e+37f ;
        salt:time = "ocean_time" ;
        salt:grid = "grid" ;
        salt:location = "face" ;
        salt:field = "salinity, scalar, series" ;

dimensions:
    longitude = 81 ;
    latitude = 48 ;
    depth = 50 ;
    time = UNLIMITED ; // (1 currently)
    bnds = 2 ;
variables:
    float longitude(longitude) ;
        longitude:standard_name = "longitude" ;
        longitude:long_name = "longitude" ;
        longitude:units = "degrees_east" ;
        longitude:axis = "X" ;
    float latitude(latitude) ;
        latitude:standard_name = "latitude" ;
        latitude:long_name = "latitude" ;
        latitude:units = "degrees_north" ;
        latitude:axis = "Y" ;
    float depth(depth) ;
        depth:long_name = "depth" ;
        depth:units = "m" ;
        depth:axis = "Z" ;
    double time(time) ;
        time:standard_name = "time" ;
        time:bounds = "time_bnds" ;
        time:units = "days since 1900-1-1 00:00:00" ;
        time:calendar = "proleptic_gregorian" ;
        time:axis = "T" ;
    double time_bnds(time, bnds) ;
    float salt(time, depth, latitude, longitude) ;
        salt:standard_name = "sea_water_salinity" ;
        salt:long_name = "Practical Salinity" ;
        salt:units = "psu" ;
        salt:_FillValue = NaNf ;
        salt:missing_value = NaNf ;
        salt:cell_methods = "time: mean" ;

I try cdo griddes file_ROMS.nc > grid.txt
cdo setgrid,grid.txt SODA_Y1985_to_2015_DJF_climatol_salt.nc outfile.nc

cdo -genbil,r488x292 -setgridtype,curvilinear atleq3_1985to2015_salt_DJF_media_climatol.nc atleq3_1985to2015_salt_DJF_media.nc

cdo remapbill,atleq3_1985to2015_salt_DJF_media_climatol.nc SODA_Y1985_to_2015_DJF_media_climatol_salt_naGRADEROMS.nc SODA_Y1985_to_2015_DJF_media_climatol_salt_INTERPOLATED.nc

Does anyone have a good idea for the SODA file to be the same as the ROMS file, so they have the same 488x292 array?

Best regards!