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

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

a.nc_ml_zonal_comparing_rezol
#a.nc_sh_zonal_comparing_rezol



cumulusid = ["eml","ias","ksc","kuo","mca","non"]
sstid = ["control","flat","peaked","1keq","3keq","3kw1","Qobs","control-5N"]

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

# rezolid = ["T39L24_non","T39L96_non", "T79L48_non", "T159L48_non"]

 rezolid = ["T159L48_non"]

# rezolid.each { |item|

cumulusid.each { |item|

#  sstid_all.each { |sst|
#  sstid.each { |sst|
#    $cumulus = item
#    $cumulus = "non"
    $rezol = "T39L48_#{item}"
#    $rezol = item
#    $expID = sst
    $expID = "3keq"
    set_dir(set_dir_id,host)
#    a.nc_ml_anm("control")

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

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

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


#    a.nc_tr
#    a.nc_mf
#    a.nc_gt
#    a.nc_sh
#    a.nc_sh_zonal
#    a.nc_ml

#    a.nc_gt_anm        unless $cumulus == "non"
#    a.nc_sh_anm 
#    a.nc_sh_zonal_anm   unless $expID == "control" 
#    a.nc_ml_anm         unless $expID == "control"
#    a.nc_tr_5deg_mean 
#    a.nc_tr_300day 
#    a.nc_tr_360day

#    }
  }



rezolid = ["T39L48_non","T39L48_eml", "T39L24_eml","T39L96_eml", "T79L48_eml", "T159L48_eml", "T39L48_ias","T39L48_ksc", "T39L48_kuo", "T39L48_mca",  "T39L24_non","T39L96_non", "T79L48_non", "T159L48_non"]
rezolid.each{ |item|
  $rezol = item
  set_dir(set_dir_id,host)
#  a.nc_gt_comparing_expid  
#  a.nc_sh_zonal_comparing_expid

}

rezolid = ["T39L48_eml","T39L48_ias","T39L48_ksc", "T39L48_kuo", "T39L48_mca", "T39L48_non"]
rezolid.each{ |item|
  $rezol = item
  set_dir(set_dir_id,host)
#  a.nc_sh_anm_comparing_expid
}

sstid = ["control","flat","peaked","1keq","3kw1","Qobs","control-5N","3keq"]
sstid = ["control"]
sstid.each{ |item|
  $expID = item
  $cumulus = "eml"
  set_dir(set_dir_id,host)
#  a.nc_sh_zonal_comparing_rezol
#  a.nc_gt_comparing_rezol
  $expID = item
  $cumulus = "non"
  set_dir(set_dir_id,host)
#  a.nc_sh_zonal_comparing_rezol
#  a.nc_gt_comparing_rezol
}

sstid = ["1keq","3kw1","3keq"]
sstid.each{ |item|
  $expID = item
  $cumulus = "eml"
  set_dir(set_dir_id,host)
#  a.nc_sh_anm_comparing_rezol
  $expID = item
  $cumulus = "non"
  set_dir(set_dir_id,host)
#  a.nc_sh_anm_comparing_rezol
}



sstid_all.each { |item|
  $expID = item 
  $cumulus = "eml"
  $rezol = "T39L48_#{$cumulus}"
  set_dir(set_dir_id,host)
#  a.nc_sh_zonal_comparing_cumulus
#  a.nc_gt_comparing_cumulus
}

["1keq","3keq","3kw1","flat3keq", "Qobs3keq", "H1998pa"].each{ |item|

  $expID = item 
  $cumulus = "eml"
  $rezol = "T39L48_#{$cumulus}"
  set_dir(set_dir_id,host)
#  a.nc_sh_anm_comparing_cumulus
}




#rezolid = ["T39L48_eml"]
#rezolid.each { |item|

cumulusid.each { |item|

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

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

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

  $expID = "flat3keq"
  set_dir(set_dir_id,host)
  #a.nc_sh_anm("flat")
#  a.nc_sh_psuv_anm("flat")
  
  $expID = "Qobs3keq"
  set_dir(set_dir_id,host)
  #a.nc_sh_anm("Qobs")
#  a.nc_sh_psuv_anm("Qobs")
  
  $expID = "H1998pa"
  set_dir(set_dir_id,host)
  #a.nc_sh_anm("H1998con")
#  a.nc_sh_psuv_anm("H1998con")

  }

