NCEPx4

(2004/04 更新)

• フォートラン用データ使用法

  データは/home9/wnabe/ncepx4の中に入っており、月ごとに一つのファイルになっています。
  読み方は、例えばnd日の4つの時間の温度データをtという数列に入れることを考えると、
  
  

  	OPEN (1,FILE='air.198601.data',STATUS='OLD',  
       +   ACCESS='SEQUENTIAL',FORM='UNFORMATTED')  
    
          do 30 id=1,nd
          do 30 it=1,4 
  	do 30 k=1,17
  	read(1)((t(i,j,k,it),i=1,lx),j=1,ly)
  30	continue
  
  	close(1)
  

  
  こんな感じです。
  
  

• データダウンロード、書き換え

  オリジナルデータはNetCDF形式で、一年ごとに配布されています。
  フォートランでNetCDFのファイルを読む場合、いくつかのライブラリを読み込む必要があり、それらはvortexとomegaにしか入っていません。
  なのでこれらのマシンにログインしてftpや解凍をしましょう。
  オリジナルのデータは/home9/wnabe/ncepx4.orgに入っています。
  なのでここに入って、
  ftp ftp.cdc.noaa.gov として、
  pub/Datasets/ncep.reanalysis/pressure/と進んでget。
  
  
  これをフォートラン用データに書き換えます。
  /home9/wnabe/ncepx4.org/readncepの中に、read_ncep.fがあります。
  この中身を欲しいデータ用に書き換えて、コンパイルして実行します。
  書き換える内容は主に3つ。read_ncep.readmeに簡単に書いてます。
  Wの時は鉛直層数17から12に減ります。
  U,V,H,T,W,OLR,毎にオフセットを変えます。
  OLRの場合はデータが空間２次元になります。
  
  コンパイルは、
  cdffrt read_ncep.f
  続いて実行して下さい。月ごとのデータが/home9/wnabe/ncepx4に書き込まれます。