#author("2019-08-26T13:49:46+09:00","default:Miyashita","Miyashita")
#author("2019-08-26T13:50:14+09:00","default:Miyashita","Miyashita")
*Julia演習問題:NetCDFデータの読み込み [#t46f4d98]
***問題 [#icd512da]
+NetCDF ファイル,wspd.mon.mean.nc を読み込みなさい.
#codeprettify(lang-julia){{
download(url, destination)
}}
+ncファイルから変数 lon, lat, time, wspd を取得しなさい.
#codeprettify(lang-julia){{
var_name = NetCDF.ncread(ncfile,"var_name")
}}
+変数 time を DateTime 型に変換しなさい.~
(補足1)変数timeの値は1800年1月1日0時0分0秒からの経過時間(単位:hour)である.~
#codeprettify(lang-julia){{
Dates.DateTime(1800,1,1)
Dates.Hour(...)
}}
+変数 wspd の最初の時刻の塗りつぶしコンター図を描きなさい.
#codeprettify(lang-julia){{
Plots.contour(lon, lat, wspd[:,:,1], fill=true, ...) # Plotsの場合
}}
+複数のステップを繰り返し描き,動画を作成しなさい.
#codeprettify(lang-julia){{
anim = Plots.@animate for k=1:n
    plot(xxx,k)
end
Plots.gif(anim, "amination_wspd.gif", fps=3)
}}
***データ [#f2f57e99]
-NCEP 再解析データの月平均高度 10m 風速データ(京大生存圏データベース)~
#ref(http://database.rish.kyoto-u.ac.jp/arch/ncep/data/ncep.reanalysis.derived/surface/wspd.mon.mean.nc)
***解答例(参考) [#xa28f506]
1〜24 ステップの 2 年間分の動画を作成してみた.~
Plotsではcolorbarの目盛り位置,目盛りラベルの調整がうまくできなかった.~
他にもPyPlot,GMTのパッケージで作成してみた.GMTでは比較的美しくできた.
Plots では colorbar の目盛り位置,目盛りラベルの調整がうまくできなかった.~
他にも PyPlot, GMT のパッケージで作成してみた.GMT では比較的美しくできた.
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/julia/MonthlyWspd.gif,716x387)

Front page   Edit Diff Attach Copy Rename Reload   New List of pages Search Recent changes   Help   RSS of recent changes