1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
|
REM GMT EXAMPLE 30
REM $Id$
REM
REM Purpose: Show graph mode and math angles
REM GMT progs: gmtmath, psbasemap, pstext and psxy
REM DOS calls: del, echo
REM
echo GMT EXAMPLE 30
set ps=example_30.ps
gmt set PS_CHAR_ENCODING Standard+
gmt psbasemap -R0/360/-1.25/1.75 -JX8i/6i -Bx90f30+u"\312" -By1g10 -BWS+t"Two Trigonometric Functions" -K --MAP_FRAME_TYPE=graph --MAP_VECTOR_SHAPE=0.5 > %ps%
REM Draw sine an cosine curves
gmt math -T0/360/0.1 T COSD = | gmt psxy -R -J -O -K -W3p >> %ps%
gmt math -T0/360/0.1 T SIND = | gmt psxy -R -J -O -K -W3p,0_6:0 --PS_LINE_CAP=round >> %ps%
REM Indicate the x-angle = 120 degrees
echo 120 -1.25 > tmp
echo 120 1.25 >> tmp
gmt psxy -R -J -O -K -W0.5p,- tmp >> %ps%
echo 360 1 18p,Times-Roman RB x = cos(@%%12%%a@%%%%) > tmp
echo 360 0 18p,Times-Roman RB y = sin(@%%12%%a@%%%%) >> tmp
echo 120 -1.25 14p,Times-Roman LB 120\312 >> tmp
echo 370 -1.35 24p,Symbol LT a >> tmp
echo -5 1.85 24p,Times-Roman RT x,y >> tmp
gmt pstext -R -J -O -K -Dj0.2c -N tmp -F+f+j >> %ps%
REM Draw a circle and indicate the 0-70 degree angle
echo 0 0 | gmt psxy -R-1/1/-1/1 -Jx1.5i -O -K -X3.625i -Y2.75i -Sc2i -W1p -N >> %ps%
echo -1 0 > tmp
echo 1 0 >> tmp
gmt psxy -R -J -O -K -W tmp >> %ps%
echo 0 -1 > tmp
echo 0 1 >> tmp
gmt psxy -R -J -O -K -W tmp >> %ps%
echo 0 0 > tmp
echo 1 0 >> tmp
gmt psxy -R -J -O -K -W1p tmp >> %ps%
echo 0 0 > tmp
echo -0.5 0.866025 >> tmp
gmt psxy -R -J -O -K -W1p tmp >> %ps%
echo -0.3333 0 > tmp
echo 0 0 >> tmp
gmt psxy -R -J -O -K -W2p tmp >> %ps%
echo -0.3333 0.57735 > tmp
echo -0.3333 0 >> tmp
gmt psxy -R -J -O -K -W2p tmp >> %ps%
echo -0.16666 0 12p,5 0 CT x > tmp
echo -0.3333 0.2888675 12p,Times-Roman 0 RM y >> tmp
echo 0.22 0.27 12p,Symbol -30 CB a >> tmp
echo -0.33333 0.6 12p,Times-Roman 30 LB 120\312 >> tmp
gmt pstext -R -J -O -K -Dj0.05i tmp -F+f+a+j >> %ps%
echo 0 0 0.5i 0 120 | gmt psxy -R -J -O -Sm0.15i+e -W1p -Gblack >> %ps%
REM Clean up
del .gmt*
del gmt.conf
|