北海道大学 大学院理学研究院・理学院宇宙理学専攻

Slides:



Advertisements
Similar presentations
ゲームプログラミング講習 第2章 関数の使い方
Advertisements

サービス管理責任者等研修テキスト 分野別講義    「アセスメントと        支援提供の基本姿勢」 <児童発達支援管理責任者> 平成27年10月1日.
ヒトの思考プロセスの解明を目的とするワーキングメモリの研究
第27講 オームの法則 電気抵抗の役割について知る オームの法則を使えるようにする 抵抗の温度変化を理解する 教科書P.223~226
コラッツ予想の変形について 東邦大学 理学部 情報科 白柳研究室 山中 陽子.
コンパイラ 第3回 字句解析 ― 決定性有限オートマトンの導出 ―
第5章 家計に関する統計 ー 経済統計 ー.
公共財 公共経済論 II no.3 麻生良文.
VTX alignment D2 浅野秀光 2011年12月15日  放射線研ミーティング.
冷却フランシウム原子を用いた 電子の永久電気双極子能率探索のための ルビジウム磁力計の研究
生命情報学 (8) スケールフリーネットワーク
前半戦 「史上最強」風 札上げクイズ.

認知症を理解し 環境の重要性について考える
フッ化ナトリウムによる洗口 2010・9・13 宮崎市郡東諸県郡薬剤師会 学校薬剤師  日高 華代子.
食品の安全性に関わる社会システム:総括 健康弱者 ハイリスク集団 HACCP (食肉処理場・食品工場) 農場でのQAP 一般的衛生管理
規制改革とは? ○規制改革の目的は、経済の活性化と雇用の創出によって、   活力ある経済社会の実現を図ることにあります。
地域保健対策検討会 に関する私見(保健所のあり方)
公共政策大学院 鈴木一人 第8回 専門化する政治 公共政策大学院 鈴木一人
医薬品ネット販売規制について 2012年5月31日 ケンコーコム株式会社.
平成26年8月27日(水) 大阪府 健康医療部 薬務課 医療機器グループ
平成26年度 呼吸器学会からの提案結果 (オレンジ色の部分が承認された提案) 新規提案 既収載の変更 免疫組織化学染色、免疫細胞化学染色
エナジードリンクの危険性 2015年6月23日 経営学部市場戦略学科MR3195稲沢珠依.
自動吸引は 在宅を変えるか 大分協和病院 院長         山本 真.
毎月レポート ビジネスの情報 (2016年7月号).
医療の歴史と将来 医療と医薬品産業 個人的経験 3. 「これからの医療を考える」 (1)医薬品の研究開発 -タクロリムスの歴史-
社会福祉調査論 第4講 2.社会調査の概要 11月2日.
2015年12月28日-2016年3月28日 掲載分.
2010度 民事訴訟法講義 補論 関西大学法学部教授 栗田 隆.
腫瘍学概論 埼玉医科大学国際医療センター 包括的がんセンター 緩和医療科/緩和ケアチーム 奈良林 至
“企業リスクへの考え方に変化を求められています。 トータルなリスクマネジメント・サービスをプロデュースします。“
情報漏えい 経済情報学科 E  西村 諭 E  釣 洋平.
金融班(ミクロ).
第11回 2009年12月16日 今日の資料=A4・4枚+解答用紙 期末試験:2月3日(水)N2教室
【ABL用語集】(あいうえお順) No 用語 解説 12 公正市場価格 13 債権 14 指名債権 15 事業収益資産 16 集合動産 17
基礎理論(3) 情報の非対称性と逆選択 公共政策論II No.3 麻生良文.
浜中 健児 昭和42年3月27日生まれ 東京都在住 株式会社ピー・アール・エフ 代表取締役 (学歴) 高 校:千葉県立東葛飾高校 卒業
COPYRIGHT(C) 2011 KYUSHU UNIVERSITY. ALL RIGHTS RESERVED
Blosxom による CMS 構築と SEO テクニック
記入例 JAWS DAYS 2015 – JOB BOARD 会社名 採用職種 営業職/技術職/その他( ) 仕事内容 待遇 募集数
ネットビジネスの 企業と特性 MR1127 まさ.
Future Technology活用による業務改革
ネットビジネス論(杉浦) 第8回 ネットビジネスと情報技術.
g741001 長谷川 嵩 g740796 迫村 光秋 g741000 西田 健太郎 g741147 小井出 真聡
自然独占 公共経済論 II no.5 麻生良文.
Autonomic Resource Provisioning for Cloud-Based Software
Webショップにおける webデザイン 12/6 08A1022 甲斐 広大.
物理的な位置情報を活用した仮想クラウドの構築
ハイブリッドクラウドを実現させるポイントと SCSKのOSSへの取組み
寺尾 敦 青山学院大学社会情報学部 第12回 情報デザイン(4) 情報の構造化と表現 寺尾 敦 青山学院大学社会情報学部
【1−1.開発計画 – 設計・開発計画】 システム開発計画にはシステム開発を効率的、効果的に実行する根拠(人員と経験、開発手順、開発・導入するシステム・アプリケーション・サービス等)を記述すること。 システム開発の開始から終了までの全体スケジュールを記載すること。 アプリケーション機能配置、ソフトウェア、インフラ構成、ネットワーク構成について概要を示すこと。
6 日本のコーポレート・ガバナンス 2008年度「企業論」 川端 望.
急成長する中国ソフトウェア産業 中国ソフトウェアと情報サービス産業の規模 総売上高は5年間で約5.3倍の成長
米国ユタ州LDS病院胸部心臓外科フェローの経験
公益社団法人日本青年会議所 関東地区埼玉ブロック協議会 JCの情熱(おもい)育成委員会 2011年度第1回全体委員会
次世代大学教育研究会のこれまでの活動 2005年度次世代大学教育研究大会 明治大学駿河台校舎リバティタワー9階1096教室
子どもの本の情報 大阪府内の協力書店の情報 こちらをクリック 大阪府内の公立図書館・図書室の情報
第2回産業調査 小島浩道.
〈起点〉を示す格助詞「を」と「から」の選択について
広東省民弁本科高校日語専業骨幹教師研修会 ①日本語の格助詞の使い分け ②動詞の自他受身の選択について   -日本語教育と中日カルチャーショックの観点から- 名古屋大学 杉村 泰.
■5Ahバッテリー使用報告 事例紹介/東【その1】 ■iphon4S(晴れの昼間/AM8-PM3) ◆約1時間で68%⇒100%
『ワタシが!!』『地域の仲間で!!』 市民が始める自然エネルギー!!
ポイントカードの未来形を形にした「MUJI Passport」
SAP NetWeaver を支える Microsoft テクノロジーの全貌 (Appendix)
ガイダンス(内業) 測量学実習 第1回.
Python超入門 久保 幹雄 東京海洋大学.
熱力学の基礎 丸山 茂夫 東京大学大学院 工学系研究科 機械工学専攻
京都民医連中央病院 CHDF学習推進委員会
資料2-④ ④下水道.
Accessによる SQLの操作 ~実際にテーブルを操作してみよう!~.
Presentation transcript:

