set range=0/25/0/1 set scale=5/5 set xga=g5a5 set yga=g0.2a0.2 set xlabel="Discharge (m@+3@+/s)" set ylabel="Height of Water surface (h/D)" set inp_dat=inp_manning_circ1.txt set fig=fig_manning_circ1.eps gawk "$1==0{printf \"0.98 0.25 8 0 0 BR D=%%.3fm\n\",$2}" %inp_dat% > _txtdat1.txt gawk "$1==0{printf \"0.98 0.15 8 0 0 BR i=%%.5f\n\",$3}" %inp_dat% >> _txtdat1.txt gawk "$1==0{printf \"0.98 0.05 8 0 0 BR n=%%.5f\n\",$4}" %inp_dat% >> _txtdat1.txt gawk "$1==-1{printf \"0.02 0.95 8 0 0 TL Qmax=%%.3fm@+3@+/s\n\",$4}" %inp_dat% > _txtdat2.txt gawk "$1==-1{printf \"0.02 0.85 8 0 0 TL (h/D=%%.3f)\n\",$2}" %inp_dat% >> _txtdat2.txt gmtset ANOT_FONT 0 gmtset ANOT_FONT_SIZE 9 gmtset LABEL_FONT 0 gmtset LABEL_FONT_SIZE 9 gmtset LABEL_OFFSET 0.1c gmtset TICK_LENGTH 0c rem psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -UManning_H-Q -P -K > %fig% psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -K > %fig% gawk "0<$1{print $4,$2}" %inp_dat% | psxy -R -JX -W3 -K -O >> %fig% pstext _txtdat1.txt -R0/1/0/1 -JX%scrtxt% -N -K -O >> %fig% pstext _txtdat2.txt -R0/1/0/1 -JX%scrtxt% -N -K -O >> %fig% echo 0.50 1.05 10 0 5 BC H-Q of circular cross section| pstext -R0/1/0/1 -JX%scrtxt% -N -O >> %fig%