以下では, 描画方法が一通り説明されていますが, エラーなく終了するかどうかの確認までできていません. ごめんなさい.
下の三つのファイルを用います.
比較のために描画する NCEP, ECMWF の再解析データ (netcdf 形式) も 準備する必要があります.
例えば NCEP のデータは, <URL:http://www.esrl.noaa.gov/psd/data/reanalysis/reanalysis.shtml> から以下のファイルをダウンロードし, draw7.5-2013-02-26-UsingU-11to20yr-pub.tgz を展開した後にできたディレクトリ draw7.5-2013-02-26-UsingU-11to20yr-pub の中に Reanal/NCEP というディレクトリ内に置いてください.
ECMWF のデータは, <URL:http://database.rish.kyoto-u.ac.jp/arch/era40/> から以下のファイルをダウンロードし, draw7.5-2013-02-26-UsingU-11to20yr-pub.tgz を展開した後にできたディレクトリ draw7.5-2013-02-26-UsingU-11to20yr-pub の中に Reanal/ECMWF というディレクトリ内に置いてください.
以下の, tar による展開は, dcpam 計算結果の *.nc ファイルのあるディレクトリに おいて行うことが想定されています.
% tar zxvf util_s2p-2012-09-12.tgz % cd util_s2p-2012-09-12 適当に Makefile を編集 (Fortran コンパイラと netcdf が必要) % make % mkdir out (結果のファイルの格納ディレクトリ) % cp s2p_Earth.nml s2p.nml (設定ファイルをコピー) s2p.nml を編集. ------------------------------------------------------------ &time ts = 2921, te = 4381 / !&time ts = 14601, te = 29201 / ------------------------------------------------------------ を ------------------------------------------------------------ !&time ts = 2921, te = 4381 / &time ts = 14601, te = 29201 / ------------------------------------------------------------ とする. % ./s2p_ncf % cd util_calcrh-2013-03-02 (比湿計算作業ディレクトリに移動) 適当に Makefile を編集 (Fortran コンパイラと netcdf が必要) % make % ./calcrh (比湿を計算. 計算結果は ../out ディレクトリに入る.)
util_s2p-2012-09-12/out の下に *.nc ファイルができます.
なお, util_s2p-2012-09-12 内のプログラムの使い方の詳細は, util_s2p-2012-09-12/README.rd を参照してください.
注意: 以下の, tar による展開は, dcpam 計算結果の *.nc ファイルのあるディレクトリにおいて行うことが想定されています.
% tar zxvf util_timemean-2012-11-02.tgz % cd util_timemean-2012-11-02 適当に Makefile を編集 (Fortran コンパイラと netcdf が必要) % make % mkdir out (結果のファイルの格納ディレクトリ) 下の内容を含む timemean.cntl を作成. ------------------------------------------------------------ &mean NDataPerDay = 4, NDataDays = 31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31 &end &file ncfn_in = '../OLRA.nc', varname = 'OLRA', ncfn_out = './out/OLRA.nc' &end &file ncfn_in = '../SurfTemp.nc', varname = 'SurfTemp', ncfn_out = './out/SurfTemp.nc' &end &file ncfn_in = '../util_s2p-2012-09-12/out/U.nc', varname = 'U', ncfn_out = './out/U.nc' &end &file ncfn_in = '../util_s2p-2012-09-12/out/Temp.nc', varname = 'Temp', ncfn_out = './out/Temp.nc' &end ------------------------------------------------------------ % ./timemean
util_timemean-2012-11-02/out の下に *.nc ファイルができます.
以下の, tar による展開は, dcpam 計算結果の *.nc ファイルのあるディレクトリに おいて行うことが想定されています.
また, 上に説明したように, 比較用のデータ置いてください.
% tar zxvf draw7.5-2013-02-26-UsingU-11to20yr-pub.tgz % cd draw7.5-2013-02-26-UsingU-11to20yr-pub % ruby drawgm.rb (annual mean, global mean の値を計算) % ruby draw_ctl.rb (作図) % cd texdir (TeX 作業ディレクトリへ移動) % make (TeX コンパイル)
main.pdf ができます.