[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[dennou-ruby:002418] Global attribute



竹広です. 

以前 Global attribute の扱いを質問しまして, それに沿って
以下のようなコード(抜粋)をかいて gp* コマンドに
とりいれてみました. 

入力ファイルが NetCDF でない場合にも対処したつもりなんですが
いかがなものでしょうか. 

それから入力が複数ファイルの場合には gphys.data.file の返り値が
NetCDF Class でなくファイル名の配列になってしまうので, 
おなじ方法がとれません. その場合には別途一つのファイルを開いて
Global attribute を読み取るぐらいでいいのでしょうか. 
もっと洒落た方法があったら教えて下さい. 

                    Takepiro(竹広真一)@数理解析研究所. 京都大学
                      	E-mail:takepiro@xxxxxxxxxxxxxx
			       takepiro@xxxxxxxxxxxxxxxxxxxx

======================================================================
gphys = GPhys::IO.open_gturl(gturl)    # 入力

outncfile=NetCDF.create(outfilename)   # 出力

# Copy global attributes (only for NetCDF file)
ncfile = gphys.data.file 
if /NetCDF/  =~ ncfile.class.to_s
  ncfile.each_att{|att|
    att.copy(outncfile)
  }
end

...

GPhys::IO.write( outncfile, gphys )
NetCDF_Conventions.add_history(outncfile, command) # Exec command written as a history

outncfile.close
======================================================================