[gradsusr] Maskout/Datafile Issue

[Justin Hicks]

GrADS users,

I am using 1 degree gridded GLDAS data along with 1 degree land and ocean
masks to output variables such as precipitation and ET. GLDAS only displays
land data, but when I apply an ocean mask (for example, the South Pacific)
and use the 'maskout' function, I get data values for the ocean.

Here is the code I am using:

'sdfopen /data2/control/GLDAS_NOAH10_M.020.ctl_1Deg_RainfallRateSfc.rmp.nc'
'set t 265'
'open /data2/control/Region10_SPC_mask_MERRA2.nc'
'd
aave(maskout(rainfsfc.1*86400,lterp(mask.2(time=00Z01JAN1980),rainfsfc.1)),lon=-179.5,lon=179.5,lat=-89.5,lat=89.5)'


This code outputs a value for the southern pacific ocean despite GLDAS only
containing land data. So, this code should be outputting an undefined value
rather than actual data. I am using lterp to interpolate any missing land
grid points...but even without using lterp I have the same issue. How could
I go about getting this sea mask to output the undefined value that I
expect?

-Justin


-- 
*Justin Hicks*

*UndergraduateAtmospheric Sciences Major*

*GIS MinorUniversity of Maryland-College Park*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20170627/146902a2/attachment.html