[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[dennou-ruby:002418] Global attribute
- To: dennou-ruby@xxxxxxxxxxx
- Subject: [dennou-ruby:002418] Global attribute
- From: Shin-ichi Takehiro <takepiro@xxxxxxxxxxxxxx>
- Date: Mon, 22 Aug 2005 07:34:11 +0900
竹広です. 
以前 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
======================================================================