[ncl-talk] can't figure out why could not plot time axis properly

[Ruksana Abedin]

Hi Rick,

Thank you for your suggestion. Yes, you are right, labels had only two
values

Variable: labels
Type: string
Total Size: 16 bytes
                  2 values
Number of Dimensions: 1
Dimensions and sizes:   [2]
Coordinates:
(0)     12
(1)     01

I have fixed this and now the x-axis shows 1-12 for the months. Thank you
so much for your help.
Can I change this to Jan, Feb... Dec using ncl? Any suggestion will be
highly appreciated.

With many thanks and regards,
Ruksana

On Thu, Nov 26, 2015 at 5:37 PM, Rick Brownrigg <brownrig at ucar.edu> wrote:

> Hi,
>
> I'm not certain, and don't have a good way to test at the moment. However,
> a hint might be the variable "incr = 1"
>
> Maybe print out the variable "labels" ... I' guess that it only contains
> two values. Likewise for for PPT_obs&time(:,1).
>
> Hope that helps, but please write back to the list if not.
>
> Rick
>
>
> On Wednesday, November 25, 2015, Ruksana Abedin <ruksana.abedin at gmail.com>
> wrote:
>
>> Hi.
>>
>> Can you please help me to find out why I could not plot the time axis
>> properly?
>>
>> Below given the script anf the ouput. The axis should show 2-10 but only
>> showing 2-3. How can I change this to Feb-Nov?
>>
>> Thank you for your kind support.
>>
>> Best regards,
>> Ruksana
>>
>>
>> ;*****************************************************
>> 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"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/time_axis_labels.ncl"
>> ;*****************************************************
>> begin
>> ;*****************************************************
>> ;read in data
>> ;*****************************************************
>> f0=addfile ("/ouce-home/students/rege1218/WAH2_SAS_Results_1985_2011/
>> wah2_sas_1985_2007.part.Bdesh.ymonmean.nc", "r");
>> f =
>> addfile("/ouce-home/students/rege1218/Observations/APHRO_Asian_Monsoon_1963-2007/APHRO_Bdesh_ppt_1985_2007_ymonmean.nc",
>> "r")
>> ;*****************************************************
>> ; parameters
>> ;*****************************************************
>>
>> PPT_mod = f0->field90(:,0,0,0)
>> PPT_mod = PPT_mod*86400
>>
>> PPT_obs = f->precip
>> time = f->time
>> x = PPT_obs&time
>> timax = dimsizes(time)-1
>>
>> fldmean = wgt_areaave_Wrap(PPT_obs,1.0,1.0,1)
>>
>> ;--convert the time proleptic_gregorian calendar to UTC date
>> utc_date   = cd_calendar(time,0)
>>
>> ;-- set date variable names
>>
>> year     = tointeger(utc_date(:,0))
>> month     = tointeger(utc_date(:,1))
>> day     = tointeger(utc_date(:,2))
>> hour      = tointeger(utc_date(:,3))
>> minute      = tointeger(utc_date(:,4))
>> second      = utc_date(:,5)
>>
>> ;-- write date as string (MM)
>> date_str_i = sprinti("%0.2i",month)
>>
>> ;-- create the time strings, plot every axis annotation
>> incr     = 1
>> labels      = (/date_str_i(0:incr)/)
>>
>> ;-- set t resources
>> res = True
>> res at trXMinF           = time(0)     ; time minimum on axis
>> res at trXMaxF           = time(timax)     ; time maximum on axis
>> res at tmXBMode         = "Explicit"           ; explicit time setting
>> res at tmXBValues         = PPT_obs&time(:incr)          ; axis tick
>> position
>> res at tmXBLabels         = labels
>>
>>
>> colors = (/"blue","violet"/)
>> labels = (/"RCM_ppt","OBS_ppt"/)
>>
>> wks = gsn_open_wks("eps","annpptcycle")
>>
>> res =True
>> res at gsnDraw    = False
>> res at gsnFrame   = False
>> res at gsnPaperOrientation = "portrait"
>>
>> res= True   ; plot mods desired
>> res =True
>> res at gsnMaximize= True
>> res at gsnDraw    = False
>> res at gsnFrame   = False
>> res at gsnPaperOrientation = "portrait"
>>
>>
>> res= True   ; plot mods desired
>> res at tiMainString = "Bangladesh with Larger Region"
>> res at tiYAxisString   = "mm/day"
>> res at xyLineThicknessF = (/3,3/)
>> res at trYMinF                 = 0.0
>> res at trYMaxF                 = 20.0
>>
>>
>> res at xyLineColor          = colors(0)
>> plot0 = gsn_csm_xy(wks,x,PPT_mod,res)
>>
>> res at xyLineColor          = colors(1)
>> plot1 = gsn_csm_xy(wks,x,fldmean,res)
>>
>> overlay(plot0,plot1)
>>
>>
>> ; Attach a legend
>>
>> lgres = True
>> lgres at lgLineColors              = colors
>> lgres at lgItemType        = "Lines"
>> lgres at lgLabelFontHeightF     = 0.2
>> lgres at lgBoxMinorExtentF      = 0.2
>> lgres at vpWidthF            = 0.2
>> lgres at vpHeightF         = 0.2
>> lgres at pmLegendDisplayMode     = "Always"
>> lgres at pmLegendWidthF          = 0.05
>> lgres at pmLegendHeightF          = 0.04
>> lgres at xyExplicitLabels         = labels
>> lgres at lgBoxMinorExtentF     = 0.3
>> legend = gsn_create_legend(wks, 2, labels, lgres)
>>
>> amres                 = True
>> amres at amJust             = "BottomRight"     ; Use bottom right corner
>> of box for determining its location.
>> amres at amOrthogonalPosF        = -0.22
>> amres at amParallelPosF         = 0.4
>>
>> annoid  =  gsn_add_annotation(plot0, legend, amres)  ; add legend to plot
>>
>> draw(plot0)
>>
>> frame(wks)
>>
>> end
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151127/1796390e/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: LRB_annpptcycle(1).jpg
Type: image/jpeg
Size: 355781 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151127/1796390e/attachment-0001.jpg