#a.nc_ml_tuom_anm
#a.nc_sh_psuv_anm
#a.nc_sh_tuv500_anm
#a.nc_sh_phiuv250_anm
#a.nc_sh_phiuv850_anm


rezolid =["T39L24_non","T39L96_non", "T79L48_non", "T159L48_non"]
rezolid.each { |rezol|
$cumulus = "non"
$rezol = rezol
$expID = "3keq" 
set_dir(set_dir_id,host)
#  a.nc_sh_phiuv850_anm 

  $expID = "3kw1"
  set_dir(set_dir_id,host)
#  a.nc_sh_phiuv850_anm 
  
  $expID = "1keq"
  set_dir(set_dir_id,host)
#  a.nc_sh_phiuv850_anm 

}

#a.nc_ml_tuom_anm
#a.nc_sh_psuv_anm
#a.nc_sh_tuv500_anm
#a.nc_sh_phiuv250_anm
#a.nc_sh_phiuv850_anm


["eml","non"].each { |item|

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

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

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

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

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

}




sstid.each{ |item|
  $cumulus = "non"
  $rezol = "T159L48_#{$cumulus}"
  $expID = item
  set_dir(set_dir_id,host)
#  a.nc_tr_300day 
#  a.nc_tr_5deg_mean 
#  a.nc_tr_360day
#  a.nc_sh_zonal_anm unless $expID == "control" 

if $expID == "1keq" || $expID == "3keq" || $expID == "3kw1"
#a.nc_sh_psuv_anm
#  a.nc_ml_tuom_anm
#a.nc_sh_tuv500_anm
#a.nc_sh_psuv_anm
#a.nc_sh_phiuv250_anm
#a.nc_sh_phiuv850_anm
end

}


sstid.each{ |item|
  $cumulus = "non"
  $rezol = "T79L48_#{$cumulus}"
  $expID = item
  set_dir(set_dir_id,host)
#  a.nc_tr_300day 
#  a.nc_tr_5deg_mean 
#  a.nc_tr_360day
#  a.nc_sh_zonal_anm unless $expID == "control" 

if $expID == "1keq" || $expID == "3keq" || $expID == "3kw1"
#a.nc_sh_psuv_anm
#a.nc_sh_tuv500_anm
#a.nc_sh_psuv_anm
#a.nc_sh_phiuv250_anm
#a.nc_sh_phiuv850_anm
end

}

sstid.each{ |item|
  $cumulus = "non"
  $rezol = "T39L96_#{$cumulus}"
  $expID = item
  set_dir(set_dir_id,host)
#  a.nc_tr_300day 
#  a.nc_tr_5deg_mean 
#  a.nc_tr_360day
#  a.nc_sh_zonal_anm unless $expID == "control" 

if $expID == "1keq" || $expID == "3keq" || $expID == "3kw1"
# a.nc_sh_psuv_anm
#  a.nc_ml_tuom_anm
#a.nc_sh_tuv500_anm
#a.nc_sh_psuv_anm
#a.nc_sh_phiuv250_anm
#a.nc_sh_phiuv850_anm
end

}


sstid.each{ |item|
  $cumulus = "non"
  $rezol = "T39L24_#{$cumulus}"
  $expID = item
  set_dir(set_dir_id,host)
#  a.nc_tr_300day 
#  a.nc_tr_5deg_mean 
#  a.nc_tr_360day
#  a.nc_sh_zonal_anm unless $expID == "control" 

if $expID == "1keq" || $expID == "3keq" || $expID == "3kw1"
# a.nc_sh_psuv_anm
#  a.nc_ml_tuom_anm
#a.nc_sh_tuv500_anm
#a.nc_sh_psuv_anm
#a.nc_sh_phiuv250_anm
#a.nc_sh_phiuv850_anm
end

}

