[Thread Prev][Thread Next][Index]

Re: with Descriptor file: ** unknown netCDF error code: -31



On Fri, 14 May 2004, Ansley Manke wrote:

> Your time unit of 2629746 corresponds to monthly data.  (1/12 of 365.2425
> days).  You have
> 
>  S_START =   1,
>  S_END   =   239,
>  S_DELTA =   1,
> 
> I assume this corresponds to 20 years of data.  But the time axis is 
> computed
> using (S_END-S_START-1)/S_DELTA = (239-1)  months.  This is less than
> 20 years; and that looks to me like it's why you're seeing the error 
> message,
> 
>  WARNING - NON-INTEGER NUMBER OF STEPS PER STEPFILE
>  INDICATED IN DESCRIPTOR - TRYING TO CONTINUE ANYWAY
> 
> Following the suggestion of using a DEFINE AXIS command
> to make the axis that you want, and using its steps to
> define your axis, if you want:
> 
> yes? def
axis/t=15-jan-1980:15-dec-1999:1/units=months/t0="15-dec-1979:00:00" tax
> yes? list t[gt=tax]
> 
> Then this suggests thatyou want S_END to be 240.


Hi Ansley, 

I gave it a try using 240 and it works!  It is a bit non-intuitive (for me
anyway) though.  Rather than using (S_END-S_START-1)/S_DELTA I would have
thought it would be (S_END-S_START+1)/S_DELTA.  I have 239 time slices so
I would (and did) guess that S_END would be 1 and S_END would be 239.
At least it works and it is a very easy fix.

Thanks very much,

Steve



> 
> Ansley Manke
> 
> 
> Steve Cousins wrote:
> 
> >>I have a series of 20 netCDF files, each with 12 time slices (except the
> >>last file only has 11 time slices). I have a descriptor file to connect
> >>these:
> >>
> >>&FORMAT_RECORD
> >>D_TYPE	= '  MC',
> >>D_FORMAT	= '  1C',
> >>/
> >>&BACKGROUND_RECORD
> >>D_TITLE	= 'BIOLOGY',
> >>D_T0TIME	= '15-dec-1979 00:00:00',
> >>D_TIME_UNIT	= 2629746.0,
> >>D_TIME_MODULO	= .false.,
> >>/
> >>&MESSAGE_RECORD
> >>D_MESSAGE     = ' ',
> >>D_ALERT_ON_OPEN      = F,
> >>D_ALERT_ON_OUTPUT    = F,
> >>/
> >>&EXTRA_RECORD
> >>/
> >>&STEPFILE_RECORD
> >> S_FILENAME    = '/neptune/usr2/pacific-g-46-99/bio.cdf.001',
> >> S_AUX_SET_NUM   =           0,
> >> S_START =   1,
> >> S_END   =   239,
> >> S_DELTA =   1,
> >> S_NUM_OF_FILES  =           20,
> >> S_REGVARFLAG    = '     ',
> >> S_INDEX = '                       '
> >>/ 
> >>&STEPFILE_RECORD
> >>  S_FILENAME    = '**END OF STEPFILES**'
> >>/
> >>
> >>In ferret I periodically (every 12 slices) get:
> >>
> >>	** unknown netCDF error code: -31
> >>
> >>Here is an example:
> >>
> >>[ferret@rocky server]# ferret
> >> 	NOAA/PMEL TMAP
> >> 	FERRET v5.53  
> >> 	Linux(g77) 2.4.x - 09/26/03
> >> 	14-May-04 13:02     
> >>
> >>yes?  set data "/usr2/fer_dsets/descr/bio-g-80-99.des"
> >> WARNING - NON-INTEGER NUMBER OF STEPS PER STEPFILE
> >> INDICATED IN DESCRIPTOR - TRYING TO CONTINUE ANYWAY
> >>yes?  shade/L=1/Z=5/levels=30 co2flx
> >>yes?  shade/L=12/Z=5/levels=30 co2flx
> >> ** unknown netCDF error code: -31
> >>yes?  shade/L=13/Z=5/levels=30 co2flx
> >>yes?  shade/L=11/Z=5/levels=30 co2flx
> >>yes?  shade/L=24/Z=5/levels=30 co2flx
> >> ** unknown netCDF error code: -31
> >>yes?  shade/L=36/Z=5/levels=30 co2flx
> >> ** unknown netCDF error code: -31
> >>yes?  shade/L=48/Z=5/levels=30 co2flx
> >> ** unknown netCDF error code: -31
> >>yes?  shade/L=47/Z=5/levels=30 co2flx
> >>yes?  shade/L=49/Z=5/levels=30 co2flx
> >>
> >>
> >>If I access the file directly (not using the Descriptor file) that
> >>corresponds to slice 13 I can plot that slice fine.
> >>
> >>I can't find what error -31 is.  Anybody have any ideas?
> >>
> >>Thanks,
> >>
> >>Steve
> >>
> >>______________________________________________________________________
> >> Steve Cousins, Ocean Modeling Group    Email: cousins@umit.maine.edu
> >> Marine Sciences, 208 Libby Hall        http://rocky.umeoce.maine.edu
> >> Univ. of Maine, Orono, ME 04469        Phone: (207) 581-4302
> >>
> >>
> >>  
> >>
> >
> >
> 



[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement