fig=fig_model.eps
gmt psxy -R-2/4/-1/4 -JX6/5 -G#c0c0c0 -P -K << EOT > $fig
 -1.0 -0.5
 -1.0  3
 -1.5  3
 -1.5 -0.5
EOT
gmt psxy -R -JX -G#cccccc -K -O << EOT >> $fig
 -0.5  0
  0.5  0
  0.5 -0.3
 -0.5 -0.3
EOT
gmt psxy -R -JX -W1 -K -O << EOT >> $fig
 -1.0 -0.5
 -1.0  3
EOT
gmt psxy -R -JX -W1 -K -O << EOT >> $fig
 -0.5  0
  0.5  0
EOT
awk 'BEGIN{a=1/4;for(y=0;y<=2;y=y+0.1){print a*y*y,y}}' | gmt psxy -R -JX -W2 -K -O >> $fig
echo ' 1 2' | gmt psxy -R -JX -SC0.4 -G0 -K -O >> $fig
echo ' 0.0  0.0  0.0  3.0' | gmt psxy -R -J -Sv0.25+s -W0.25 -G0 -O -K >> $fig
echo ' 1.0  2.0  1.0  3.0' | gmt psxy -R -J -Sv0.25+s -W0.25 -G0 -O -K >> $fig
echo '-1.0  2.4  0.0  2.4' | gmt psxy -R -J -Sv0.2+s+b+e+a60 -W0.25 -G0 -O -K >> $fig
echo ' 0.0  2.4  1.0  2.4' | gmt psxy -R -J -Sv0.2+s+b+e+a60 -W0.25 -G0 -O -K >> $fig
gmt pstext -R -JX -F+f+a+j -N -K -O << EOT >> $fig
-0.5 2.6 12p 0 BC @~f@~
 0.5 2.6 12p 0 BC u
EOT
echo '0 0' | gmt psxy -R -J -Sp -O >> $fig