Re: help

Also, it is recommended to use the new dim_***_n functions where possible

eg, change

    olrSeaClm = dim_avg_Wrap( olrClm(lat|:,lon|:,month|i_season(ns,:)) )

to
    olrSeaClm = dim_avg_n_Wrap( olrClm(month|i_season(ns,:),:,:), 0 )

===
This had nothing to do with your PLCHHQ error.

Good lick

kalim ullah wrote:
> Dear Imran, thanks for your useful guidance and information.
> regrads,
> kalim
> --- On *Thu, 12/17/09, IMRAN NADEEM /<qphoton@gmail.com>/* wrote:
>
>
> From: IMRAN NADEEM <qphoton@gmail.com>
> Subject: Re: help
> To: "kalim ullah" <kalimrid@yahoo.com>
> Date: Thursday, December 17, 2009, 6:23 AM
>
> Dear Kalim,
> I found on the NCL webpage at:
>
> http://www.ncl.ucar.edu/Document/Graphics/error_msg.shtml
>
>
> that error is caused by missing ".hluresfile". Here you can find
> the ".hluresfile"
>
> http://www.ncl.ucar.edu/Document/Graphics/hlures.shtml
>
>
> copy it and put it in your home directory.
>
> Regards
> Imran Nadeem
>
>
>
>
>
> On Thu, Dec 17, 2009 at 3:17 PM, kalim ullah <kalimrid@yahoo.com
> <http://us.mc537.mail.yahoo.com/mc/compose?to=kalimrid@yahoo.com>>
> wrote:
>
> Dear Imran, first thanks for your quick response, so the error
> is like this:
> PLCHHQ - CHARACTER NUMBER 19 (a) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 20 (y) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 21 ()) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 26 (9) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 28 (9) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 29 (-) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 33 (8) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 13 (c) IS NOT A LEGAL FUNCTION CODE
> .................
> ......
> regards,
> kalim
>
>
> --- On *Thu, 12/17/09, IMRAN NADEEM /<qphoton@gmail.com
> <http://us.mc537.mail.yahoo.com/mc/compose?to=qphoton@gmail.com>>/*
> wrote:
>
>
> From: IMRAN NADEEM <qphoton@gmail.com
> <http://us.mc537.mail.yahoo.com/mc/compose?to=qphoton@gmail.com>>
> Subject: Re: help
> To: "kalim ullah" <kalimrid@yahoo.com
> <http://us.mc537.mail.yahoo.com/mc/compose?to=kalimrid@yahoo.com>>
> Date: Thursday, December 17, 2009, 6:04 AM
>
> post the error you are getting while running the script?
>
> On Thu, Dec 17, 2009 at 2:32 PM, kalim ullah
> <kalimrid@yahoo.com
> <http://us.mc537.mail.yahoo.com/mc/compose?to=kalimrid@yahoo.com>>
> wrote:
>
> Dear NCL Users,
> I am trying to find the climatology and using the
> example (climo_6.ncl
> <http://www.ncl.ucar.edu/Applications/Scripts/climo_6.ncl>)
> but the script is not running properly, please anyone
> have a look on the script given below:
> thanks and regards,
> kalim
> ************************
> 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"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
> begin
> plev = 850 ; 850hPa
> ymStrt = 197901 ; start yyyymm
> ymLast = 200812 ; last
> yrStrt = ymStrt/100
> yrLast = ymLast/100
> wksType = "x11"
> wksName = "climo" ; ": "+yrStrt+"_"+yrLast
> diri = "/disk1/home/bwd/" ; input directory
> ;***********************************************************
> fili = "sst.mnmean.nc <http://sst.mnmean.nc/>"
> f = addfile (diri+fili , "r")
> time = f->time ; yyyymmdd (all times)
> ;yyyymm = time/100 ; yyyymm
> yyyymm = ut_calendar(time, -1) ; convert to yyyymm
> ntStrt = ind(yyyymm.eq.ymStrt) ; index start
> ntLast = ind(yyyymm.eq.ymLast) ; index last
> ;sst = f->sst(ntStrt:ntLast,:,:) ; year-month
> ; Climatology
> sst = short2flt( f->sst(ntStrt:ntLast,:,:))
> sstClm = clmMonTLL( sst) ; (month,lat,lon)
> delete(yyyymm)
> delete(time)
> delete(sst)
> ; land-sea mask distributed with NCL
> fmsk =
> addfile("$NCARG_ROOT/lib/ncarg/data/cdf/landsea.nc
> <http://landsea.nc/>","r")
> lsmask = landsea_mask(fmsk->LSMASK,sstClm&lat,sstClm&lon)
> sstClm = mask(sstClm,lsmask.eq.1,False) ; mask out land
> points
> printVarSummary( sstClm )
> ;printMinMax( sstClm, True )
> ;***********************************************************
> ; Read year-month PRC; Create climatology for desired period
> ;***********************************************************
> fili = "precip.mon.mean.nc <http://precip.mon.mean.nc/>"
> f = addfile(diri+fili , "r")
> time = f->time ; days since ...
> yyyymm = ut_calendar(time, -1) ; convert to yyyymm
> ntStrt = ind(yyyymm.eq.ymStrt) ; index start
> ntLast = ind(yyyymm.eq.ymLast) ; index last
> prc = f->precip(ntStrt:ntLast,:,:) ; year-month
> prcClm = clmMonTLL( prc) ; Climatology
> ; (12,lat,lon)
> delete(yyyymm)
> delete(time)
> delete(prc)
> printVarSummary( prcClm )
> printMinMax( prcClm, True )
> ;***********************************************************
> ; Read year-month U ; Create climatology for desired period
> ; Convert type short to float
> ;************************************************************
> fili = "uwnd.mon.mean.nc <http://uwnd.mon.mean.nc/>"
> f = addfile (diri+fili , "r")
> time = f->time ; days since ...
> yyyymm = ut_calendar(time, -1) ; convert to yyyymm
> ntStrt = ind(yyyymm.eq.ymStrt) ; index start
> ntLast = ind(yyyymm.eq.ymLast) ; index last
> uwnd = short2flt( f->uwnd(ntStrt:ntLast,{plev},:,:) )
> uwndClm = clmMonTLL( uwnd ) ; Climatology
> ; (12,lat,lon)
> delete(time)
> delete(yyyymm)
> delete(uwnd)
> ;printVarSummary( uwndClm )
> ;printMinMax( uwndClm, True )
> ;***********************************************************
> ; Read year-month OLR; Create climatology for desired period
> ; Convert type short to float
> ;***********************************************************
> fili = "olr.mon.mean.nc <http://olr.mon.mean.nc/>"
> f = addfile (diri+fili , "r")
> time = f->time ; days since ...
> yyyymm = ut_calendar(time, -1) ; convert to yyyymm
> ntStrt = ind(yyyymm.eq.ymStrt) ; index start
> ntLast = ind(yyyymm.eq.ymLast) ; index last
> olr = short2flt( f->olr(ntStrt:ntLast,:,:) ) ; year-month
> olrClm = clmMonTLL( olr ) ; Climatology
> ; (12,lat,lon)
> delete(yyyymm)
> delete(time)
> delete(olr)
> ;printVarSummary( olrClm )
> ;printMinMax( olrClm, True )
> ;************************************************
> ; plots ; Compute 6-month climatologies
> ;************************************************
> ; (2,6)
> season =(/(/ 5, 6, 7, 8, 9,10/)\ ; May-Oct [summer]
> , (/ 1, 2, 3, 4,11,12/) /) ; Nov-Apr [winter]
> i_season = season - 1 ; NCL indices
> season_label = (/ "May-Oct Climatology", "Nov-Apr
> Climatology"/)
> plot = new ( 3, "graphic")
> wks = gsn_open_wks(wksType, wksName)
> gsn_define_colormap(wks,"amwg")
> ;************************************************
> ; resource list for first (color) data array
> ;************************************************
> res1 = True
> res1@gsnDraw = False ; don't draw
> res1@gsnFrame = False ; don't advance frame
> res1@gsnStringFontHeightF = 0.0125 ; make larger than
> default
> res1@lbOrientation = "Vertical" ; vertical label bar
> res1@lbLabelFontHeightF = 0.01 ; make labels larger
> res1@pmLabelBarOrthogonalPosF = -0.025 ; move closer to plot
> res1@lbLabelAutoStride = True ; optimal label stride
> res1@gsnSpreadColors = True ; use full range of colors
> res1@mpMinLonF = 30. ; select a subregion
> res1@mpMaxLonF = 120.
> res1@mpMinLatF = 0.
> res1@mpMaxLatF = 50.
> res1@mpLandFillColor = "background" ; color of land
> ;res1@cnFillDrawOrder = "Predraw"
> ;************************************************
> ; resource list for second (contour only) data array
> ;************************************************
> res2 = True
> res2@gsnDraw = False
> res2@gsnFrame = False
> res2@cnLevelSpacingF = 2. ; set contour spacing
> res2@cnInfoLabelOn = False ; do not want
> res2@cnLineLabelsOn = True
> res2@cnLabelMasking = True
> ;************************************************
> ; resource list for panel
> ;************************************************
> resP = True ; modify the panel plot
> resP@gsnMaximize = True ; make large
> ;resP@gsnPaperOrientation = "portrait" ; force "portrait"
> delete(sstClm@long_name) ; do not want
> res1@gsnRightString = "SST (C): contour"
> do ns =0,1 ; 2 seasons
> res1@cnLevelSpacingF = 20.0 ; set contour spacing
> res1@gsnLeftString = "OLR: ("+olrClm@units+")"
> olrSeaClm = dim_avg_Wrap(
> olrClm(lat|:,lon|:,month|i_season(ns,:)) )
> sstSeaClm = dim_avg_Wrap(
> sstClm(lat|:,lon|:,month|i_season(ns,:)) )
> plot(0) =
> gsn_csm_contour_map_overlay(wks,olrSeaClm,sstSeaClm,res1,res2)
> res1@gsnLeftString = "NCEP: U: ("+uwndClm@units+")"
> res1@cnLevelSpacingF = 2.0 ; set contour spacing
> uwndSeaClm=
> dim_avg_Wrap(uwndClm(lat|:,lon|:,month|i_season(ns,:)) )
> plot(1) =
> gsn_csm_contour_map_overlay(wks,uwndSeaClm,sstSeaClm,res1,res2)
> res1@gsnLeftString = "GPCP PRC: ("+prcClm@units+")"
> res1@cnLevelSelectionMode = "ExplicitLevels" ; set
> explicit contour levels
> ;res1@cnLevels = (/1,2,4,6,8,10,12,14/)
> res1@cnLevels = ispan(1,14,1)
> prcSeaClm = dim_avg_Wrap(
> prcClm(lat|:,lon|:,month|i_season(ns,:)) )
> plot(2) = gsn_csm_contour_map_overlay(wks,
> prcSeaClm,sstSeaClm,res1,res2)
> delete(res1@cnLevelSelectionMode)
> delete(res1@cnLevels)
> resP@txString = season_label(ns) +": "+ yrStrt+"-"+yrLast
> gsn_panel(wks,plot,(/3,1/),resP) ; now draw as one plot
> end do
> 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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Dec 17 08:20:34 2009