Shell Programming and Scripting

#data: 14

#! /bin/bash
# 
# read each block and output to a file
awk '
$1 == ">" {
        if(f) close(f)
        f="file" ++n".txt"
        next
}
{       print > f
}' UNZA2250.txt
# count the # of block to plot
nx=`awk '{c+=gsub(s,s)}END{print c}' s='>' UNZA2250.txt`
echo $nx
# --------------------------------------
gnuplot << EOF
# set size 1.7,1.3
set terminal postscript eps enhanced color "Helvetica" 30
set output "fig.eps"
set ylabel "TEC"
set xlabel "hour"
set xrange [0:24]
set yrange [0:120]
set nokey

filename(n) = sprintf("file%d.txt", n)
plot for [i=1:$nx] filename(i) using 1:4 with lines
EOF
# view the output image
gv fig.eps &
# 
# remove the text files used for ploting
rm file*.txt