[ncl-install] hello please can anyone can help me in my problem

z3303149 at unsw.edu.au z3303149 at unsw.edu.au
Thu Jun 25 23:55:08 MDT 2009



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


More information about the ncl-install mailing list