Re: Problem with dim_avg_n

Hi Louis,
I just ran a quick test here and everything
looks fine to me. I notice in your script that
for the dim_avg_n line that you are subsetting
time from 0:23, but for the dim_avg line you are
not. Could that explain the difference in the plots?
Best regards,
Adam

louis Vonder wrote:
> Dear NCL users,
>
> Here attached two figure using the NCL functions "dim_avg_n" and "dim_avg".
>
> I am wondering because the produced figures look completly different !!!!
>
> I don't why the figure using dim_avg_n is smoothed.
>
> It is a bug or there is something which is bad in my script?
>
>
> Thanks for your help
>
> Louis
>
>
> ;*****************************************************
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> ;*****************************************************
> begin
>
>
>
>
> a = addfile ("TotPrecc.nc", "r")
>
>
> lat = a->lat
> lon = a->lon
> time = a->time
> t= a->TotPrecc
> t1 = t(time|0:23,lat|:,lon|:)
> t2= t1(lat|:,lon|:,time|0)
> t2=dim_avg(t1(lat|:,lon|:,time|:))
> ; t1=dim_avg_n(t(0:23,:,:),0)
> wks = gsn_open_wks ("ps", "wrf_dim_avg")
>
>
>
> res = True ; plot mods desired
> gsn_define_colormap(wks,"gui_default") ; choose colormap*
> res_at_gsnSpreadColors = True ; use full range
> of colors
> res_at_gsnSpreadColorStart = 7
> res_at_gsnSpreadColorEnd = 101
>
>
> res_at_cnFillOn = True ; turn on color
> res_at_cnLinesOn = False ; no contour lines
> ;res_at_gsnSpreadColors = True ; use full
> color map
> res_at_lbLabelAutoStride = True ; every other
> label
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> ;;;;;;;;;;;;;;;;;;;;;;;;draw national boundery
> map;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> ; Map
> stuff
> ;
> res_at_mpDataBaseVersion = "Ncarg4_0" ; use finer
> database ;
> res_at_mpOutlineOn =
> True ;
> res_at_mpFillOn = False ; turn off gray
> fill ;
> res_at_mpOutlineBoundarySets = "National" ; turn on country
> boundaries ;
>
> res_at_mpNationalLineColor = "Black" ;"Red"
> res_at_mpNationalLineDashPattern = 0 ;1 ;for dashpattern
> res_at_mpNationalLineThicknessF = 1.5
>
> res_at_mpGeophysicalLineColor = "Black" ; color of cont.
> outlines ;
> res_at_mpGeophysicalLineThicknessF = 1.5 ; thickness of
> outlines;; ;
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
>
> ; res_at_cnLevelSelectionMode = "ManualLevels" ; manual levels
> ; res_at_cnMinLevelValF = 0.
> ;res_at_cnMaxLevelValF = 300.
> ; res_at_cnLevelSpacingF = 4.
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> res_at_cnLevelSelectionMode = "ExplicitLevels"
> res_at_cnLevels = (/ .5, 1., 2.,5., 10.,15.,20. , \
> 30., 50, 75, 100, 110./)
> ; res_at_cnFillColors = (/"White","White","DarkOliveGreen1", \
> ; "DarkOliveGreen3","Chartreuse", \
> ;
> "Chartreuse3","Green","ForestGreen", \
> ; "Yellow","Orange","Red","Violet"/)
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> res_at_pmTickMarkDisplayMode = "Always" ; turn on tick marks
>
> res_at_tiMainString = "with dim_avg"
>
> res_at_gsnAddCyclic = False ; regional data
>
> ; The following resources are Required to plot this projection correctly
>
> ;res_at_mpProjection = "Stereographic" ; projection
>
> res_at_mpLimitMode = "Corners" ; method to zoom
> ;res_at_mpLeftCornerLatF = min(lat)
> ;res_at_mpLeftCornerLonF = min(lon)
> ;res_at_mpRightCornerLatF = max(lat)
> ;res_at_mpRightCornerLonF = max(lon)
>
>
> res_at_mpLeftCornerLatF = -1.143
> res_at_mpLeftCornerLonF = -0.256
> res_at_mpRightCornerLatF = 12.983
> res_at_mpRightCornerLonF = 22.084
>
>
> ;res_at_mpRelativeCenterLon = True ; set a center lon
> ;res_at_mpCenterLonF = lon(147,172) ; center lon
> ;res_at_mpRelativeCenterLat = True ; set a center lat
> ; res_at_mpCenterLatF = 90. ; center lat
> ; res_at_tfDoNDCOverlay = True ; do not
> transform data
>
> plot = gsn_csm_contour_map (wks,t2,res)
> ;plot = gsn_csm (wks,res)
>
>
> end
>
>
>
> ------------------------------------------------------------------------
>
>
> ------------------------------------------------------------------------
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
__________________________________________________
Adam Phillips 
                 asphilli_at_ucar.edu
National Center for Atmospheric Research   tel: 
(303) 497-1726
Climate and Global Dynamics Division 
fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000 
http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk