nclでwrfのアウトプットを高度面、気圧面に内挿する時はループしなくて良い

nclはスクリプト言語で、ループの数を減らすことが実行時間の短縮につながるのだが…

参考URL
wrf_user_intrp3d https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_intrp3d.shtml

今までこの関数を使う時、馬鹿正直に鉛直層数分ループして、内挿させていたが、ループしなくていいじゃないか（公式もそう書いてあった）

つまり...

height = (/1800,2000,2200/) ;meter
n_height = dimsizes(height)
do iz = 0,n_height-1,1
vv_levs(iz,:,:)   = wrf_user_intrp3d(vv(:,:,:),zz(:,:,:),"h", height(iz),0.,False)
end do

ではなく、

height = (/1800,2000,2200/) ;meter
vv_levs(:,:,:) = wrf_user_intrp3d(vv(:,:,:),zz(:,:,:),"h", height,0.,False)

で、良い。笑
これで一つループを減らせる。

NCLでWRFの"contour on contour"のoverlay

参考URL
見失ってしまった。

NCLでWRFの関数やらを使って描画しようとして、コンターとシェイドを重ねたい時に、忘れるのが

pltres@PanelPlot = True
plot = wrf_map_overlays(a,wks,(/contour01,contour02/),pltres,mpres)

PanelPlotの記述。
本来は複数のパネルを作成する時のオプションだけど、wrf_map_overlayの類を使って、重ねる時は、Trueを宣言する必要がある。
忘れがち。