hello please can anyone can help me in my problem

when i run the program that i see this shown below please can anyone help me

Variable: p
Type: float
Total Size: 101078376 bytes
            25269594 values
Number of Dimensions: 4
Dimensions and sizes: [lev | 98] x [time | 23] x [lat | 101] x [lon | 111]
Coordinates:
Number Of Attributes: 1
  _FillValue : -999
(0) gsn_csm_pres_hgt: Fatal: The first dimension of the input data must
(0) have a coordinate variable called 'lev.'
(0) Cannot create plot.
fatal:Illegal right-hand side type for assignment
fatal:Execute: Error occurred at or near line 55 in file test.ncl~

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

;the input data directory for two dates
diri1 = "/srv/ccrc/data04/s3236814/mm5/nnrp/"
diri2 = "/srv/ccrc/data04/s3236814/mm5/nnrp/"

;the output directory
diro = "../out/"

f1 = addfile(diri1+"pp.199811.mon.nc","r")
f2 = addfile(diri2+"MM5_header.nc","r")
pp=f1->pp
ptop = f2->ptop
ps=f2->pstarcrs
lev= f2->sigma_level
lat = f1->latitcrs
lon=f1->longicrs
ntim = 98
klev = 32

time=f1->time

ppps = conform(pp,ps,(/2,3/))

pplev = conform(pp,lev,1)

p = pp - ptop - ppps*pplev

 p!1 = "time"
 p!0 = "lev"
 p!2 = "lat"
 p!3 = "lon"

printVarSummary(p)

wks = gsn_open_wks("x11",diro+"panel")

plot = gsn_csm_pres_hgt(wks,p(0,:,:,0),False)

end
_______________________________________________
ncl-install mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-install
Received on Thu Jun 25 2009 - 23:55:08 MDT