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

 $local_path = '/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","3keq","3kw1","Qobs","control-5N"]
sstid = ["control","flat","Qobs","3keq","flat3keq","Qobs3keq",
  "H1998con","H1998pa"]
groupid = ["aguforape","agcm5_adj","agcm5_kuo"]
cumulusid = ["eml","ias","ksc","kuo","mca","non"]
rezolid = ["T39L48_eml", "T39L24_eml","T39L96_eml", "T79L48_eml", "T159L48_eml"]

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

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

set_dir_id  = groupid[1]

$cumulus = "eml"
$expID = "control"
set_dir(set_dir_id,host)
ape_set
a.nc_tr


#a.nc_gt_comparing_rezol
#a.nc_gt_comparing_rezol

#a.nc_sh_zonal_comparing_rezol
#a.nc_sh_zonal_comparing_rezol


=begin

# control
$cumulus = "eml"

$rezol = "T39L48_#{$cumulus}"
$expID = "3keq"
set_dir(set_dir_id,host)
#a.nc_sh_anm 

cumulusid.each { |item|

$cumulus = item
$rezol = "T39L48_#{$cumulus}"
$expID = "3keq"
set_dir(set_dir_id,host)
#a.nc_sh_anm 

$expID = "flat3keq"
set_dir(set_dir_id,host)
#a.nc_sh_anm("flat")

$expID = "Qobs3keq"
set_dir(set_dir_id,host)
#a.nc_sh_anm("Qobs")

$expID = "H1998pa"
set_dir(set_dir_id,host)
#a.nc_sh_anm("H1998con")

}

$cumulus = "eml"
$expID = "3keq"

rezolid.each { |item|

$rezol = item
set_dir(set_dir_id,host)
#a.nc_sh_anm 

}

$cumulus = "eml"
$expID = "3keq"
$rezol = "T39L48_#{$cumulus}"
set_dir(set_dir_id,host)
#a.nc_sh_anm_comparing_rezol

cumulusid.each{ |item|
  $cumulus = item
  $rezol = "T39L48_#{$cumulus}"
  $expID = "3keq"
  set_dir(set_dir_id,host)
#  a.nc_sh_anm_comparing_expid
}

["3keq","flat3keq","Qobs3keq","H1998pa"].each{ |item|
  $expID = item 
  $cumulus = "eml"
  $rezol = "T39L48_#{$cumulus}"
  set_dir(set_dir_id,host)
  a.nc_sh_anm_comparing_cumulus

}


$cumulus = "eml"

$rezol = "T39L48_#{$cumulus}"
$expID = "3keq"
set_dir(set_dir_id,host)
a.nc_sh_anm 

=end

=begin
#a.nc_sh_zonal_comparing_expid
#a.nc_gt_comparing_expid  
#a.nc_sh_zonal_comparing_expid

sstid = ["control","flat","Qobs","3keq","H1998con" ]

sstid.each{ |item|

  $expID = item

  a.nc_sh_zonal_comparing_cumulus
  a.nc_gt_comparing_cumulus

  $cumulus = "ksc"
  $rezol = "T39L48_#{$cumulus}"
  set_dir(set_dir_id,host)
  
#  a.nc_tr_300day 
  a.nc_tr_5deg_mean 
  a.nc_tr
  a.nc_gt
  a.nc_sh
  a.nc_sh_zonal
  a.nc_ml
  a.nc_mf
  a.nc_gt_anm        unless $cumulus == "non"
  a.nc_sh_anm 
  a.nc_sh_zonal_anm   unless $expID == "control" || $cumulus == "non" 
  a.nc_ml_anm         unless $expID == "control"
  a.nc_mf_anm         unless $expID == "control"

}


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

