#! /bin/tcsh -f

set file=20011217 <----file name
set file0=$file.5h_uv <---- 実際のデータ
set fileps=$file.ps <---出力するps file の名前
set range=135/146/42/50 <---描く図の緯度経度範囲
set scale=M10 <---図の大きさ
#------------------ここでべくとるの太さなどの指定

set futosa=0.03
set tatenaga=0.1
set yokonaga=0.1
set type=-Sv$futosa/$tatenaga/$yokonaga
###topo###
set topo=map_gmt <--- 地形のデータ
set resol=0.05 <--- 緯度経度データの解像度
set color_topo=cpt_topo <---カラーパレットの名前
set gfile_topo=$topo.grd
#############################################
#山地用のカラーパレットを作る
makecpt -T-100/2500/100 -Ctopo -V > $color_topo
#地形データをGMT用のバイナリデータに落す($gfile_topoができたバイナリ)
xyz2grd $topo -G$gfile_topo -I$resol -R$range -V
#山地をプロット
grdimage $gfile_topo -J$scale -R$range -C$color_topo -Y7 -X4 -P -K >! $fileps
#ベクトルを描く $typeが上で設定したベクトル
psxy $file0 -J$scale -R$range $type -W0.1 -Ba2/a2:."$file0":neWs -G255/100/100 -P -O -K -V >> $fileps
# 右側に山地のスケールを書く
psscale -C$color_topo -Ba500f100::/:"[m]": -D5/5/7/0.4 -X6 -O >> $fileps
gs $fileps

ベクトルを書くときのデータセットは以下のようなものです。
20011217.5h_uvのデータの配列は

経度 緯度 角度(北を０度で時計回り；単位degrees) ベクトルの長さ
131.125 42.125 295.623 0.800
131.375 42.125 295.575 0.750
131.375 42.375 295.900 0.648
131.625 42.125 300.987 0.733
131.625 42.375 307.980 0.597
131.625 42.625 325.951 0.483
131.875 42.125 299.006 0.614
131.875 42.375 300.557 0.480
131.875 42.625 300.850 0.423
132.125 42.625 301.110 0.514
132.375 42.125 300.217 0.565
132.375 42.375 300.919 0.520
132.625 42.125 301.144 0.599
132.625 42.375 301.264 0.608
132.875 42.375 301.364 0.659
133.125 42.125 301.314 0.621








モドリマスカ？