# -*- coding: utf-8 -*-
require "numru/ggraph"
include NumRu

ek_file  = "./ek_analysis/ek-analysis.nc"
ens_file = "./ens_analysis/ens-analysis.nc"

#----- 変数定義 -----
hvisc = 2.017 * 10**(-19)

#----- 使用データ -----
ek_diss_rate_tro  = ((GPhys::IO.open ek_file, 'tro_epsilon5').val.abs)[0] 
ek_diss_rate_cli  = ((GPhys::IO.open ek_file, 'cli_epsilon7').val.abs)[0] 
ek_diss_rate_total = ek_diss_rate_tro + ek_diss_rate_cli

ens_diss_rate_tro = ((GPhys::IO.open ens_file, 'tro_eta5').val.abs)[0] 
ens_diss_rate_cli = ((GPhys::IO.open ens_file, 'cli_eta7').val.abs)[0] 
ens_diss_rate_total = ens_diss_rate_tro + ens_diss_rate_cli
 
#----- 下層の散逸率の計算 ---
ek_diss_rate_low = ek_diss_rate_tro - ek_diss_rate_cli

#--------- エネルギー散逸波数計算 ---------
ek_diss_tro = ( ek_diss_rate_tro**(1.0/3.0) / hvisc )**(3.0/22.0)
ek_diss_cli = ( ek_diss_rate_cli**(1.0/3.0) / hvisc )**(3.0/22.0)
ek_diss_total = ( ek_diss_rate_total**(1.0/3.0) / hvisc )**(3.0/22.0)
ek_diss_low = ( ek_diss_rate_low**(1.0/3.0) / hvisc )**(3.0/22.0)
#--------- エンストロフィー散逸波数計算 ---------
ens_diss_tro = (ens_diss_rate_tro**(1.0/3.0) / hvisc)**(1.0/8.0)
ens_diss_cli = (ens_diss_rate_cli**(1.0/3.0) / hvisc)**(1.0/8.0)
ens_diss_total = (ens_diss_rate_total**(1.0/3.0) / hvisc)**(1.0/8.0)

print "ek_diss_tro=",ek_diss_tro,"\n"
print "ek_diss_cli=",ek_diss_cli,"\n"
print "ek_diss_total=",ek_diss_total,"\n"
print "ens_diss_tro=",ens_diss_tro,"\n"
print "ens_diss_cli=",ens_diss_cli,"\n"
print "ens_diss_total=",ens_diss_total,"\n"