[gradsusr] Masking using 'maskout' command

[Justin Hicks]

GrADS users,

I am planning on using masks created using ENVI to mask oceanic data points
so that they display as NaN. These points have been assigned the value of
zero, and land values are given the value of 1.

To use the maskout command, I know the masked points must be assigned a
value less than 0, so I assume the command would be 'd
maskout(variable,mask-0.5)'.

My mask is in a binary .img file, and I was wondering if I need to write a
control file for this, as it is not working as is.

A .hdr file shows the contents of the mask:

ENVI
description = {
  File Imported into ENVI.}
samples = 360
lines   = 180
bands   = 1
header offset = 0
file type = ENVI Standard
data type = 4
interleave = bsq
sensor type = Unknown
byte order = 0
map info = {Geographic Lat/Lon, 1.5000, 1.5000, -179.50000000, 89.50000000,
1.00
00000000e+00, 1.0000000000e+00, WGS-84, units=Degrees}
coordinate system string =
{GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["W
GS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.01745
32925199433]]}
wavelength units = Unknown


The sequence of bytes are also stored as little-endian.

Not sure if I should approach this like any other data file or not.

Any help is appreciated,

-Justin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160622/a4c5de5e/attachment.html