绘制岩性
### Example B5: A geological cross-section
### Step 1: Remove output-to-be
rm tmp.ps
### Step 2: Plot the bottom rock layer. The depth is from -100 to -90 m
# and color this layer using a sort of yellowish color
psxy -R0/3/-100/0 -JX1/7 -P -X2 -Y3 -K \
-B/a20f10 -G230/200/40 -L -W1 <<END>tmp.ps
0 -100
3 -100
3 -90
0 -90
0 -100
END
### Step 3: Write a label on the RHS of this layer
pstext -R -JX -O -K -N <<END>>tmp.ps
3.6 -95 15 0 5 ML Sandstone
END
### Step 4: Plot the 2nd rock layer. The depth is from -90 m to -70 m
# and fill this layer using a pattern. See Page 60 of
# the "GMT Technical References and Cookbook" in Room 201.
psxy -R -JX -O -K -Gp300/29 -L -W1 <<END>>tmp.ps
0 -90
3 -90
3 -70
0 -70
0 -90
END
### Step 5: Write a label on the RHS of this layer
pstext -R -JX -O -K -N <<END>>tmp.ps
3.6 -80 15 0 5 ML Granite
END
### Step 6: Plot the 3rd rock layer. The depth is from -70 m to -30 m
# and fill this layer using a pattern.
psxy -R -JX -O -K -GP300/32 -L -W1 <<END>>tmp.ps
0 -70
3 -70
3 -30
0 -30
0 -70
END
### Step 7: Write a label on the RHS of this layer
pstext -R -JX -O -K -N <<END>>tmp.ps
3.6 -50 15 0 5 ML Volcanic ash
END
### Step 8: Plot the 4th rock layer. The depth is from -30 m to 0 m
# and fill this layer using the blue color.
psxy -R -JX -O -K -G0/0/255 -L -W1 <<END>>tmp.ps
0 -30
3 -30
3 0
0 0
0 -30
END
### Step 9: Write a label on the RHS of this layer
pstext -R -JX -O -K -N <<END>>tmp.ps
3.6 -15 15 0 5 ML Sea-water
END
### Step 10: Title and label
pstext -R -JX -O -N <<END>>tmp.ps
1.5 4 16 0 0 MC A Cross-Section
-2.0 -50 15 90 5 MC Depth (m)
END
ps2jpeg150 tmp.ps
页:
[1]