NCEP CSFR grb2 errors

Hello--

I'm trying to process the NCEP CSFR GRB2 temperature output for 1979-2009 and am running into some problems.

The attached script opens the files in a loop and tries to read in a specified variable (ALBDO_P8_L1_GGA0) from each file. The script runs for the first fifteen months with the usual NCL grib warnings, e.g.:
warning: Entry (205) not found in code table file /usr/share/ncarg/grib2_codetables/ncep/4/4.10.table

However, when it reaches the 16th file,the grib error changes slightly to

warning: NclGRIB2: codetable file "/usr/share/ncarg/grib2_codetables/ncep/1/1.3.table" not a valid GRIB2 code table.

It then crashes on the subsequent file (the 17th) with the following error:
fatal: NclGRIB2: codetable directory "/usr/share/ncarg/grib2_codetables" invalid
fatal:Could not open (../Ncep/flxf03.gdas.198006.grb2)
fatal:Either file (f) isn't defined or variable (ALBDO_P8_L1_GGA0) is not a variable in the file
fatal:Execute: Error occurred at or near line 19 in file testncep2.ncl

I would suspect this was something with the libraries, but the script runs for multiple months before getting this error. Also, the crash always occurs on the 17th file, even if I change the decade of data I'm processing, indicating that the error is not due to any one specific file. I'm running NCL v. 5.1.1. Two years of monthly data are loaded to (ftp.cgd.ucar.edu) with the beginning file extension of flxf03.gdas.YYYYMM.grb2.

Any ideas?

Thank you,

Allison

   ;************************************
   ; Script to read in NCEP CSFR data 1979-2009
   ;************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin

; list all files
fili = systemfunc("ls ../Ncep/*.grb2")
nfil = dimsizes(fili)

; file loop
do nf=0, nfil-1
  print(nf+" "+fili(nf))
  f = addfile(fili(nf),"r") ; grib
  ntest=f->ALBDO_P8_L1_GGA0

  printVarSummary(ntest)

  delete(f)
  delete(ntest)
end do

end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jul 20 13:10:11 2011