sdfopen /grove3/iesuw/ccm1/CCM1-CNTROLPI-MONE0510.nc
enable print oceanmask.prt
set grid off
set grads off
set t 1 12


*maskout(expression, mask)
*when the "mask" value is < 0, values in expression are set to the missing data value
*since oro = 0 over oceans, oro-1 is negative over the oceans and oceans are masked out

define tmask = maskout(ts,oro-1)
set t 1
set cint 5
d tmask
draw title TS with OCEAN MASKED OUT
print
disable print