[ncl-talk] Time Labels -
Prasad Kunjeer
pskunjeer at gmail.com
Thu Sep 26 09:18:33 MDT 2019
Sir,
I am trying to plot a simple graph from a WRF output. The X axis of the
plot is time and i want to print the time table for the same. I am using
the function 'time_axis_labels'. I have attached the script:
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/time_axis_labels.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
;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"
; This file still has to be loaded manually
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/time_axis_labels.ncl"
srcDirName = "/home/prasad/res/"
srcFileName = "newdo3.nc"
srcFilePath = srcDirName + srcFileName
;---Retrieve either one level, or all levels. Use '-1' for all.
sfile = addfile(srcFilePath,"r")
ua = wrf_user_getvar(sfile,"ua",-1) ; On mass grid
va = wrf_user_getvar(sfile,"va",-1)
printVarSummary(ua) ;
(Time,bottom_top,south_north,west_east)
printVarSummary(va) ;
(Time,bottom_top,south_north,west_east)
Times = sfile->Times
Time_0 = wrf_times_c(Times, 0)
printVarSummary(Times)
printVarSummary(Time_0)
print(Time_0)
wks = gsn_open_wks("png","time_labels") ; send graphics to PNG file
;--------------------------------------------------
; The time_axis_label function adds additional
; resources to "res" to produce nicely-formatted
; time labels on X axis. This function only works
; if you have a time "units" recognized by the
; cd_calendar function.
;--------------------------------------------------
pltres = True
restime = True
restick = True
restick at ttmFormat = "%N/%D %H:%M"
time_axis_labels(Time_0,restime,restick)
contour_vor = wrf_contour(sfile,wks,ua(:,:,100,100),restime)
plot = wrf_overlays(sfile,wks,(/contour_vor/),pltres)
The WRF file contains the 25 time steps at 3 hour interval. I have attached
the output of printvarSummary and print of the time varibale.
Variable: Time_0
Type: double
Total Size: 200 bytes
25 values
Number of Dimensions: 1
Dimensions and sizes: [Time | 25]
Coordinates:
Time: [ 0.. 72]
Number Of Attributes: 4
units : hours since 2018-09-22 12:00:00
calendar : standard
long_name : Time
description : Time
Variable: Time_0
Type: double
Total Size: 200 bytes
25 values
Number of Dimensions: 1
Dimensions and sizes: [Time | 25]
Coordinates:
Time: [ 0.. 72]
Number Of Attributes: 4
units : hours since 2018-09-22 12:00:00
calendar : standard
long_name : Time
description : Time
(0) 0
(1) 3
(2) 6
(3) 9
(4) 12 etc.
But when I plot the figure, the output takes 1 hour interval rather than 3
hour interval.
Please guide in correcting the output of graph
Regards
प्रसाद कुंजीर / Prasad Kunjeer
वैज्ञानिक - सी / Scientist - C
केन्द्रीय जल और विद्युत अनुसंधान शाला / Central Water and Power Research
Station
पुणे 411 024/ Pune 411 024
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190926/a20a3356/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: time_label.png
Type: image/png
Size: 150924 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190926/a20a3356/attachment.png>
More information about the ncl-talk
mailing list