Re: gsn_csm_xy (3 lines on one plot)

Hi Debasish et al.
commenting these lines out has no effect.
It was precisely xy_2.ncl that I was trying to mimick, but to no avail. It
must be something
simple that I am just not seeing.

thanks,
Dave

On Mon, Oct 11, 2010 at 10:43 AM, <debasish@sca.uqam.ca> wrote:

> Dear Dave,
> Could please try to run the script by commenting out following linemm
>
> ;res@gsnDraw = False
> ;res@gsnFrame = False
>
> ;draw(plot)
> ;frame(wks)
>
> Also check the example xy_2.ncl
> (http://www.ncl.ucar.edu/Applications/Scripts/xy_2.ncl)
>
> Thanks
>
> Debasish
>
> > Dear All,
> > should be a simple question, I am just a bit rusty with NCL.
> > I want to plot 3 timeseries of data points on a single xy graph.
> >
> > I can plot each one individually fine, but when I try to plot all three
> > only
> > one shows up.
> >
> > thanks in advance,
> > Dave
> >
> > -------------------------------------ncl
> > code------------------------------------------
> > load "/usr/local/lib/ncarg/nclex/gsun/gsn_code.ncl"
> > load "/usr/local/lib/ncarg/nclex/gsun/gsn_csm.ncl"
> > load "/usr/local/lib/ncarg/nclscripts/csm/contributed.ncl"
> > load "/usr/local/lib/ncarg/nclscripts/csm/shea_util.ncl"
> >
> > begin
> >
> > ; Read in Data
> >
> > ncol = 14
> > nlvl = 525599
> >
> > diri = "/home/dadams/Desktop/CPTEC_data/ZF2_swaying_tower/"
> > fili = "GPS_met_CTT_GIPSY_INPA_NOAA_2009.dat"
> > ZF2data = asciiread(diri+fili,(/nlvl,ncol/), "float")
> >
> > julian_day = ZF2data(:,1)
> > pw_RBDC = ZF2data(:,9)*10.0
> > pw_INPA = ZF2data(:,12)*10.0
> > pw_NOAA = ZF2data(:,13)*10.0
> >
> > julian_day!0 = "julian_day"
> > pw_RBDC!0 = "julian_day"
> > pw_INPA!0 = "julian_day"
> > pw_NOAA!0 = "julian_day"
> >
> > julian_day&julian_day = julian_day
> > pw_RBDC&julian_day = julian_day
> > pw_INPA&julian_day = julian_day
> > pw_NOAA&julian_day = julian_day
> >
> > pw_RBDC@units = "mm"
> > pw_INPA@units = "mm"
> > pw_NOAA@units = "mm"
> >
> >
> > ; create new data set for plotting the 3 lines on one figure
> > data1 = new((/3,dimsizes(julian_day&julian_day)/),float)
> > data1(0,:) = pw_INPA
> > data1(1,:) = pw_RBDC
> > data1(2,:) = pw_NOAA
> > data1@units = "mm"
> > data1@_FillValue = -9999.99
> >
> > wks = gsn_open_wks ("pdf","Figure_1") ; open ps or pdf
> > file
> >
> > ; plot figure
> > res = True ; plot mods desired
> > res@gsnDraw = False
> > res@gsnFrame = False
> >
> > res@xyMarkLineModes = "markers"
> > res@xyMarkers = (/3,3,3/)
> > res@xyMarkerColor = (/"blue","black","red"/)
> > res@tiYAxisString = "(PWV mm)" ; axis string
> >
> > res@trXMaxF = 197.0
> > res@trXMinF = 193.0
> > res@trYMaxF = 56.0 ; pwv in mm
> > res@trYMinF = 32.0
> >
> > plot = gsn_csm_xy(wks,julian_day({193:197}),data1(:,{193:197}),res)
> >
> > draw(plot)
> > frame(wks)
> > end
> >
> ---------------------------------------------------------------------------------------------------------------------------------------------
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Oct 11 09:11:37 2010