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
This archive was generated by hypermail 2.2.0 : Fri Oct 09 2009 - 08:29:22 MDT