sdfopen /grove3/iesuw/ccm1/CCM1-QUPPLTPL-SEAE1520.nc
sdfopen /grove3/iesuw/ccm1/CCM1-PPLTPL-SEAE1520.nc
sdfopen /grove3/iesuw/ccm1/CCM1-QUPMNTMN-SEAE1520.nc
sdfopen /grove3/iesuw/ccm1/CCM1-PMNTMN-SEAE1520.nc

*enable output file
enable print four_panel.prt 

*turn grads grid off
set grid off

*"scaled" mproj does not preserve aspect/ratio
*plot will fill whole area
set mproj scaled

*set map line attributes: color, style, thickness
set map 1 1 6

*set xlopts, ylopts, clopts, control axis label display 
*set xlopts color thickness size
set xlopts 1 4 0.15
set ylopts 1 4 0.15
set clopts -1 3 0.14

*set map dimensions
set lat 0 60
set lon 40 180

*set time
set t 1

*set dimensions on page of first map
*upper right corner
*set vpage xmin xmax ymin ymax
set vpage 0.25 5.75 4.12 8.5 

*turn grads label off
*need to do this before every plot
set grads off

*set contour properties

*set clevs sets explicit contour levels
set clevs 0 0.5 1 2 4 8 16 32

*set cthick sets thickness of contour lines
set cthick 1

*set substitution string for x label
*if set only to %, no x axis labels are displayed
set xlab %

*draw prectmm from the first file (.1 specifies first file)
d prectmm.1

*draw plot title
draw title DJF QUPPLTPL

*second plot
*must reset vpage
set vpage 5.5 11 4.12 8.5 

set grads off
set clevs 0 0.5 1 2 4 8 16 32
set cthick 1
*turn off x and y axis labels
set xlab %
set ylab %
d prectmm.2
draw title PPLTPL
set vpage 0.25 5.75 0 4.37 
set grads off
set clevs 0 0.5 1 2 4 8 16 32
set cthick 1
d prectmm.3
draw title QUPMNTMN
set vpage 5.5 11 0 4.37 
set grads off
set clevs 0 0.5 1 2 4 8 16 32
set cthick 1
*turn off y axis labels
set ylab %
d prectmm.4
draw title PMNTMN
print
disable print