#!/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")
sstid.push("H1998con")
sstid.push("H1998pa")

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]


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


def set_directry_agcm5
  $groupid = $groupid_hash[$rezol]
  $ncfile_path = "/work11/ape/NetCDF/#{$rezol}/"
  $file_label = "#{$rezol}_#{$expID}"
end

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

#  a.nc_sh_zonal_comparing_rezol
  $expID = sstid[num]
#  set_dir(set_dir_id,host)
#  a.nc_gt_comparing_rezol
#  a.nc_ml_zonal_comparing_rezol

}


["1keq", "3keq", "3kw1", "flat3keq", "Qobs3keq", "H1998pa"].each{ |expid|
  $expID = expid
  set_dir(set_dir_id,host)
#   a.nc_sh_anm_comparing_rezol

}

#["AGCM5_adj", "AGCM5_kuo", "AGCM5_kuo-non"].each{ |rezolid| 
["AGCM5_kuo-nosc"].each{ |rezolid| 

  $expID = "control"
  $rezol = rezolid
  set_directry_agcm5

  puts $groupid 

  a.nc_sh_anm_comparing_expid
  a.nc_ml_zonal_comparing_expid
  a.nc_gt_comparing_expid  
  a.nc_sh_zonal_comparing_expid 

}




=begin

# 0--5
drowid = 5

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

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


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


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

  a.nc_sh_zonal_comparing_rezol
  $expID = sstid[num]
  set_dir(set_dir_id,host)
  a.nc_gt_comparing_rezol
  a.nc_ml_zonal_comparing_rezol

}


["1keq", "3keq", "3kw1", "flat3keq", "Qobs3keq", "H1998pa"].each{ |expid|
  $expID = expid
  set_dir(set_dir_id,host)
   a.nc_sh_anm_comparing_rezol

}

["AGCM5_adj", "AGCM5_kuo", "AGCM5_kuo-non"].each{ |rezolid| 

  set_dir_id  = rezolid
  $expID = "control"
  set_dir(set_dir_id,host)

  a.nc_sh_anm_comparing_expid
  a.nc_ml_zonal_comparing_expid
 a.nc_gt_comparing_expid  
 a.nc_sh_zonal_comparing_expid 

}

=end