Takuya Miyashita

Julia演習問題:NetCDFデータの読み込み †

問題 †

1. NetCDF ファイル，wspd.mon.mean.nc を読み込みなさい．

   download(url, destination)

2. ncファイルから変数 lon, lat, time, wspd を取得しなさい．

   var_name = NetCDF.ncread(ncfile,"var_name")

3. 変数 time を DateTime 型に変換しなさい．
   (補足1)変数timeの値は1800年1月1日0時0分0秒からの経過時間(単位:hour)である．
   

   Dates.DateTime(1800,1,1)
   Dates.Hour(...)

4. 変数 wspd の最初の時刻の塗りつぶしコンター図を描きなさい．

   Plots.contour(lon, lat, wspd[:,:,1], fill=true, ...)　# Plotsの場合

5. 複数のステップを繰り返し描き,動画を作成しなさい.

   anim = Plots.@animate for k=1:n
       plot(xxx,k)
   end
   Plots.gif(anim, "amination_wspd.gif", fps=3)

データ †

• NCEP 再解析データの月平均高度 10m 風速データ(京大生存圏データベース)
  
  wspd.mon.mean.nc

解答例(参考) †

1〜24 ステップの 2 年間分の動画を作成してみた．
Plots では colorbar の目盛り位置，目盛りラベルの調整がうまくできなかった．
他にも PyPlot, GMT のパッケージで作成してみた．GMT では比較的美しくできた．