Project

General

Profile

mergetime error

Added by Felix Emlinger over 3 years ago

Hello,

I do have 360 files with climate projection data, containing precipitation values and a 0.11° grid over whole Europe. I want to cut out a specific grid cell and merge all files together (in total around 25 Gb). When I am trying to use mergetime on a small amount of 5 files, I always get the error:

cdo mergetime (Abort): Grid size of the input parameter pr do not match!
terminate called without an active exception
Aborted

Here my command:

cdo -z zip -mergetime -remapnn,lon=8.730004_lat=51.491288 pr_EUR-11_CCCma-CanESM2_rcp85_r1-r1i1p1_OURANOS-CRCM5_kba_1h_20*.nc Merge.nc

This does not make any sense in my point of view, because when I do check the sinfo all the files have the exact same coordinate range.

Output of: cdo sinfov pr_EUR-11_CCCma-CanESM2_rcp85_r1-r1i1p1_OURANOS-CRCM5_kba_1h_20750*.nc

felixemlinger@DESKTOP-0S2QC8F:~/Test2$ cdo sinfov pr_EUR-11_CCCma-CanESM2_rcp85_r1-r1i1p1_OURANOS-CRCM5_kba_1h_20750*.nc
File format : NetCDF4 classic zip
-1 : Institut Source T Steptype Levels Num Points Num Dtype : Parameter name
1 : unknown unknown v instant 1 1 78400 1 F32z : pr
Grid coordinates :
1 : curvilinear : points=78400 (280x280)
lon : -24.98154 to 37.91363 degrees_east
lat : 30.42796 to 64.665 degrees_north
mapping : rotated_latitude_longitude
rlon : -21.66499 to 9.025009 by 0.11 degrees
rlat : -16.775 to 13.915 by 0.11 degrees
Vertical coordinates :
1 : surface : levels=1
Time coordinate : 744 steps
RefTime = 2075-01-01 00:00:00 Units = days Calendar = 365_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
2075-01-01 00:30:00 2075-01-01 01:30:00 2075-01-01 02:30:00 2075-01-01 03:30:00

File format : NetCDF4 classic zip
-2 : Institut Source T Steptype Levels Num Points Num Dtype : Parameter name
1 : unknown unknown v instant 1 1 78400 1 F32z : pr
Grid coordinates :
1 : curvilinear : points=78400 (280x280)
lon : -24.98154 to 37.91363 degrees_east
lat : 30.42796 to 64.665 degrees_north
mapping : rotated_latitude_longitude
rlon : -21.66499 to 9.025009 by 0.11 degrees
rlat : -16.775 to 13.915 by 0.11 degrees
Vertical coordinates :
1 : surface : levels=1
Time coordinate : 672 steps
RefTime = 2075-01-01 00:00:00 Units = days Calendar = 365_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
2075-02-01 00:30:00 2075-02-01 01:30:00 2075-02-01 02:30:00 2075-02-01 03:30:00
File format : NetCDF4 classic zip
-4 : Institut Source T Steptype Levels Num Points Num Dtype : Parameter name
1 : unknown unknown v instant 1 1 78400 1 F32z : pr
Grid coordinates :
1 : curvilinear : points=78400 (280x280)
lon : -24.98154 to 37.91363 degrees_east
lat : 30.42796 to 64.665 degrees_north
mapping : rotated_latitude_longitude
rlon : -21.66499 to 9.025009 by 0.11 degrees
rlat : -16.775 to 13.915 by 0.11 degrees
Vertical coordinates :
1 : surface : levels=1
Time coordinate : 720 steps
RefTime = 2075-01-01 00:00:00 Units = days Calendar = 365_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
2075-04-01 00:30:00 2075-04-01 01:30:00 2075-04-01 02:30:00 2075-04-01 03:30:00
2075-04-30 20:30:00 2075-04-30 21:30:00 2075-04-30 22:30:00 2075-04-30 23:30:00
cdo sinfon: Processed 4 variables over 2880 timesteps [0.20s 57MB].
I did remove the timesteps.

Another question I do have is if it is normal, that using the remapnn command on a single 50 MB files takes 250 seconds.

I would be really happy if I could get some help with my problem.

Greetings:
Felix


Replies (2)

RE: mergetime error - Added by Ralf Mueller over 3 years ago

in your command

cdo -z zip -mergetime -remapnn,lon=8.730004_lat=51.491288 pr_EUR-11_CCCma-CanESM2_rcp85_r1-r1i1p1_OURANOS-CRCM5_kba_1h_20*.nc Merge.nc
you use a shell wildcard which possibly expands to multiple input files for the remapnn operator. BUT this operators does not work with multiple inputs, but only with a single one. That's why you call does not work as expected.

Try

cdo -z zip -mergetime -apply,-remapnn,lon=8.730004_lat=51.491288 [ pr_EUR-11_CCCma-CanESM2_rcp85_r1-r1i1p1_OURANOS-CRCM5_kba_1h_20*.nc ] Merge.nc 
with the current 1.9.9 release of CDO instead. 1.9.8 might also work.

hth
ralf

RE: mergetime error - Added by Felix Emlinger over 3 years ago

Worked perfectly, thanks alot :))

    (1-2/2)