CDO

Hello,

I am having some trouble changing projections from Lambert Conformal Conic to WGS 84 lat/lon. I have tried gdalwarp which works, but it converts my NetCDF into multiple bands, as if a geoTIFF (but with .nc extension). I have read a better option is using CDO combined with sellonlatbox? However, I am having trouble determining the right information/format either for creating a template.nc or grid.txt. Please advise appropriate detailed steps. Also, further guidance is appreciated should I want to apply procedure to other NetCDF files (i.e. EURO-Cordex).

NetCDF file size is too large to be uploaded here, but available at https://drive.google.com/file/d/1K0DaNHT_ovGbiGtA8FmlAbo1F5F-Xo3P/view?usp=sharing.

Some details:

dimensions:
    time = UNLIMITED;   // (1827 currently)
    x = 373;
    y = 513;
    bnds = 2;
  variables:
    double height;
      :axis = "Z";
      :long_name = "height";
      :positive = "up";
      :standard_name = "height";
      :units = "m";

int Lambert_Conformal;
      :grid_mapping_name = "lambert_conformal_conic";
      :standard_parallel = 24.0; // double
      :longitude_of_central_meridian = 45.99999999999999; // double
      :latitude_of_projection_origin = 24.0; // double
      :false_easting = 1939997.2595554162; // double
      :false_northing = 2639996.270735208; // double
      :earth_radius = 6371220.0; // double
      :proj4 = "+proj=lcca +lat_1=24.000000 +lat_0=24.000000 +lon_0=46.000000 +k_0=1.0 +x_0=1939997.259555 +y_0=2639996.270735 +a=6371220.000000 +b=6371220.000000";

double lat(y=513, x=373);
      :long_name = "Latitude";
      :standard_name = "latitude";
      :units = "degrees_north";
      :_ChunkSizes = 513U, 373U; // uint

double lon(y=513, x=373);
      :long_name = "Longitude";
      :standard_name = "longitude";
      :units = "degrees_east";
      :_ChunkSizes = 513U, 373U; // uint

float tas(time=1827, y=513, x=373);
      :standard_name = "air_temperature";
      :long_name = "Near-Surface Air Temperature";
      :units = "K";
      :grid_mapping = "Lambert_Conformal";
      :coordinates = "height lat lon";
      :_FillValue = 1.0E20f; // float
      :missing_value = 1.0E20f; // float
      :cell_methods = "time: mean";
      :_ChunkSizes = 1U, 513U, 373U; // uint

double time(time=1827);
      :standard_name = "time";
      :long_name = "time";
      :bounds = "time_bnds";
      :units = "days since 1949-12-1 00:00:00";
      :calendar = "standard";
      :axis = "T";
      :_ChunkSizes = 524288U; // uint

double time_bnds(time=1827, bnds=2);
      :_ChunkSizes = 1U, 2U; // uint

double x(x=373);
      :axis = "X";
      :long_name = "X Coordinate Of Projection";
      :standard_name = "projection_x_coordinate";
      :units = "m";
      :_ChunkSizes = 373U; // uint

double y(y=513);
      :axis = "Y";
      :long_name = "Y Coordinate Of Projection";
      :standard_name = "projection_y_coordinate";
      :units = "m";
      :_ChunkSizes = 513U; // uint

Thank you kindly.