#author("2023-11-23T16:07:24+09:00","external:moriat","moriat") #author("2023-11-23T16:08:08+09:00","external:moriat","moriat") #topicpath * 専攻演習(ゼミ) [#af46c217] ** 春学期 [#a0cdef54] → [[講義のページ/Pythonのお勉強]] ** 秋学期 [#i70796e2] *** JRA-3Q [#xfbd606d] - JRA-3Qとは~ 気象庁による新しい長期再解析データ (cf. [[気象庁第3次長期再解析(JRA-3Q)>https://jra.kishou.go.jp/JRA-3Q/index_ja.html]]) - データの入手 -- [[DIAS アカウント申請>https://auth.diasjp.net/account/public/ja/guest/]] -- [[DIAS によるダウンロードページ>https://data.diasjp.net/dl/storages/filelist/dataset:645]] - ファイル名 (cf. [[気象業務支援センターによる説明>http://www.jmbsc.or.jp/jp/offline/hd0025.html]])~ 基本フルセットの場合 -- Const 不変データ -- Clim9120 平年値 -- Hist/Monthly 月平均 -- Hist/Daily 日別~ |ファイル名 | 物理量 | 時間間隔(時刻)|h |anl_surf125| 地表面解析値| 6 時間毎(00, 06, 12, 18UTC) | |anl_snow125| 積雪深解析値| 日毎(18UTC) | |anl_p125| 等圧面解析値 |6 時間毎(00, 06, 12, 18UTC) | |anl_isentrop125| 等温位面解析値 |6 時間毎(00, 06, 12, 18UTC) | |anl_land125| 陸面解析値 |6 時間毎(00, 06, 12, 18UTC) | |fcst_surf125| 2 次元物理量瞬間値 |毎時(00, 01, ..., 23UTC)| |fcst_land125| 陸面予報値 |毎時(00, 01, ..., 23UTC) | |bnd_ocean125| 海面境界条件 |毎時(00, 01, ..., 23UTC) | |fcst_phy2m125| 2 次元物理量平均値 |毎時(00(00〜01UTC の平均値), 01, ..., 23UTC)| |fcst_phyland125| 陸面物理量平均値 |毎時(00(00〜01UTC の平均値), 01, ..., 23UTC) | *** 気象庁予報データ [#z31c54ba] - データについて~ -- 気象庁予報データとは~ 気象庁が各種気象情報を提供するために数値計算を実行している。その実行結果を、気象予報業務を行うものに提供するためのデータ。~ 基本的にはデータ提供に費用がかかりるため、購入するべきもの。しかし、教育・研究目的で公開されている。~ -- MSM~ 全球予報モデル(GSM)に対して、日本近辺を計算したメソスケールモデル(MSM)とがある。 - Pythonでの扱い~ ※ Rで扱う場合には、[[2022年度専攻演習のページ>../Seminar2022]]を参照のこと。 *** Python と Jupyter Notebook [#s62fef16] - 背景~ Google Colaboratory は便利でパワフルだ。しかし、例えば長い処理を行うときにはタイムアウトを心配しなければならない。~ そこで、似たような環境を手元のノートパソコンで実現する。~ ~ - インストール -- Python(アップグレードも同じ手続きでできる。)~ --- [[Python.org>https://www.python.org/]] から、Downloads > Windows を選択し、Latest Python 3 Release - Python 3.X.X をクリックする。 --- 環境に合わせて、インストーラーをダウンロードする。64bit なら、Windows installer (64-bit) --- ダウンロードされたファイルを実行する。~ □ Add python.exe to PATH にチェックを入れておく。~ そのうえで、"Install Now" をクリックする。 -- Jupyter Notebook~ --- Python のインストールを終えたあとで、コマンドプロンプトを開く。管理者である必要はない。 --- インストールのコマンドを入力。 pip install notebook --- "pip" の新しいバージョンがある、と、(英語で)表示されたら支持に従って更新しておく。 -- wgrib2~ 気象関係のデータを扱うために wgrib2 を導入する。~ wgrib2 は、grib2 形式のファイルを netCDF に変換するために用いる。~ pygrib というパッケージでも使えそうだが、インストールに失敗した。~ そこで、安心して作業するために、原因追求ではなく、wgrib2 を使うことに。 --- ダウンロード~ [[NOAAのサイト>https://www.ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/]] から、Windows10 > v3.1.3 (最新版) へ移動し、wgrib2.exe をダウンロードする。 --- 移動~ コマンドプロンプトで、 path を実行すると、コマンド名だけでプログラムを実行できる場所が表示される。例えば、~ C:\Users\ユーザー名\AppData\Local\Programs\Python\Python312\Scripts\ があるだろう。ここに、wgrib2.exe を移動する。エクスプローラーを開いて、上のパス(ディレクトリ名)をコピペしてそこを開き、貼り付ける。 -- Cygwin~ wgrib2.exe は、cygwin で提供されるライブラリ(DLL)を利用している。 --- ダウンロード~ [[Cygwinのサイト>https://cygwin.com/install.html]]からsetup-x86_64.exe をダウンロードして実行~ 途中、どのサイトからダウンロードするかを尋ねられる。jp で終わるものを選ぶ。 --- パッケージの選択~ パッケージの選択画面(Select Packages)になったら、View を Full にして、 gcc-fortran gcc-core を skip から新しめのバージョンを選ぶ。次へ を繰り返してインストールする。 --- パスの設定~ コントロールパネル > システムの詳細設定の表示 > 環境変数 > Path を選んで編集~ 新規で C:\cygwin64\bin を追加 -- Visual Studio~ 通常必要ないが、Cartopy のようなパッケージを導入する場合に必要になる。 --- ダウンロード~ [[Microsoftの公式サイト>https://visualstudio.microsoft.com/ja/downloads/]]からVisual Studioのインストーラーをダウンロードする。2022 が最新? --- インストーラーの実行~ ダウンロードしたインストーラーを実行。ワークロードの選択では、「C++によるデスクトップ開発」を選択。~ インストールボタンをクリックする。 ~ -- パッケージのインストール --- pip install pickleshare requests cartopy metpy netcdf4 ~ - 実行 -- Jupyter Notebook の実行 --- コマンドプロンプトで次のコマンドを実行する。 jupyter notebook ~ --- ブラウザを選択するように指示されたら、適当なブラウザを選ぶ。~ ユーザーインターフェースは、ブラウザで実現される。 --- Google Colaboratory と同様に、プログラムはセルごとに表示される。Shift+Enterで実行する。~ 実行中は、ブラウザの画面右上の表示が Python 3 (ipykernel) ○ から Python 3 (ipykernel) ● に変化する。実行が終わると、再び○になる。 ~ -- ipynb ファイルをダブルクリックして開く~ Google Colaboratory のファイルをダウンロードすると、ipynb という拡張子がついたファイルがダウンロードされる。~ これをダブルクリップして開くようにしたい。~ ・ jupyter-notebook.exe を探す。~ --- jupyter-notebook.exe を探す。~ 例: C:\Users\mori\AppData\Local\Programs\Python\Python312\Scripts\ ・ 当該ファイルがあるフォルダーをエクスプローラで開き、ファイルにポインタを合わせて右クリック~ ・ プログラムから開く > 別のプログラムを選択 > その他のアプリ で、上の jupyter-notebook.exe を指定。~ --- 当該ファイルがあるフォルダーをエクスプローラで開き、ファイルにポインタを合わせて右クリック~ --- プログラムから開く > 別のプログラムを選択 > その他のアプリ で、上の jupyter-notebook.exe を指定。~ このとき、□ 常にこのアプリを使って… にチェックを入れておく。~ ~ *** データを扱うサンプルプログラム [#c3fde59d] | 目的 | Google Colaboratory | Jupyter Notebook |h |レーダーエコー強度を取得して図にする|CENTER: [[●>https://colab.research.google.com/drive/1hlsyg7GOO8_SgRKndCmkwHe4LxqY-EF5?usp=sharing]] |CENTER: &ref(20231115_レーダーエコー強度の取得と表示.ipynb,,●);| | 気象庁予報データの風速分布を図にする |CENTER: [[●>https://colab.research.google.com/drive/19sKiuSWS8ktugy2EQB_H3gX_hwP2RZui?usp=sharing]] |CENTER: &ref(20231117_気象庁の予報データから水平風速場を描く.ipynb,,●);| | 静止画(png)を動画(mp4)に変換する|CENTER:[[●>https://colab.research.google.com/drive/1HUiftB7Z8dHd_ZvMSbZgvhuDrJVehNIU?usp=sharing]] |CENTER: &ref(20231122_静止画pngファイルから動画mp4ファイルを作る.ipynb,,●);| ** 過去の専攻演習(ゼミ)のページ [#vb16f810] -- [[2022年度>../Seminar2022]] -- [[2021年度>../Seminar2021]] -- [[2015年度>../Seminar2015]] -- [[2014年度>../Seminar2014]] -- [[2013年度>../Seminar2013]] -- [[2012年度>../Seminar2012]] -- [[2009年度>../Seminar2009]]