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

$local_path = '/home/yukiko/work/ape/yukiko/lib'
$: << $local_path

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

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

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

sstid = ["control","flat","peaked","1keq","3kw1","Qobs","control-5N","3keq"]
sstid_all = ["control","flat","peaked","1keq","3keq","3kw1","Qobs","control-5N","flat3keq","Qobs3keq","H1998con","H1998pa"]
sstid_nocon = ["flat","peaked","1keq","3keq","3kw1","Qobs","control-5N"]
sstid_tmp = ["control","flat","Qobs","3keq","flat3keq","Qobs3keq",
  "H1998con","H1998pa"]
groupid = ["aguforape","agcm5_adj","agcm5_kuo"]
cumulusid = ["eml","ias","ksc","kuo","mca","non"]

 rezolid = ["T39L48_non", "T39L24_non","T39L96_non", "T79L48_non", "T159L48_non"]
# rezolid = ["T39L48_non", "T39L24_non","T39L96_non", "T79L48_non"]
#rezolid = ["T39L48_eml", "T39L24_eml","T39L96_eml", "T79L48_eml", "T159L48_eml"]
#rezolid = ["T39L48_eml", "T39L24_eml","T39L96_eml", "T79L48_eml"]

#rezolid = ["T39L48_eml", "T39L24_eml","T39L96_eml", "T79L48_eml", "T159L48_eml", "T39L48_ias","T39L48_ksc", "T39L48_kuo", "T39L48_mca", "T39L48_non"]
# rezolid = ["T39L48_ias","T39L48_ksc", "T39L48_kuo", "T39L48_mca", "T39L48_non"]


host = "eva01"
 a = Ape_mkfig.new 3
# a = Ape_mkfig.new 2
# a = Ape_mkfig.new 1

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

set_dir_id  = groupid[0]


backwindid_hash = {
  "T39L48_eml" => -7.5, 
  "T39L24_eml" => -8.5, 
  "T39L96_eml" => -7.5, 
  "T79L48_eml" => -11.5, 
  "T159L48_eml"=> -11.5, 
  "T39L48_non"=> -6.0, 
  "T39L24_non"=> -5.0, 
  "T39L96_non"=> -6.0, 
  "T79L48_non"=> -6.5, 
  "T159L48_non"=> -6.5, 
  "T39L48_ias"=> -6.0, 
  "T39L48_ksc"=> -5.5, 
  "T39L48_kuo"=> -5.0, 
  "T39L48_mca"=> -5.5
}

END{


rezolid = [ 
    "T39L48_eml",
    "T39L24_eml","T39L96_eml", "T79L48_eml", "T159L48_eml", 
    "T39L48_non",
    "T39L24_non","T39L96_non", "T79L48_non","T159L48_non",
    "T39L48_ias","T39L48_ksc", "T39L48_kuo", "T39L48_mca"
]
  rezolid.each{ |item|
    $expID = "control"
    $cumulus = "non"
    $rezol = item
    
    $iryu = backwindid_hash[$rezol]

    set_dir(set_dir_id,host)
    a.nc_tr_wkplot_iryu
  }


}


class Ape_mkfig

  def nc_tr_wkplot_iryu
    
    if $rezol == "T159L48_eml"
      @data = Ape.new("/home/yukiko/work/ape/yukiko/data/wk-stspct-qtuz/#{$rezol}_#{$expID}_wk_40smooth_qtuz.nc")		
    else
      @data = Ape.new("/home/yukiko/work/ape/yukiko/data/wk-stspct/#{$rezol}_#{$expID}_wk_40smooth.nc")		
    end

    # 時空間スペクトル
    gphys = @data.gphys_open("tr_tppn_sym_spct").
      rename("tr_tppn_sym_bgw_spct").
      set_lost_axes("background wind = #{$iryu}")
    dim1 = gphys.coord(0).shape_current.to_s.to_i
    dim2 = gphys.coord(1).shape_current.to_s.to_i
    
    mkfig_plot(gphys[((dim1+1)/2-31)..((dim1+1)/2+29),0..80])
    
  end

end