Ammu,
What about using histograms?
http://www.ncl.ucar.edu/Applications/histo.shtml
This function creates bar charts that contain the frequencies of ranges of datas, or exact values of data.
--Mary
On Oct 3, 2010, at 10:57 PM, Ammu Priya wrote:
> Hai!,,
> i read dennis reply.i am really sorry.I just removed all my commented lines.Actually this script doesn't give me any error.the only thing is am not sure about the calculation part .I want the frequency count for all the three data(june,july,Aug,Sep) sets.in the range of 0-0.001 ,0.001-2,.....till 40. precipitation counts for all the three data sets.I hope graphical resources is fine.i am not sure whether pdfx function gives the number counts.that s the only thing has to be clarified.here nothing is unappropriate..
> Thank you.
>
>
> load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$HOME/ncl/lib/scripts/my_scripts.ncl"
>
> begin
>
> eps = 0
>
> fps = "Precbar_GPCP-TRMM-IMD_june98_frequency_col"
> fps@OverWrite = True
> fps@View = True
> diragpcp="/user1/akila/plot/frequncydistribution/gpcp1deg/"
> diraimd="/user1/akila/plot/frequncydistribution/imd1deg/"
> diratrmm="/user1/akila/plot/frequncydistribution/trmm1deg/"
> filesgpcp = systemfunc("ls "+diragpcp+"prec.gpcp.1dd.p1d.1998*")
> filesimd = systemfunc("ls "+diraimd+"prec.imd.mo.360x180.1998*")
> filestrmm = systemfunc("ls "+diratrmm+"prec.trmm.3b42.360x180.1998*")
> print(filesgpcp)
> print(filesimd)
> print(filestrmm)
>
> fgpcp = addfiles(filesgpcp,"r")
> ListSetType (fgpcp,"cat")
>
> fimd = addfiles(filesimd,"r")
> ListSetType (fimd,"cat")
>
> ftrmm = addfiles(filestrmm,"r")
> ListSetType (ftrmm,"cat")
>
> precgpcp = addfiles_GetVar(fgpcp,filesgpcp,"prec")
> precimd = addfiles_GetVar(fimd,filesimd,"prec")
> prectrmm = addfiles_GetVar(ftrmm,filestrmm,"prec")
>
> printVarSummary(precgpcp)
> printVarSummary(precimd)
> printVarSummary(prectrmm)
>
> gpcpjjas= precgpcp(time|151:272,lat|:,lon|:)
> gpcpjune= precgpcp(time|151:180,lat|:,lon|:)
> gpcpjuly= precgpcp(time|181:211,lat|:,lon|:) gpcpsep= precgpcp(time|243:272,lat|:,lon|:)
> printVarSummary(gpcpjune)
> t=gpcpjune(:,{8:28},{70:90})
> printVarSummary(t)
> print(t)
> printVarSummary(gpcpjuly)
> printVarSummary(gpcpaug)
> printVarSummary(gpcpsep)
> imdjjas =precimd(time|151:272,lat|:,lon|:)
> imdjune =precimd(time|151:180,lat|:,lon|:)
> imdjuly=precimd(time|181:211,lat|:,lon|:)
> imdaug= precimd(time|212:242,lat|:,lon|:)
> imdsep=precimd(time|243:272,lat|:,lon|:)
> printVarSummary(imdjune)
> printVarSummary(imdjuly)
> printVarSummary(imdaug)
> printVarSummary(imdsep)
> trmmjjas =prectrmm(time|151:272,lat|:,lon|:)
> trmmjune =prectrmm(time|151:180,lat|:,lon|:)
> trmmjuly=prectrmm(time|181:211,lat|:,lon|:)
> trmmaug=prectrmm(time|212:242,lat|:,lon|:)
> trmmsep=prectrmm(time|243:272,lat|:,lon|:)
> printVarSummary(trmmjune)
> printVarSummary(trmmjuly)
> printVarSummary(trmmaug)
> printVarSummary(trmmsep)
> x1=ndtooned(gpcpjune(:,{8:28},{70:90}))
> printMinMax(x1,0)
> y1=ndtooned(imdjune(:,{8:28},{70:90}))
> z1=ndtooned(trmmjune(:,{8:28},{70:90}))
> printVarSummary(x1)
> printVarSummary(y1)
> printVarSummary(z1)
> printMinMax(x1,0)
> printMinMax(y1,0)
> printMinMax(z1,0)
> print(x1)
> printVarSummary(y1)
> printMinMax(y1,0)
> printVarSummary(z1)
> printMinMax(z1,0)
> opt=True
> opt@bin_min =0
> opt@bin_max =50
> ry=pdfx(gpcpjune(:,{8:28},{70:90}),0,opt)
> opt1 = True
> opt1@bin_min =0.4
> opt1@bin_max =50
> ry1=pdfx(imdjune(:,{8:28},{70:90}),0,opt1)
> opt2 = True
> opt2@bin_min =0.7
> opt2@bin_max =50
> ry2=pdfx(imdjune(:,{8:28},{70:90}),0,opt2)
> print(ry@bin_center+" "+ry)
> printVarSummary(ry)
> printVarSummary(ry1)
> print(ry1@bin_center+" "+ry1)
> printVarSummary(ry2)
> print(ry2@bin_center+" "+ry2)
>
> if(eps .eq. 1)then
> wks = gsn_open_wks("eps",fps);
> else
> wks = gsn_open_wks("x11",fps);
> end if
>
> sres = True
> sres@gsnMaximize = True
> sres@vpWidthF = 0.7
> sres@vpHeightF = 0.5
> sres@vpXF = .10
> sres@trXMinF = 0
> sres@trXMaxF = 40
> sres@trYMinF = 0
> sres@trYMaxF = 61
> sres@gsnDraw = True
> sres@gsnFrame = False
> sres@gsnXYBarChart = True
> sres@gsnXYBarChartBarWidth = 0.55
>
> sres@tmXBMode = "Explicit"
> ; sres@tmXBValues = ispan(2,40,2)
> sres@tmXBValues = fspan(1.45,40.55,20)
> sres@tmXBLabelAngleF = 315
> sres@tmXBLabels = (/"0-0.001","0.001-2","2-4","4-6","6-8","8-10","10-12","12-14","14-16","16-18","18-20","20-22","22-24","24-26","26-28","28-30","30-32","32-34","34-36","36-38","38-40"/)
> sres@tmXBLabelFont=0.4
> sres@tmXBLabelFontHeightF = 0.0100
> sres@tmXTLabelFontHeightF = 0.0110
> sres@tmYLLabelFontAspectF =1.0
> sres@tmYLLabelFontHeightF = 0.0105
> sres@tiMainFontHeightF = 0.015
> sres@tiMainFont = "helvetica"
> sres@tiMainString = "Frequency distribution,JUNE-98"
> sres@tiYAxisFont = "helvetica"
> sres@tiYAxisString = "count(%)"
>
> sres@gsnXYBarChartColors = (/"red"/)
> plot1=gsn_csm_xy(wks,ry@bin_center,ry,sres)
> sres@gsnXYBarChartColors = (/"blue"/)
> plot2=gsn_csm_xy(wks,ry1@bin_center,ry1,sres)
> sres@gsnXYBarChartColors = (/"green"/)
> plot2=gsn_csm_xy(wks,ry2@bin_center,ry2,sres)
>
> lbres = True ; labelbar only resources
> lbres@vpWidthF = 0.25 ; labelbar width
> lbres@vpHeightF = 0.05 ; labelbar height
> lbres@lbBoxMajorExtentF = 0.36 ; puts space between color boxes
> lbres@lbFillColors = (/"red"/)
> lbres@lbMonoFillPattern = True ; Solid fill pattern
> lbres@lbLabelFontHeightF = 0.015 ; font height. default is small
> lbres@lbLabelJust = "CenterLeft" ; left justify labels
> lbres@lbPerimOn = False
>
> labels = (/"GPCP"/)
> gsn_labelbar_ndc(wks,1,labels,0.66,0.80,lbres)
> lbres@lbFillColors = (/"blue"/)
> labels = (/"IMD"/)
> gsn_labelbar_ndc(wks,1,labels,0.65,0.76,lbres) ; draw right labelbar column
> lbres@lbFillColors = (/"green"/)
> labels = (/"TRMM"/)
> gsn_labelbar_ndc(wks,1,labels,0.665,0.72,lbres) ; draw right labelbar column
>
> frame(wks)
> if(eps .eq. 1)then
> mkpdf(fps)
> end if
>
> 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 4 08:51:17 2010
This archive was generated by hypermail 2.1.8 : Mon Oct 04 2010 - 08:55:54 MDT