IT pass HikiWiki - [Exp2022]GPhys スクリプト課題案 Diff

  • Added parts are displayed like this.
  • Deleted parts are displayed like this.

{{toc}}


* ((<2009 年のレポート課題|[Exp2009]シェルスクリプト課題>)), ((<2010 年のレポート課題|[Exp2010]シェルスクリプト課題>)), ((<2011 年のレポート課題|[Exp2011]シェルスクリプト課題>)), ((<2012 年のレポート課題|[Exp2012]シェルスクリプト課題>)), ((<2013 年のレポート課題|[Exp2013]シェルスクリプト課題>)), ((<2014 年のレポート課題|[Exp2014]シェルスクリプト課題>)), ((<2015 年のレポート課題|[Exp2015]シェルスクリプト課題>)), ((<2016 年のレポート課題|[Exp2016]シェルスクリプト課題>)), ((<2017 年のレポート課題|[Exp2017]シェルスクリプト課題>)), ((<2018 年のレポート課題|[Exp2018]シェルスクリプト課題>)), ((<2019 年のレポート課題|[Exp2019]シェルスクリプト課題>)), ((<2020 年のレポート課題|[Exp2020]GPhysスクリプト課題>)), ((<2021 年のレポート課題|[Exp2021]GPhysスクリプト課題>))  も参考にしてみてください.




== 問題 1

