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

Dennis Shea shea at ucar.edu
Fri Jun 26 06:16:49 MDT 2009


You posted to the wrong mail list.
ncl-install is for installation questions.
ncl-talk is for question like this one.

[1]
My *guess* is it should be
   p!0 = "time"
   p!1 = "lev"


[2]
Read the documentation
http://www.ncl.ucar.edu/Document/Graphics/Interfaces/gsn_csm_pres_hgt.shtml

The 'lev' and 'lat' dimension should have coordinates
and units associated with them.

z3303149 at unsw.edu.au wrote:
> 
> 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


More information about the ncl-install mailing list