#!/usr/bin/env ruby
# ----------------------------------------------
# local load path

# $local_path = '/work11/ape/yukiko/lib'
# $local_path = '/home/yukiko/tmp/ape-data/lib'
 $local_path = '/home/yukiko/work/ape/yukiko/lib'
$: << $local_path

# ----------------------------------------------
# 必要なライブラリ, モジュールの読み込み

load "#{$local_path}/ape-view.rb"

# ----------------------------------------------

class  Ape 
  def __show_line_index(str_ary,line_index_ary, line_type_ary, x=0.15, y=0.1)
    charsize = 0.8 * DCL.uzpget('rsizec1')
    len = 0.04     # 線の長さ
    dvx = 0.03
    dvy = charsize*1.6
    
    raise TypeError,"Array expected" if ! str_ary.is_a?(Array)
    vxmin,vxmax,vymin,vymax = DCL.sgqvpt
    vx = 0.15
    vy = 0.13 - charsize/2
#    vx = x
#    vy = y - charsize/2
    str_ary.size.times{|num|
      DCL::sgplzv([vx, vx + len],[vy, vy], line_type_ary[num],
		  line_index_ary[num])
      DCL::sgtxzv(vx + len + 0.01 ,vy,str_ary[num],
                  charsize, 0, -1, 13 )
      vy -= dvy
      if num == 2
        vx = 0.34
        vy = 0.13 - charsize/2
      elsif num == 5
        vx = 0.53
        vy = 0.13 - charsize/2
      end
    }
    nil
  end

end


# ---------

dirpath = "/home/yukiko/work/yukiko/eva01-home/work/aqua-mradl-ronbun/composite_nc/"
$file_label= ""
$fig_path = "/home/yukiko/work/yukiko/eva01-home/work/aqua-mradl-ronbun/figs/tmp/"
apename = "sea surface temperature"
lost_axis = ["(mean) lon:0..357.188"]


a = ape_new("#{dirpath}/sstcon.nc")
gphys = a.go("con")[true,true,0].
  set_att("units","K").
  set_att("ape_name",apename). 
  rename("sst-H1998").
  set_lost_axes("")

gphys_mean = gphys.mean(0).rename("sst-H1998_zonal").
  set_lost_axes(lost_axis)


#a.plot(gphys) 
a.mkdump(gphys)  ; a.mkdump(gphys_mean) 
#a.mkps(gphys) ; a.mkps(gphys_mean) 
a.grcls