Module: NumRu::GGraph

Defined in:

lib/numru/ggraph.rb,
lib/numru/ganalysis/histogram.rb,
lib/numru/ggraph_on_merdional_section.rb

Constant Summary

TONE_TONF_THRES =

10000

@@nannot =

thresfold to swith uetonf and uetone under the auto=true option in the tone method

0

@@fig =

Misc::KeywordOptAutoHelp.new(
  ['new_frame',true,      'whether to define a new frame by DCL.grfrm (otherwise, DCL.grfig is called)'],
  ['no_new_fig',false,      'If true, neither DCL.grfrm nor DCL.grfig is called (overrides new_frame) -- Then, you need to call one of them in advance. Convenient to set DCL parameters that are reset by grfrm or grfig.'],
  ['itr',     1,     'coordinate transformation number'],
  ['viewport',[0.2,0.8,0.2,0.8], '[vxmin, vxmax, vymin, vymax]'],
  ['eqdistvpt',false, 'modify viewport to equidistant for x and y (only for itr=1--4)'],
  ['window',  nil, '(for itr<10,>50) [uxmin, uxmax, uymin, uymax]. each element allowd nil (only for itr<5,>50)'],
  ['xreverse','positive:down,units:hPa', '(for itr<10,>50) Assign max value to UXMIN and min value to UXMAX if condition is satisfied (nil:never, true:always, String: when an attibute has the value specified ("key:value,key:value,..")'],
  ['yreverse','positive:down,units:hPa', '(for itr<10,>50) Assign max value to UYMIN and min value to UYMAX if condition is satisfied (nil:never, true:always, String: when an attibute has the value specified ("key:value,key:value,..")'],
  ['round0',     false,     'expand window range to good numbers (effective only to internal window settings)'],
  ['round1',     false,     'expand window range to good numbers (effective even when "window" is explicitly specified)'],
  # for 5<=itr<=7 (Rectangular Curvilinear Coordinates)
  ['similar', nil, '3-element float array for similar transformation in a rectangular curvilinear coordinate, which is fed in DCL:grssim:[simfac,vxoff,vyoff],  where simfac and [vxoff,vyoff] represent scaling factor and origin shift, respectively.'],
  # for 10<=itr<=50 (map projection):
  ['map_axis', nil, '(for all map projections) 3-element float array to be fed in DCL::umscnt: [uxc, uxy, rot], where [uxc, uyc] represents the tangential point (or the pole at top side for cylindrical projections), and rot represents the rotation angle. If nil, internally determined. (units: degrees)'],
  ['map_radius', nil, '(for itr>=20: conical/azimuhal map projections) raidus around the tangential point. (units: degrees)'],
  ['map_fit',nil,'(Only for itr=10(cylindrical) and 11 (Mercator)) true: fit the plot to the data window (overrides map_window and map_axis); false: do not fit (then map_window and map_axis are used); nil: true if itr==10, false if itr==11'],
  ['map_rsat',nil,'(Only for itr=30) satellite distance from the earth\'s center (Parameter "RSAT" for sgpack)'],
  ['map_window', [-180,180,-75,75], '(for itr<20: cylindrical map projections) lon-lat window [lon_min, lon_max, lat_min, lat_max ] to draw the map (units: degres)']
)

@@next_fig =

nil

@@axes =

Misc::KeywordOptAutoHelp.new(
  ['xside',  'tb',  'Where to draw xaxes (combination of t, b and u)'],
  ['yside',  'lr',  'Where to draw yaxes (combination of l, r and u)'],
  ['xtitle',  nil,  'Title of x axis (if nil, internally determined)'],
  ['ytitle',  nil,  'Title of y axis (if nil, internally determined)'],
  ['title',   nil,  'Title of the figure'],
  ['xunits',  nil,  'Units of x axis (if nil, internally determined)'],
  ['yunits',  nil,  'Units of y axis (if nil, internally determined)'],
  ['xtickint',  nil,
        'Interval of x axis tickmark (if nil, internally determined)'],
  ['ytickint',  nil,
        'Interval of y axis tickmark (if nil, internally determined)'],
  ['xlabelint',  nil,
           'Interval of x axis label (if nil, internally determined)'],
  ['ylabelint',  nil,
           'Interval of y axis label (if nil, internally determined)'],
  ['xloglabelall',  false,
           'Show lavels for all log-level tick marks (x-axes) (e.g.,1000,900,800,... inseatd of 1000,500,200,...)'],
  ['yloglabelall',  false,
           'Show lavels for all log-level tick marks (y-axes) (e.g.,1000,900,800,... inseatd of 1000,500,200,...)'],
  ['xmaplabel',  nil,
           'If "lon"("lat"), use DCLExt::lon_ax(DCLExt::lat_ax) to draw xaxes; otherwise, DCL::usxaxs is used.'],
  ['ymaplabel',  nil,
           'If "lon"("lat"), use DCLExt::lon_ax(DCLExt::lat_ax) to draw yaxes; otherwise, DCL::usyaxs is used.'],
  ['time_ax',  nil,  'Type of calendar-type time axis: nil (=> auto slection), false (do not use the time axis even if the units of the axis is a time one with since field), "h" (=> like nil, but always use the hour-resolving datetime_ax method in dclext_datetime_ax.rb), or "ymd" (=> like "h" but for y-m-d type using DCL.uc[xy]acl)']
)

@@next_axes =

nil

@@map =

Misc::KeywordOptAutoHelp.new(
  ['lim',          true,  'draw map lim (t or f)'],
  ['grid',         true,  'draw map grid (t or f)'],
  ['vpt_boundary', false,  'draw viewport boundaries (f, t or 1,2,3.., representing the line width)'],
  ['wwd_boundary', false,  'draw worksation window boundaries (f, t or 1,2,3.., representing the line width)'],
  ['fill',  false,  'fill the map if coast_world or coast_japan is true (t or f)'],
  ['coast_world',  false,  'draw world coast lines (t or f)'],
  ['border_world', false,  'draw nation borders (t or f)'],
  ['plate_world',  false,  'draw plate boundaries (t or f)'],
  ['state_usa',  false,    'draw state boundaries of US (t or f)'],
  ['coast_japan',  false,  'draw japanese coast lines (t or f)'],
  ['pref_japan',   false,  'draw japanese prefecture boundaries (t or f)'],
  ['dgridmj', nil, 'the interval between the major lines of latitudes and longitudes. If nil, internally determined. (units: degrees) (this is a UMPACK parameter, which is nullified when uminit or grfrm is called)'],
  ['dgridmn', nil, 'the interval between the minor lines of latitudes and longitudes. If nil, internally determined. (units: degrees) (this is a UMPACK parameter, which is nullified when uminit or grfrm is called)']
)

@@next_map =

nil

@@data_prep_options =

Misc::KeywordOptAutoHelp.new(
  ['exchange', false, 'whether to exchange x and y axes'],
  ['transpose', false, 'if true, exchange x and y axes'],
  ['xintv', 1, 'interval of data sampling in x'],
  ['yintv', 1, 'interval of data sampling in y'],
  ['xcoord', nil, 'Name of the coordinate variable for x-axis'],
  ['ycoord', nil, 'Name of the coordinate variable for y-axis'],
  ['slice', nil, 'An Array to be pathed to the GPhys#[] method to subset the data before plotting (order applied: slice -> cut -> mean)'],
  ['cut', nil, 'An Array or Hash to be pathed to the GPhys#cut method to subset the data before plotting (order applied: slice -> cut -> mean)'],
  ['mean', nil, 'An Array to be pathed to the GPhys#mean method to take mean of the data before plotting (order applied: slice -> cut -> mean)']
)

@@line_options =

Misc::KeywordOptAutoHelp.new(
  ['title', nil, 'Title of the figure(if nil, internally determined)'],
  ['annotate', true, 'if false, do not put texts on the right margin even when newframe==true'],
  ['index', 1, 'line/mark index'],
  ['type', 1, 'line type'],
  ['label', nil, 'if a String is given, it is shown as the label'],
  ['max', nil, 'maximam data value'],
  ['min', nil, 'minimam data value'],
  ['legend', nil, 'legend to annotate the line type and index. nil (defalut -- do not show); a String as the legend; true to use the name of the GPhys as the legend'],
  ['legend_vx', nil, '(effective if legend) viewport x values of the lhs of the legend line (positive float); or nil for automatic settting (shown to the right of vpt); or move it to the left relatively (negtive float)'],
  ['legend_dx', nil, '(effective if legend) length of the legend line'],
  ['legend_vy', nil, '(effective if legend) viewport y value of the legend (Float; or nil for automatic settting)'],
  ['legend_size', nil, '(effective if legend) character size of the legend'],
  ['map_axes', false, '[USE IT ONLY WHEN itr=10 (cylindrical)] If true, draws axes by temprarilly switching to itr=1 and calling GGraph::axes.'],
  @@data_prep_options
)

@@mark_options =

---

Misc::KeywordOptAutoHelp.new(
      ['title', nil, 'Title of the figure(if nil, internally determined)'],
      ['annotate', true, 'if false, do not put texts on the right margin even when newframe==true'],
      ['index', 1, 'mark index'],
      ['type', 2, 'mark type'],
      ['size', 0.01, 'marks size'],
      ['max', nil, 'maximam data value'],
      ['min', nil, 'minimam data value'],
      ['legend', nil, 'legend to annotate the mark type, index, and size. nil (defalut -- do not to show); a String as the legend; true to use the name of the GPhys as the legend'],
      ['legend_vx', nil, '(effective if legend) viewport x values of the lhs of the legend line (positive float); or nil for automatic settting (shown to the right of vpt); or move it to the left relatively (negtive float)'],
      ['legend_vy', nil, '(effective if legend) viewport y value of the legend (Float; or nil for automatic settting)'],
      ['legend_size', nil, '(effective if legend) character size of the legend'],
      ['map_axes', false, '[USE IT ONLY WHEN itr=10 (cylindrical)] If true, draws axes by temprarilly switching to itr=1 and calling GGraph::axes.'],
      @@data_prep_options
)

@@linear_contour_options =

Misc::KeywordOptAutoHelp.new(
  ['min',      nil,       'minimum contour level'],
  ['max',      nil,       'maximum contour level'],
  ['nlev',     nil,       'number of levels'],
  ['interval', nil,       'contour interval'],
  ['nozero',   nil,       'delete zero contour'],
  ['coloring', false,     'set color contours with ud_coloring'],
  ['clr_min',  13,        '(if coloring) minimum color number for the minimum data values'],
  ['clr_max',  99,       '(if coloring) maximum color number for the maximum data values']
)

@@next_linear_contour_options =

nil

@@contour_levels =

Misc::KeywordOptAutoHelp.new(
  ['log',   nil,   'approximately log-scaled levels (by using DCLExt::quasi_log_levels)'],
  ['log_cycle',  3,   '(if log) number of levels in one-order (1 or 2 or 3)'],
  ['levels',   nil,   'contour levels (Array/NArray of Numeric)'],
  ['index',    nil,   '(if levels) line index(es) (Array/NArray of integers, Integer, or nil)'],
  ['line_type',nil,   '(if levels) line type(s) (Array/NArray of integers, Integer, or nil)'],
  ['label',    nil,   '(if levels) contour label(s) (Array/NArray of String, String, true, false, nil). nil is recommended.'],
  ['label_height',nil,'(if levels) label height(s) (Array/NArray of Numeric, Numeric, or nil). nil is recommended.']
)

@@set_contour_levels =

false

@@next_contour =

nil

@@contour_options =

Misc::KeywordOptAutoHelp.new(
  ['title', nil, 'Title of the figure(if nil, internally determined)'],
  ['annotate', true, 'if false, do not put texts on the right margin even when newframe==true'],
  ['map_axes', false, '[USE IT ONLY WHEN itr=10 (cylindrical)] If true, draws axes by temprarilly switching to itr=1 and calling GGraph::axes.'],
  ['keep', false, 'Use the contour levels used previously'],
  @@linear_contour_options,
  @@contour_levels,
  @@data_prep_options
)

@@linear_tone_options =

Misc::KeywordOptAutoHelp.new(
  ['min',      nil,       'minimum tone level'],
  ['max',      nil,       'maximum tone level'],
  ['nlev',     nil,       'number of levels'],
  ['interval', nil,       'contour interval']
)

@@next_linear_tone_options =

nil

@@tone_levels =

Misc::KeywordOptAutoHelp.new(
  ['log',   nil,   'approximately log-scaled levels (by using DCLExt::quasi_log_levels)'],
  ['log_cycle',  3,   '(if log) number of levels in one-order (1 or 2 or 3)'],
  ['levels',   nil,   'tone levels  (Array/NArray of Numeric). Works together with patterns'],
  ['patterns', nil,   'tone patters (Array/NArray of Numeric). Works together with levels']
)

