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 (313d)