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

$local_path = '/work11/ape/yukiko/lib'
$: << $local_path

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

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

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

sstid = ["control","flat","peaked","1keq","3keq","3kw1","Qobs","control-5N"]
sstid.push("HS1986")
sstid.push("flat3keq")
sstid.push("Qobs3keq")

groupid = ["aguforape","agcm5_adj","agcm5_kuo","agcm5_kuo-nosc"]

host = "eva01"
a = Ape_mkfig_agcm.new 3
# a = Ape_mkfig_agcm.new 1

set_dir_id  = groupid[1]

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


# 0--5
drowid = 4

sstid = ["1keq","3keq","3kw1","flat3keq","Qobs3keq","H1998pa"]
sstid_control =   { 
  "1keq" => "control", 
  "3keq" => "control", 
  "3kw1" => "control", 
  "flat3keq" => "flat", 
  "Qobs3keq" => "Qobs", 
  "H1998pa" => "H1998con"
}

# control
$expID = sstid[0]
set_dir(set_dir_id,host)
#a.nc_sh_phiuv850_anm(sstid_control[$expID])


["agcm5_adj","agcm5_kuo","agcm5_kuo-nosc"].each{ |groupid|

set_dir_id  = groupid

sstid.size.times{ |num|
  $expID = sstid[num]
  set_dir(set_dir_id,host)
  if drowid == 0
    a.nc_sh_phiuv850_anm(sstid_control[$expID])
  elsif drowid == 1
    a.nc_ml_tuom_anm(sstid_control[$expID])
  elsif drowid == 2
    a.nc_sh_tuv500_anm(sstid_control[$expID])
  elsif drowid == 3
    a.nc_sh_psuv_anm(sstid_control[$expID])
  elsif drowid == 4
    a.nc_sh_phiuv250_anm(sstid_control[$expID])
  elsif drowid == 5
    a.nc_sh_anm(sstid_control[$expID])
  end

}

}

# control
$expID = sstid[0]
set_dir(set_dir_id,host)
#a.nc_sh_phiuv850_anm(sstid_control[$expID])




=begin

# 0--4
drowid = 4

sstid = ["1keq","3keq","3kw1","flat3keq","Qobs3keq","H1998pa"]
sstid_control =   { 
  "1keq" => "control", 
  "3keq" => "control", 
  "3kw1" => "control", 
  "flat3keq" => "flat", 
  "Qobs3keq" => "Qobs", 
  "H1998pa" => "H1998con"
}

# control
$expID = sstid[0]
set_dir(set_dir_id,host)
a.nc_sh_phiuv850_anm(sstid_control[$expID])

sstid.size.times{ |num|
  $expID = sstid[num]
  set_dir(set_dir_id,host)
  if drowid == 0
    a.nc_sh_phiuv850_anm(sstid_control[$expID])
  elsif drowid == 1
    a.nc_ml_tuom_anm(sstid_control[$expID])
  elsif drowid == 2
    a.nc_sh_tuv500_anm(sstid_control[$expID])
  elsif drowid == 3
    a.nc_sh_psuv_anm(sstid_control[$expID])
  elsif drowid == 4
    a.nc_sh_phiuv250_anm(sstid_control[$expID])
  end

}


# control
$expID = sstid[0]
set_dir(set_dir_id,host)
a.nc_sh_phiuv850_anm(sstid_control[$expID])


=end

=begin
test.nc_ml_tuom_anm
test.nc_sh_tuv500_anm
test.nc_sh_psuv_anm
test.nc_sh_phiuv250_anm
test.nc_sh_phiuv850_anm

=end




=begin

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

set_dir_id  = groupid[3]

# control
$expID = sstid[0]
set_dir(set_dir_id,host)
a.nc_gt

# control, 3keq, flat
sstid = ["control","flat","peaked","1keq","3keq","3kw1","Qobs","control-5N","H1998con","HS1986"]

sstid.size.times{ |num|
  $expID = sstid[num]
  set_dir(set_dir_id,host)

  a.nc_gt
  a.nc_tr
  a.nc_sh_zonal
  a.nc_sh
  a.nc_ml
  
  a.nc_sh_zonal_anm unless num == 0
  a.nc_gt_anm
  a.nc_sh_anm
  a.nc_ml_anm  unless num == 0

}

a.nc_gt


# flat3keq
$expID = "flat3keq"
set_dir(set_dir_id,host)

a.nc_gt
a.nc_tr
a.nc_sh_zonal
a.nc_sh
a.nc_ml

a.nc_sh_zonal_anm("flat")
a.nc_gt_anm("flat")
a.nc_sh_anm("flat")
a.nc_ml_anm("flat")


# Qobs3keq
$expID = "Qobs3keq"
set_dir(set_dir_id,host)

a.nc_gt
a.nc_tr
a.nc_sh_zonal
a.nc_sh
a.nc_ml

a.nc_sh_zonal_anm("Qobs")
a.nc_gt_anm("Qobs")
a.nc_sh_anm("Qobs")
a.nc_ml_anm("Qobs")

# H1998pa
$expID = "H1998pa"
set_dir(set_dir_id,host)

a.nc_gt
a.nc_tr
a.nc_sh_zonal
a.nc_sh
a.nc_ml

a.nc_sh_zonal_anm("H1998con")
a.nc_gt_anm("H1998con")
a.nc_sh_anm("H1998con")
a.nc_ml_anm("H1998con")


# おしまい
a.nc_gt
a.nc_sh_anm("H1998con")

=end