[ncl-talk] WRF Time labels - repost
Rashed Mahmood
rashidcomsis at gmail.com
Sun Oct 6 01:31:28 MDT 2019
Hi Prasad,
>From the figure you had attached, it seems that the NCL is plotting 3
hourly data correctly. For example, count the minor tic mark lines, i.e
small vertical grey lines on x-axis. The problem is that it is not showing
results to the last value of your data.
Try setting the following two resources for both/one of your resources set,
i.e. restime and restic:
restime at trXMinF = 0 ; OR min(Time_0)
restime at trXMaxF = 72; OR 25; OR max(Time_0)
restic at trXMinF = 0 ; OR min(Time_0)
restic at trXMaxF = 72; OR 25; OR max(Time_0)
PS: I have not used WRF plotting routines, so I am not quite sure if these
would work. Anyway, worth a try.
Cheers,
Rashed
On Sat, Oct 5, 2019 at 2:43 AM Prasad Kunjeer via ncl-talk <
ncl-talk at ucar.edu> wrote:
> Hi Adam,
>
> I have modified the script as per your instructions. The output of the
> PrintVarSummary (ua) is given below ;
>
> Variable: ua (subsection)
> Type: float
> Total Size: 3200 bytes
> 800 values
> Number of Dimensions: 2
> Dimensions and sizes: [Time | 25] x [bottom_top | 32]
> Coordinates:
> * Time: [0..72]*
> Number Of Attributes: 6
> FieldType : 104
> MemoryOrder : XYZ
> description : x-wind component
> units : m s-1
> stagger :
> coordinates : XLONG XLAT
>
> But still, there is no change in the graph.
>
> Regards
> प्रसाद कुंजीर/ Prasad Kunjeer
> वैज्ञानिक - सी / Scientist - C
> केन्द्रीय जल और विद्युत अनुसंधान शाला / Central Water and Power Research
> Station
> पुणे 411 024/ Pune 411 024
>
>
> On Sat, Oct 5, 2019 at 3:59 AM Adam Phillips <asphilli at ucar.edu> wrote:
>
>> Hi Prasad,
>> Your ur array does not have a time coordinate variable as shown here:
>> Type: float
>> Total Size: 3200 bytes
>> 800 values
>> Number of Dimensions: 2
>> Dimensions and sizes: [Time | 25] x [bottom_top | 32]
>> *Coordinates:*
>> Number Of Attributes: 6
>>
>> Thus, NCL defaults to setting your times to go from 0 to 24 (as you have
>> 25 timesteps). This is why you are only seeing two labels at tmXBValues 0
>> and 18, as the range that is shown on the plot goes from 0->24. To fix
>> this, I would recommend assigning a coordinate variable to ua as follows:
>> ua&Time = ispan(0,72,3)
>>
>> That way the Times associated with ua run from 0->72, and your tick mark
>> labels cover the same span.
>>
>> If you have any further questions please respond to the ncl-talk email
>> list.
>> Adam
>>
>>
>> On Thu, Oct 3, 2019 at 3:17 AM Prasad Kunjeer <pskunjeer at gmail.com>
>> wrote:
>>
>>> SIr,
>>>
>>> The outputs of:
>>>
>>>
>>> *printVarSummary(restime) is*
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> *Variable: restimeType: logicalTotal Size: 4 bytes 1
>>> valuesNumber of Dimensions: 1Dimensions and sizes: [1]Coordinates: Number
>>> Of Attributes: 4 tmXBMode : Explicit tmXBValues : ( 0, 18, 36,
>>> 54, 72 ) tmXBMinorValues : <ARRAY of 25 elements> tmXBLabels : ( 09/22
>>> 12:00, 09/23 06:00, 09/24 00:00, 09/24 18:00, 09/25 12:00 )*
>>>
>>>
>>> *printVarSummary(ua(:,:,100,100)) is *
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> *Type: floatTotal Size: 3200 bytes 800 valuesNumber of
>>> Dimensions: 2Dimensions and sizes: [Time | 25] x [bottom_top |
>>> 32]Coordinates: Number Of Attributes: 6 FieldType : 104 MemoryOrder :
>>> XYZ description : x-wind component units : m s-1 stagger : coordinates
>>> : XLONG XLAT*
>>>
>>> In restime it is showing labels from 09/22 12:00, 09/23 06:00, 09/24
>>> 00:00, 09/24 18:00, 09/25 12:00 but during plotting it is not showing the
>>> same.
>>>
>>> Regards
>>> प्रसाद कुंजीर/ Prasad Kunjeer
>>> वैज्ञानिक - सी / Scientist - C
>>> केन्द्रीय जल और विद्युत अनुसंधान शाला / Central Water and Power
>>> Research Station
>>> पुणे 411 024/ Pune 411 024
>>>
>>>
>>> On Wed, Oct 2, 2019 at 11:09 PM Adam Phillips <asphilli at ucar.edu> wrote:
>>>
>>>> Hi Prasad,
>>>> Can you send ncl-talk the output from the following printVarSummary
>>>> call:
>>>> time_axis_labels(Time_0,restime,restick)
>>>> *printVarSummary(restime)*
>>>> *printVarSummary(ua(:,:,100,100))*
>>>> contour_vor = wrf_contour(sfile,wks,ua(:,:,100,100),restime)
>>>> Adam
>>>>
>>>>
>>>> On Wed, Oct 2, 2019 at 9:43 AM Prasad Kunjeer via ncl-talk <
>>>> ncl-talk at ucar.edu> wrote:
>>>>
>>>>> Sorry I just forgot to attach the file.
>>>>>
>>>>> On Wed, Oct 2, 2019, 9:04 PM Prasad Kunjeer <pskunjeer at gmail.com>
>>>>> wrote:
>>>>>
>>>>>> 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.
>>>>>>
>>>>>> The time interval for WRF file is 3 hour and there are 25 time steps.
>>>>>> The first time steps is 22/09/2018 1200. The time should start from
>>>>>> 22/09/2018 and end at 25/09/2018 with 25 time steps (22/09/2018 1200,
>>>>>> 22/09/2018 1500 etc upto 25/09/2018 1200). But the graph ends at
>>>>>> 23/09/2018 1200 which correspond to 25 time steps at 1 hour interval. I
>>>>>> want graph to end at 25/09/2018 1200.
>>>>>>
>>>>>> Please guide in correcting the output of graph
>>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at ucar.edu
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>>
>>>> --
>>>> Adam Phillips
>>>> Associate Scientist, Climate and Global Dynamics Laboratory, NCAR
>>>> www.cgd.ucar.edu/staff/asphilli/ 303-497-1726
>>>>
>>>> <http://www.cgd.ucar.edu/staff/asphilli>
>>>>
>>>
>>
>> --
>> Adam Phillips
>> Associate Scientist, Climate and Global Dynamics Laboratory, NCAR
>> www.cgd.ucar.edu/staff/asphilli/ 303-497-1726
>>
>> <http://www.cgd.ucar.edu/staff/asphilli>
>>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191006/563bb359/attachment.html>
More information about the ncl-talk
mailing list