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