#!/bin/csh

### Example B5: A geological cross-section

### Step 1: Remove output-to-be

/bin/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