#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)