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-talkReceived on Thu Mar 09 2006 - 18:38:04 MST
This archive was generated by hypermail 2.2.0 : Tue Mar 14 2006 - 08:36:41 MST