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

 $local_path = '/home/yukiko/work/ape/yukiko/lib'
$: << $local_path

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

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

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

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

# ----------------------------------------------
groupid = ["aguforape","agcm5_adj","agcm5_kuo"]
host = "eva01"
set_dir_id  = groupid[0]

END{

sstid = ["control"]

rezolid = [
#  "T39L48_non", "T39L24_non","T39L96_non", "T79L48_non", "T159L48_non",
#  "T39L48_eml", "T39L24_eml","T39L96_eml", "T79L48_eml", 
"T159L48_eml"#,
#  "T39L48_ias"#, "T39L48_ksc", "T39L48_kuo", "T39L48_mca" 
]

  rezolid.each { |item|
    $expID = "control"
    $rezol = item
    set_dir(set_dir_id,host)
    #  a.nc_tr_wkplot_mono
    a.nc_tr_wkplot

  }


}

class Ape_mkfig

  def nc_tr_wkplot
    
    @data = Ape.new("/home/yukiko/work/ape/yukiko/data/wk-stspct-qtuz/#{$rezol}_#{$expID}_wk_40smooth_qtuz.nc")

    @data.netcdf_open.var_names.each { |name|
      unless name == "wvn" || name == "freq" || name == "freq0" || name == "noname" 
	if name =~ /lw_toa|mslp|ps/
	  lost_axis = []
	elsif name =~ /250/
	  lost_axis = ["sigma=0.25"]
	elsif name =~ /500/
	  lost_axis = ["sigma=0.50"]
	elsif name =~ /170/
	  lost_axis = ["sigma=0.17"]
	elsif name =~ /570/
	  lost_axis = ["sigma=0.57"]
	elsif name =~ /950/
	  lost_axis = ["sigma=0.95"]
	elsif name =~ /850/
	  lost_axis = ["sigma=0.85"]
	elsif name =~ /760/
	  lost_axis = ["sigma=0.76"]
	end
	
	# 時空間スペクトル
	gphys = @data.gphys_open(name). 
	  add_lost_axes(lost_axis)
	dim1 = gphys.coord(0).shape_current.to_s.to_i
	dim2 = gphys.coord(1).shape_current.to_s.to_i

#	if name =~ /_asym_org_spct/
#	  gphys = gphys.rename(name.sub("asym", "sym"))
#	elsif name =~ /_sym_org_spct/
#	  gphys = gphys.rename(name.sub("sym", "asym"))
#	end
#	gphys = gphys.rename(name.sub("ps", "mslp") )

        puts name
	if name =~ /bg_spct/ 
	  mkfig_plot(gphys[((dim1+1)/2-31)..((dim1+1)/2+29),0..80].log.
                       add_lost_axes("logarithm of data") )
	else
	  mkfig_plot(gphys[((dim1+1)/2-31)..((dim1+1)/2+29),0..80])
	end

      end
    }
    
  end

  def nc_tr_wkplot_mono
    
#    @data = Ape.new("/work11/ape/yukiko/data/wk-stspct/#{$rezol}_#{$expID}_wk_40smooth.nc")
    @data = Ape.new("/home/yukiko/eva01/work11/ape/yukiko/data/wk-stspct/#{$rezol}_#{$expID}_wk_40smooth.nc")
	
    # 時空間スペクトル
    gphys = @data.gphys_open("tr_tppn_sym_spct").
      rename("tr_tppn_sym_spct_mono")
    dim1 = gphys.coord(0).shape_current.to_s.to_i
    dim2 = gphys.coord(1).shape_current.to_s.to_i

    mkfig_plot(gphys[((dim1+1)/2-31)..((dim1+1)/2+29),0..80])
    
  end

end