NCLで一番ありがちな「次元が違います」「メタデータがありません」エラーをチェックする

ncl-dimension-check.ncl

begin
  f1 = addfile("sample1.nc","r")
  f2 = addfile("sample2.nc","r")

  v1 = f1->v(0,:,:)             ; 0 にした次元は落ちるので、2次元になる
  v2 = f2->v(0,:,:)

  printVarSummary(v1)           ; 次元を表示してチェック

  diff = (v2 - v1) - 273.15     ; 算数をするとメタデータが消えたりするので
  copy_VarCoords(v1, diff)      ; 他の変数からコピーしてあげる
end

・とにかく悩んだら printVarSummary してチェックする
・メタデータが欠落したらてきとうなところから copy_VarCoords してあげる