#!/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"


# gif, netcdf の格納場所
$fig_path = "/work11/yukiko/waveCISK/figs/tmp/"
$ncfile_path = "/work11/yukiko/spmodel-yukiko/data/"


#----------------------------------------------
# ユーザ設定

# 鉛直層 
$layer = 4

# 加熱型 0:linear, 1:positive
$heat = 0

# 加熱幅 0: control, 1: 4times X half-width
$width = 0

# 変数 height, u, v
$variable = "height"



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

END{

  # --
  $layer = 4; $heat = 0; $width = 0
  set_filename
  t = ape_new($filename)
  g = t.go($variable)
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  File.rename("#{$fig_path}_#{g.name}.gif", "#{$fig_path}#{$gifname}")  

  # --
  $layer = 4; $heat = 1; $width = 0
  set_filename
  t = ape_new($filename)
  g = t.go($variable)
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  File.rename("#{$fig_path}_#{g.name}.gif", "#{$fig_path}#{$gifname}")
  
  # --
  $layer = 10; $heat = 0; $width = 0
  set_filename
  t = ape_new($filename)
  g = t.go($variable)
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  File.rename("#{$fig_path}_#{g.name}.gif", "#{$fig_path}#{$gifname}")

  # --
  $layer = 10; $heat = 1; $width = 0
  set_filename
  t = ape_new($filename)
  g = t.go($variable)
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  File.rename("#{$fig_path}_#{g.name}.gif", "#{$fig_path}#{$gifname}")

  # --
  $layer = 4; $heat = 0; $width = 1
  set_filename
  t = ape_new($filename)
  g = t.go($variable)
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  File.rename("#{$fig_path}_#{g.name}.gif", "#{$fig_path}#{$gifname}")

  # --
  $layer = 4; $heat = 1; $width = 1
  set_filename
  t = ape_new($filename)
  g = t.go($variable)
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  File.rename("#{$fig_path}_#{g.name}.gif", "#{$fig_path}#{$gifname}")

  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  t.mkdump(g.cut(true,0,true,true)[true,0,true].
	   set_att("ape_name","#{g.data.get_att("long_name")}").
	   set_lost_axes($lost_axis))
  File.rename("#{$fig_path}_#{g.name}.gif", "#{$fig_path}#{$gifname}")

  t.grcls


}

#----------------------------------------------
# 設定

def set_filename

  heat_ary = ["linear","positive"]
  width_ary = ["","-4width"]

  $filename = 
    "#{$ncfile}wavecisk-#{$layer}layer-#{heat_ary[$heat]}#{width_ary[$width]}.nc"

  print "#{$filename}@#{$variable}\n"

  $gifname = 
    "wavecisk-#{$layer}layer-#{heat_ary[$heat]}#{width_ary[$width]}-xteqz1-#{$variable}.gif"

  $file_label = ""
  
  heat_ary = ["linear","positive_only"]
  
  sign = "YU.YAMADA #{Time.now.strftime("%Y/%m/%d %H:%M:%S JST")}"
  
  $lost_axis = [
    sign,
    "y=0, z=1", 
    "rezol: k=l=64, z=#{$layer}", 
    "#{heat_ary[$heat]}_heat, linear_equation", 
    "beta=1, eta1 = 1.5, eta2 = -1.5", 
  ]

  $lost_axis.push("init-heat: 4 times X-width") if $width == 1

end