sdfopen /grove3/iesuw/ccm1/CCM1-CNTROLPI-MONE0510.nc enable print landmask.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 = 1 over land, -oro is negative over the land and land are masked out define tmask = maskout(ts,-oro) set t 1 set cint 5 d tmask draw title TS with LAND MASKED OUT print disable print