#  a.nc_sh_zonal_comparing_cumulus
#  a.nc_gt_comparing_cumulus
#  a.nc_sh_zonal_comparing_cumulus

  set_dir(set_dir_id,host)
  
  a.nc_tr
  a.nc_gt
  a.nc_sh
  a.nc_sh_zonal
  a.nc_ml
  a.nc_mf
  a.nc_gt_anm("flat")        unless $cumulus == "non"
  a.nc_sh_anm("flat") 
  a.nc_sh_zonal_anm("flat")   unless $expID == "control" || $cumulus == "non" 
  a.nc_ml_anm("flat")         unless $expID == "control"
  a.nc_mf_anm("flat")         unless $expID == "control"


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

  a.nc_sh_zonal_comparing_cumulus
  a.nc_gt_comparing_cumulus
  a.nc_sh_zonal_comparing_cumulus

  $expID = "Qobs3keq"
  $cumulus = "ksc"
  $rezol = "T39L48_#{$cumulus}"
  set_dir(set_dir_id,host)
  
  a.nc_tr
  a.nc_gt
  a.nc_sh
  a.nc_sh_zonal
  a.nc_ml
  a.nc_mf
  a.nc_gt_anm("Qobs")        unless $cumulus == "non"
  a.nc_sh_anm("Qobs") 
  a.nc_sh_zonal_anm("Qobs")   unless $expID == "control" || $cumulus == "non" 
  a.nc_ml_anm("Qobs")         unless $expID == "control"
  a.nc_mf_anm("Qobs")         unless $expID == "control"


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

  a.nc_sh_zonal_comparing_cumulus
  a.nc_gt_comparing_cumulus
  a.nc_sh_zonal_comparing_cumulus


  $expID = "H1998pa"
  $cumulus = "ksc"
  $rezol = "T39L48_#{$cumulus}"
  set_dir(set_dir_id,host)
  
  a.nc_tr
  a.nc_gt
  a.nc_sh
  a.nc_sh_zonal
  a.nc_ml
  a.nc_mf
  a.nc_gt_anm("H1998con")        unless $cumulus == "non"
  a.nc_sh_anm("H1998con") 
  a.nc_sh_zonal_anm("H1998con") unless $expID == "control" || $cumulus == "non"
  a.nc_ml_anm("H1998con")         unless $expID == "control"
  a.nc_mf_anm("H1998con")         unless $expID == "control"


set_dir(set_dir_id,host)

a.nc_sh_zonal_comparing_cumulus
a.nc_gt_comparing_cumulus
a.nc_sh_zonal_comparing_cumulus


cumulusid = ["ksc"]

cumulusid.each{ |item|

  $cumulus = item

  set_dir(set_dir_id,host)
  
  a.nc_tr
  a.nc_gt
  a.nc_sh
  a.nc_sh_zonal
  a.nc_ml
  a.nc_mf
  a.nc_gt_anm        unless $cumulus == "non"
  a.nc_sh_anm 
  a.nc_sh_zonal_anm   unless $expID == "control" || $cumulus == "non" 
  a.nc_ml_anm         unless $expID == "control"
  a.nc_mf_anm         unless $expID == "control"

}

# ------

$expID = sstid[1]

set_dir(set_dir_id,host)

a.nc_sh_zonal_comparing_cumulus
a.nc_gt_comparing_cumulus
a.nc_sh_zonal_comparing_cumulus

cumulusid = ["ksc"]

cumulusid.each{ |item|

  $cumulus = item

  set_dir(set_dir_id,host)
  
  a.nc_tr
  a.nc_gt
  a.nc_sh
  a.nc_sh_zonal
  a.nc_ml
  a.nc_mf
  a.nc_gt_anm        unless $cumulus == "non"
  a.nc_sh_anm 
  a.nc_sh_zonal_anm   unless $expID == "control" || $cumulus == "non" 
  a.nc_ml_anm         unless $expID == "control"
  a.nc_mf_anm         unless $expID == "control"

}

a.nc_tr

=end

=begin


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

a.nc_sh
a.nc_sh_anm("flat")
a.nc_sh_zonal
a.nc_sh_zonal_anm("flat")
a.nc_gt
a.nc_gt_anm("flat")

$expID = "Qobs3keq"
set_dir(set_dir_id,host)
a.nc_sh
a.nc_sh_anm("Qobs")
a.nc_sh_zonal
a.nc_sh_zonal_anm("Qobs")
a.nc_gt
a.nc_gt_anm("Qobs")

$expID = "H1998pa"
set_dir(set_dir_id,host)
a.nc_sh
a.nc_sh_anm("H1998con")
a.nc_sh_zonal
a.nc_sh_zonal_anm("H1998con")
a.nc_gt
a.nc_gt_anm("H1998con")


$expID = "H1998con"
set_dir(set_dir_id,host)
a.nc_sh
a.nc_sh_anm
a.nc_sh_zonal
a.nc_sh_zonal_anm
a.nc_gt
a.nc_gt_anm
a.nc_gt_anm



# control, 3keq, flat
sstid = ["control","flat","3keq"]

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")

a.nc_gt

=end


# ------------------------------------------------
=begin

set_dir_id  = groupid[1]
$expID = sstid[0]
set_dir(set_dir_id,host)

a.nc_gt


8.times { |num|

$expID = sstid[num]
set_dir(set_dir_id,host)

a.nc_sh_zonal
a.nc_sh_zonal_anm unless num == 0
a.nc_gt_anm
a.nc_sh_anm
a.nc_ml_anm  unless num == 0

}

set_dir_id  = groupid[2]
8.times { |num|

$expID = sstid[num]
set_dir(set_dir_id,host)

a.nc_sh_zonal
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

8.times { |num|

$expID = sstid[num]
set_dir(set_dir_id,host)

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

}

set_dir_id  = groupid[2]
8.times { |num|

$expID = sstid[num]
set_dir(set_dir_id,host)

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

}

a.nc_gt


=end