consects
Added by Grace Redmond over 11 years ago
Hi,
I am trying to calculate the number of consecutive dry days in a 30 year netCDF file. I tried the eca_cdd command but I kept getting 'segmentation fault'
I am now trying to use consects as follows:
cdo consects -gtc,1 30_year_precip.nc out.nc
The output file appears to have a field for every day, what exactly do these fields mean? Also if I try to view the file I get:
ncvarget: ncid 5: Invalid argument
I want to find the maximum number of consecutive dry days at each grid point (I'm assuming I would have to use cdo timmax -consects -gtc,1
Thanks!
Replies (9)
RE: consects - Added by Ralf Mueller over 11 years ago
Hi Grace!
There is a ticket for the development of the consects/consecsum operators: #14. I hope it is some how illustrative. Could you upload the input file or at least the lets say the first year?
regards
ralf
RE: consects - Added by Jaison-Thomas Ambadan over 11 years ago
I've tested the "eca_cdd" (CDO v1.6.x) on "GPCP Version 2.2 Combined Precipitation Dataset (Final)" NetCDF file from NOAA. It works fine, so I guess the 'segmentation fault' is file specific - may be you could upload the input file as Ralf suggested
RE: consects - Added by Grace Redmond over 11 years ago
Hi, I have uploaded the first year of the file
RE: consects - Added by Grace Redmond over 11 years ago
Sorry I though I attached it, I think it was too big. Here is a 3 month file
| aphrodite_1971_jja.nc (35.4 MB) aphrodite_1971_jja.nc | 3 months of daily precip |
RE: consects - Added by Jaison-Thomas Ambadan over 11 years ago
I've tested both "eca_cdd" and the other method you mentioned. It works fine and the output (attached) look same in ncview. "ncdump -h" will show the details of the commad sequence and the cdo version that I used. Could you please upload the output of "cdo -V"
| eca_cdd_file.nc (794 KB) eca_cdd_file.nc | |||
| other_way_file.nc (400 KB) other_way_file.nc |
RE: consects - Added by Grace Redmond over 11 years ago
Thanks, I can view the files in ncview ok, but if I try to xconv -i them I get a segmentation fault with the eca_cdd file, and not with the other_way. This isn't necessarily a problem, I just wanted to make sure eca_cdd is working properly.
cdo -V:
Climate Data Operators version 1.6.1 (http://code.zmaw.de/projects/cdo)
Compiler: gcc -std=gnu99 -g -O2 -pthread
version: gcc (GCC) 4.4.6 20120305 (Red Hat 4.4.6-4)
Compiled: by itmh on eld026 (x86_64-unknown-linux-gnu) Jul 10 2013 15:06:17
with: PTHREADS NC4 OPeNDAP Z
filetype: srv ext ieg grb nc nc2 nc4 nc4c
CDI library version : 1.6.1 of Jul 10 2013 15:06:11
CGRIBEX library version : 1.6.1 of Jun 27 2013 15:38:33
netCDF library version : 4.1.3 of Jun 20 2013 16:23:57 $
HDF5 library version : 1.8.11
SERVICE library version : 1.3.1 of Jul 10 2013 15:06:00
EXTRA library version : 1.3.1 of Jul 10 2013 15:05:56
IEG library version : 1.3.1 of Jul 10 2013 15:05:58
FILE library version : 1.8.2 of Jul 10 2013 15:05:56
RE: consects - Added by Ralf Mueller over 11 years ago
eca_cdd creates a very long variable name, maybe it's to long for xconv