((<NCEP/NCAR Reanalysis|URL:https://psl.noaa.gov/data/gridded/data.ncep.reanalysis.html>))のデータを用いて,
以下の # つの中からどれか一つの図を描画するスクリプトを作成しなさい.
そして読み取れることを説明しなさい.


* 赤道における東西平均東西風の高度分布の時間変化
  * 現象の説明
    * 成層圏の赤道域では, 東西風が約二年の周期で変化することが知られ
      ており, 成層圏準二年周期振動 (QBO: Quasi-biennial oscillation)
      と呼ばれています. これは, その特異な周期と構造から, 直観的には
      その原因を理解し難い興味深い現象です.
  * 補足・ヒント
    * 横軸を時間, 縦軸を圧力 (高度) として東西平均 (経度方向に平均)
      した東西風の分布を描画するとよいでしょう.

===ここから問題文を書く

* 東太平洋赤道域における 1997 年 11 月と 1988 年 12 月の月平均海面水温の平年値との偏差の分布
  * 現象の説明
    * 東太平洋の赤道域では, 海面水温が平年より高くなるエルニーニョ現象と, 平年より低くなるラニーニャ現象が数年おきに発生します. これらの現象は, 世界中の天候に影響を与えると考えられています.
  * 補足・ヒント
    * 1997 年 11 月は 1997 年春に発生し, 1998 年春に終息したエルニーニョ現象の最盛期に当たります.
    * 1988 年 12 月は 1988 年春に発生し, 1989 年春に終息したラニーニャ現象の最盛期に当たります.
    * エルニーニョ現象とラニーニャ現象の 2 枚の図を出力してください.
    * 偏差は, (指定した期間の海面水温の平均) - (指定した年の前年までの 30 年間の各月の平均値) で計算してください.
      * 例えば, (1997 年 11 月の海面水温の平均) - (1967 年から 1996 年の 11 月の海面水温の平均値).
    * ちなみに気象庁では、エルニーニョ監視海域(5ºS-5ºN, 150º-90ºW)の月平均海面水温の基準値(その年の前年までの 30 年間の各月の平均値)との差の 5 か月移動平均値(その月および前後 2 か月を含めた 5 か月の平均をとった値)が 6 か月以上連続して +0.5 度以上になった状態をエルニーニョ現象, −0.5 度以下となった場合をラニーニャ現象と定義しています.

* 赤道における東西平均東西風の高度分布の時間変化
  * 現象の説明
    * 成層圏の赤道域では, 東西風が約二年の周期で変化することが知られ
      ており, 成層圏準二年周期振動 (QBO: Quasi-biennial oscillation)
      と呼ばれています. これは, その特異な周期と構造から, 直観的には
      その原因を理解し難い興味深い現象です.
  * 補足・ヒント
    * 横軸を時間, 縦軸を圧力 (高度) として東西平均 (経度方向に平均)
      した東西風の分布を描画するとよいでしょう.

* ある一年におけるアジアモンスーンの描画
  * 現象の説明
    * モンスーンとは, 季節によってその向きが交替する季節風のことです. この現象が特に顕著に現れる地域として, ユーラシア大陸の南部から東部にかけての地域が挙げられます. 季節風の向きの変化に伴い, これらの地域には雨季と乾季が訪れます.
  * 補足・ヒント
    * アジア地域が収まるような領域について, 好きな年を選んで, 一ヶ月ごとに時間平均した東西風 (level: 10m), 南北風 (level: 10m), 降水量(level: Surface, 物理量の名前は Precipitation Rate)を重ねて描画し, 12 枚の図をまとめて一つの動画を作成すると良いでしょう.

* 太陽放射と気温の分布
  * 現象の説明
    * 地軸の傾きによって受け取る太陽放射の量は日平均太陽放射フラックスは年変化する. この太陽放射の実効値は日太陽放射フラックスは気温を決定する要素の一つである.
  * 補足・ヒント
    * ある地点での年平均気温から日平均気温の差を横軸に, 年間の太陽放射の実効値とその年間平均太陽放射フラックスとその日の平均太陽放射の実効値の日平均太陽放射フラックスの差を縦軸に分布図を一年分プロットすればいいでしょう.
  * また別の地点の同じ年のデータを使い, 比較考察もあればなお良いです.

  * ある地点である年の値が異常な値かもしれないと疑り深い方は, 平年気温(直近30年間の平均値)と複数年のデータとの標準偏差をとるとよいかもしれません.
* 火星大気の東西平均した子午面温度分布の季節変化
  * 注意 : この問題のみ例外的に NCEP/NCAR Reanalysis データを使わず, 惑星大気大循環モデル DCPAM を用いて行われた, 火星の実験の結果 ((<(こちら)|URL:http://itpass.scitec.kobe-u.ac.jp/~itpass/exp/fy2022/220811/practice_kadai/dcpam_mars/>)) を用います.
      * それぞれのファイルには一つの物理量の値が保存されています.
        * Temp.nc     : 大気温度
        * U.nc        : 東西風
        * V.nc        : 南北風
        * Ps.nc       : 地面気圧
        * SurfTemp.nc : 地面温度
      * ファイルに保存されたデータの時間間隔は 1/4 火星日であり, データの期間は 1 火星年分 (669 火星日分) です. 時間軸の最初は火星の北半球の春分点に対応しています.
      * データの詳細はメタデータを確認しましょう.
  * 現象の説明
    * 火星の大気の構造は様々な要因によって地球の大気とは異なっており,
      温度の季節変化も地球とは異なる特徴を示します.
  * 補足・ヒント
    * 東西平均 (経度方向に平均) した温度を 1/12 火星年ごとに時間平均して 12 枚の図を描画し, それらをまとめることで季節変化の動画を作成すると良いでしょう.
      * 1 火星年が 669 火星日ですから, 1/12 火星年は約 56 火星日に対応します.



==== レポートに書く内容

レポート内には以下を記述してください.
ただしリンクを貼る際は((*相対パスで指定するようにしてください*)).

* どの図を描画するスクリプトを作成したのか

* スクリプトファイル quiz1.rb.txt へのリンク

  上記で作成した quiz1.rb を quiz1.rb.txt という名前に変更してください. ブラウザ上から読めるようにするためです.
* スクリプトの解説

  スクリプト内の各行が何をしているのか, ((*スクリプト内にコメント文として書く*)).
* 描画に使用するデータファイルへのリンク

  リンク先はダウンロード元の URL を記してください.

* スクリプトの使い方の解説

  これをもとに TA が実際にテストします. 
  実行する際のスクリプト, データの置き場を明示してください.

#* 実行例 result1.txt, result1.png へのリンク
* 作成した動画ファイル (result1.gif) へのリンク


* 作成した図から読み取れること

  自分で作成した図から読み取れることを説明してください.
  必要に応じて参考文献, 参考ウェブページを引用しても構いません.
  その際には下の「参考資料」の欄に記入してください.

* 参考資料

  URL をベタ書きにするのではなく, そのページの作者やタイトル, および何を参考にしたのかも明記する.
* 共同作業した人の名前
* 工夫したこと, 感想など

#=== 実行結果の例
#
# (作成中)
#
# $ ruby quiz1.rb


== 問題 2

=== 興味のある大気現象の描画スクリプトの作成

自分が興味のある大気構造について, 図を描く GPhys を用いたスクリプトを作成しなさい. 動画でも構わない.
また, 作成した図について説明しなさい (何の図なのか). そして読み取れることを説明しなさい.

選ぶ題材は, 問題 1 で選んだもの以外のいずれかでも構わない.
問題 1 に挙げたもの以外の図を描くときには,
((<ここ|URL:https://psl.noaa.gov/data/gridded/data.ncep.reanalysis.html>))に ある NCEP NCAR Reanalysis データを使うとよいでしょう.

再解析データは実習中に説明した方法によってダウンロードしてください.

例えば, 2000 年における圧力面上の大気温度の日平均値のデータは下のようにたどることで見つけることができます.
* ((<NCEP/NCAR Reanalysis 1: Summary|URL:https://psl.noaa.gov/data/gridded/data.ncep.reanalysis.html>))
  の "We have separated the data documentation into seven sections" にある, "Pressure level" をクリック.
  * "Pressure level" : 圧力面上の値
* "Download/Plot Data" の "Air Temperature", "Daily", "Pressure" の行の "see list" をクリック.
  * "Air Temperature" : 大気温度
  * "Daily" : 日平均値
* "Level" の項目が "Multiple levels" の行の "See lists" をクリック.
  * "Multiple levels" : 複数面
* "air.2000.nc" をクリック.
  * ファイルの URL をコピーするには, リンク (air.2000.nc) を右クリックして「リンクのアドレスをコピー」のような項目をクリックする (名称はブラウザに依存する) と良いでしょう

物理量の名称は下のようなものがあるでしょう.
# RT
英語, 日本語, 備考
Air Temperature, 大気温度,
Geopotential Height, ジオポテンシャル高度, 重力ポテンシャル÷重力加速度
Relative Humidity, 相対湿度, 飽和水蒸気圧に対する水蒸気圧の割合
Specific Humidity, 比湿, 大気密度に対する水蒸気密度の割合
Omega (Vertical Velocity), 鉛直風, dp/dt (圧力の時間微分; 圧力座標での鉛直風)
U-Wind, 東西風,
V-Wind, 南北風,

なお, 上では圧力面上の物理量を保存したファイルの場所を書いていますが, その他の物理量を用いて図を描いても構いません.


==== レポートに書く内容

レポート内には以下を記述してください.
ただしリンクを貼る際は((*相対パスで指定するようにしてください*)).

* スクリプトファイル quiz2.rb.txt へのリンク

  上記で作成した quiz2.rb を quiz2.rb.txt という名前に変更してください. ブラウザ上から読めるようにするためです.
* スクリプトの解説

  スクリプト内の各行が何をしているのか, ((*スクリプト内にコメント文として書く*)).
* 描画に使用するデータファイルへのリンク

  リンク先はダウンロード元の URL を記してください.

* スクリプトの使い方の解説

  これをもとに TA が実際にテストします. 
  実行する際のスクリプト, データの置き場を明示してください.

#* 実行例 result2.txt, result2.png へのリンク
* 作成した図 (動画ファイル) (result2-1.png, result2-2.png, ... 等) へのリンク


* 作成した図の説明

  自分で作成した図が何の図なのかを説明してください.

  図を説明するには, 「いつ」「どこ」「物理量」「平均方法 (瞬間値? 平均値?)」「図の軸」のような情報が必要です.
  例えば, 下のような記述になるでしょう.
    * 2000 年 10 月 1 日 0 UT の 200 hPa 気圧面における東西風の瞬間値です.
    * 2000 年 10 月の 500 hPa 気圧面における温度の 1 カ月平均値です.
    * 2000 年の緯度 45°における海面更生気圧の経度分布の時間変化です.


* 作成した図から読み取れること

  自分で作成した図から読み取れることを説明してください.
  必要に応じて参考文献, 参考ウェブページを引用しても構いません.
  その際には下の「参考資料」の欄に記入してください.

* 参考資料

  URL をベタ書きにするのではなく, そのページの作者やタイトル, および何を参考にしたのかも明記する.
* 共同作業した人の名前
* 工夫したこと, 感想など

#=== 実行結果の例
#
# (作成中)
#
# $ ruby quiz1.rb


==== 作成する図の候補

作成する図はそれぞれの興味に従って選んで構いません.
もし思いつくことがない人は下の項目を参考にするとよいでしょう.

* 東西平均温度/東西風/比湿の緯度-圧力分布 (の季節変化の動画)
* 神戸市 (に最も近い点) における地表 (1000 hPa における) 温度/気圧の時間変化
* 神戸市 (に最も近い点) における温度/東西風/南北風の時間-圧力分布
* YYYY 年 MM 月に発生した台風の温度/風速/比湿分布 (の動画)