LBLRTM メモ
本家ウェブページ
メモ
- gfortran ではうまく動かない
- pgi fortran では動く
準備
<TOPDIR> = ~/work/AER
使用ファイル (バージョン)
- aer_v_3.5.tar.gz
- aer_lblrtm_v12.7.tar.gz
- aerlnfl_v3.0.tar.gz
- aerradsum_v2.7.tar.gz
ファイル展開
$ tar xvf aer_v_3.5.tar.gz
aer_v_3.5 ができた
$ tar xvf aer_lblrtm_v12.7.tar.gz
lblrtm ができた
$ tar xvf aerlnfl_v3.0.tar.gz
lnfl_v3.0.1 ができた
$ tar xvf aerradsum_v2.7.tar.gz
radsum ができた
コンパイル
lnfl
$ cd <TOPDIR> $ cd lnfl_v3.0.1/build $ lv README.build_instructions $ make -f make_lnfl linuxPGIsgl
lblrtm
$ cd <TOPDIR> $ cd lblrtm/build $ lv README.build_instructions $ make -f make_lblrtm linuxPGIsgl
radsum
$ cd <TOPDIR> $ cd radsum/build $ lv README.build_instructions $ make -f make_radsum linuxPGIsgl
実行
lnfl
$ cd <TOPDIR> $ cd lnfl_v3.0.1 $ mkdir -p run/all_00000-50025 $ cd run/all_00000-50025 $ cp ../../run_examples/run_example_infrared/run_lnfl . $ cp ../../run_examples/run_example_infrared/TAPE5_ex . $ ln -s ../../*_param . $ ln -s ../../../aer_v_3.5/line_file/aer_v_3.5 .
run_lnfl を変更
time lnfl $line_file -> time ./lnfl $line_file
TAPE5_ex を変更
300. 3500.
-> 0. 50025.
実行
$ ./run_lnfl
lblrtm & radsum
$ cd <TOPDIR> $ mkdir -p run/trp $ cd run/trp $ ln -s ../../lblrtm/lblrtm_v12.7_linux_pgi_sgl lblrtm $ ln -s ../../radsum/radsum_v2.7_linux_pgi_sgl radsum $ ln -s ../../lnfl_v3.0.1/run/all_00000-50025/TAPE3_aer_v_3.5 TAPE3 $ cp ../../radsum/run_examples/script.run_radsum . $ cp ../../radsum/run_examples/tape5.lbl.sample . $ cp ../../radsum/run_examples/tape6.lbl.sample . $ cp ../../radsum/run_examples/IN_RADSUM . $ cp tape5.lbl.sample tape5.lbl.trp
script.run_radsum を変更 (変更後ファイル)
lblrtm -> ./lblrtm \cp tape5.lbl.sample TAPE5 -> \cp tape5.lbl.trp TAPE5 \mv TAPE6 tape6.lbl.sample -> #\mv TAPE6 tape6.lbl.sample radsum -> ./radsum \rm ODdef* TAPE3? TAPE6? TAPE?? TAPE7 TAPE9 TAPE5 -> #\rm ODdef* TAPE3? TAPE6? TAPE?? TAPE7 TAPE9 TAPE5
- TRP プロファイルを設定したファイル
- 計算結果消去用 Makefile
- OUTPUT_RADSUM netcdf 変換プログラム
- TAPE5 用プロファイルデータ表示プログラム
しかし, 上向きフラックスが特に何だか変なんだよな. 地面付近でフラックスが減りすぎる.