fig=fig_model_1.eps
gmt psxy -R-2/5/-1/2 -JX7/3 -W0.5 -P -K << EOT > $fig
>
0.0 0.0
4.5 0.0
>
0.0 0.0
0.0 1.5
EOT
gmt psxy -R -JX -G#cccccc -W1 -K -O << EOT >> $fig
0 -0.1
0  0.1
3  0.1
3 -0.1
EOT
awk 'BEGIN{pi=3.141592654;for(i=90;i<=270;i++){print 0.5*cos(i/180*pi)-0.7,0.5*sin(i/180*pi)}}' | gmt psxy -R -J -W2 -K -O >> $fig
awk 'BEGIN{pi=3.141592654;for(i=-90;i<=90;i++){print 0.5*cos(i/180*pi)+3.7,0.5*sin(i/180*pi)}}' | gmt psxy -R -J -W2 -K -O >> $fig
gmt set MAP_VECTOR_SHAPE 0
gmt psxy -R -J -Sv0.2+s+e+a45  -G0 -W1.5 -O -K << EOT >> $fig
 0.0  0.0  0.0  1.0
 3.0  0.0  3.0  1.0
-0.7 -0.5 -0.5 -0.5
 3.7  0.5  3.5  0.5
-0.7  0.0  0.0  0.0
 3.0  0.0  3.7  0.0
EOT
gmt pstext -R -J -F+f+a+j -N -K -O << EOT >> $fig
 4.7  0.0 10p,Helvetica 0 MC x
 0.0  1.8 10p,Helvetica 0 MC y
 0.2  1.0 10p,Helvetica 0 ML S@-i@-,v@-i@-
 2.8  1.0 10p,Helvetica 0 MR S@-j@-,v@-j@-
-1.2  0.5 10p,Helvetica 0 MR M@-i@-,@~q@~@-i@-
 4.1 -0.5 10p,Helvetica 0 ML M@-J@-,@~q@~@-j@-
-0.15  0.3 10p,Helvetica 0 MR N@-i@-,u@-i@-
 3.1 -0.2 10p,Helvetica 0 ML N@-j@-,u@-j@-
 1.5  -0.4 10p,Helvetica 0 MC L
EOT
echo '0 0' | gmt psxy -R -J -Sp -O >> $fig


fig=fig_model_2.eps
gmt psxy -R-2.5/3.5/-1/4 -JX6/5 -SC0.2 -G0 -P -K << EOT > $fig
1.5 2.0
EOT
gmt psxy -R -J -Sw -W0.5 -K -O << EOT >> $fig
0 0 1.7 0 30
0 0 1.6 0 30
EOT
gmt psxy -R -J -Sv0.2+s -G0 -W0.5 -K -O << EOT >> $fig
 0     0      3     0
 0     0      0     3
 1.5   0      1.5   2
 0     2      1.5   2
 0     0      3     1.732
 0     0     -1.732 3
 1.991 1.150  1.5   2
-0.491 0.850  1.5   2
EOT
gmt set MAP_VECTOR_SHAPE 0
gmt psxy -R -J -Sv0.2+s+e+a45  -G0 -W1.5 -O -K << EOT >> $fig
0  0   1.5   0.0
0  0   0.0   2.0
0  0   1.991 1.150
0  0  -0.491 0.850
EOT
gmt pstext -R -J -F+f+a+j -N -K -O << EOT >> $fig
 3.3  0.0 10p,Helvetica  0 MC X
 0.0  3.3 10p,Helvetica  0 MC Y
 1.5 -0.3 10p,Helvetica  0 MC U
-0.3  2.0 10p,Helvetica  0 MC V
 3.2  1.8 10p,Helvetica 30 MC x
-1.8  3.2 10p,Helvetica 30 MC y
 2.2  1.0 10p,Helvetica 30 MC u
-0.7  0.8 10p,Helvetica 30 MC v
 1.0  0.3 10p,Helvetica  0 MC @~f@~
EOT
echo '0 0' | gmt psxy -R -J -Sp -O >> $fig