#!/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 = [ "AGUforAPE", "ECMWF" , "GSFC", "LASG", "NCAR"]
#rezolid = [ "GSFC", "LASG", "NCAR"]
#rezolid = [ "ECMWF_07"]
rezolid = [ "ECMWF_07", "CSIRO"]

  rezolid.each { |item|
    $expID = "control"
    $rezol = item
    set_directry
    #  a.nc_tr_wkplot_mono
    a.nc_tr_wkplot

  }


}

def set_directry
  $groupid = $groupid_hash[$rezol]
  $ncfile_path = "/home/yukiko/tr_request/#{$rezol}/"
  $file_label = "#{$rezol}_#{$expID}"
end

class Ape_mkfig

  def nc_tr_wkplot
    
    @data = Ape.new("/home/yukiko/tr_request/netcdf/wk1999/#{$rezol}_#{$expID}_wk_40smooth_qtuz.nc")

    @data.netcdf_open.var_names.each { |name|
      unless name == "wvn" || name == "freq" || name == "freq0" || name == "noname" 

#|| name == "tr_u600_sym_org_spct" || name == "tr_u600_asym_org_spct" || name == "tr_u925_bg_spct" || name == "tr_u925_sym_org_spct" || name == "tr_u925_asym_org_spct" || name == "tr_v250_bg_spct" || name == "tr_v250_sym_org_spct" || name == "tr_v250_asym_org_spct" 
	if name =~ /lw_toa|mslp|ps/
	  lost_axis = []
	elsif name =~ /150/
	  lost_axis = ["plev=15000 Pa"]
	elsif name =~ /250/
	  lost_axis = ["plev=25000 Pa"]
	elsif name =~ /600/
	  lost_axis = ["plev=60000 Pa"]
	elsif name =~ /700/
	  lost_axis = ["plev=70000 Pa"]
	elsif name =~ /850/
	  lost_axis = ["plev=85000 Pa"]
	elsif name =~ /925/
	  lost_axis = ["plev=92500 Pa"]
	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