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

file =  "./ek_analysis/ek-analysis.nc"
# ----- 使用データ -----
tro_term1 = (GPhys::IO.open file, 'tro_epsilon1').val
tro_term2 = (GPhys::IO.open file, 'tro_epsilon2').val 
tro_term3 = (GPhys::IO.open file, 'tro_epsilon3').val 
tro_term4 = (GPhys::IO.open file, 'tro_epsilon4').val 
tro_term5 = (GPhys::IO.open file, 'tro_epsilon5').val 
cli_term1 = (GPhys::IO.open file, 'cli_epsilon1').val 
cli_term2 = (GPhys::IO.open file, 'cli_epsilon2').val 
cli_term3 = (GPhys::IO.open file, 'cli_epsilon3').val 
cli_term4 = (GPhys::IO.open file, 'cli_epsilon4').val 
cli_term5 = (GPhys::IO.open file, 'cli_epsilon5').val 
cli_term6 = (GPhys::IO.open file, 'cli_epsilon6').val 
cli_term7 = (GPhys::IO.open file, 'cli_epsilon7').val 

print "Energy dissipation rate \n"
print "--- barotropic mode --- \n"
print "term I   = ",tro_term1[0],"\n"
print "term II  = ",tro_term2[0],"\n"
print "term III = ",tro_term3[0],"\n"
print "term IV  = ",tro_term4[0],"\n"
print "term V   = ",tro_term5[0],"\n"
print "--- baroclinic mode --- \n"
print "term I   = ",cli_term1[0],"\n"
print "term II  = ",cli_term2[0],"\n"
print "term III = ",cli_term3[0],"\n"
print "term IV  = ",cli_term4[0],"\n"
print "term V   = ",cli_term5[0],"\n"
print "term VI  = ",cli_term6[0],"\n"
print "term VII = ",cli_term7[0],"\n"