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 風速データ(京大生存圏データベース)

解答例(参考)

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

MonthlyWspd.gif

Front page   Edit Diff Attach Copy Rename Reload   New List of pages Search Recent changes   Help   RSS of recent changes
Last-modified: 2018-07-12 (Thu) 12:50:51 (693d)