研究のページ/領域気象モデル/WRFの出力結果を描画のバックアップ(No.7)

バックアップ一覧
差分を表示
現在との差分を表示
ソースを表示
研究のページ/領域気象モデル/WRFの出力結果を描画へ行く。
- 1 (2013-04-16 (火) 19:03:30)
- 2 (2013-04-16 (火) 19:53:21)
- 3 (2013-04-17 (水) 15:11:08)
- 4 (2013-04-19 (金) 11:38:32)
- 5 (2013-05-14 (火) 19:07:58)
- 6 (2013-05-15 (水) 19:15:53)
- 7 (2013-05-19 (日) 23:07:32)
- 8 (2013-06-01 (土) 20:00:12)

更新： 2013-05-15 (水) 19:15:16

WRF の出力結果を描画することについて †

出力データの概要
出力データは、"netCDF" (Network Common Deta Form) という形式で出力される。
そこで、netCDF 形式のデータが扱えるソフトであれば処理や描画ができる。

ARWの User's Guide (Chapter 9) には、いくつかの処理ソフトが紹介されている。
- NCL 　(User's Guide 9-2～)
  NCAR Graphics Command Laguage。NCAR で開発されているので、WRF とも相性がいいと予想される。独自のコマンド体系で描画する。
  NCAR Graphics (NCARG)は以前使ったことがあるけど、全然違う印象を受ける。
- RIP4　(User's Guide 9-20～)
  これも NCAR Graphics に基づいた対話型のソフトウエア。WRF が MM5 として開発されていたころから描画ように使われていたようである。
- ARWpost　(User's Guide 9-29～)
  気象でよく用いられる GrADS でデータを扱うことが出きるようにするためのコンバータ。
- UUP　(User's Guide 9-36～)
  NCEP Unified Post Processor。以前は WPP (WRF Post Processor)。GEMPAK, GrADS などで描画する必要がある？
- VAPOR　(User's Guide 9-59～)
  Visualization and Analysis Platform for Ocean, Atmosphere, and Solar Researchers。3次元表示が得意なようである。

ARWpost + GrADS †

導入目的(余談)

netCDF データは自己記述型のデータである。そこで、netCDF データを入手し、適当な描画ツールを用いれば、どの
ような変数が、どのような配列として格納されているか、知ることができ、そうした情報を元に、描画すればよい。

…と思うと、そんなことはなく、WRF ARW の出力ファイルは、netCDF ファイルでありながら、どの変数が何を意味
しているのか、正確なところはソースを読み込む必要があることが多い。具体的に、個人的には、変数 T で困った。
こうしたことを、WRF の達人に聞いたところ、ARWpost を使って、GrADS で描画するのがいいと勧められた。GrADS は
プログラミング環境として貧弱なので、できるだけ使いたくなかった。しかし、このような状況では致し方ない。
使ってみることにした。

環境整備
- 導入パッケージ
  - grads

インストール
入力した部分を -> で示す。

-> $ ./configure
   ** WARNING: No path to NETCDF and environment variable NETCDF not set.
   ** would you like me to try to fix? [y]
-> y
   Enter full path to NetCDF include directory on your system
-> /usr/include
   Enter full path to NetCDF library directory on your system
-> /usr/lib

   …
   ------------------------------------------------------------------------
   …
      3.    PC Linux i486 i586 i686 x86_64, gfortran compiler
   …
   Enter selection [1-3] : 3
   …

-> $ ./compile
   …

実行
適切な namelist.ARWpost を作成してから実行する。
注意
どうやら、gfortran ではうまく動かないらしい。少なくとも私の環境ではうまく動かなかった。
WRFのフォーラムにも記述があり、どうやら解決していないらしい。
残念！

↑

R + netCDF †

環境整備
- 導入パッケージ
  - r-base
    統計パッケージ R の基本的なパッケージがまとめて入る。
  - ncdf
    netCDF のためのパッケージ。一般的に、導入方法は二つある。root で次の作業のどちらかを行う。
    CRAN サイトからパッケージのソースをダウンロードして、root で
```
R CMD INSTALL ncdf_1.6.6.tar.gz
```
    とするか、root で R を立ち上げて、
```
> install.packages("ncdf")
```
    とする。
    ※ 実行後は、/usr/local/lib/R/site-library/ncdf/ をチェック！パーミッションが正しくないと使えない。
簡単な描画
典型的には、R の起動後、次のようにする。
```
> library(ncdf)
> nc <- open.ncdf("wrfout_d01_2007-09-01_00:00:00")
> print(nc)
> filled.contour( get.var.ncdf(nc, "P") [,,1,3])
```
※ print(nc) で、格納されている変数名がわかる。

↑

ncview †

　　とりあえず使うには重宝する。

debian で配布されるパッケージの説明

説明: NetCDF 形式のファイル用の視覚的な X11 向けブラウザ
ncview を使うと、素早く簡単にボタンを押すだけで NetCDF ファイルを見ることが
できます。データの単純な動画の表示、さまざまな次元での表示、実際のデータの
値の表示、カラーマップの変更、データの反転や他の単純な視覚的操作ができま
す。

関連ページ
- ncview ホームページ

研究のページ/領域気象モデル/WRFの出力結果を描画 のバックアップ(No.7)

WRF の出力結果を描画することについて †

ARWpost + GrADS †

R + netCDF †

ncview †

研究のページ/領域気象モデル/WRFの出力結果を描画のバックアップ(No.7)