###################################################################
under construction(まだ落書き程度です…。)
###################################################################
gmtsetを結構使いますね。
PLOT_CLOCK_FORMATやらPLOT_DATE_FORMATを使って文字の表示をいろいろ
選べます。変更した変数を-Bオプションから読みこみます、これは従来の
GMTと一緒の概念ですが、やたらと多い…。
たとえば時間、7時19分を07:19と表示させたい場合には
# gmtset PLOT_CLOCK_FORMAT hh:mmとする。
日にち8月23日をAug 23と表示したい場合には
# gmtset PLOT_DATE_FORMAT "o dd" として、-Bオプションでaを使う。
以下を参照。
Table 4.3:
I and i interval annotation modifiers.
| Modifier |
Description |
| f (or F) |
Full name of interval, e.g., January (or JANUARY) |
| a (or A) |
Short abbreviation of interval, e.g., Jan (or JAN) |
| c (or C) |
1-char abbreviation of interval, e.g., J (or J) |
|
とりあえずデータのフォーマットは
# cat data.txt
2001-08-31T07:19:16 2.2546 34.3636 0.9372 27.4465
2001-08-31T14:28:21 2.3100 34.3061 0.8411 27.3960
2001-08-31T20:39:11 2.3335 34.2101 0.9488 27.3172
2001-09-1T03:23:59 2.3308 34.2230 0.9452 27.3278
2001-09-1T10:53:00 2.2313 34.3656 1.0161 27.4500
2001-09-1T15:51:23 2.3109 34.3102 0.8381 27.3992
こんな感じにする。
#!/bin/zsh -f
mintime=`sed -n '1p' data.txt |awk '{print $1}'`
maxtime=`sed -n '$p' data.txt |awk '{print $1}'`
ranged=$mintime/$maxtime/26.6/27.8
gmtset PLOT_DATE_FORMAT "o dd" PLOT_CLOCK_FORMAT hh:mm
gmtset HEADER_FONT 5 ANOT_FONT 5 LABEL_FONT 5 ANOT_FONT_SIZE 15
psbasemap -R$ranged -JX15T/-10 -B":.B2 @~s@-q@-@~(time series):" -P -K > $psfile
awk '{print $1,$5}' data.txt |\
psxy -R$ranged -JX -Ba6Hf1hIa1D/a0.2nSEw -O -Wto5/255/0/0 >>$psfile
-Ba6Hf1hIa1D の説明(不充分)
a6:6時間ごとに数字を付ける
Hf1:gmtsetで指定したPLOT_CLOCK_FORMATを読みこみ1時間毎にticを付ける
h:Plot as 2-digit interger(0-24)
I:lower interval annotation spacing
a1:Short abbreviationi of interval e.g. Augそれを１か月後とに表示
D:Format annotation using PLOT_DATE_FORMAT
