Hi all,
I am try to create histogram of number of convective events for each month
of the year.
The input file is simply the integer value of the day of the year for the
reported convective event. However, the plots are incorrect.
My version of ncl is NCAR Command Language Version 6.0.0
Any suggestions?
thanks,
Dave Adams
; test.ncl
; This ncl program reads in integer days of the year
; for for convective events to make a
; histogram of monthly convective events.
load "/usr/local/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/contributed.ncl"
load "/$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
ncol1 = 1
nlvl = 278
seasonal_data =
asciiread("Seasonal_dates_convective_events.dat",(/nlvl,ncol1/), "integer")
day = seasonal_data(:,0)
wks = gsn_open_wks("pdf", "Time_local_hour_event_histogram")
res = True
res@gsnDraw = False
res@gsnFrame = False
res@gsnHistogramDiscreteBinValues = ispan(30,360,30)
plot = gsn_histogram(wks,day,res)
delete(res@gsnDraw)
delete(res@gsnFrame)
res@tmXBMode = "Explicit" ; Define tick mark labels.
res@tmXBValues = plot@MidBarLocs
res@tmXBLabels = (/"J","F","M","A","M","J","J","A","S","O","N","D"/)
plot=gsn_histogram(wks,day,res) ; create histogram
draw(plot)
end
~
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.1.8 : Thu Feb 28 2013 - 14:47:31 MST