@@set_tone_levels =

false

@@next_tone =

nil

@@tone_options =

Misc::KeywordOptAutoHelp.new(
  ['title', nil, 'Title of the figure(if nil, internally determined)'],
  ['annotate', true, 'if false, do not put texts on the right margin even when newframe==true'],
  ['ltone', true, 'Same as udpack parameter ltone'],
  ['auto', false, 'Swith DCL.uetone and DCL.uetonf depending on the data size'],
  ['tonf', false, 'Use DCL.uetonf instead of DCL.uetone'],
  ['tonb', false, 'Use DCL.uetonb instead of DCL.uetone'],
  ['tonc', false, 'Use DCL.uetonc instead of DCL.uetone'],
  ['clr_min', nil,  'if an integer (in 10..99) is specified, used as the color number for the minimum data values. (the same can be done by setting the uepack parameter "icolor1")'],
  ['clr_max', nil,  'if an integer (in 10..99) is specified, used as the color number for the maximum data values. (the same can be done by setting the uepack parameter "icolor2")'],
  ['map_axes', false, '[USE IT ONLY WHEN itr=10 (cylindrical)] If true, draws axes by temprarilly switching to itr=1 and calling GGraph::axes.'],
  ['keep', false, 'Use the tone levels and patterns used previously'],
  ['color_bar', false, 'Add a color bar: THIS IS ONLY FOR QUICK LOOK. Use the GGraph::color_bar method explicitly for full option control'],
  @@linear_tone_options,
  @@tone_levels,
  @@data_prep_options
)

@@scatter_options =

< scatter diagram >

Misc::KeywordOptAutoHelp.new(
      ['title', "", 'Title of the figure(if nil, internally determined)'],
      ['annotate', true, 'if false, do not put texts on the right margin even when newframe==true'],
      ['index', 1, 'mark index'],
      ['type', 2, 'mark type'],
      ['size', 0.01, 'marks size'],
      ['map_axes', false, '[USE IT ONLY WHEN itr=10 (cylindrical)] If true, draws axes by temprarilly switching to itr=1 and calling GGraph::axes.'],
      ['correlation', false, 'calculate Pearson\'s correlation coefficient and show it on the bottom margin'],
      @@data_prep_options
)

@@next_scatter =

nil

@@color_scatter_options =

< scatter diagram colored by values >

Misc::KeywordOptAutoHelp.new(
      ['title', nil, 'Title of the figure(if nil, internally determined)'],
      ['annotate', true, 'if false, do not put texts on the right margin even when newframe==true'],
      ['index', 3, 'mark index (1-9)'],
      ['type', 10, 'mark type'],
      ['size', 0.01, 'marks size'],
      ['map_axes', false, '[USE IT ONLY WHEN itr=10 (cylindrical)] If true, draws axes by temprarilly switching to itr=1 and calling GGraph::axes.'],
      ['clr_min', nil,  'if an integer (in 10..99) is specified, used as the color number for the minimum data values. (the same can be done by setting the uepack parameter "icolor1")'],
      ['clr_max', nil,  'if an integer (in 10..99) is specified, used as the color number for the maximum data values. (the same can be done by setting the uepack parameter "icolor2")'],
      ['keep', false, 'Use the tone levels and patterns used previously'],
      ['correlation', false, 'calculate Pearson\'s correlation coefficient and show it on the bottom margin'],
      @@linear_tone_options,
      @@tone_levels,
      @@data_prep_options
)

@@next_color_scatter =

nil

@@add_mode_vectors_options =

Misc::KeywordOptAutoHelp.new(
      ['lineindex', 1, 'line index'],
      ['linetype', 1, 'line type'],
      ['fact', 2, 'scaling factor (line length = stddev * fact for each side)'],
      ['style', 'line', 'style of displaying modes (line, arrow, ellipse)']
)

@@vxfxratio =

for flow_vect

nil

@@vyfyratio =

for flow_vect

nil

@@vfratio =

for flow_vect_anyproj

nil

@@vector_options =

Misc::KeywordOptAutoHelp.new(
  ['title', nil, 'Title of the figure(if nil, internally determined)'],
  ['annotate', true, 'if false, do not put texts on the right margin even when newframe==true'],
  ['map_axes', false, '[USE IT ONLY WHEN itr=10 (cylindrical)] If true, draws axes by temprarilly switching to itr=1 and calling GGraph::axes.'],
  ['flow_vect', true, 'If true, use DCLExt::flow_vect to draw vectors; otherwise, DCL::ugvect is used.'],
  ['flow_vect_anyproj', nil, 'Whether to use flow_vect_anyproj. If nil, up to the current projection number (when >=2); if true, always; if false, never. (precedence if higher than the flow_vect parameter)'],
  ['keep', false, 'Use the same vector scaling as in the previous call. -- Currently, works only when "flow_vect" is true'],
  ['factor', 1.0, '(Effective only if flow_vect) scaling factor to strech/reduce the arrow lengths'],
  ['unit_vect', false, 'Show the unit vector'],
  ['max_unit_vect', false, '(Effective only if flow_vect && unit_vect) If true, use the maximum arrows to scale the unit vector; otherwise, normalize in V coordinate.'],
  ['ux_unit', nil, '(If Numeric) length of the x direction unit vector (precedence of this option is lower than max_unit_vect)'],
  ['uy_unit', nil, '(If Numeric) length of the y direction unit vector (precedence of this option is lower than max_unit_vect)'],
  ['len_unit', nil, '(Effective only when flow_vect_anyproj is used) If Numeric, length of the the unit vector in terms of (fx,fy); if nil (defualt), the unitvecto length is set to be the maximum one'],
  ['flow_itr5', false, 'If true, use DclExt::flow_itr5 to draw vectors on 2-dim polar coordinate. Should be set DCL.sgstrn(5)'],
  ['polar_thinning', nil, '(Effective only when flow_vect_anyproj is used) If Numeric, specifies the maximum grid isotrpy, with which thinning is made near the poles (recommended value: 0.3 ~ 1)'],
  ['distvect_map', true, '(effective only for flow_vect_anyproj when map projection) by default (true) it is assumed that the vector (fx,fy) is based on lengths (such as wind velocities in m/s and fluxes in which wind velocities are incorporated [q*u, q*v]); set it to false if the vector is based on angles (such as the time derivatives of longitude and latitude).'],
   @@data_prep_options
)

@@regression_line_options =

Misc::KeywordOptAutoHelp.new(
  ['x_given_y', false, 'If false (default), regress y for given x. If true, x given y.'],
  ['annot_slope', nil, 'Whether to show the slope on the right margin. nil: show when itr=1,4, do not show when itr=2,3; true: show; false: not show'],
  ['annot_intercept', false, '[Meaningful only when itr=1 (linear scaling plot)] if true, show the intercept on the right margin'],
  ['limit', false, 'If true, regression line is shown only over the range between max and min of the independent variable. If false, regression line is written over the whole range of the independent variable'],
  ['clip', true, "If true, use DCL's clipping not to show the line outside the viewport"],
  ['index', 1, 'line index'],
  ['type', 1, 'line type']
)

@@next_regression_line =

nil

@@histogram_options =

Misc::KeywordOptAutoHelp.new(
            ['window', [nil,nil,nil,nil], "window bounds"],
            ['title', "histogram", "window title"],
            ['exchange', false, "exchange x and y"],
            ['fill', false, "fill bars"]
)

Class Method Summary

• .__calendar_ax(axtype, xax, side, sunits, calendar, ttl, tickint = nil, labelint = nil) ⇒ Object

  axtype: nil (–>auto), ‘h’ (to resolve hours), or ‘ymd’.

• .__good_range(x0, x1) ⇒ Object
• .__round_window(window) ⇒ Object
• .__shorten_path(str, maxlen) ⇒ Object

  private methods in the module.

• ._get_axinfo(vary) ⇒ Object
• ._scatter_display_corr(gpx, gpy, opts) ⇒ Object
• ._set_contour_levels_(opts) ⇒ Object
• ._set_tone_levels_(opts) ⇒ Object
• .add_mode_vectors(mean, modes, options = nil) ⇒ Object
• .annotate(str_ary, noff = nil) ⇒ Object
• .annotate_at(str_ary, vx, vy, charsize = nil, align = -1.0,, dvyfact = -1.5,, noff = nil) ⇒ Object
• .axes(xax = nil, yax = nil, options = nil) ⇒ Object
• .axes_or_map_and_ttl(gp, opts, xax, yax, gp2 = nil) ⇒ Object
• .clear_contour_levels ⇒ Object
• .clear_tone_levels ⇒ Object
• .color_bar(*args) ⇒ Object
• .color_scatter(fx, fy, fz, newframe = true, options = nil) ⇒ Object
• .contour(gphys, newframe = true, options = nil) ⇒ Object
• .current_itr(fig_yet_to_be_called = false) ⇒ Object
• .data_prep_1D(gphys, opts) ⇒ Object
• .data_prep_1D2(fx, fy, opts) ⇒ Object
• .data_prep_2D(gps, newframe, opts) ⇒ Object
• .data_prep_multidim(gphys, opts) ⇒ Object
• .fig(xax = nil, yax = nil, options = nil) ⇒ Object
• .gropn_1_if_not_yet ⇒ Object
• .histogram(gphys, newframe = true, options = nil) ⇒ Object (also: histogram1D)
• .itr_is?(itr, fig_yet_to_be_called = false) ⇒ Boolean
• .line(gphys, newframe = true, options = nil) ⇒ Object
• .map(options = nil) ⇒ Object
• .map_trn?(fig_yet_to_be_called = false) ⇒ Boolean
• .margin_info(program = nil, data_source = nil, char_height = nil, date = nil, xl = 0.0, xr = 0.0, yb = nil, yt = 0.0) ⇒ Object
• .mark(gphys, newframe = true, options = nil) ⇒ Object
• .name_selector(gp, limit) ⇒ Object
• .next_axes(options) ⇒ Object
• .next_color_scatter(options) ⇒ Object
• .next_contour(options) ⇒ Object
• .next_fig(options) ⇒ Object
• .next_linear_contour_options(options) ⇒ Object
• .next_linear_tone_options(options) ⇒ Object
• .next_map(options) ⇒ Object
• .next_regression_line(options) ⇒ Object
• .next_scatter(options) ⇒ Object
• .next_tone(options) ⇒ Object
• .next_unit_vect_options(options) ⇒ Object
• .polar_coordinate_boundaries(xax = nil, yax = nil) ⇒ Object
• .profile(gphys, newframe = true, options = nil) ⇒ Object

  ———————————————————- Added at 080207 –.

• .regression_line(x, y, options = nil) ⇒ Object
• .scatter(fx, fy, newframe = true, options = nil) ⇒ Object
• .set_axes(options) ⇒ Object
• .set_color_scatter(options) ⇒ Object

  @@color_scatter_options.delete(“transpose”) # delete is yet to be supported @@color_scatter_options.delete(“exchange”) # (2010-01).

• .set_contour(options) ⇒ Object
• .set_contour_levels(options) ⇒ Object
• .set_fig(options) ⇒ Object
• .set_linear_contour_options(options) ⇒ Object
• .set_linear_tone_options(options) ⇒ Object
• .set_map(options) ⇒ Object
• .set_regression_line(options) ⇒ Object
• .set_scatter(options) ⇒ Object

  @@scatter_options.delete(“transpose”) # delete is yet to be supported @@scatter_options.delete(“exchange”) # (2010-01).

• .set_tone(options) ⇒ Object
• .set_tone_levels(options) ⇒ Object
• .set_unit_vect_options(options) ⇒ Object

  < vector >.

• .sim_trn? ⇒ Boolean
• .timeplot(gphys, newframe = true, options = nil) ⇒ Object

  timeplot 0802011 ######.

• .title(string) ⇒ Object
• .tone(gphys, newframe = true, options = nil) ⇒ Object
• .tone_and_contour(gphys, newframe = true, options = nil) ⇒ Object
• .vector(fx, fy, newframe = true, options = nil) ⇒ Object
• .vector_on_merdional_section(fx, fy, newframe = true, options = nil) ⇒ Object

Class Method Details

.__calendar_ax(axtype, xax, side, sunits, calendar, ttl, tickint = nil, labelint = nil) ⇒ Object

axtype: nil (–>auto), ‘h’ (to resolve hours), or ‘ymd’

# File 'lib/numru/ggraph.rb', line 1547

