Thanks Mary,
Every thing is perfect now.
Best wishes
--- 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?
À: "louis Vonder" <appopson_at_yahoo.fr>
Cc: "ncl forum" <ncl-talk_at_ucar.edu>
Date: Mercredi 17 Septembre 2008, 17h27
Hi Louis,
Thanks for sending your files. You are absolutely right, and the problem
is that NCL is picking "nice" values for the XY plot, and hence
the Y axis is going from 0 to 16.
Try adding these two lines to your XY plot resources:
xyres_at_trYMinF = min(u&latitude)
xyres_at_trYMaxF = max(u&latitude)
I thought there was a resource you could set telling NCL not to
pick nice values, but I can't seem to find it. Meanwhile, the above
method should work.
--Mary
On Wed, 17 Sep 2008, louis Vonder wrote:
>
> 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
>
>
Received on Wed Sep 17 2008 - 10:01:24 MDT
This archive was generated by hypermail 2.2.0 : Wed Sep 17 2008 - 10:32:36 MDT