大量の図をbashスクリプトで作図するときのテンプレート

draw-ncl.sh

#!/bin/bash

# draw a figure
function draw_figure(){
  local ifile=$1                                # local にしておくと、うっかり同じ変数名を
  local mon=$2                                  # 使ってしまったときも安心
  local ofile=$3
  echo "draw: $ifile, $ofile, $mon"

cat > tmp.ncl << EOF                            # ここからEOFまでをtmp.nclというファイルに入れる
begin
  f = addfile=("${ifile}","r")                  # ここで ifile で指定したファイルを開いている
  u = f->v(${mon},:,:)                          # mon で指定した時間を切り出し
  wks  = gsn_open_wks("pdf","figure.${mon}")    # 出力ファイルの名前も mon をいれて定義
  plot = gsn_contour(wks,u,False)
end
EOF
  ncl tmp.ncl                                   # ncl を起動
  rm tmp.ncl                                    # tmp.ncl を削除
}


# main
for m in `seq 1 12`; do                         # 1から12のループ
  draw_figure inputdata.nc $m fig-monthly       # 出力は fig-monthly.1.pdf...みたいになる
done