NCEP reanalysis data in grib format

I am working with NCEP reanalysis data in grib format. My NCL script
reads wrong names.

I have two files: prs.grib.mean.y1957, and prs.grib.mean.y1958.
"getfilevarnames" gives me two different lists:

prs.grib.mean.y1958 prs.grib.mean.y1957
(0) initial_time0 initial_time0
(1) lv_ISBL1 lv_ISBL1
(2) lon_2 lon_2
(3) lat_2 lat_2
(4) lv_ISBL4 lv_ISBL4
(5) lv_ISBL5 lv_ISBL5
(6) initial_time0_encoded initial_time0_encoded
(7) initial_time0_hours initial_time0_hours
(8) NCIP_2_ISBL_113 CBTMW_2_ISBL_113
(9) VSTM_2_ISBL_113 CBTZW_2_ISBL_113
(10) USTM_2_ISBL_113 CBMZW_2_ISBL_113
(11) CPOZP_2_ISBL_113 CBQW_2_ISBL_113
(12) CPOFP_2_ISBL_113 CBTW_2_ISBL_113
(13) POP_2_ISBL_113 CBVQ_2_ISBL_113
(14) PROBN_2_ISBL_113 CBUQ_2_ISBL_113
(15) PROB_2_ISBL_113 CBVW_2_ISBL_113
(16) HLCY_2_ISBL_113 CBUW_2_ISBL_113
(17) R_H_2_ISBL_113 RH_2_ISBL_113
(18) SPF_H_2_ISBL_113 SPFH_2_ISBL_113
(19) REL_D_2_ISBL_113 REL_D_2_ISBL_113
(20) REL_V_2_ISBL_113 REL_V_2_ISBL_113
(21) V_VEL_2_ISBL_113 V_VEL_2_ISBL_113
(22) V_POT_2_ISBL_113 V_POT_2_ISBL_113
(23) STRM_2_ISBL_113 STRM_2_ISBL_113
(24) V_GRD_2_ISBL_113 V_GRD_2_ISBL_113
(25) U_GRD_2_ISBL_113 U_GRD_2_ISBL_113
(26) VTMP_2_ISBL_113 VTMP_2_ISBL_113
(27) TMP_2_ISBL_113 TMP_2_ISBL_113
(28) HGT_2_ISBL_113 HGT_2_ISBL_113

The list in the left (prs.grib.mean.y1958) is not correct, because
when I plot "VSTM_2_ISBL_113", for instance, it gives me reasonable
values for "CBMZW_2_ISBL_113". So, I know that the name is switched
somehow. A colleague of mine tried to upload the two data on Matlab,
and the two data have the same name convention.

My NCL script file is very simple:
----------------------------------------------------------------------
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
;***********************************************
; to use the environmental variable
;***********************************************
grib_in = addfile("prs.grib.mean.y1957.grib", "r")
names = getfilevarnames(grib_in)
print(names)
----------------------------------------------------------------------

I don't have a clue......

Thank you.

Sang-ki

-- 
----------------------------------------------------------------------
Dr. Sang-ki Lee
RSMAS/CIMAS
University of Miami
4600 Rickenbacker Causeway
Miami, FL 33149 USA
Tel) 305-361-4521
Fax) 305-361-4412
E-mail) sang-ki.lee_at_noaa.gov
Homepage) http://www.aoml.noaa.gov/phod/people/sklee.html
----------------------------------------------------------------------
Mailing-address:
Dr. Sang-ki Lee
AOML/PHOD
4301 Rickenbacker Causeway
Miami, FL, 33149 USA
----------------------------------------------------------------------
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk