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

 $local_path = '/home/yukiko/work/ape/yukiko/lib'
# $local_path = '/home/yukiko/tmp/ape-data/lib'
$: << $local_path

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

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

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


END{

  mknc_file_marge

}

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


def mknc_file_marge

  $gr2ncfile_path = "./tmp/"
  file =  
    NetCDF.create("ECMWF_control_wk_40smooth_qtuz.nc")
#    NetCDF.create("CSIRO_control_wk_40smooth_qtuz.nc")
  puts "ECMWF_control_wk_40smooth_qtuz.nc create start"

  file_name = Array.new
  Dir.foreach($gr2ncfile_path) { |item|
    file_name.push("#{$gr2ncfile_path}/#{item}") if item =~ /ECMWF_cont/
#    file_name.push("#{$gr2ncfile_path}/#{item}") if item =~ /CSIRO_/
  }

  file_name.each{ |moto|
    
    puts moto
    $t = Ape.new(moto)

    $t.netcdf_open.var_names.each { |name| 
      if name =~ /tr_/
        puts name
        gphys = $t.go(name)
        
        # netCDF 初期化
        GPhys::NetCDF_IO.write(file, gphys )
      
      end
    }      

  }

  # 大域属性
  file.put_att("Conventions", "CF-1.0")
  file.put_att("title","Aqua Planet: ECMWF_control Experiment, Wheeler and Kiladis, 1999 Plot")
#  file.put_att("title","Aqua Planet: CSIRO_control Experiment, Wheeler and Kiladis, 1999 Plot")
  file.put_att("history", "Original data produced: #{Time.now} from TR netCDF by Yukiko YAMADA (AGU for APE)")
  file.put_att("institution", "ECMWF")
#  file.put_att("institution", "CSIRO")
  file.close
  
end