(例) p(x, y) : x とy の2次元データ
GPhys使用
pbar = p.avg(0)
b = (p – pbar) / pbar
# ↑ 次元数が異なる
b.copy(NetCDF.create(filename))
GPhys不使用
物理量同士の計算
pval = p.get
pbarval = pval.avg(0)
bval = Narray.float(*pval.shape[1..-1])
for i in 0..(pval.shape[1]-1)
bval = (pval[i, true] – pbarval)/ pbarval
newfile = NetCDF.create(newfilename)
newfile.def_var(“x”, “float”, [xd])
xd = newfile.def_dim(“x”, p.shape[0])
yd = newfile.def_dim(“y”, p.shape[1])
自己記述性的ファイルを手作業で構成すると相当な行数になる
end
ループが必要
(まだまだつづく)