GMT で地震動の加速度応答スペクトルをグラフ化した事例です.
入力データファイルは以下の書式となっています.
1〜2行目は情報,3行目は列の解説,4行目以降が数値データです.
2011/03/11 14:46:00 TCGH16 EW2 Max.acc=1196.697 nt,200 period(sec),acc(gal),vel(kine),dis(cm) 2.000000e-002,1.196764e+003,2.924900e-001,1.213241e-002 2.063122e-002,1.198466e+003,3.113090e-001,1.292889e-002 2.128236e-002,1.200244e+003,3.314467e-001,1.377853e-002 ・・・・・ |
1列目は固有周期,2列目は加速度応答値,2列目は速度応答値,3列目は変位応答値です.
rem ============================================ rem Drawing by GMT (response spectrum) rem ============================================ set range=0.01/10/10/10000 set scale=12l/10l set xga=g3a1 set yga=g3a1 set xlabel="Period (sec)" set ylabel="Response acceleration (gal)" rem ************************************************************************** rem * set range=0.01/10/10/10000 rem * x軸範囲を0.01〜10,y軸範囲を10〜19999に設定 rem * set scale=12l/10l rem * x軸を長さ12cmの対数軸,y軸を長さ10cmの対数軸に設定 rem * set xga=g3a1 x軸目盛を設定 rem * set yga=g3a1 y軸目盛を設定 rem * set xlabel="Period (sec)" x軸ラベルを設定 rem * set ylabel="Response acceleration (gal)" y軸ラベルを設定 rem ************************************************************************** gmtset ANOT_FONT_SIZE 12 gmtset LABEL_FONT_SIZE 12 gmtset TICK_LENGTH 0c rem ============================================ rem Drawing of TCGH16 rem ============================================ set inpl=inp_legend_t.txt set inp_1=dat_spc_TCGH16_EW2.csv set inp_2=dat_spc_TCGH16_NS2.csv set inp_3=dat_spc_TCGH16_UD2.csv set fig_out=fig_spc_t.eps rem ************************************************************************** rem * set inpl=inp_legend_t.txt 凡例作成用入力ファイル rem * set inp_1=dat_spc_TCGH16_EW2.csv 入力ファイル名1 rem * set inp_2=dat_spc_TCGH16_NS2.csv 入力ファイル名2 rem * set inp_3=dat_spc_TCGH16_UD2.csv 入力ファイル名3 rem * ここでは3本の加速度応答スペクトルを描画するので3入力ファイルを指定 rem * set fig_out=fig_spc_t.eps 出力画像ファイル名 rem ************************************************************************** rem ------------------------------- rem Making of batch file for legend rem ------------------------------- gawk -f awk_klegend.awk %inpl% > legend.bat rem ************************************************************************** rem * gawk -f awk_klegend.awk %inpl% > legend.bat rem * 凡例出力用バッチファイルの作成(awkスクリプト実行) rem ************************************************************************** rem ------------------------------- rem Plotting rem ------------------------------- psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -X5 -Y6 -K > %fig_out% gawk "BEGIN{FS=\",\"}{if(3 |
ファイル名 | 概要 |
---|---|
bat_spc.txt | 作図用バッチファイル |
inp_legend_t.txt | 凡例作成用入力ファイル |
legend.txt | awk_klegendで作成した凡例描画用バッチファイル |
fig_spc_t.png | 出力画像 |