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

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

sstid = ["control","flat","peaked","1keq","3keq","3kw1","Qobs","control-5N"]
#sstid = ["control","flat","Qobs","3keq","flat3keq","Qobs3keq",
#  "H1998con","H1998pa"]
cumulusid = ["eml","ias","kuo","mca","non"]

$rezol = "T39L48_eml"
$ncfile_path = "/home/yukiko/tmp/ape-data/NetCDF/#{$rezol}/"
$ncfile_write_path = "/home/yukiko/tmp/ape-data/NetCDF/ohfuchi/#{$rezol}/"
$expID = "control"

sstid.each{ |item|

$expID = item

t = ape_new("#{$ncfile_path}AGUforAPE-03a_ML_#{$expID}.nc")
gphys = t.go("ml_v")
gphys = gphys.mean(0).strm_function("ml_psi").
  set_att("ape_name","mass_stream_function").
  set_att("units","Kg s-1").
  set_att("cell_methods","time: mean (interval: 6 hours comment: averaged over 3-year period following 6-month spin-up) longitude: mean")

f=NetCDF.create("#{$ncfile_write_path}AGUforAPE_#{$rezol}_#{$expID}_ml_psi.nc")
GPhys::NetCDF_IO.write(f,gphys)
f.put_att("Conventions", "CF-1.0")
f.put_att("title", "Aqua Planet: Multiple-Level 3-D Means from \'control\' Experiment")
f.put_att("history","Original data produced: 2004/01/22 on ES; area-weighted global means calculated for APE ML standard output\n#{Time.now} yukiko> calculated from ml_v on #{$rezol}/AGUforAPE-03a_ML_#{$expID}.nc")
f.put_att("institution","AGU for APE; AFES Working Team, GFD Dennou Club, University of Tokyo")
f.put_att("source","AFES v1.15")
f.close

}