Titles on Spectral Plot ?

From: Helen Parish <hparish_at_nyahnyahspammersnyahnyah>
Date: Wed, 18 Feb 2009 05:40:37 -0800

Does anyone know how to put titles on the axes of a spectral plot
which uses "specx_anal" ?. I have tried using

res_at_tiXaxisString = "Days"

analogous to the example given in Application Examples, however I get
the message :

warning:tiXaxisString is not a valid resource in
zeusranddaily7608.60km.spec.test_xy at this time

Does anyone know what I need to use here to label the axes ?.

I enclose my current script below.

Thanks,
Helen.

;***********************
; fourtest.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

; Keep all files open when have a lot of files :
   setfileoption("nc","SuppressClose",False)

   diri = "./"
; fili = systemfunc("cd "+diri+" ; ls zeusorigspongej*.cam2.h0*nc")
; fili = systemfunc("cd "+diri+" ; ls
zeusrand09125daily*.cam2.h0*nc")
   fili = systemfunc("cd "+diri+" ; ls
zeusrand09125daily*.cam2.h0*76-08*nc")

   f = addfiles (fili,"r")
   U = addfiles_GetVar(f,fili,"U")
   lat = addfiles_GetVar(f,fili,"lat")
   lon = addfiles_GetVar(f,fili,"lon")
   lev = addfiles_GetVar(f,fili,"lev")
   printVarSummary( U )

   dimu = dimsizes( U )
   ntim = dimu(0)
   klvl = dimu(1)
   nlon = dimu(2)
   mlon = dimu(3)

  uavg = dim_avg_Wrap(U)

;***********************
; Create Plot
;***********************
   d = 1
; sm = 3
   sm = 1
   pct = 0.1

; kl = 10
    kl = 21
    nl = 96
    ml = 40
    nt = 0
; spec = specx_anal(u(lev_p|kl, lat|nl, lon|:, time|nt),d,sm,pct)
   spec = specx_anal(U(lev|kl, lat|nl, lon|ml, time|:),d,sm,pct)

   wks = gsn_open_wks ("pdf",
"zeusranddaily7608.60km.spec.test" ) ; open workstation
; wks = gsn_open_wks ("ps", "rand09125k260117speca" ) ;
open workstation
   gsn_define_colormap(wks,"rainbow") ; choose colormap

   res = True
; res_at_cnFillOn = True
; res_at_lbLabelAutoStride = True
; res_at_gsnMaximize = True ; if [ps, eps, pdf] make large
; res_at_gsnSpreadColors = True ; span color map

      res_at_gsnCenterString = "zeusranddaily7608.60km"
      res_at_tiXaxisString = "Days"
; res_at_tiYaxisString = "Variance"

; plot = gsn_csm_xy(wks,spec_at_frq,spec_at_spcx,res)
    plot = gsn_csm_xy(wks,1/spec_at_frq,spec_at_spcx,res)

end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Feb 18 2009 - 06:40:37 MST

This archive was generated by hypermail 2.2.0 : Thu Feb 19 2009 - 09:54:51 MST