はじめてのNCL NCLの実例 NCLのインストール NCLのチュートリアル プログラムソースの解説 いざ実践 計算機セミナー

Presentation on theme: "はじめてのNCL NCLの実例 NCLのインストール NCLのチュートリアル プログラムソースの解説 いざ実践 計算機セミナー"— Presentation transcript:

1 はじめてのNCL NCLの実例 NCLのインストール NCLのチュートリアル プログラムソースの解説 いざ実践 計算機セミナー
はじめてのNCL 近本 喜光　(東大・気候システム) 本日の発表内容 NCLの実例 NCLのインストール NCLのチュートリアル プログラムソースの解説 いざ実践

2 1. NCLの実例

3 NCLでできること 絵を描く データ解析をする 特徴 オブジェクト指向 多種類の統計解析が簡単にできる Freeで利用でき、ユーザー数が多い
欠点 マニュアルがすべて英語 日本では、まだまだユーザーが少ない

4 NCLの実例: その１ BGM 観測 同化 CTL予測 LAF予測 BGM予測 p1 m1 LAF

5 NCLの実例: その２ ASSIM GOOD BAD
上段: SST, U_sfc, V_sfc における偏差 下段: 赤道域のTO,UO,WOにおける偏差

6 NCLを使った論文の紹介 Chikamoto et al. 2007: Evidence of Growing Bred Vector Associated with the Tropical Intraseasonal Oscillation Geophys. Res. Lett., Vol. 34, L04806, doi: /2006GL028450 Chikamoto and Tanimoto 2006: Air-sea Humidity Effects on the Generation of Tropical Atlantic SST Anomalies during the ENSO Events Geophys. Res. Lett., Vol. 33, L19702, doi: /2006GL02723 Chikamoto and Tanimoto 2005: Role of specific humidity anomalies in Caribbean SST response to ENSO J. Meteor. Soc. Japan, Vol. 83, No. 6, pp

7 ２. NCLのインストール

8 インストールの主な手順 ESGにアカウント登録 NCLファイルをダウンロード NCLファイルを展開 パスを通す インストールできたか確認

9 NCLのホームページ

10 NCLのインストール手順

11 ESGの top page

12 NCLのコードを選択

13 自分のシステムに合ったNCLを選択 Vine Linux 4.2 の場合 (32bit, i686, gcc3.x)

14 NCLのインストール $ rpm –qa | grep gcc gcc-3.3.6-0vl7 $ mkdir ncl $ cd ncl
$ tar zxvf ncl_ncarg Linux_i686_gcc3.tar.gz $ pwd /home/chika44/soft/ncl $ vi ~/.bashrc export NCARG_ROOT=/home/chika44/soft/ncl export PATH=${PATH}:${NCARG_ROOT}/bin $ source ~/.bashrc gccのバージョン

15 インストールできたか確認 こうなれば成功 $ ncl
Copyright (C) All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 5.0.0 The use of this software is governed by a License Agreement. See for more details. ncl 0> こうなれば成功

16 3. NCLのチュートリアル

17 NCLを使い始める２つの方法 Bottom up 式 ： 文法からしっかり知りたい人向け Top down 式
　　　：　文法からしっかり知りたい人向け Top down 式 　　　：　習うより慣れろタイプの人向け

18 Top down 式チュートリアル

19 最初に抑えるべき項目

20 Anomalies のリンクをたどってみる
プログラムの内容 プログラムソース 実行結果 使用した関数の説明

21 Can we run?

22 4. プログラムソースの解説

23 木を見て森を見る おまじない begin データ解析 描画 end

24 おまじないについて 基本： なにも考えずにこの４つを load しておけばよい. 応用： 自分で作成したサブルーチンを集めて、１つの
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl" 基本： なにも考えずにこの４つを load しておけばよい. 応用： 自分で作成したサブルーチンを集めて、１つの nclファイルとし、それを load することもできる

25 データ解析について ファイル fn 内にある変数 TS を読み込み、変数 ts に代入 コメント r: 読み込み c:作成 コメント
東西平均を削除 printVarSummary(ts) の１行を加えて実行するとtsの説明がでる anom(lat|:,lon|:,time|:)=dim_rmvmean(ts(lat|:,lon|:,time|:)) とすれば、時間平均を引くことができる

26 描画について wks = gsn_open_wks(“ps”,”dev”) plot = new(2,graphic) res = True
= False = False plot(0) = gsn_csm_contour_map_ce(wks,ts(0,:,:),res) plot(1) = gsn_csm_contour_map_ce(wks,anom(0,:,:),res) gsn_panel(wks,plot,(/2,1/),False) 出力ファイル名 ファイル形式

27 5. いざ実践