Project

General

Profile

consects

Added by Grace Redmond almost 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 Grace Redmond almost 11 years ago

Sorry, I meant ltc,1 not gtc,1

RE: consects - Added by Ralf Mueller almost 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 almost 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 almost 11 years ago

Hi, I have uploaded the first year of the file

RE: consects - Added by Ralf Mueller almost 11 years ago

where?

RE: consects - Added by Grace Redmond almost 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 almost 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"

RE: consects - Added by Grace Redmond almost 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 almost 11 years ago

eca_cdd creates a very long variable name, maybe it's to long for xconv

    (1-9/9)