北海道大学 大学院理学研究院・理学院宇宙理学専攻 地球惑星情報学I 大気大循環モデル 石渡正樹 北海道大学 大学院理学研究院・理学院宇宙理学専攻 情報実験第11回 2018年7月13日

目次 地球惑星科学分野における数値モデル 大気大循環モデルとは DCPAM

地球惑星科学分野 における数値モデル

数値モデル 数値計算の重要性はどんどん増加 観測できない量を知る 定量的な情報を得る 数値モデルは知見の集積場(これについては次回に) 科学の第4の手法と呼ばれることもある 観測できない量を知る 観測困難な場所の状態 未来の状態 定量的な情報を得る どのようなプロセスが最も重要か?などを考えるヒントを得る 数値モデルは知見の集積場(これについては次回に)

数値モデルの例: 地球化学分野 MELTS 鉱物組成を与えると 熱力学方程式を 用いた第一原理計算 により相図を作成 http://melts.ofm-research.org 鉱物組成を与えると 熱力学方程式を 用いた第一原理計算 により相図を作成 開発者は Mark Ghiorso氏 (OFM Research, ワシントン大学)らしい 灰長石-透輝石系の相図 1570℃ の場合

数値モデルの応用例: 地球内部物理学分野 マントル対流のシミュレーション 時間 色: 温度 等値線: マグマ量 Ogawa (2014) J. Geophys. Res., 119, 867-883