sstid_all.each{ |item|
  $cumulus = "non"
  $rezol = "T39L48_#{$cumulus}"
  $expID = item
  set_dir(set_dir_id,host)
#  a.nc_tr_300day 
#  a.nc_tr_5deg_mean 
#  a.nc_tr_360day
#  a.nc_sh_zonal_anm unless $expID == "control" 

if $expID == "1keq" || $expID == "3keq" || $expID == "3kw1"
# a.nc_sh_psuv_anm
#  a.nc_ml_tuom_anm
#a.nc_sh_tuv500_anm
#a.nc_sh_psuv_anm
#a.nc_sh_phiuv250_anm
#a.nc_sh_phiuv850_anm
end

}


sstid.each{ |item|

  $expID = item
  $cumulus = "non"
  $rezol = "T159L48_#{$cumulus}"
  set_dir(set_dir_id,host)



#  a.nc_ml_zonal_comparing_rezol
#  a.nc_sh_zonal_comparing_rezol
#  a.nc_gt_comparing_rezol


#  a.nc_tr
#  a.nc_sh     unless  $expID == "3keq"
#  a.nc_gt
#  a.nc_sh_zonal unless  $expID == "3keq"
#  a.nc_ml
#  a.nc_sh_anm  unless  $expID == "3keq"

}



=begin

$cumulus = "non"
$rezol = "T39L24_#{$cumulus}"
$expID = "control"
set_dir(set_dir_id,host)
a.nc_tr_300day_spct

$cumulus = "non"
$rezol = "T39L96_#{$cumulus}"
$expID = "control"
set_dir(set_dir_id,host)
a.nc_tr_300day_spct

$cumulus = "non"
$rezol = "T79L48_#{$cumulus}"
$expID = "control"
set_dir(set_dir_id,host)
a.nc_tr_300day_spct

$cumulus = "non"
$rezol = "T159L48_#{$cumulus}"
$expID = "control"
set_dir(set_dir_id,host)
#a.nc_tr_300day_spct


DCL.grcls


#a.nc_ml_zonal_comparing_rezol
#a.nc_sh_zonal_comparing_rezol
#a.nc_gt_comparing_rezol



$cumulus = "non"
$rezol = "T39L24_#{$cumulus}"
$expID = "control"
set_dir(set_dir_id,host)
  a.nc_tr
  a.nc_gt
  a.nc_sh
  a.nc_sh_zonal
  a.nc_ml
  a.nc_sh_anm 



sstid = ["control","peaked","flat","Qobs" ]

sstid.each{ |item|

  $expID = item
  $cumulus = "non"
  $rezol = "T39L96_#{$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_sh_anm 
}

sstid = ["control","peaked","flat" ]

sstid.each{ |item|

  $expID = item
  $cumulus = "non"
  $rezol = "T79L48_#{$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_sh_anm 
}


sstid = ["control","peaked" ]

sstid.each{ |item|

  $expID = item
  $cumulus = "non"
  $rezol = "T159L48_#{$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_sh_anm 
}



sstid_all.each{ |item|

  $expID = item
  $cumulus = "eml"
  $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"

}

sstid_nocon.each { |item|

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

a.nc_sh_zonal_comparing_rezol
a.nc_gt_comparing_rezol

}




sstid.each { |item|

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

a.nc_sh_zonal_comparing_rezol
a.nc_gt_comparing_rezol

}


rezolid.each{ |item|
  $rezol = item
  set_dir(set_dir_id,host)
  a.nc_sh_zonal_comparing_expid
  a.nc_gt_comparing_expid  
}

puts "hoge"

sstid_tmp.each { |item|
  $expID = item 
  $cumulus = "eml"
  $rezol = "T39L48_#{$cumulus}"
  set_dir(set_dir_id,host)
  a.nc_sh_zonal_comparing_cumulus
  a.nc_gt_comparing_cumulus
}


sstid.each{ |item|

  $expID = item
  $cumulus = "eml"
  $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"

}


# 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 

}


sstid_nocon.each { |item|

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

a.nc_sh_zonal_comparing_rezol
a.nc_gt_comparing_rezol

}

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


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







