用GMT画M-t图

vfory

1  时间类型的数据格式必须为     yyyy-mm-ddT[HH]:[MM]:[SS.xxx]，中间的那个T表示为时间格式的数据，必需的；
2  在利用-R参数指定绘图区间时，同样也要用到这样的时间格式；
3 在利用-B参数指定边框刻度时，注意需要根据所绘的时间长度来指定边框刻度的标识，常用的有Y(年)、O（月）、d（日）、H（时）、M（分）、C（秒）
4 一般利用JX[x]投影方式
 
#!/usr/bin/perl -w
open(SRC,"wenchuan.eqt");
@aa=<SRC>;
chomp @aa;
close SRC;

$start_t="2008-05-01T";  #开始时间
$end_t="2011-10-30T";  #结束时间
$start_m=0;          #最小震级
$end_m=8;          #最大震级
$region="$start_t/$end_t/$start_m/$end_m";
`psbasemap -JX7i/4i -R$region -Ba1Yf4O:"time":/1:"magnitude":WSne -K -P > mt.ps`;  #每年标注一次，每4个月画一个刻度
open(REC,">tmp");
foreach $a(@aa){
  $yr=substr($a,1,4);
  $mm=substr($a,5,2);
  $dd=substr($a,7,2);
  $HH=substr($a,9,2);
  $MM=substr($a,11,2);
  $SS=substr($a,13,2);
  $mag=substr($a,28,3);
  $line1=$yr."-".$mm."-".$dd."T".$HH.":".$MM.":".$SS."   $mag";
  $line2=$yr."-".$mm."-".$dd."T".$HH.":".$MM.":".$SS."   $start_m";
  print REC "$line1\n$line2\n>\n";
}
close REC;
`psxy tmp -R -JX -m -K -P -O >> mt.ps`;