require "numru/netcdf"
include NumRu

file = NetCDF.create("test.nc")      # 新規作成
nx, ny = 10, 5
file.def_dim("x",nx)                 # 次元の定義
file.def_dim("y",ny)
file.def_dim("t",0)                  # tはUNLIMITED
require "date"
file.put_att("history","created by #{$0}  #{Date.today}")
                                     # グローバル属性の設定
x = file.def_var("x","sfloat",["x"]) # 変数の定義
y = file.def_var("y","sfloat",["y"])
t = file.def_var("t","sfloat",["t"])
v1 = file.def_var("v1","sfloat",["x","y"])
v1.put_att("long_name","test 1")     # 属性の設定
v1.put_att("units","1")
v2 = file.def_var("v2","sfloat",["x","y","t"])
v2.put_att("long_name","test 2")
v2.put_att("units","1")
file.enddef                          # defineモード終わり


x.put( NArray.float(nx).indgen! )    # 値を入れる
y.put( NArray.float(ny).indgen! )

z = NArray.float(nx,ny).indgen!*0.1
v1.put(z)
v1.put( NArray.float(nx).fill!(20), "start"=>[0,2],"end"=>[-1,2])
                                     # startからendまで値を入れる
v2.put(z, "start"=>[0,0,0],"end"=>[-1,-1,0])
t.put( 0, "index"=>[0])              # indexの場所に値を入れる
v2.put(-z, "start"=>[0,0,1],"end"=>[-1,-1,1])
t.put( 1, "index"=>[1])

file.close
print `ncdump test.nc`               # できたファイルを見る