def __calendar_ax(axtype, xax, side, sunits, calendar, ttl, \
                 tickint=nil, labelint=nil)
  window = DCL.sgqwnd
  viewport = DCL.sgqvpt
  /(.*) *since *(.*)/ =~ sunits
  if (!$1 or !$2)
    raise("Units mismatch. Requires time units that includes 'since'")
  end
  tun = Units[$1]
  dayun = Units['days']
  if xax
    t0 = window[0]
    t1 = window[1]
  else
    t0 = window[2]
    t1 = window[3]
  end
  time = UNumeric.new(t0,sunits)
  tstr = time.to_datetime(0.1,calendar)
  jd0 = tstr.strftime('%Y%m%d').to_i
  tlen = tun.convert( t1-t0, dayun )
  if !axtype
    if tlen < 5
      axtype = 'h'
    else
      axtype = 'ymd'
    end
  end

  if xax
    DCL.grswnd(0.0, tlen, window[2], window[3] )
    DCL.grstrf
    if axtype == 'h'
      opts = {'cside'=>side}
      opts['dtick1'] = tickint if tickint
      opts['dtick2'] = labelint if labelint
      DCLExt.datetime_ax(tstr, tstr+tlen, opts )
    else
      DCL.ucxacl(side,jd0,tlen)
    end
    DCL.grswnd(*window)
    DCL.grstrf
  else
    DCL.grswnd(window[0], window[1], 0.0, tlen)
    DCL.grstrf
    if axtype == 'h'
      opts = {'yax'=>true, 'cside'=>side}
      opts['dtick1'] = tickint if tickint
      opts['dtick2'] = labelint if labelint
      DCLExt.datetime_ax(tstr, tstr+tlen, opts )
    else
      DCL.ucyacl(side,jd0,tlen)
    end
    DCL.grswnd(*window)
    DCL.grstrf
  end
end

.__good_range(x0, x1) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1490

def __good_range(x0,x1)
  x0 = x0.to_f
  x1 = x1.to_f
  xw = (x1-x0).abs
  e10, lx  = (Math::log10(xw)).divmod(1.0)
  ef = 5.0 * 10**(-e10)
  [(x0*ef).floor/ef,(x1*ef).ceil/ef]
end

.__round_window(window) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1499

def __round_window(window)
  x0,x1 = __good_range(window[0],window[1])
  y0,y1 = __good_range(window[2],window[3])
  [ x0, x1, y0, y1 ]
end

.__shorten_path(str, maxlen) ⇒ Object

private methods in the module

# File 'lib/numru/ggraph.rb', line 1082

