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

# $local_path = '/work11/ape/yukiko/lib'
 $local_path = '/GFD_Dennou_Work2/yukiko/eva01/work11/ape/yukiko/lib'
# $local_path = '/home/yukiko/tmp/ape-data/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", "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]

$cumulus = "eml"
$rezol = "T39L48_#{$cumulus}"
$expID = "control"
set_dir(set_dir_id,host)


class Ape_mkfig
  # #{$groupid}_TR_control.nc
  def nc_tr_mono
    @data = Ape.new("#{$ncfile_path}#{$groupid}_TR_#{$expID}.nc")
    
    # x-t ダイヤグラム
    gphys = @data.gphys_open("tr_tppn")
    gphys = gphys.cut(true,0,true)[true,-401..-1]
    
    grid_0 = gphys.grid_copy.axis(0)
    grid_1 = Axis.new().set_pos(VArray.new(NArray.sfloat(401).indgen!/4).rename("time").put_att("units","days"))
    grid = Grid.new(grid_0,grid_1)
    gphys = GPhys.new(grid, gphys.data )
    mkfig_plot(gphys.lon_lotate.rename("tr_tppn_mono").set_lost_axes(""))
  end


  def nc_sst_mono
#    $ncfile_path = "/work11/ape/NetCDF/sstncfiles/"
    $ncfile_path = "/GFD_Dennou_Work2/yukiko/eva01/work11/ape/NetCDF/sstncfiles/"

    $expID = "control"
    $file_label = ""
    gphys_ary = []
    $rezol = "T39"
    @data = Ape.new("#{$ncfile_path}#{$rezol}_SST.nc")
    gphys = @data.gphys_open("sst_#{$expID}")
    gphys = gphys.
      set_att("ape_name","sea surface temperature").
      set_att("line_name","").mean(0).
      rename("sst_zonal").set_lost_axes("")
    mkfig_plot(gphys)

  end

end


$cumulus = "eml"
$rezol = "T39L48_#{$cumulus}"
$expID = "control"
set_dir(set_dir_id,host)

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

rezolid.each { |item|
$cumulus = "eml"
$rezol = item
$expID = "control"
set_dir(set_dir_id,host)
#a.nc_tr_mono
a.nc_sst_mono
}