Thanks Mary,
As you can see in the attached (which I loaded for plotting) there is no missing value.
This behaviour wonder me !
Looking again the Example 10 (Panel), you can see that the XY graph (in the right panel) presents the same discrepancy as mine.
--- En date de : Mer 17.9.08, Mary Haley <haley_at_ucar.edu> a écrit :
De: Mary Haley <haley_at_ucar.edu>
Objet: Re: [ncl-talk] Whats wrong in this panel attach plot?
À: appopson_at_yahoo.fr
Cc: "Caspar Ammann" <ammann_at_ucar.edu>, "ncl forum" <ncl-talk_at_ucar.edu>
Date: Mercredi 17 Septembre 2008, 15h24
Louis,
I agree with Caspar, I think you may have some missing values in your XY plot. Try printing "mean_std" and see if the beginning and ending set of points havebeen set to missing.
--Mary
On Sep 17, 2008, at 6:33 AM, louis Vonder wrote:
Sorry Caspar,
I forgot to add my script in the previous message.
Thanks
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
xlat = fspan(1, 15.25, 58)
nlat = 58
xtime = fspan(1, 48, 48)
ntime = dimsizes(xtime)
xlat_at_units = "degrees_north"
xlat_at_long_name = "Latitude"
xtime_at_units = "LST"
xtime_at_long_name = "Time"
M_std = asciiread("Mean_Std_mat_land", (/nlat, 2/), "float")
mean_std = transpose(M_std)
fili = systemfunc("ls Ano_mat_land")
nfil = dimsizes( fili )
do n = 0, nfil - 1
name = stringtocharacter(fili(n))
; name1 = name(0:18)
name2 = chartostring(name)
dummy_data = - asciiread(fili(n), (/nlat, ntime/), "float")
; Set up coordinate arrays so we can overlay on map.
dummy_data!0 = "latitude"
dummy_data!1 = "time"
dummy_data&latitude = xlat
dummy_data&time = xtime
dummy_data_at_units = "%"
dummy_data@_FillValue = -9999.
dummy_data_at_missing_value = dummy_data@_FillValue
dummy_data_at_long_name = "Fractional coverage"
diro = "./" ; Output directory
filo = name+".nc" ; Output file
system("/bin/rm -f " + diro + filo) ; remove if exists
u = dummy_data
;=================================================;
; PLOT 1
;=================================================;
wks = gsn_open_wks("png", name2) ; open a ps file
gsn_define_colormap(wks,"wh-bl-gr-ye-re")
res = True
res_at_gsnDraw = False ; don't draw yet
res_at_gsnFrame = False ; don't advance frame yet
res_at_cnFillOn = True ; turn on color
res_at_gsnSpreadColors = True ; use full range of colormap
res_at_cnInfoLabelOn = False
res_at_cnInfoLabelPerimOn = True
res_at_cnLineLabelsOn = True ; turn on line labels
res_at_cnLevelSelectionMode = "AutomaticLevels" ; set manual contour levels
res_at_gsnSpreadColors = True ; use full color map
res_at_gsnSpreadColorEnd = -3 ; don't use added gray
res_at_lbOrientation ="horizontal" ; vertical label bar
res_at_lbLabelStride = 2
res_at_tmXUseBottom = False ; Keep top axis independent of bottom.
res_at_tmXBLabelsOn = True ; have tick mark labels
res_at_tmXTOn = False ; have tick marks
res_at_tmXBMode = "Explicit" ; label independently
res_at_tmXBValues = (/2, 6, 10, 14, 18, 22, 26, 30, 34, 38, 42, 46/)
res_at_tmXBLabels = (/"01", "03", "05", "07", "09", "11", "13", "15", "17", "19", "21", "23"/)
res_at_tmXBMinorValues = (/2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28 , 30, 32, 34, 36, 40, 42, 44, 46, 48/)
res_at_tiMainString = " " ; add titles
res_at_tiXAxisString = " "
res_at_tiYAxisString = "Latitude"
res_at_gsnLeftString = " " ; add the gsn titles
res_at_gsnCenterString = " "
res_at_gsnRightString = " "
plot = gsn_csm_contour(wks, u, res) ; contour the variable
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;; PLOT 2
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
xyres = True ; xy plot mods desired
xyres_at_vpWidthF = .20 ; set width of second plot
xyres_at_tmXBMinorOn = False ; no minor tickmarks
xyres_at_tmXBLabelStride = 2 ; label stride
xyres_at_gsnDraw = False ; don't draw yet
xyres_at_gsnFrame = False ; don't advance frame yet
xyres_at_txFontHeightF = .015 ; change font height
plot2 = gsn_csm_xy(wks, mean_std, u&latitude, xyres)
newplot = gsn_attach_plots(plot, (/plot2/), res, xyres)
draw(plot)
frame(wks)
end do
end
--- En date de : Mer 17.9.08, Caspar Ammann <ammann_at_ucar.edu> a écrit :
De: Caspar Ammann <ammann_at_ucar.edu>
Objet: Re: [ncl-talk] Whats wrong in this panel attach plot?
À: "louis Vonder" <appopson_at_yahoo.fr>
Cc: "ncl forum" <ncl-talk_at_ucar.edu>
Date: Mercredi 17 Septembre 2008, 13h17
Louis, I can't see your script, but you appear to have applied a slight smoother to the selected time slices. If so, then the end points where the smoother window starts to extend outside the data range will be set to missing. You could use data from further north or south to compute the smooth meridional profile, but only plot in the window you selected.Caspar
On Sep 17, 2008, at 3:01 AM, louis Vonder wrote:
Dear NCL users,
Here attached script (Inspired from www.ncl.ucar.edu/Applications/panel.html Example 10) that I am using to plot attached panel.
I am also sending data I used and figure produced by this script.
I don't know why the second figure (the one using gsn_csm_xy) is not covering the all Y-axis?
There is something wrong in my script?
Thanks for your help.
<Ano_mat_land><Mean_Std_mat_land><Ano_mat_land.png>_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Caspar M. AmmannNational Center for Atmospheric ResearchClimate and Global Dynamics Division - Paleoclimatology1850 Table Mesa DriveBoulder, CO 80307-3000email: ammann_at_ucar.edu tel: 303-497-1705 fax: 303-497-1348
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.2.0 : Wed Sep 17 2008 - 10:32:36 MDT