数値モデルの応用例: 天文分野 渦巻き銀河の形成のシミュレーション 青色:ガス,黄色:恒星 (1) ガス集積 (2) 小銀河の による恒星 の形成 (2) 小銀河の 合体による 成長 (3) 渦巻き 構造の形成 (4) 銀河の 出来上がり http://4d2u.nao.ac.jp/t/var/download/spiral2.html Saito and Wada (2004) ApJ, 615, L93-L96

数値モデルに共通した特徴 数値モデルの「実体」は巨大なプログラム 通常は多数のファイルの集合体 フリーソフトウエアであるものも、そうでないものも存在する 使用する場合にはライセンスの確認を! (ライセンスについてはINEX第9回参照)

大気大循環モデルとは

大気大循環モデルとは 大気の循環・熱力学量・物質分布の時間発展を計算するソフトウエア 英語では Atmospheric General Circulation Model (AGCM) ちなみに海洋 大循環モデルは OGCM 高橋他(2012)

AGCM使用例(1):天気予報 予報天気図 実況天気図 気象庁 天気図 (https://www.jma.go.jp/jp/g3/)

AGCM使用例(2):温暖化予測 年平均表面気温変化 (モデル平均、1986-2005 平均からの差) 2016-203 2081-2100 2016-203 2081-2100 RCP85 温度の時間変化の図(Fig. SPM7, Fig. TS14:103枚目)と平面分布(Fig. SPM8, TS15)を入れたい TS はTechnical summary. RCP60 シナリオ の違い RCP45 RCP26 IPCC(2013)

AGCM使用例(3):惑星気候研究 火星、金星、木星、系外惑星など 例:同期回転惑星(系外惑星の1種) 同期回転惑星:固定された 夜半球と昼半球を持つ 計算結果:365日平均場の視点を                変えた動画 A B Star Ω*=0.5, S=1600W/m2 色: 表面温度 矢印: 水平風 等値線: 降水 黄色点: 恒星直下点

基礎方程式 運動方程式 静水圧の式 質量保存則 エネルギー保存則 水蒸気の式 状態方程式 流体力学の方程式 𝑑𝑢 𝑑𝑡 − 𝑓+ 𝑢 tan 𝜑 𝑎 𝑣= − 1 𝜌 𝜕𝑝 𝑎 cos 𝜑 𝜕𝜆 + 𝐹 λ 運動方程式 𝑑𝑣 𝑑𝑡 + 𝑓+ 𝑢 tan 𝜑 𝑎 𝑢= − 1 𝜌 𝜕𝑝 𝑎𝜕𝜑 + 𝐹 φ 0=− 1 𝜌 𝜕𝑝 𝜕𝑧 −𝑔 静水圧の式 𝑑𝜌 𝑑𝑡 +𝜌 1 𝑎 cos 𝜑 𝜕𝑢 𝜕𝜆 + 𝜕 𝑣 cos 𝜑 𝜕𝜑 + 𝜕𝑤 𝜕𝑧 =0 質量保存則 𝐶 𝑣 𝑑𝑇 𝑑𝑡 − 𝑝 𝜌 2 𝑑𝜌 𝑑𝑡 =𝑄 エネルギー保存則 𝑑𝑞 𝑑𝑡 =𝑆 水蒸気の式 𝑝=𝜌𝑅𝑇 状態方程式

空間分布の表現(1):格子点法 離散的な点(格子点)上の値のみを考える 𝑑𝑢 𝑑𝑧 𝑢 𝑘+1 − 𝑢 𝑘 ∆𝑧 微分は格子点上の値の差であらわす 多くのGCMでは鉛直方向には格子点法を採用 格子点間隔より小さいスケールの現象は表現できない z 𝑢 𝑘𝑚𝑎𝑥 モデル上端 ・・・・ 𝑑𝑢 𝑑𝑧 𝑢 𝑘+1 − 𝑢 𝑘 ∆𝑧 𝑢 𝑘+1 ∆𝑧 𝑢 𝑘 格子点数が多いほど モデル解像度は高くなる ・・・・ 𝑢 1 地表面

空間分布の表現(2):スペクトル法 物理量を直交関数系で展開 例:1次元問題、三角関数で展開(フーリエ級数) 使う関数系によっては、微分の計算が楽になる、 空間積分が高精度でできる 用いる関数(項)の数が多いほど解像度が高くなる 例:1次元問題、三角関数で展開(フーリエ級数) 𝑢 𝑥 = 𝑎 0 + 𝑘=1 𝐾 𝑎 𝑘 sin 𝑘𝑥 + 𝑘=1 𝐾 𝑏 𝑘 cos 𝑘𝑥 波数3まで こんな関数も フーリエ展開 できる 波数30まで 波数50まで

空間分布の表現(2):球面のスペクトル法 物理量を球面調和関数で展開 水平積分を高精度で行うことができる  気象庁モデル、DCPAMなどで利用されている 物理量A(λ φ)の展開式 𝐴 𝑖𝑗 = 𝑚=−𝑀 𝑀 𝑛=|𝑚| 𝑁 𝐴 𝑛 𝑚 𝑌 𝑛 𝑚 ( λ 𝑖 , φ 𝑗 ) λ:経度, φ:緯度 m:東西波数, n:全波数 球面調和関数: 𝑌 𝑛 𝑚 λ 𝑖 , φ 𝑗 = 𝑃 𝑛 𝑚 ( sin φ ) 𝑒 𝑖𝑚λ ルジャンドル陪関数 𝑃 𝑛 𝑚 μ = 𝑛+1 𝑛− 𝑚 ! 𝑛+ 𝑚 ! 1−𝜇 𝒎 𝟐 2 𝑛 𝑛! 𝑑 𝑛+ 𝑚 𝑑𝜇 𝑛+ 𝑚 ( 𝜇 2 −1) 𝑛

ここまでのまとめ 大気大循環モデルは風速・気温などの空間分布・時間変化を求める やっていることは微分方程式の数値積分 数値積分をおこなうには数値計算独特のやり方をしないといけない 𝑑𝑢 𝑑𝑡 − 𝑓+ 𝑢 tan 𝜑 𝑎 𝑣= − 1 𝜌 𝜕𝑝 𝑎 cos 𝜑 𝜕𝜆 + 𝐹 λ 𝑑𝑣 𝑑𝑡 + 𝑓+ 𝑢 tan 𝜑 𝑎 𝑢= − 1 𝜌 𝜕𝑝 𝑎𝜕𝜑 + 𝐹 φ 数値積分 0=− 1 𝜌 𝜕𝑝 𝜕𝑧 −𝑔 差分化 格子点法 スペクトル法 ・・・・・・

DCPAM

DCPAMとは 地球流体電脳倶楽部で開発しているAGCM https://www.gfd-dennou.org/library/dcpam 他の惑星への拡張を意識 可読性を重視 プログラムが読みやすく なるように数式の書き方を 工夫 フリーソフトウエア 日本のAGCMで フリーソフトウエアとして 公開されているものは 非常に少ない 可読性・可変性はわからない(小高)

地球流体電脳倶楽部 知識の情報化, 知見プラットフォームの構築を目指した有志集団 https://www.gfd-dennou.org 地球惑星に関する知見の ネットワーク上への蓄積, そのための道具作り ネットワーク上の「教科書」 地球流体室内実験集 知見の集積場としての地球 流体計算ソフトウェア群 数値データの可視化ツール開発 それらのためのサーバ運営・管理(全国3か所) dennou-k はDebian のインストールで使用

DCPAMで計算するためには 実行ファイルを作る(ビルド作業) プログラムファイルをコンパイル(機械語に翻訳)する オブジェクトファイル:コンパイルで作られる機械語ファイル 情報実験機で使うコンパイラはGNU Fortran (GFortran) 使用するライブラリと結合して実行ファイルを作る ライブラリ:複数のソフトウエアで共通して使用されるプログラムを「外から使える形で」まとめたもの 実行ファイル:計算機が命令を実行できるファイル 結合をするためにはリンカと呼ばれるコマンドを使う コンパイル 結合 プログラムファイル:*.f90 オブジェクトファイル*.o 実行ファイル makeコマンド で一発! 外部のライブラリlib****.a

DCPAMがおこなう処理の流れ DCPAMは数値積分の実行前に以下を読み込む 設定ファイル 初期値データファイル 境界条件ファイル 積分の進行に応じて結果をデータファイルに出力 設定 ファイル 初期値 データ 境界条件 データ 数値積分 (反復計算) 終了 処理 初期設定 出力データ 温度 東西風 降水量 ・・・・ 温度

出力データ の例 詳しくは次回に! $ ncdump –h SurfTemp.nc netcdf SurfTemp { dimensions: lon = 128 ; (略) variables: float lon(lon) ; lon:long_name = “longitude” ; (略) float SurfTemp(time, lat, lon) ; SurfTemp:long_name = “surface temperature” ; (略) // global attributes: :title = “….." ; (略) data: lon = 0, ……….   (略) SurfTemp = 237.6823, …. } 出力データ の例 詳しくは次回に!

まとめ 大気大循環モデルは微分方程式の数値積分をおこない風速・気温などの空間分布・時間変化を求める 大気大循環モデルを実行するには、コンパイル・結合の作業が必要(多くの科学計算ソフトウエアも同様) 𝑑𝑢 𝑑𝑡 − 𝑓+ 𝑢 tan 𝜑 𝑎 𝑣= − 1 𝜌 𝜕𝑝 𝑎 cos 𝜑 𝜕𝜆 + 𝐹 λ 数値積分 ・・・・ コンパイル 結合 プログラムファイル オブジェクトファイル*.o 実行ファイル ライブラリlib****.a

実技編では DCPAMのインストール作業 地球設定実験 $ make

参考書, 参考文献 Oram, A., Talbott, S. 著, 矢吹道郎監訳, 1997: make 改訂版, オライリー・ジャパン IPCC, 2013: Climate change 2013 The physical science basis. http://www.ipcc.ch/ Saito, T., Wada, K., 2004: Coevolution of galactic cores and spiral galaxies. Astrophysical J., 615, L93-L96 Wirth, K., Advanced Modeling Programs: MELTS. https://serc.Carleton.edu/research_education/equilibria/melts.html Ogawa (2014) A positive feedback between magmatism and mantle upwelling in terrestrial planets: Implications for the Moon. J. Geophys. Res., 119, 867-883 Satoh, M., 2004: Atmospheric circulation dynamics and general circulation models. Springer Mecklenburg, R. 著、矢吹道郎監訳, 菊池彰訳, 2005: GNU Make 第3版, オライリー・ジャパン 坂井弘亮, 2010: リンカ・ローダ実践開発テクニック, CQ出版社 伊理正夫, 藤野和建, 1985:数値計算の常識, 共立出版 地球流体電脳倶楽部dcmodelプロジェクト, 2017: 大気大循環モデルDCPAM, http://www.gfd-dennou.org/library/dcpam ビジネスアイコン無料素材http://business-icon.com/highresolution/l_005.png