rem ******************************************** rem Fourier spectrum ratio (normal-log chart) rem ******************************************** rem Data file format (CSV data) rem freq : Frequency rem fsp1 : Smoothing spectrum of numerator rem fsp2 : Smoothing spectrum of denominator rem fsr : Spectrum ratio (fsp1/fsp2) rem ******************************************** rem Drawing by GMT (response spectrum) rem ******************************************** set rangeW=0/20/0.1/1000 set rangeE=0/20/1/10000 set scale=12/10l set xga=g5a5 set yga=g1a1 set xlabel="Frequency (Hz)" set ylabelW="Spectrum ratio" set ylabelE="Fourier spectrum (gal*sec)" gmtset ANOT_FONT_SIZE 12 gmtset LABEL_FONT_SIZE 12 gmtset TICK_LENGTH 0c rem =============================== rem Setting of files (EW) rem =============================== set inpl=inp_legend_fsr_EW.txt set inp_1=inp_fsr_t.csv set fig=fig_fsr_EW1.eps rem =============================== rem Prerapation for legend drawing rem =============================== echo 12 $ 10 $ tr $ 0.3 $ 0.3 > %inpl% echo 10 $ 0.5 $ 5 >> %inpl% echo TCGH16_EW2 $ 5t30_10_10_10:0 $ --- $ --- $ --- >> %inpl% echo TCGH16_EW1 $ 5t10_10:0 $ --- $ --- $ --- >> %inpl% echo Spectrum ratio $ 5 $ --- $ --- $ --- >> %inpl% gawk -f \WANtaroHP_plot\gmt_basic\awk_klegend.awk %inpl% > legend.bat rem =============================== rem Plotting rem =============================== psbasemap -R%rangeW% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabelW%:WSen -P -X5 -Y6 -K > %fig% gawk "BEGIN{FS=\",\"}{if(2<=NR)print $2,$5}" %inp_1% | psxy -R -J -B -W5 -P -O -K >> %fig% psbasemap -R%rangeE% -JX%scale% -B::/%yga%:%ylabelE%:E -P -O -K >> %fig% gawk "BEGIN{FS=\",\"}{if(2<=NR)print $2,$4}" %inp_1% | psxy -R -J -B -W5t30_10_10_10:0 -P -O -K >> %fig% gawk "BEGIN{FS=\",\"}{if(2<=NR)print $2,$3}" %inp_1% | psxy -R -J -B -W5t10_10:0 -P -O -K >> %fig% rem ================================ rem Drawing of Peak frequency (no -N option) rem ================================ gawk -f awk_ssmv0.awk %inp_1% > _temp.txt pstext _temp.txt -R%rangeW% -J -P -O -K >> %fig% gawk -f awk_ssmv1.awk %inp_1% > _temp.txt pstext _temp.txt -R%rangeW% -J -P -O -K >> %fig% rem *** marginal note *** echo 0 1200 10 0 0 BL Band=1.0Hz | pstext -R%rangeW% -J -N -P -O -K >> %fig% rem =============================== rem Drawing of legend rem =============================== call legend.bat rem ******************************************** rem Drawing by GMT (response spectrum) rem ******************************************** set rangeW=0/20/0.1/1000 set rangeE=0/20/0.1/1000 set scale=12/10l set xga=g5a5 set yga=g1a1p set xlabel="Frequency (Hz)" set ylabelW="Spectrum ratio" set ylabelE="Fourier spectrum (gal*sec)" gmtset ANOT_FONT_SIZE 12 gmtset LABEL_FONT_SIZE 12 gmtset TICK_LENGTH 0c rem =============================== rem Setting of files (EW) rem =============================== set inpl=inp_legend_fsr_EW.txt set inp_1=inp_fsr_t.csv set fig=fig_fsr_EW2.eps rem =============================== rem Plotting rem =============================== gawk -f \WANtaroHP_plot\gmt_basic\awk_klegend.awk %inpl% > legend.bat psbasemap -R%rangeW% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabelW%:WSen -P -X5 -Y6 -K > %fig% gawk "BEGIN{FS=\",\"}{if(2<=NR)print $2,$5}" %inp_1% | psxy -R -J -B -W5 -P -O -K >> %fig% psbasemap -R%rangeE% -JX%scale% -B::/%yga%:%ylabelE%:E -P -O -K >> %fig% gawk "BEGIN{FS=\",\"}{if(2<=NR)print $2,$4}" %inp_1% | psxy -R -J -B -W5t30_10_10_10:0 -P -O -K >> %fig% gawk "BEGIN{FS=\",\"}{if(2<=NR)print $2,$3}" %inp_1% | psxy -R -J -B -W5t10_10:0 -P -O -K >> %fig% rem ================================ rem Drawing of Peak frequency (no -N option) rem ================================ gawk -f awk_ssmv0.awk %inp_1% > _temp.txt pstext _temp.txt -R%rangeW% -J -P -O -K >> %fig% gawk -f awk_ssmv1.awk %inp_1% > _temp.txt pstext _temp.txt -R%rangeW% -J -P -O -K >> %fig% rem *** marginal note *** echo 0 1200 10 0 0 BL Band=1.0Hz | pstext -R%rangeW% -J -N -P -O -K >> %fig% rem =============================== rem Drawing of legend rem =============================== call legend.bat set range= set scale= set xga= set yga= set xlabel= set ylabel= set inpl= set inp_1= set fig= del .gmt* del _*