def __shorten_path(str,maxlen)
  return str if str.length <= maxlen
  astr = str.split(/\//)
  while(str.length > maxlen)
    astr.length!=2 ? at = (astr.length)/2 : at = 0
    astr.delete_at( at )
    (a = astr.dup)[at, 0] = '...'    # insert(at,'...') if Ruby>=1.8
    str = a.join('/')
  end
  str
end

._get_axinfo(vary) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1709

def _get_axinfo(vary)
  if vary.nil?
    hash = {
            'title'=>'',
            'units'=>''
           }
  else
    #raise "Not 1D" if vary.rank!=1
    hash = {
            'title'=>(vary.long_name || vary.name), #.gsub('_','' )
            'units'=>(vary.get_att('units') || '')             #.gsub('_',' ')
           }
  end
  hash
end

._scatter_display_corr(gpx, gpy, opts) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2887

def _scatter_display_corr(gpx, gpy, opts)
  require "numru/ganalysis/covariance" 
  r, n = gpx.correlation(gpy,0)
  r = r.val if r.is_a?(GPhys)
  DCL::uxsttl("B","r=#{DCL.chval('A',r)}", 0)
  return [r, n]
end

._set_contour_levels_(opts) ⇒ Object

Raises:

• (ArgumentError)

# File 'lib/numru/ggraph.rb', line 2360

def _set_contour_levels_(opts)
  raise ArgumentError, "'levels' must be explicitly set" if !opts['levels']
  DCLExt.ud_set_contour(opts['levels'],opts['index'],opts['line_type'],
              opts['label'],opts['label_height'])
  @@set_contour_levels=true
  nil
end

._set_tone_levels_(opts) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2556

def _set_tone_levels_(opts)
  if !opts['levels'] && !opts['patterns']
  end
  if opts['levels']
    levels = opts['levels']
  else
    raise ArgumentError,"'levels' must be explicitly set"
  end
  if opts['patterns']
    patterns = opts['patterns']
  else
    nlev = levels.length - 1
    itpat = DCL.ueiget('itpat')
    iclr1 = DCL.ueiget('icolor1')
    iclr2 = DCL.ueiget('icolor2')
    patterns = (0...nlev).collect{|i|
      (((iclr2-iclr1)/(nlev-1.0)*i+iclr1).round) * 1000 + itpat
    }
  end
  DCLExt.ue_set_tone(levels,patterns)
  @@set_tone_levels=true
  nil
end

.add_mode_vectors(mean, modes, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2902

def add_mode_vectors(mean, modes, options=nil)
  opts = @@add_mode_vectors_options.interpret(options)
  unless modes.shape == [2,2] || modes.shape[2,1] || modes.shape == [2]
    raise ArgumentError, "the second arg must be a NArray of the shape [2,2], [2,1], or [2]."
  end
  unless (mean.length == 2) && (mean.is_a?(GPhys) || mean.is_a?(NArray) || mean.is_a?(Array))
    raise ArgumentError, "the 1st arg must be a GPhys, NArray, or Array with length of 2" 
  end
  mean = mean.val if mean.is_a?(GPhys)
  modes = modes.val * opts['fact']
  index = opts['lineindex']
  type = opts['linetype']
  modes = modes.newdim(-1) if modes.rank == 1
  for i in 0...(modes.shape[1])
    case opts['style']
    when "line"
      DCL::sgplzu([mean[0]-modes[0,i],mean[0]+modes[0,i]],
                  [mean[1]-modes[1,i],mean[1]+modes[1,i]],
                  type,index)
    when "arrow"
      DCL::sglazu(mean[0]-modes[0,i],mean[1]-modes[1,i],
                  mean[0]+modes[0,i],mean[1]+modes[1,i],
                  type,index)
    when "ellipse"
      raise "Drawing an ellipse needs two modes" unless modes.shape[1]==2
      arg = NArray.int(361).indgen!.to_f / 180.0 * NMath::PI
      s = NMath.sin(arg)
      c = NMath.cos(arg)
      x = mean[0] + modes[0,0] * c + modes[0,1] * s
      y = mean[1] + modes[1,0] * c + modes[1,1] * s
      DCL::sgplzu(x,y,type,index)
    else
      raise "invalid value for the option 'style' in add_mode_vectors: #{kind}"
    end
  end
end

.annotate(str_ary, noff = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1152

def annotate(str_ary, noff=nil)
  charsize = 0.7 * DCL.uzpget('rsizec1')
  vxmin,vxmax,vymin,vymax = DCL.sgqvpt
  vx = vxmax + 0.01
  vy = vymax - charsize/2
  nannot = annotate_at(str_ary,vx,vy,charsize,-1.0,-1.5,noff)
  @@nannot = nannot
end

.annotate_at(str_ary, vx, vy, charsize = nil, align = -1.0,, dvyfact = -1.5,, noff = nil) ⇒ Object

Raises:

• (TypeError)

# File 'lib/numru/ggraph.rb', line 1161

def annotate_at(str_ary,vx,vy,charsize=nil, align=-1.0, dvyfact=-1.5, noff=nil)
  raise TypeError,"Array expected" if ! str_ary.is_a?(Array)
  charsize = 0.7 * DCL.uzpget('rsizec1') unless(charsize)
  dvy = charsize*dvyfact
  vy += noff*dvy if noff
  str_ary.each{|str|
    DCL::sgtxzr(vx,vy,str,charsize,0,align,1)
    vy += dvy
  }
  nannot = str_ary.length
  nannot += noff if noff
  nannot
end

.axes(xax = nil, yax = nil, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1606

def axes(xax=nil, yax=nil, options=nil)
  if @@next_axes
    options = ( options ? @@next_axes.update(options) : @@next_axes )
    @@next_axes = nil
  end
  opts = @@axes.interpret(options)
  if itr51 = (DCL.sgqtrn == 51)
	wnd = DCL::sgqwnd
	DCL::sgswnd( *(cwd = DCL::sgqcwd) )
	DCL::sgstrn(1)
	DCL::sgstrf
  end
  if opts['xside'].length > 0
	xai = _get_axinfo(xax)
    if opts['xmaplabel'] == "lon"
      opts['xside'].split('').each{|s|     # scan('.') also works
        DCLExt.lon_ax('cside'=>s, 'dtick1'=>opts['xtickint'], 'dtick2'=>opts['xlabelint'])
      }
      DCL.uxsttl('b', (opts['xtitle'] || xai['title']), 0.0)
    elsif opts['xmaplabel'] == "lat"
      opts['xside'].split('').each{|s|     # scan('.') also works
        DCLExt.lat_ax('xax'=>true, 'cside'=>s, 'dtick1'=>opts['xtickint'], 'dtick2'=>opts['xlabelint'])
      }
      DCL.uxsttl('b', (opts['xtitle'] || xai['title']), 0.0)
    else
      sunits = opts['xunits'] || xai['units']
      units = Units[sunits]
      if opts['time_ax'] != false && \
         /since/ =~ sunits && units =~ Units['days since 0001-01-01'] && \
         ( calendar = ( xax.get_att('calendar') || nil ) ; \
           UNumeric::supported_calendar?(calendar) )
        opts['xside'].split('').each do |s|   # scan('.') also works
          ttl = opts['xtitle'] || xai['title']
          __calendar_ax(opts['time_ax'], true, s, sunits, calendar, ttl,
            	   opts['xtickint'], opts['xlabelint'])
        end
      else
        DCL.uscset('cxttl', (opts['xtitle'] || xai['title']) )
        DCL.uscset('cxunit', sunits)
        DCL.uspset('dxt', opts['xtickint'])  if(opts['xtickint'])
        DCL.uspset('dxl', opts['xlabelint']) if(opts['xlabelint'])
        if(opts['xloglabelall'])
          DCL.uspset('nlblx', 4)
          DCL.ulsxbl([1,2,3,4,5,6,7,8,9])
        end
        opts['xside'].split('').each{|s|     # scan('.') also works
          DCL.usxaxs(s)
        }
      end
    end
  end
  if opts['yside'].length > 0
    yai = _get_axinfo(yax)
    if opts['ymaplabel'] == "lon"
      opts['yside'].split('').each{|s|   # scan('.') also works
        DCLExt.lon_ax('yax'=>true, 'cside'=>s, 'dtick1'=>opts['ytickint'], 'dtick2'=>opts['ylabelint'])
      }
      DCL.uysttl('l', (opts['ytitle'] || yai['title']), 0.0)
    elsif opts['ymaplabel'] == "lat"
      opts['yside'].split('').each{|s|   # scan('.') also works
        DCLExt.lat_ax('cside'=>s, 'dtick1'=>opts['ytickint'], 'dtick2'=>opts['ylabelint'])
      }
      DCL.uysttl('l', (opts['ytitle'] || yai['title']), 0.0)
    else
      sunits=(opts['yunits'] || yai['units'])
      units = Units[sunits]
      if opts['time_ax'] != false && \
         /since/ =~ sunits && units =~ Units['days since 0001-01-01'] && \
         ( calendar = ( yax.get_att('calendar') || nil ) ; \
           UNumeric::supported_calendar?(calendar) )
        opts['yside'].split('').each do |s|   # scan('.') also works
          ttl = opts['ytitle'] || yai['title']
          __calendar_ax(opts['time_ax'], false, s, sunits, calendar, ttl,
            	   opts['ytickint'], opts['ylabelint'])
        end
      else
        DCL.uscset('cyttl', (opts['ytitle'] || yai['title']) )
        DCL.uscset('cyunit', sunits )
        DCL.uspset('dyt', opts['ytickint'])  if(opts['ytickint'])
        DCL.uspset('dyl', opts['ylabelint']) if(opts['ylabelint'])
        if(opts['yloglabelall'])
          DCL.uspset('nlbly', 4)
          DCL.ulsybl([1,2,3,4,5,6,7,8,9])
        end
        opts['yside'].split('').each{|s|   # scan('.') also works
          DCL.usyaxs(s)
          if s=='l' && sunits && sunits != ''
            DCL.uzpset('roffxt', DCL.uzpget('rsizec1')*1.5 )
          end
        }
      end
    end
	title(opts['title']) if opts['title']
  end
  if itr51
	DCL::sgstrn(51)
	DCL::sgswnd( *wnd )
	DCL::sgscwd( *cwd )
	DCL::sgstrf
  end
  nil
end

.axes_or_map_and_ttl(gp, opts, xax, yax, gp2 = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2411

def axes_or_map_and_ttl(gp, opts, xax, yax, gp2=nil)
  if !gp2
    ttl = opts['title'] || name_selector(gp, 30)
  else
    ttl = opts['title'] || '(' << name_selector(gp, 15) << ',' << name_selector(gp2, 15) << ')'
  end
  if map_trn?
    map
    if opts['map_axes'] && itr_is?(10)
      vpt = DCL.sgqvpt
      wnd = DCL.sgqwnd
      if wnd.include?(DCL.glrget('rundef'))
        map_fit = false
        wnd = DCL.sgqtxy
        cnt = DCL.umqcnt
      else
        map_fit = true
      end
      vrat = ( (vpt[3]-vpt[2]) / (vpt[1]-vpt[0]) ).abs
      wrat = ( (wnd[3]-wnd[2]) / (wnd[1]-wnd[0]) ).abs
      if wrat < vrat
        vyoff = (vpt[3]+vpt[2])/2
        dvy2 =  (vpt[1]-vpt[0])*wrat/2
        vpt[2] = vyoff - dvy2
        vpt[3] = vyoff + dvy2
      else
        vxoff = (vpt[1]+vpt[0])/2
        dvx2 =  (vpt[3]-vpt[2])/wrat/2
        vpt[0] = vxoff - dvx2
        vpt[1] = vxoff + dvx2
      end
      trn = DCL.sgqtrn
      if map_fit
        fig(xax,yax, {'itr'=>1, 'new_frame'=>false, 'viewport'=>vpt})
        axes(xax, yax, 'title'=>ttl)
        fig(xax,yax,{'itr'=>trn, 'new_frame'=>false, 'viewport'=>vpt, 'map_fit'=>true})
      else
        xax_map = xax[0..1].copy
        xax_map[0] = cnt[0] + wnd[0]
        xax_map[1] = cnt[0] + wnd[1]
        yax_map = yax[0..1].copy
        yax_map[0] = wnd[2]
        yax_map[1] = wnd[3]
        fig(xax_map,yax_map,{'itr'=>1, 'new_frame'=>false, 'viewport'=>vpt})
        axes(xax_map, yax_map, 'title'=>ttl)
        fig(xax_map,yax_map,{'itr'=>trn, 'new_frame'=>false, 'viewport'=>vpt, 'map_fit'=>false, 'map_axis'=>cnt, 'map_window'=>wnd})
      end
    else
      title( ttl )
    end
  elsif itr_is?(5)
    polar_coordinate_boundaries(xax,yax)
    title( ttl )
  else
    axes(xax, yax, 'title'=>ttl)
  end
  annotate(gp.lost_axes) if opts['annotate']
end

.clear_contour_levels ⇒ Object

# File 'lib/numru/ggraph.rb', line 2368

def clear_contour_levels
  @@set_contour_levels=false
  DCL.udiclv   # this may not be needed
  nil
end

.clear_tone_levels ⇒ Object

# File 'lib/numru/ggraph.rb', line 2579

def clear_tone_levels
  @@set_tone_levels=false
  nil
end

.color_bar(*args) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1133

def color_bar(*args)
  DCLExt.color_bar(*args)
end

.color_scatter(fx, fy, fz, newframe = true, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2800

def color_scatter(fx, fy, fz, newframe=true, options=nil)
  gropn_1_if_not_yet
  if newframe!=true && newframe!=false
    raise ArgumentError, "4th arg (newframe) must be true or false"
  end
  if @@next_color_scatter
    options = ( options ? @@next_color_scatter.update(options) : @@next_color_scatter )
    @@next_color_scatter = nil
  end
  opts = @@color_scatter_options.interpret(options)
  # interpret cut, slice, etc
  gpx, closerx = data_prep_multidim(fx, opts)
  gpy, closery = data_prep_multidim(fy, opts)
  gpz, closerz = data_prep_multidim(fz, opts)
  # flatten gphys data -----
  fakeaxis = Axis.new.set_pos(VArray.new(NArray.int(gpx.length).indgen!,
                                         {"long_name"=>"fake",
                                           "units"=>""},
                                         "fake"))
  fakegrid = Grid.new(fakeaxis)
  gpx, gpy, gpz = [gpx, gpy, gpz].collect{|gp|
    na = gp.val.reshape!(gp.length)
    attr = gp.data.attr_copy
    newgp = GPhys.new(fakegrid,VArray.new(na,attr,gp.name))
    newgp.set_lost_axes(gp.lost_axes)
    newgp
  }
  # ------------------------
  xax = gpx.data
  yax = gpy.data
  if newframe
    fig(xax, yax)
    axes_or_map_and_ttl(gpz, opts, xax, yax)
  end

  zval = gpz.data.val

  # < levels setting --- taken from tone >

  if opts['clr_min']
    icolor1_bak = DCL.uepget('icolor1')
    DCL.uepset('icolor1', opts['clr_min'])
  end
  if opts['clr_max']
    icolor2_bak = DCL.uepget('icolor2')
    DCL.uepset('icolor2', opts['clr_max'])
  end

  rmiss = DCL.glrget('rmiss')

  if !opts['keep'] or DCL.ueqntl==0
    if !opts['levels'] && opts['log']
      levels, = DCLExt.quasi_log_levels_z(zval, opts['nlev'],
                              opts['max'], opts['min'], opts['log_cycle'])
      levels.unshift(rmiss) if !opts['max'] && levels[0]<0
      levels.push(rmiss)    if !opts['max'] && levels[-1]>0
      opts['levels'] = levels
    end
    if opts['levels']
      backup = @@set_tone_levels
      _set_tone_levels_(opts)
      @@set_tone_levels = backup
    elsif opts['patterns']
      # Here only the colors are used from patterns
      raise "option 'patterns' is not effective unless 'levels' is specified"
    elsif !@@set_tone_levels
      DCL.ueitlv
      DCLExt.ue_set_linear_levs(zval, opts)
    else
      DCL.ueitlv if !@@set_tone_levels
    end
  end
  DCL::uemrkz(xax.val, yax.val, zval, opts['type'], opts['index'], opts['size'])
  if opts['correlation']
    return _scatter_display_corr(gpx, gpy, opts) 
  else
    return nil
  end
  #< finish >
ensure
  DCL.uepset('icolor1', icolor1_bak) if opts['clr_min']
  DCL.uepset('icolor2', icolor2_bak) if opts['clr_max']
  closerx.call if closerx
  closery.call if closery
  closerz.call if closerz
end

.contour(gphys, newframe = true, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2470

def contour(gphys, newframe=true, options=nil)
  gropn_1_if_not_yet
  if newframe!=true && newframe!=false
    raise ArgumentError, "2nd arg (newframe) must be true or false"
  end
  if @@next_contour
    options = ( options ? @@next_contour.update(options) : @@next_contour )
    @@next_contour = nil
  end
  if @@next_linear_contour_options
    options = ( options ? @@next_linear_contour_options.update(options) :
                          @@next_linear_contour_options )
    @@next_linear_contour_options = nil
  end
  opts = @@contour_options.interpret(options)

  xax, yax, closer, gp = data_prep_2D([gphys], newframe, opts)
  if newframe
    fig(xax, yax)
  end
  if newframe
    axes_or_map_and_ttl(gp,opts, xax, yax)
  end
  if !opts['keep'] or DCL.udqcln==0
    if !opts['levels'] && opts['log']
      levels, mjmn = DCLExt.quasi_log_levels_z(gp.data.val, opts['nlev'],
                              opts['max'], opts['min'], opts['log_cycle'])
      opts['levels'] = levels
      opts['index'] = NArray.to_na(mjmn)*2 + 1 if !opts['index']
      opts['label'] = true if opts['label'].nil?
    end
    if (opts['levels'])
      backup = @@set_contour_levels
      _set_contour_levels_(opts)
      @@set_contour_levels = backup
    elsif !@@set_contour_levels
      DCLExt.ud_set_linear_levs(gp.data.val, opts)
    end
  end
  saved_udparams = nil
  if DCL.udqcln >= 3
    lv1, = DCL.udqclv(1)
    lv2, = DCL.udqclv(2)
    lv3, = DCL.udqclv(3)
    if lv2-lv1 != lv3-lv2
      saved_udparams = DCLExt.ud_set_params({'lmsg'=>false})
    end
  end
  DCL.uwsgxa(xax.val) if xax.rank==1
  DCL.uwsgya(yax.val) if yax.rank==1
  DCL.udcntz(gp.data.val)
  DCLExt.ud_set_params(saved_udparams) if saved_udparams
  closer.call if closer
  nil
end

.current_itr(fig_yet_to_be_called = false) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1228

def current_itr( fig_yet_to_be_called=false )
  if fig_yet_to_be_called
    current_itr = ( @@next_fig && @@next_fig['itr'] ) || @@fig['itr']
  else
    current_itr = DCL.sgqtrn
  end
  current_itr
end

.data_prep_1D(gphys, opts) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1903

def data_prep_1D(gphys, opts)
  gp = gphys
  if gp.is_a?(String)
    gp = GPhys::IO.str2gphys(gp) 
    closer = Proc.new{gp.close}
  else
    closer = nil
  end
  gp = gp[*opts['slice']] if opts['slice']
  if (cut = opts['cut'])  # substitution, not ==
    if cut.is_a?(Array) 
      gp = gp.cut(*cut) 
    else
      gp = gp.cut(cut) 
    end
  end
  gp = gp.mean(*opts['mean']) if opts['mean']
  gp = gp.first1D
  window = ( @@next_fig['window'] if @@next_fig ) || @@fig['window']
  window = ( window ? window.dup : [nil, nil, nil, nil])
  max = opts['max']; min = opts['min']
  if !opts['exchange']
    xd = gp.coord(0)
    yd = gp.data
    window[2] = min if min
    window[3] = max if max
  else
    yd = gp.coord(0)
    xd = gp.data
    window[0] = min if min
    window[1] = max if max
  end
  [xd, yd, closer, gp, window]
end

.data_prep_1D2(fx, fy, opts) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1938

def data_prep_1D2(fx, fy, opts)
  gpx = fx
  gpy = fy
  to_close = Array.new
  if gpx.is_a?(String)
    gpx = GPhys::IO.str2gphys(gpx) 
    to_close.push(gpx)
  end
  if gpy.is_a?(String)
    gpy = GPhys::IO.str2gphys(gpy) 
    to_close.push(gpy)
  end
  if to_close.length > 0
    closer = Proc.new{to_close.each{|gp| gp.close}}
  else
    closer = nil
  end

  gpx = gpx[*opts['slice']] if opts['slice']
  if (cut = opts['cut'])  # substitution, not ==
    if cut.is_a?(Array) 
      gpx = gpx.cut(*cut) 
    else
      gpx = gpx.cut(cut) 
    end
  end
  gpx = gpx.mean(*opts['mean']) if opts['mean']
  gpx = gpx.first1D

  gpy = gpy[*opts['slice']] if opts['slice']
  if (cut = opts['cut'])  # substitution, not ==
    if cut.is_a?(Array) 
      gpy = gpy.cut(*cut) 
    else
      gpy = gpy.cut(cut) 
    end
  end
  gpy = gpy.mean(*opts['mean']) if opts['mean']
  gpy = gpy.first1D

  [gpx, gpy, closer]
end

.data_prep_2D(gps, newframe, opts) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1827

def data_prep_2D(gps, newframe, opts)
  opend_here = []
  gps.collect!{|gp| 
    if gp.is_a?(String)
      gp = GPhys::IO.str2gphys(gp)
      opend_here.push(gp) 
    end
    gp = gp[*opts['slice']] if opts['slice']
    if (cut = opts['cut'])  # substitution, not ==
      if cut.is_a?(Array) 
        gp = gp.cut(*cut) 
      else
        gp = gp.cut(cut) 
      end
    end
    gp = gp.mean(*opts['mean']) if opts['mean']
    gp = gp.first2D
    gp = gp.transpose(1,0) if opts['transpose'] || opts['exchange']
    if opts['xintv'] || opts['yintv']
      gp = gp.copy
      sh = gp.shape
      if opts['xintv'] && ((xi=opts['xintv']) >= 2)
        idx = NArray.int((sh[0]/xi.to_f).ceil).indgen!*xi  # 0,xi,2*xi,..
        gp = gp[idx, true]
      end
      if opts['yintv'] && ((yi=opts['yintv']) >= 2)
        idx = NArray.int((sh[1]/yi.to_f).ceil).indgen!*yi  # 0,yi,2*yi,..
        gp = gp[true, idx]
      end
    end
    gp
  }
  if opend_here.length > 0
    closer = Proc.new{opend_here.each{|gp| gp.data.file.close}}
  else
    closer = nil
  end
  xd = gps[0].coord(0)
  yd = gps[0].coord(1)
  if xc=opts['xcoord']
    xa = gps[0].coord(xc)
    if xa.rank == 1 && xa.length != xd.length
	  raise("Len of #{xc} do not agree with the #{xd.name} dim")
	elsif xa.rank == 2 && xa.shape != gps[0].shape
	  raise("Invalid 2D coord shape. #{xc} : #{xa.shape.inspect}" +
		" vs #{gps[0].shape.inspect}")
    end
	xd = xa
  end
  if yc=opts['ycoord']
    ya = gps[0].coord(yc)
    if ya.rank == 1 && ya.length != yd.length
	  raise("Len of #{yc} do not agree with the #{yd.name} dim")
	elsif ya.rank == 2 && ya.shape != gps[0].shape
	  raise("Invalid 2D coord shape. #{yc} : #{ya.shape.inspect}" +
		" vs #{gps[0].shape.inspect}")
    end
	yd = ya
  end
  if map_trn?(newframe)
    # // convert to the anti-clockwise coordinate in case for map filling
    #    (Because of a bug in UMFMAP in DCL 5.3)  -->
    if ( (xd[-1].val-xd[0].val)*(yd[-1].val-yd[0].val) < 0 )
      gps.collect!{|gp| gp.copy[true,-1..0]}
      yd = gps[0].coord(1)
    end
    # <-- ///
    gps.collect!{|gp| gp.cyclic_ext(0, 360.0)} #cyclic extentn along lon
    xd = gps[0].coord(0)
  elsif itr_is?(5, newframe)
    gps.collect!{|gp| gp.cyclic_ext(1, 360.0)} #cyclic ext along azimuth
    yd = gps[0].coord(1)
  end
  [xd, yd, closer, *gps]
end

.data_prep_multidim(gphys, opts) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1981

def data_prep_multidim(gphys, opts)
  gp = gphys
  if gp.is_a?(String)
    gp = GPhys::IO.str2gphys(gp) 
    closer = Proc.new{gp.close}
  else
    closer = nil
  end
  gp = gp[*opts['slice']] if opts['slice']
  if (cut = opts['cut'])  # substitution, not ==
    if cut.is_a?(Array) 
      gp = gp.cut(*cut) 
    else
      gp = gp.cut(cut) 
    end
  end
  gp = gp.mean(*opts['mean']) if opts['mean']
  [gp, closer]
end

.fig(xax = nil, yax = nil, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1247

def fig(xax=nil, yax=nil, options=nil)

  # xax and yax are needed (i.e. Axis objects) if not map projection

  if @@next_fig
    options = ( options ? @@next_fig.update(options) : @@next_fig )
    @@next_fig = nil
  end
  opts = @@fig.interpret(options)

  if opts['no_new_fig']
    # do nothing
  elsif opts['new_frame']
    DCL.grfrm
  else
    DCL.grfig
  end
  raise "viewport's length must be 4" if opts['viewport'].length != 4
  DCL.grsvpt(*opts['viewport'])
  @@nannot = 0

  if opts['window']
    if !opts['window'].is_a?(Array) || opts['window'].length!=4
      raise "Option 'window' must be an Array of length==4"
    end
  end

  itr = opts['itr']
  DCL.grstrn(itr)

  map_fit = ( (itr==10 and opts['map_fit']!=false) or
              (itr==11 and opts['map_fit']==true)  )
  DCL.sgpset("RSAT",opts['map_rsat']) if (opts['map_rsat'])

  if ( (1<=itr and itr<=4) or (51<=itr and itr<=99) or map_fit )
    window = opts['window']
    window = ( window ? window.dup : [nil, nil, nil, nil])
    if window.include?(nil)
      raise(ArgumentError, "xax and yax must be provided") if !xax or !yax
      if (xreverse=opts['xreverse']).is_a?(String)
        atts = opts['xreverse'].split(',').collect{|v| v.split(':')}
        xreverse = false
        atts.each{|key,val|
          xreverse = ( key.downcase=="units" ? (Units.new(xax.get_att(key)) =~ Units.new(val)) : (xax.get_att(key) == val) )
          break if xreverse
        }
      end
      if (yreverse=opts['yreverse']).is_a?(String)
        atts = opts['yreverse'].split(',').collect{|v| v.split(':')}
        yreverse = false
        atts.each{|key,val|
          yreverse = ( key.downcase=="units" ? (Units.new(yax.get_att(key)) =~ Units.new(val)) : (yax.get_att(key) == val) )
          break if yreverse
        }
      end
      if xreverse
        xrange = [ xax.max, xax.min ]
      else
        xrange = [ xax.min, xax.max ]
      end
      if yreverse
        yrange = [ yax.max, yax.min ]
      else
        yrange = [ yax.min, yax.max ]
      end
      [ xrange, yrange ].each do |range|
        r0 = range[0]
        r1 = range[1]
        if r0 == r1
          if r0.nil?
            raise("Cannot set a window. Maybe plotting all-missing data?")
          end
          # if the max and min of the range is the same, shift them a bit
          range[0] = DCL.rgnlt(r0)
          range[1] = DCL.rgngt(r0)
          if range[0] == range[1] # if min and max are still the same
            range[0] -= 1
            range[1] += 1
          end
        end
      end
      default_window=[xrange[0], xrange[1], yrange[0], yrange[1]]
      default_window = __round_window(default_window) if opts['round0']
      if (if [].respond_to?("nitems")  
            window.nitems == 0
          else
            window.none?
          end) # if all elements is nil
        window = default_window
      else
        window.each_index do |i|
          window[i] = default_window[i] if window[i] == nil
        end
      end
      window = __round_window(window) if opts['round1']
    end
  end

  case itr
  when 1..4,51..99
    # all but for map projections and curvilinear coordinates
	if xax.rank==2 || yax.rank==2
	  if xax.rank==2
 nx,ny = xax.shape
	  else
 nx,ny = yax.shape
	  end
	  if x2d = (xax.rank==2)
 ux = NArray.sfloat(nx).indgen!
 cx = xax.val
	  else
 ux = xax.val
 cx = NArray.sfloat(nx,ny)
 cx[0,0] = DCL::glrget('rundef')
	  end
	  if y2d = (yax.rank==2)
 uy = NArray.sfloat(ny).indgen!
 cy = yax.val
	  else
 uy = yax.val
 cy = NArray.sfloat(nx,ny)
 cy[0,0] = DCL::glrget('rundef')
	  end
	  DCL::grswnd(0.0,ux.max,0.0,uy.max)
	  DCL::grscwd(*window)
	  DCL::g2sctr(ux,uy,cx,cy)
	  DCL::grstrn(51)
	  DCL.uwsgxa(ux)
	  DCL.uwsgya(uy)
	else
	  DCL.grswnd(*window)
	end
    if itr <= 4 and opts['eqdistvpt']
      # modify viewport to equidistant
      vx1,vx2,vy1,vy2 = opts['viewport']
      wx1,wx2,wy1,wy2 = window
      vxw = (vx2-vx1).abs
      vyw = (vy2-vy1).abs
      wxw = (wx2-wx1).abs
      wyw = (wy2-wy1).abs
      if ( [vxw,vyw,wxw,wyw].min > 0 )
        rx = wxw / vxw
        ry = wyw / vyw
        if rx > ry
          rc = ry/rx
          vym = (vy1+vy2)/2.0
          vy1 = vym - rc*vyw/2.0
          vy2 = vym + rc*vyw/2.0
        elsif rx < ry
          rc = rx/ry
          vxm = (vx1+vx2)/2.0
          vx1 = vxm - rc*vxw/2.0
          vx2 = vxm + rc*vxw/2.0
        end
        DCL.grsvpt(vx1,vx2,vy1,vy2)
      end
    end
  when 5,6
      if opts['similar']
        similar=opts['similar']
        DCL.grssim(*similar)
      elsif opts['window']
        if defined?(DCL::DCLVERSION) && DCL::DCLVERSION >= '5.3'
          DCL.sgscwd(*opts['window'])
        else
          raise "You need DCL 5.3 or later to use the 'window' parameter in this transform (#{itr}). Use 'similar' instead."
        end
      else
        case itr
        when 5
          vxmin,vxmax,vymin,vymax=DCL.sgqvpt
          raise(ArgumentError, "xax must be provided") if !xax
          simfac = (vymax-vymin)/(xax.max*2)    # only xax is used
          similar = [simfac,0.0,0.0]
          DCL.grssim(*similar)
        else
          raise NotImplementedError, "Sorry, automatic window setting is yet to be available for the transform #{itr}. Please specify the 'similar' or 'window' parameter"
        end
      end

  when 10..19
    sv = DCL.umpget('lglobe')
    if !map_fit
      map_axis = opts['map_axis'] || [180.0, 0.0, 0.0]
      DCL::umscnt( *map_axis )
      map_window = opts['map_window']
      raise "map_window's length must be 4" if map_window.length != 4
      map_lat_range = [ map_window[2], map_window[3] ]
      if map_lat_range.is_a?(Range)
        map_lat_range = [map_lat_range.first, map_lat_range.end]
      end
      case ( map_axis[2] - map_axis[1] ) % 360
      when 0
        tyrange = map_lat_range
      when 180
        tyrange = [ -map_lat_range[1], -map_lat_range[0] ]
      else
        tyrange = [ -75, 75 ]      # latange is ignored
      end
      DCL::grstxy(map_window[0], map_window[1], tyrange[0], tyrange[1])
      # 'lglobe' becomes true when itr >= 12 or data coverage is almost global for Mercator
      DCL.umpset('lglobe', (itr >= 12) || ((itr == 11) && (tyrange[0].abs >= 75) && (tyrange[1].abs >= 75) && (tyrange[0] * tyrange[1] < 0)))
    else
      lon_cent =( window[0] + window[1] ) / 2.0
      dlon2 = ( window[1] - window[0] ) / 2.0
      map_axis = [ lon_cent, 0.0, 0.0 ]
      DCL::sgswnd(*window)
      DCL::umscnt( *map_axis )
      DCL::grstxy( -dlon2, dlon2, window[2], window[3] )
      # 'lglobe' becomes true when data coverage is almost global for Mercator
      DCL.umpset('lglobe', (itr == 11) && (window[2].abs >= 75) && (window[3].abs >= 75) && (window[2] * window[3] < 0))
    end
    DCL::umpfit
  when 20..34
    map_axis = opts['map_axis'] || [180.0, 90.0, 0.0]
    map_radius = opts['map_radius'] || 70.0
    DCL::umscnt( *map_axis )
    case itr
    when 31
      v = opts['viewport']
      vptsize = [ v[1]-v[0], v[3]-v[2] ].min
      DCL::grssim( vptsize*0.25/Math::tan(0.5*Math::PI*map_radius/180.0), 0.0, 0.0 )
    when 22
      v = opts['viewport']
      vptsize = [ v[1]-v[0], v[3]-v[2] ].min
      DCL::grssim( vptsize*0.35*(90.0/map_radius), 0.0, 0.0 )
      DCL::grstxy(-180.0, 180.0, 90.0-map_radius, 90.0)
    else
      DCL::grstxy(-180.0, 180.0, 90.0-map_radius, 90.0)
    end
    sv = DCL.umpget('lglobe')
    DCL.umpset('lglobe', true)
    DCL::umpfit
    DCL.umpset('lglobe', sv)
  else
    raise "unsupported transformation number: #{itr}"
  end

  DCL.grstrf
  DCL.umpset('lglobe', sv)
  nil
end

.gropn_1_if_not_yet ⇒ Object

# File 'lib/numru/ggraph.rb', line 1098

def gropn_1_if_not_yet
  begin
    DCL.stqwrc
  rescue
    DCL.gropn(1)
  end
end

.histogram(gphys, newframe = true, options = nil) ⇒ Object Also known as: histogram1D

# File 'lib/numru/ganalysis/histogram.rb', line 162

def histogram(gphys, newframe=true, options=nil)
  gropn_1_if_not_yet
  if newframe!=true && newframe!=false
    raise ArgumentError, "2nd arg (newframe) must be true of false"
  end
  unless gphys.rank == 1
    raise ArgumentError, "rank of gphys must be 1"
  end
  unless gphys.axis(0).cell?
    raise ArgumentError, "axis must be cell type"
  end
  # if window is specified via GGraph#fig or GGraph#set_fig, use it.
  opts = @@histogram_options.interpret(options)
  exchange = opts["exchange"]
  window = opts["window"] || @@fig['window']
  unless exchange
    x = gphys.axis(0).cell_bounds
    y = gphys
    window[2] ||= 0
  else
    y = gphys.axis(0).cell_bounds
    x = gphys
    window[0] ||= 0
  end
  itr = @@fig['itr'] || DCL::sgpget("itr")
  if (itr==2 || itr==4)
    tmp = y.val
    if tmp.min * tmp.max < 0
      if tmp.min.abs < tmp.max
        mask = tmp.lt(0)
      else
        mask = tmp.gt(0)
      end
    else
      mask = tmp.ne(0)
    end
    if tmp.is_a?(NArrayMiss)
      tmp.set_mask(mask * tmp.get_mask)
    else
      tmp = NArrayMiss.to_nam_no_dup(tmp, mask)
    end
    y.replace_val(tmp)
    window[2] = tmp.abs.min * (tmp.min < 0 ? -1 : 1)
  end
  if (itr==3 || itr==4)
    tmp = x.val
    if tmp.min * tmp.max < 0
      if tmp.min.abs < tmp.max
        mask = tmp.lt(0)
      else
        mask = tmp.gt(0)
      end
    else
      mask = tmp.ne(0)
    end
    if tmp.is_a?(NArrayMiss)
      tmp.set_mask(mask * tmp.get_mask)
    else
      tmp = NArrayMiss.to_nam_no_dup(tmp, mask)
    end
    x.replace_val(tmp)
    window[0] = tmp.abs.min * (tmp.min < 0 ? -1 : 1)
  end
  opts["window"] = window

  if newframe
    fig(x, y, "window"=>opts["window"])
    axes(x, y, "title"=>opts["title"])
  end
  lmiss = DCL::gllget("lmiss")
  DCL::gllset("lmiss", true)
  unless exchange
    if opts["fill"]
      DCL::uvbxa(x.val, [window[2]] * y.length, y.val)
    end
    if (itr == 2) || (itr == 4)
      bottom = [y.val.min] * y.length
    else
      bottom = [window[2]] * y.length
    end
    DCL::uvbxf(x.val, bottom, y.val)
  else
    if opts["fill"]
      DCL::uhbxa([window[0]] * x.length, x.val, y.val)
    end
    if (itr == 3) || (itr == 4)
      bottom = [x.val.min] * x.length
    else
      bottom = [window[0]] * x.length
    end
    DCL::uhbxf(bottom, x.val, y.val)
  end
  DCL::gllset("lmiss", lmiss)
  return nil
end

.itr_is?(itr, fig_yet_to_be_called = false) ⇒ Boolean

Returns:

• (Boolean)

# File 'lib/numru/ggraph.rb', line 1224

def itr_is?( itr, fig_yet_to_be_called=false )
  itr == current_itr( fig_yet_to_be_called )
end

.line(gphys, newframe = true, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2018

def line(gphys, newframe=true, options=nil)
  gropn_1_if_not_yet
  if newframe!=true && newframe!=false
    raise ArgumentError, "2nd arg (newframe) must be true or false"
  end
  opts = @@line_options.interpret(options)
  x, y, closer, gp, window = data_prep_1D(gphys, opts)
  if newframe
    fig(x, y, {'window'=>window})
    axes_or_map_and_ttl(gp, opts, x, y)
  end

  if opts['label']
    lcharbk = DCL.sgpget('lchar')
    DCL.sgpset('lchar',true)
    DCL.sgsplc(opts['label'])
  end
  DCL.uulinz(x.val, y.val, opts['type'], opts['index'])
  DCL.sgpset('lchar',lcharbk) if opts['label']

  legend = opts['legend']
  if legend
    legend = gp.name if legend==true
    DCLExt::legend(legend,opts['type'],opts['index'],true,
                   opts['legend_size'],opts['legend_vx'],opts['legend_dx'],
                   opts['legend_vy'],newframe)
  end
  closer.call if closer
  nil
end

.map(options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1776

def map(options=nil)
  if @@next_map
    options = ( options ? @@next_map.update(options) : @@next_map )
    @@next_map = nil
  end
  opts = @@map.interpret(options)
  DCL.umplim if opts['lim']
  DCL.umpset('dgridmj',opts['dgridmj']) if opts['dgridmj']
  DCL.umpset('dgridmn',opts['dgridmn']) if opts['dgridmn']
  DCL.umpgrd if opts['grid']
  if opts['vpt_boundary']
    if opts['vpt_boundary'].is_a?(Integer)
      idx = opts['vpt_boundary']
    else
      idx = 1
    end
    DCL.slpvpr( idx )
  end
  if opts['wwd_boundary']
    if opts['wwd_boundary'].is_a?(Integer)
      idx = opts['wwd_boundary']
    else
      idx = 1
    end
    DCL.slpwwr( idx )
  end
  DCL.umpmap('coast_world') if opts['coast_world']
  DCL.umpmap('coast_japan') if opts['coast_japan']
  if opts['fill']
    DCL.umfmap('coast_world') if opts['coast_world']
    DCL.umfmap('coast_japan') if opts['coast_japan']
  end
  DCL.umpmap('border_world') if opts['border_world']
  DCL.umpmap('plate_world') if opts['plate_world']
  DCL.umpmap('state_usa') if opts['state_usa']
  DCL.umpmap('pref_japan') if opts['pref_japan']
  nil
end

.map_trn?(fig_yet_to_be_called = false) ⇒ Boolean

Returns:

• (Boolean)

# File 'lib/numru/ggraph.rb', line 1210

def map_trn?( fig_yet_to_be_called=false )
  if fig_yet_to_be_called
    itr = ( @@next_fig && @@next_fig['itr'] ) || @@fig['itr']
  else
    itr = DCL.sgqtrn
  end
  case itr
  when 1..9,51..99
    false
  else
    true
  end
end

.margin_info(program = nil, data_source = nil, char_height = nil, date = nil, xl = 0.0, xr = 0.0, yb = nil, yt = 0.0) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1106

def margin_info(program=nil, data_source=nil, char_height=nil, date=nil,
                xl=0.0, xr=0.0, yb=nil, yt=0.0)

  program = File.expand_path($0)       if !program
  if date
    program = __shorten_path(program,77) + '  ' + Date.today.to_s
  else
    program = __shorten_path(program,99)
    if date.nil?
      program = program + '  ' + Date.today.to_s if program.length < 77
    end
  end
  if !data_source
    data_source = Dir.pwd
    data_source = '' if data_source == program.sub(/\/[^\/]+$/,'')
  end
  data_source = __shorten_path(data_source,99)
  sz = 1.0/( program.length + data_source.length )
  char_height = [0.008, sz].min        if !char_height
  yb = 2.0 * char_height               if !yb

  DCL.slmgn(xl, xr, yb, yt)
  DCL.slsttl(program,     'b', -1.0, -1.0, char_height, 1)
  DCL.slsttl(data_source, 'b',  1.0, -1.0, char_height, 2)
  nil
end

.mark(gphys, newframe = true, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2110

def mark(gphys, newframe=true, options=nil)
  gropn_1_if_not_yet
  if newframe!=true && newframe!=false
    raise ArgumentError, "2nd arg (newframe) must be true or false"
  end
  opts = @@mark_options.interpret(options)
  x, y, closer, gp, window = data_prep_1D(gphys, opts)
  if newframe
    fig(x, y, {'window'=>window})
    axes_or_map_and_ttl(gp, opts, x, y)
  end
  DCL.uumrkz(x.val, y.val, opts['type'], opts['index'], opts['size'])

  legend = opts['legend']
  if legend
    legend = gp.name if legend==true
    DCLExt::legend(legend,opts['type'],opts['index'],false,
                   opts['legend_size'],opts['legend_vx'],nil,
                   opts['legend_vy'],newframe,opts['size'])
  end
  closer.call if closer
  nil
end

.name_selector(gp, limit) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2398

def name_selector(gp, limit)
  if gp.long_name && gp.name
    if gp.long_name.length > limit && gp.long_name.length > gp.name.length && !gp.name.empty?
      gp.name 
    else 
      gp.long_name
    end
  else
    gp.long_name || gp.name
  end
end

.next_axes(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1537

def next_axes(options)
  if options.is_a?(Hash)
    @@next_axes = options
  else
    raise TypeError,"Hash expected"
  end
  nil
end

.next_color_scatter(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2791

def next_color_scatter(options)
  if options.is_a?(Hash)
    @@next_color_scatter = options
  else
    raise TypeError,"Hash expected"
  end
  nil
end

.next_contour(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2379

def next_contour(options)
  if options.is_a?(Hash)
    @@next_contour = options
  else
    raise TypeError,"Hash expected"
  end
  nil
end

.next_fig(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1201

def next_fig(options)
  if options.is_a?(Hash)
    @@next_fig = options
  else
    raise TypeError,"Hash expected"
  end
  nil
end

.next_linear_contour_options(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2149

def next_linear_contour_options(options)
  if options.is_a?(Hash)
    @@next_linear_contour_options = options
  else
    raise TypeError,"Hash expected"
  end
  nil
end

.next_linear_tone_options(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2536

def next_linear_tone_options(options)
  if options.is_a?(Hash)
    @@next_linear_tone_options = options
  else
    raise TypeError,"Hash expected"
  end
  nil
end

.next_map(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1767

def next_map(options)
  if options.is_a?(Hash)
    @@next_map = options
  else
    raise TypeError,"Hash expected"
  end
  nil
end

.next_regression_line(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 3040

def next_regression_line(options)
  if options.is_a?(Hash)
    @@next_regression_line = options
  else
    raise TypeError,"Hash expected"
  end
  nil
end

.next_scatter(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2713

def next_scatter(options)
  if options.is_a?(Hash)
    @@next_scatter = options
  else
    raise TypeError,"Hash expected"
  end
  nil
end

.next_tone(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2589

def next_tone(options)
  if options.is_a?(Hash)
    @@next_tone = options
  else
    raise TypeError,"Hash expected"
  end
  nil
end

.next_unit_vect_options(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2944

def next_unit_vect_options(options)
  DCLExt.next_unit_vect_options(options)
end

.polar_coordinate_boundaries(xax = nil, yax = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1725

def polar_coordinate_boundaries(xax=nil,yax=nil)

  slln=DCL.sgpget('LLNINT') ; slgc=DCL.sgpget('LGCINT')
  DCL.sgpset('LLNINT',true) ; DCL.sgpset('LGCINT',true)

  xmin=xax.min ; xmax=xax.max ; ymin=yax.min ; ymax=yax.max
  azimuth_len = ymax-ymin
  azimuth_len = azimuth_len.val if azimuth_len.is_a?(UNumeric)
  unless( (azimuth_len - 360.0).abs < 1e-2 )
    # not a full circle
    DCL.sgplzu([xmin,xmax],[ymin,ymin],1,5)
    DCL.sgplzu([xmin,xmax],[ymax,ymax],1,5)
  end

  DCL.sgplzu([xmin,xmin],[ymin,ymax],1,5)
  DCL.sgplzu([xmax,xmax],[ymin,ymax],1,5)

  DCL.sgpset('LLNINT',slln) ; DCL.sgpset('LGCINT',slgc)
  nil
end

.profile(gphys, newframe = true, options = nil) ⇒ Object

———————————————————- Added at 080207 –

# File 'lib/numru/ggraph.rb', line 2049

def profile(gphys, newframe=true, options=nil)
  gropn_1_if_not_yet
  if newframe!=true && newframe!=false
    raise ArgumentError, "2nd arg (newframe) must be true or false"
  end
  
  opts = @@line_options.interpret(options)
  check_max = opts['max']
  opts['max'] = nil
  opts['exchange']= true
  x, y, closer, gp, window = data_prep_1D(gphys, opts)

  if newframe
    fig(x, y, {'window'=>window})
    axes_or_map_and_ttl(gp, opts, x, y)
  end
  if opts['label']
    lcharbk = DCL.sgpget('lchar')
    DCL.sgpset('lchar',true)
    DCL.sgsplc(opts['label'])
  end

  if opts['index']
     iat=opts['index']
  end
  if check_max
     if x.val.max > check_max
        iat=21
     end
  end
  DCL.uulinz(x.val, y.val, opts['type'], iat )
  DCL.sgpset('lchar',lcharbk) if opts['label']


  legend = opts['legend']
  if legend
    legend = gp.name if legend==true
    DCLExt::legend(legend,opts['type'],opts['index'],true,
                   opts['legend_size'],opts['legend_vx'],opts['legend_dx'],
                   opts['legend_vy'],newframe)
  end
  nil
end

.regression_line(x, y, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 3049

def regression_line(x,y, options=nil)
  if @@next_regression_line
    options = ( options ? @@next_regression_line.update(options) : @@next_regression_line )
    @@next_regression_line = nil
  end
  opts = @@regression_line_options.interpret(options)

  x = x.val if x.respond_to?(:val)  # GPhys or VArray like
  y = y.val if y.respond_to?(:val)  # GPhys or VArray like
  x = NArray.to_na(x) if x.is_a?(Array)
  y = NArray.to_na(y) if y.is_a?(Array)
   

  xs = x   # scaled x:  x or log(x)  (see below)
  ys = y   # scaled y:  y or log(y)  (see below)

  xlog = ylog = false  # can be overwritten below
  itr = current_itr()
  case itr
  when 1
    # nothing to do
  when 2
    ys = NMMath.log(ys)
    ylog = true
  when 3
    xs = NMMath.log(xs)
    xlog = true
  when 4
    xs = NMMath.log(xs)
    ys = NMMath.log(ys)
    xlog = ylog = true
  else
    raise "Drawing a linear regression line is limited to itr = 1~4 (linear or log plots), while the current itr is #{itr}. (itr is a coordinate transformation parameter set by fig (or set_fig or next_fig).)"
  end

  y_given_x = !opts['x_given_y']

  xr = NArray.sfloat(2)  # --> x ends of the regression line
  yr = NArray.sfloat(2)  # --> y ends of the regression line
  if y_given_x
    a,b = linear_regression(xs,ys)
    if opts["limit"]
      xr[0] = x.min
      xr[1] = x.max
    else
      xr[0],xr[1], = DCL.sgqwnd
    end
    xrs = !xlog ? xr : NMath::log(xr)
    yrs = a + b*xrs
    yr = !ylog ? yrs : NMath::exp(yrs)
  else
    a,b = linear_regression(ys,xs)
    if opts["limit"]
      yr[0] = y.min
      yr[1] = y.max
    else
      dum0,dum1, yr[0], yr[1] = DCL.sgqwnd
    end
    yrs = !ylog ? yr : Math::log(yr)
    xrs = a + b*yrs
    xr = !xlog ? xrs : Math::exp(xrs)
  end

  before = DCLExt.sg_set_params({'lclip'=>opts['clip']})
  DCL.sgplzu(xr,yr,opts["type"],opts["index"]) 
  DCLExt.sg_set_params(before)

  ai = opts['annot_intercept'] || opts['annot_intercept'].nil? && itr==1
  as = opts['annot_slope'] || opts['annot_slope'].nil? && (itr==1 or itr==4)
  if as && ai 
    if y_given_x
      annot = ["regr  y=a+bx:", 
               "  a=#{DCL.chval('A',a)}", "  b=#{DCL.chval('A',b)}"]
    else
      annot = ["regr  x=a+by:", 
               "  a=#{DCL.chval('A',a)}", "  b=#{DCL.chval('A',b)}"]
    end
    annotate(annot, @@nannot)
  elsif as  
    if y_given_x
      annot = ["slope: #{DCL.chval('A',b)}"]
    else
      annot = ["regress x given y:", " x(y) slope: #{DCL.chval('A',b)}"]
    end
    annotate(annot, @@nannot)
  end

  [a,b]
end

.scatter(fx, fy, newframe = true, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2722

def scatter(fx, fy, newframe=true, options=nil)
  gropn_1_if_not_yet
  if newframe!=true && newframe!=false
    raise ArgumentError, "3rd arg (newframe) must be true or false"
  end
  if @@next_scatter
    options = ( options ? @@next_scatter.update(options) : @@next_scatter )
    @@next_scatter = nil
  end
  opts = @@scatter_options.interpret(options)
  # interpret cut, slice, etc
  gpx, closerx = data_prep_multidim(fx, opts)
  gpy, closery = data_prep_multidim(fy, opts)
  # flatten gphys data -----
  fakeaxis = Axis.new.set_pos(VArray.new(NArray.int(gpx.length).indgen!,
                                         {"long_name"=>"fake",
                                           "units"=>""},
                                         "fake"))
  fakegrid = Grid.new(fakeaxis)
  gpx, gpy = [gpx, gpy].collect{|gp|
    na = gp.val.reshape!(gp.length)
    attr = gp.data.attr_copy
    newgp = GPhys.new(fakegrid,VArray.new(na,attr,gp.name))
    newgp.set_lost_axes(gp.lost_axes)
    newgp
  }
  # ------------------------
  x = gpx.data
  y = gpy.data
  if newframe
    fig(x, y)
    axes_or_map_and_ttl(gpx, opts, x, y, gpy)
  end
  DCL.uumrkz(x.val, y.val, opts['type'], opts['index'], opts['size'])
  if opts['correlation']
    return _scatter_display_corr(gpx, gpy, opts) 
  else
    return nil
  end
  # close if necessary -----
  closerx.call if closerx
  closery.call if closery
  nil
end

.set_axes(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1532

def set_axes(options)
  @@axes.set(options)
end

.set_color_scatter(options) ⇒ Object

@@color_scatter_options.delete(“transpose”) # delete is yet to be supported @@color_scatter_options.delete(“exchange”) # (2010-01)

# File 'lib/numru/ggraph.rb', line 2787

def set_color_scatter(options)
  @@color_scatter_options.set(options)
end

.set_contour(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2374

def set_contour(options)
  @@contour_options.set(options)
end

.set_contour_levels(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2356

def set_contour_levels(options)
  opts = @@contour_levels.interpret(options)
  _set_contour_levels_(opts)
end

.set_fig(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1196

def set_fig(options)
  @@fig.set(options)
end

.set_linear_contour_options(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2145

def set_linear_contour_options(options)
  @@linear_contour_options.set(options)
end

.set_linear_tone_options(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2532

def set_linear_tone_options(options)
  @@linear_tone_options.set(options)
end

.set_map(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1762

def set_map(options)
  @@map.set(options)
end

.set_regression_line(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 3036

def set_regression_line(options)
  @@regression_line_options.set(options)
end

.set_scatter(options) ⇒ Object

@@scatter_options.delete(“transpose”) # delete is yet to be supported @@scatter_options.delete(“exchange”) # (2010-01)

# File 'lib/numru/ggraph.rb', line 2709

def set_scatter(options)
  @@scatter_options.set(options)
end

.set_tone(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2584

def set_tone(options)
  @@tone_options.set(options)
end

.set_tone_levels(options) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2552

def set_tone_levels(options)
  opts = @@tone_levels.interpret(options)
  _set_tone_levels_(opts)
end

.set_unit_vect_options(options) ⇒ Object

< vector >

# File 'lib/numru/ggraph.rb', line 2941

def set_unit_vect_options(options)
  DCLExt.set_unit_vect_options(options)
end

.sim_trn? ⇒ Boolean

Returns:

• (Boolean)

# File 'lib/numru/ggraph.rb', line 1237

def sim_trn?
  itr = DCL.sgqtrn
  case itr
  when 5..7
    true
  else
    false
  end
end

.timeplot(gphys, newframe = true, options = nil) ⇒ Object

timeplot 0802011 ######

# File 'lib/numru/ggraph.rb', line 2159

def timeplot(gphys, newframe=true, options=nil)
  gropn_1_if_not_yet
  if newframe!=true && newframe!=false
    raise ArgumentError, "2nd arg (newframe) must be true or false"
  end
  if ! defined?(@@timeplot_options)
    @@timeplot_options = Misc::KeywordOptAutoHelp.new(
      ['size', 0.01, 'marks size'],
      ['lamin', nil, 'minimam latitude value'],
      ['lamax', nil, 'maximam latitude value'],
      ['lomin', nil, 'minimam longitude value'],
      ['lomax', nil, 'maximam longitude value'],
      ['ltone', true, 'Same as udpack parameter ltone'],
      ['tonf', false, 'Use DCL.uetonf instead of DCL.uetone'],
      ['tonb', false, 'Use DCL.uetonb instead of DCL.uetone'],
      ['tonc', false, 'Use DCL.uetonc instead of DCL.uetone'],
      ['clr_min', 13,  'if an integer (in 10..99) is specified, used as the color number for the minimum data values. (the same can be done by setting the uepack parameter "icolor1")'],
      ['clr_max', 99,  'if an integer (in 10..99) is specified, used as the color number for the maximum data values. (the same can be done by setting the uepack parameter "icolor2")'],
      ['keep', false, 'Use the tone levels and patterns used previously'],
      @@line_options
    )
  end

  opts = @@timeplot_options.interpret(options)
  opts['type'] = 2
  if opts['clr_min']
    icolor1_bak = DCL.uepget('icolor1')
    DCL.uepset('icolor1', opts['clr_min'])
  end
  if opts['clr_max']
    icolor2_bak = DCL.uepget('icolor2')
    DCL.uepset('icolor2', opts['clr_max'])
  end
  cmn=opts['clr_min']+1
  cmx=opts['clr_max']-1

  if opts['lamin'] && opts['lamax']
     latmn= opts['lamin']
     latmx= opts['lamax']
     if latmn > latmx
        raise ArgumentError,"Specific Settings parameter is invalid.(1)"
     end
  elsif opts['lamin'] || opts['lamax']
     raise ArgumentError,"Specific Settings parameter is invalid.(2)"
  else
     latmn=false
     latmx=false
  end

  if opts['lomin'] && opts['lomax']
     lonmn= opts['lomin']
     lonmx= opts['lomax']
     if lonmn > lonmx
        raise ArgumentError,"Specific Settings parameter is invalid.(3)"
     end
  elsif opts['lomin'] || opts['lomax']
     raise ArgumentError,"Specific Settings parameter is invalid.(4)"
  else
     lonmn=false
     lonmx=false
  end

  gp = gphys.first2D
  gp = gp.transpose(1,0) if opts['transpose'] || opts['exchange']
#--
  x = gp.coord(0)
  y = gp.coord(1)
  zv = gp.data.val
       
  idim=0
  case x.name.downcase
  when "time"
     idim=0
  else
     case y.name.downcase
     when "time"
       idim=1
     end
  end

  latp=VArray.new
  lonp=VArray.new
  lat_flag = false
  lon_flag = false
  if gp.has_assoccoord?
    gp.assoc_coords.coordnames.each do |auxname|
      if auxname == "Latitude"
        latp=gp.coord(auxname)
        lat_flag = true
      elsif auxname == "Longitude"
        lonp=gp.coord(auxname)
        lon_flag = true
      end
    end
  end

  #------ for color bar
  DCL.uwsgxa(x.val)
  DCL.uwsgya(y.val)
  if !opts['keep'] or DCL.ueqntl==0
    if !opts['levels'] && opts['log']
      levels, = DCLExt.quasi_log_levels_z(gp.data.val, opts['nlev'],
                              opts['max'], opts['min'], opts['log_cycle'])
      rmiss = DCL.glrget('rmiss')
      levels.unshift(rmiss) if !opts['max'] && levels[0]<0
      levels.push(rmiss)    if !opts['max'] && levels[-1]>0
      opts['levels'] = levels
    end
    if opts['levels']
      backup = @@set_tone_levels
      _set_tone_levels_(opts)
      @@set_tone_levels = backup
    elsif opts['patterns']
      raise "option 'patterns' is not effective unless 'levels' is specified"
    elsif !@@set_tone_levels && opts['ltone']
      DCL.ueitlv
      DCLExt.ue_set_linear_levs(gp.data.val, opts)
    else
      DCL.ueitlv if !@@set_tone_levels
    end
  end
  #-------
  if newframe
    fig(x, y)
    axes_or_map_and_ttl(gp, opts, x, y)
  end
  zmn=gp.data.val.min
  zmx=gp.data.val.max

  zmax = ((cmx-cmn)+cmn).round*10+3
  zmin = cmn.round*10+3
  #-------
  xx=x.val.to_a
  yy=y.val.to_a
  plat=latp.val.to_a if lat_flag
  plon=lonp.val.to_a if lon_flag
  zz=zv

  k=0
  for i in 0..yy.length-1
    for j in 0..xx.length-1
       tflg=0
       idim == 0 ? h=j : h=i

       tflg=1 if latmn && latmx && ( plat[h] < latmn || plat[h] > latmx )
       tflg=2 if lonmn && lonmx && ( plon[h] < lonmn || plon[h] > lonmx )
       if tflg == 0
          xx1=NArray.to_na(Array.new(1,xx[j]*1.0))
          yy1=NArray.to_na(Array.new(1,yy[i]*1.0))
          z = zz[j,i] if zz[j,i].is_a?(Numeric)
          z = zmn if zz[j,i].is_a?(Array)
          iz=((z-zmn)/(zmx-zmn)*(cmx-cmn)+cmn).round*10+3
          DCL.uumrkz(xx1, yy1, 230, iz, 0.005 )
       end
       k+=1
    end
  end

  legend = opts['legend']
  if legend
    legend = gp.name if legend==true
    DCLExt::legend(legend,opts['type'],opts['index'],false,
                   opts['legend_size'],opts['legend_vx'],nil,
                   opts['legend_vy'],newframe,opts['size'])
  end
  DCL.uepset('icolor1', icolor1_bak) if opts['clr_min']
  DCL.uepset('icolor2', icolor2_bak) if opts['clr_max']
  #----------------------------- color bar unit -- 08.01.31 ----
  vx1, vx2, vy1, vy2 = DCL.sgqvpt
  x1=vx2+0.02
  y1=vy1-0.02
  if  "#{gp.data.units}".length >0
      DCL.sgtxzv(x1, y1, "#{gp.data.units}", 0.02, 0, -1 ,11 )
  end
  #-----------------------------------------
  nil
end

.title(string) ⇒ Object

# File 'lib/numru/ggraph.rb', line 1137

def title(string)
  if string
    if map_trn? || itr_is?(5)
      v = DCL::sgqvpt
      txhgt = DCL.uzpget('rsizec2')     # rsizec2: larger text
      vx = (v[0]+v[1])/2
      vy = v[3] + (0.5 + DCL.uzpget('pad1')) * txhgt
      DCL::sgtxzr(vx , vy, string, txhgt, 0, 0, 3)
    else
      DCL.uxmttl('t',string,0.0)
    end
  end
  nil
end

.tone(gphys, newframe = true, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2616

def tone(gphys, newframe=true, options=nil)
  gropn_1_if_not_yet
  if newframe!=true && newframe!=false
    raise ArgumentError, "2nd arg (newframe) must be true or false"
  end
  if @@next_tone
    options = ( options ? @@next_tone.update(options) : @@next_tone )
    @@next_tone = nil
  end
  if @@next_linear_tone_options
    options = ( options ? @@next_linear_tone_options.update(options) :
                          @@next_linear_tone_options )
    @@next_linear_tone_options = nil
  end
  opts = @@tone_options.interpret(options)
  xax, yax, closer, gp = data_prep_2D([gphys], newframe, opts)

  if opts['clr_min']
    icolor1_bak = DCL.uepget('icolor1')
    DCL.uepset('icolor1', opts['clr_min'])
  end
  if opts['clr_max']
    icolor2_bak = DCL.uepget('icolor2')
    DCL.uepset('icolor2', opts['clr_max'])
  end
  if newframe
    fig(xax, yax)
  end
  DCL.uwsgxa(xax.val) if xax.rank==1
  DCL.uwsgya(yax.val) if yax.rank==1
  if !opts['keep'] or DCL.ueqntl==0
    if !opts['levels'] && opts['log']
      levels, = DCLExt.quasi_log_levels_z(gp.data.val, opts['nlev'],
                              opts['max'], opts['min'], opts['log_cycle'])
      rmiss = DCL.glrget('rmiss')
      levels.unshift(rmiss) if !opts['max'] && levels[0]<0
      levels.push(rmiss)    if !opts['max'] && levels[-1]>0
      opts['levels'] = levels
    end
    if opts['levels']
      backup = @@set_tone_levels
      _set_tone_levels_(opts)
      @@set_tone_levels = backup
    elsif opts['patterns']
      raise "option 'patterns' is not effective unless 'levels' is specified"
    elsif !@@set_tone_levels && opts['ltone']
      DCL.ueitlv
      DCLExt.ue_set_linear_levs(gp.data.val, opts)
    else
      DCL.ueitlv if !@@set_tone_levels
    end
  end

  val = gp.data.val
  if ( opts['ltone'] && 
      ( opts['tonf'] || opts['auto'] && val.length > TONE_TONF_THRES ) )
    DCL.uetonf(val)
  elsif opts['tonb'] && opts['ltone']
    DCL.uetonb(val)
  elsif opts['tonc'] && opts['ltone']
    DCL.uetonc(val)
  else
    DCL.uetone(val)
  end

  if newframe
    axes_or_map_and_ttl(gp,opts, xax, yax)
  end
  DCL.uepset('icolor1', icolor1_bak) if opts['clr_min']
  DCL.uepset('icolor2', icolor2_bak) if opts['clr_max']

  if opts['color_bar']
    color_bar "log"=>opts['log']
  end
  closer.call if closer
  nil
end

.tone_and_contour(gphys, newframe = true, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2339

def tone_and_contour(gphys, newframe=true, options=nil)
  options ? copt=@@contour_options.select_existent(options) : copt=nil
  options ? topt=@@tone_options.select_existent(options) : topt=nil
  tone(gphys, newframe, topt)
  contour(gphys, false, copt)
end

.vector(fx, fy, newframe = true, options = nil) ⇒ Object

# File 'lib/numru/ggraph.rb', line 2971

def vector(fx, fy, newframe=true, options=nil)
  gropn_1_if_not_yet
  if newframe!=true && newframe!=false
    raise ArgumentError, "2nd arg (newframe) must be true or false"
  end
  opts = @@vector_options.interpret(options)
  xax, yax, closer, fx, fy = data_prep_2D([fx,fy], newframe, opts)
  if newframe
    fig(xax, yax)
  end
  if newframe
    axes_or_map_and_ttl(fx, opts, xax, yax, fy)
  end
  itr = current_itr
  xaxv = xax.val
  yaxv = yax.val
  DCL.uwsgxa(xaxv)
  DCL.uwsgya(yaxv)
  if opts['flow_itr5']
    if itr_is?(5)
      DCLExt.flow_itr5( fx, fy, opts['factor'], opts['unit_vect'] )
    else
      raise "flow_itr5 option should use with itr=5."
    end
  elsif opts['flow_vect_anyproj'] ||
        ( opts['flow_vect_anyproj'].nil? && opts['flow_vect'] && itr >=2 )
    @@vfratio, flenmax = 
        DCLExt.flow_vect_anyproj(fx.val, fy.val, xaxv, yaxv, 
                                 opts['factor'], 1, 1, opts['distvect_map'],
                                 (opts['keep'] && @@vfratio),
                                 opts['polar_thinning'] )
    if opts['unit_vect']
      len_unit = opts['len_unit'] || flenmax
      DCLExt.unit_vect_single(@@vfratio, len_unit)
    end
  elsif opts['flow_vect']
    uninfo = DCLExt.flow_vect(fx.val, fy.val, opts['factor'], 1, 1,
              (opts['keep']&& @@vxfxratio), (opts['keep'] && @@vyfyratio) )
    @@vxfxratio, @@vyfyratio, = uninfo
    if opts['unit_vect']
      unless opts['max_unit_vect']
        uninfo[2] = opts['ux_unit']  # nil or a Numeric to specify the len 
        uninfo[3] = opts['uy_unit']  # nil or a Numeric to specify the len 
      end
      DCLExt.unit_vect(*uninfo)
    end
  else
    before=DCLExt.ug_set_params({'lunit'=>true}) if opts['unit_vect']
    DCL.ugvect(fx.val, fy.val)
    DCLExt.ug_set_params(before) if opts['unit_vect']
  end
  closer.call if closer
  nil
end

.vector_on_merdional_section(fx, fy, newframe = true, options = nil) ⇒ Object

# File 'lib/numru/ggraph_on_merdional_section.rb', line 72

def vector_on_merdional_section(fx, fy, newframe=true, options=nil)  
  if ! defined?(@@vector_on_merdional_section_options)
    @@vector_on_merdional_section_options = Misc::KeywordOptAutoHelp.new(
      ['newfig', true, 'if false, do not cleared before figure setting.'],
      ['title', nil, 'Title of the figure(if nil, internally determined)'],
      ['annotate', true, 'if false, do not put texts on the right margin even when newframe==true'],
      ['transpose', false, 'if true, exchange x and y axes'],
      ['flow_vect', true, 'If true, use DCLExt::flow_vect to draw vectors; otherwise, DCL::ugvect is used.'],
      ['xintv', 1, '(Effective only if flow_vect) interval of data sampling in x'],
      ['yintv', 1, '(Effective only if flow_vect) interval of data sampling in y'],
      ['factor', 1.0, '(Effective only if flow_vect) scaling factor to strech/reduce the arrow lengths'],
      ['unit_vect', false, 'Show the unit vector'],
      ['use_before_scale', false, '(Effective only unless flow_vect) If true, use scale factor before vector.'],
      ['max_unit_vect', false, '(Effective only if flow_vect && unit_vect) 
        If true, use the maximum arrows to scale the unit vector; otherwise, normalize in V coordinate.']
    )
  end
  opts = @@vector_on_merdional_section_options.interpret(options)
  fx = fx.first2D.copy
  fy = fy.first2D.copy
  sh = fx.shape
  if sh != fy.shape
    raise ArgumentError, "shapes of fx and fy do not agree with each other"
  end
  fx = fx.transpose(1,0) if opts['transpose']
  fy = fy.transpose(1,0) if opts['transpose']
  if ((xi=opts['xintv']) >= 2)
    idx = NArray.int(sh[0]/xi).indgen!*xi     # [0,xi,2*xi,..]
    fx = fx[idx, true]
    fy = fy[idx, true]
  end
  if ((yi=opts['xintv']) >= 2)
    idx = NArray.int(sh[1]/yi).indgen!*yi     # [0,yi,2*yi,..]
    fx = fx[true, idx]
    fy = fy[true, idx]
  end
  xax = fx.coord(0)
  yax = fy.coord(1)
  aphi_ax, z_ax, was_proportional_to_p = \
        GPhys::EP_Flux::preparate_for_vector_on_merdional_section(xax, yax)
  if was_proportional_to_p
	itr = 2 
  else
	itr = 1
  end
  if newframe
	nextfig = @@next_fig.dup  if ( @@next_fig != nil ) # backup next_fig
    fig(xax, yax, {'itr'=>itr}) 
    axes(xax, yax)
    if opts['title']
      ttl = opts['title']
    else
      fxnm = fx.data.get_att('long_name') || fx.name
      fynm = fy.data.get_att('long_name') || fy.name
      ttl =   '('+fxnm+','+fynm+')'
    end
    title( ttl )
    annotate(fx.lost_axes) if opts['annotate']
	@@next_fig = nextfig      if ( @@next_fig != nil )
  end
  fig(aphi_ax, z_ax, {'new_frame'=>false, 'itr'=>1, 'yreverse'=>false}) \
                                                      if (opts['newfig'])
  DCL.uwsgxa(aphi_ax.val) 
  DCL.uwsgya(z_ax.val)
  if opts['flow_vect']
	if opts['use_before_scale']
	  vxfxratio = @@uninfo[0]; vxfyratio = @@uninfo[1]
	  before=DCLExt.ug_set_params( {'LNRMAL'=>false, 'LMSG'=>false,
	 'XFACT1'=>1.0, 'YFACT1'=>1.0} )
	  DCL.ugvect(vxfxratio*fx.val, vxfyratio*fy.val)
	  DCLExt.ug_set_params(before) 
	  if opts['unit_vect']
 if opts['max_unit_vect']
   DCLExt.unit_vect(*@@uninfo)
 else
   DCLExt.unit_vect(*@@uninfo[0..1])
 end
	  end
	else
	  @@uninfo = DCLExt.flow_vect(fx.val, fy.val, opts['factor'] )
	  if opts['unit_vect']
 if opts['max_unit_vect']
   DCLExt.unit_vect(*@@uninfo)
 else
   DCLExt.unit_vect(*@@uninfo[0..1])
 end
	  end
	end
  else
	if opts['use_before_scale']
	  vxfxratio = @@uninfo[0]; vxfyratio = @@uninfo[1]
	  before1=DCLExt.ug_set_params({'lunit'=>true}) if opts['unit_vect']
	  before2=DCLExt.ug_set_params( {'LNRMAL'=>false, 'LMSG'=>false,
	 'XFACT1'=>1.0, 'YFACT1'=>1.0} )
	  DCL.ugvect(vxfxratio*fx.val, vxfyratio*fy.val)
	  DCLExt.ug_set_params(before1) if opts['unit_vect']
	  DCLExt.ug_set_params(before2) 
	else
	  before=DCLExt.ug_set_params({'lunit'=>true}) if opts['unit_vect']
	  DCL.ugvect(fx.val, fy.val)
	  DCLExt.ug_set_params(before) if opts['unit_vect']
	end
	nil
  end
end