Project

General

Profile

Refer CORDEX data to non-rotated coordinates

Added by Cecilia Campero almost 3 years ago

I have a CORDEX CAM-44 data in rotated coordinates, but I need the file in non-rotated lat/long. By using cdo sinfo pr_rcp26_2021_2025.nc, I can see both, the rotated and the non-rotated coordinate variable values.

File format : NetCDF4 classic zip
-1 : Institut Source T Steptype Levels Num Points Num Dtype : Parameter ID
1 : unknown unknown v instant 1 1 29190 1 F32z : -1
Grid coordinates :
1 : curvilinear : points=29190 (210x139)
lon : -127.5838 to -19.40505 degrees_east
lat : -19.46236 to 46.38 degrees_north
mapping : rotated_latitude_longitude
rlon : -52.8 to 39.16 by 0.44 degrees
rlat : -28.6 to 32.12 by 0.44 degrees
Vertical coordinates :
1 : surface : levels=1
Time coordinate : 1800 steps
RefTime = 1949-12-01 00:00:00 Units = days Calendar = 360_day Bounds = true
YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss
2021-01-01 12:00:00 2021-01-02 12:00:00 2021-01-03 12:00:00 2021-01-04 12:00:00
2021-01-05 12:00:00 2021-01-06 12:00:00 2021-01-07 12:00:00 2021-01-08 12:00:00
2021-01-09 12:00:00 2021-01-10 12:00:00 2021-01-11 12:00:00 2021-01-12 12:00:00

But, with ncdump -h pr_rcp26_2021_2025.nc, the precipitation variable is refered to the rotated coordinates (rlat,rlon). My question is, how can I refer the precipitation data to the non-rotated coordinates? Thanks in advance

dimensions:
        rlat = 139 ;
        rlon = 210 ;
        time = UNLIMITED ; // (1800 currently)
        bnds = 2 ;
variables:
        double lat(rlat, rlon) ;
                lat:standard_name = "latitude" ;
                lat:long_name = "latitude" ;
                lat:units = "degrees_north" ;
        double lon(rlat, rlon) ;
                lon:standard_name = "longitude" ;
                lon:long_name = "longitude" ;
                lon:units = "degrees_east" ;
        float pr(time, rlat, rlon) ;
                pr:grid_mapping = "rotated_pole" ;
                pr:_FillValue = 1.e+20f ;
                pr:missing_value = 1.e+20f ;
                pr:standard_name = "precipitation_flux" ;
                pr:long_name = "Precipitation" ;
                pr:units = "kg m-2 s-1" ;
                pr:coordinates = "lon lat" ;
                pr:cell_methods = "time: mean" ;
        double rlat(rlat) ;
                rlat:standard_name = "grid_latitude" ;
                rlat:long_name = "latitude in rotated pole grid" ;
                rlat:units = "degrees" ;
                rlat:axis = "Y" ;
        double rlon(rlon) ;
                rlon:standard_name = "grid_longitude" ;
                rlon:long_name = "longitude in rotated pole grid" ;
                rlon:units = "degrees" ;
                rlon:axis = "X" ;
        char rotated_pole ;
                rotated_pole:grid_mapping_name = "rotated_latitude_longitude" ;
                rotated_pole:grid_north_pole_latitude = 75.74 ;
                rotated_pole:grid_north_pole_longitude = -246.02 ;
        double time(time) ;
                time:standard_name = "time" ;
                time:units = "days since 1949-12-01 00:00:00" ;
                time:calendar = "360_day" ;
                time:long_name = "time" ;
                time:bounds = "time_bnds" ;
                time:axis = "T" ;
        double time_bnds(time, bnds) ;

// global attributes:
                :Conventions = "CF-1.4" ;
                :contact = "rossby.cordex@smhi.se" ;
                :creation_date = "2016-06-09-T18:15:09Z" ;
                :experiment = "RCP2.6" ;
                :experiment_id = "rcp26" ;
                :driving_experiment = "MOHC-HadGEM2-ES, rcp26, r1i1p1" ;
                :driving_model_id = "MOHC-HadGEM2-ES" ;
                :driving_model_ensemble_member = "r1i1p1" ;
                :driving_experiment_name = "rcp26" ;
                :frequency = "day" ;
                :institution = "Swedish Meteorological and Hydrological Institute, Rossby Centre" ;
                :institute_id = "SMHI" ;
                :model_id = "SMHI-RCA4" ;
                :rcm_version_id = "v1" ;
                :project_id = "CORDEX" ;
                :CORDEX_domain = "CAM-44" ;
                :product = "output" ;
                :references = "http://www.smhi.se/en/Research/Research-departments/climate-research-rossby-centre" ;
                :tracking_id = "3488aee8-871d-42d4-b2e2-316cae45e58f" ;
                :rossby_comment = "201601: CORDEX Central America 0.44 deg | RCA4 v1 | MOHC-HadGEM2-ES | r1i1p1 | rcp26 | L40" ;
                :rossby_run_id = "201601" ;
                :rossby_grib_path = "/home/rossby/prod/201601/raw/" ;

Replies (1)

RE: Refer CORDEX data to non-rotated coordinates - Added by Ralf Mueller almost 3 years ago

hi Cecilia!

I think CDO correctly represents to coordinate. it just displays both systems with sinfo. If you want to put the precip onto another grid you should use normal remapping I think.

best wishes
ralf

    (1-1/1)