band_pass_hovmueller_plot

Generate a plot using the output from band_pass_hovmueller.

Available in version 5.1.0 and later.

Prototype

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"        ; These three libraries are automatically
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"         ; loaded from NCL V6.2.0 onward.
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"     ; No need for user to explicitly load.

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl" ; Still need to load this library.

	procedure band_pass_hovmueller_plot (
		x    [*][*] : numeric,  
		pltDir  [1] : string,   
		pltType [1] : string,   
		pltName [1] : string,   
		opt     [1] : logical   
	)

Arguments

The variable returned from band_pass_hovmueller. It has dimensions (time,lon).

The directory to which the output will be directed. diro="./" means the current directory.

Graphical type to be generated: "png", "x11", "ps", "eps", "pdf", "ncgm".

The name associated with the generated graphic.

If opt=False, the function will operate under default mode regardless of any attributes associated with the variable.

If opt=True, then attributes associated with opt may be used to alter the default behavior. Currently, minimal options are available:

• opt@yearFraction=True will result in the left axis using year and fraction-of-year rather than the textual month and year. Further, if this option is set to True, the opt@yearFractionSpacingF can be used to specify the tick mark spacing (default=0.25).

Description

Using the variable returned from band_pass_hovmueller, and its attributes, generate a Hovmueller plot. Options for the style of the time axis are provided.

See Also

band_pass_hovmueller_cam (not yet developed),
band_pass_hovmueller,
band_pass_area_time,
band_pass_latlon_time,
filwgts_lanczos,
wgt_areaave,
wgt_areaave_Wrap,
wkSpaceTime,
wkSpaceTime_cam

Examples

Example 1

Consider daily "olr". Extract MJO (Madden-Julian Oscillation) information using the band pass 20-to-100 period and 201 weights as suggested by the WMO CLIVAR Working Group.

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/diagnostics_cam.ncl"


  f     = addfile("...", "r") 
  x     = f->OLR                 ; (time,lat,lon) , ntim=1461

  latS  = -10.              ; select region
  latN  =  10.
  lonL  =  60.              ; Indian Ocean
  lonR  = 280.              ; East Pacific

  spd   = 1                 ; daily means
  bpf   = (/20,100,201/)    ; MJO ... suggested by MJO WG [WMO] 
  wgty  = latRegWgt(x&lat, "double", 0) 
  opt      = False
                            ; pass in area of interest
  mjo_area = band_pass_hovmueller (x(:,{latS:latN},{lonL:lonR} \
                                ,spd, bpf, wgty({latS:latN}), opt)

  printVarSummary( mjo_area )
  printMinMax( mjo_area, True)

The edited output (currently) looks like:

       Variable: mjo_area
       Number of Dimensions: 1
       Dimensions and sizes:  [var | 3] x [time | 1461]
       Coordinates:
                  time: [17479080..17514120]
       Number Of Attributes: 18
        units :       W/m2
        long_name :   filtered: weighted area avg: Daily OLR
        band_pass_start :     20
        band_pass_last :      100
        band_pass_Nwgts :     201
        var_0:                band pass
        var_1:                raw areal means
        var_2:                local variances
       
        filtered: weighted area avg: Daily OLR: min=-6.58817   max=5.37212