鉛直方向は分解能が150mより細かい（はず）。 1次元スパイラルで補完できるとよいかもしれない。

○PointデータからGridデータへ変換 GMTではレンジビンを1~534で計算。 札幌までの距離が13㎞なので、札幌で座標０㎞となるように 距離＝13000ー距離 でX座標を設定。 出力範囲（xmin:xmax）を±6㎞で指定。 ○GradsでNetCDFを読む。 ctlではXの初期値を-6000と…