def land_unmasked(res='4x5', debug=False):
"""
Get land mask from GEOS-Chem LWI ( at given resolution)
"""
from funcs4GEOSC import get_land_map # Kludge, use GEOS-Chem LWI
# Create a np.ma mask
if debug:
print 'land_mask called for: ', res
m = np.ma.masked_not_equal( get_land_map(res=res),1 )
if debug:
print mask, mask.shape,
return m.mask
def land_unmasked(res='4x5', debug=False):
"""
Get land mask from GEOS-Chem LWI ( at given resolution)
"""
from funcs4GEOSC import get_land_map # Kludge, use GEOS-Chem LWI
# Create a np.ma mask
if debug:
print 'land_mask called for: ', res
m = np.ma.masked_not_equal(get_land_map(res=res), 1)
if debug:
print mask, mask.shape,
return m.mask
def ocean_unmasked(res='4x5', debug=False):
"""
Get ocean mask from GEOS-Chem LWI ( at given resolution)
NOTES:
- this currently returns mask with all except ocean masked
- NEEDS UPDATE
"""
from funcs4GEOSC import get_land_map
if debug:
print 'ocean_mask called for: ', res
# Create a mask from land/water/ice indices
m = np.ma.masked_not_equal( get_land_map(res=res),0 )
if debug:
print mask, mask.shape,
return m.mask
def ocean_unmasked(res='4x5', debug=False):
"""
Get ocean mask from GEOS-Chem LWI ( at given resolution)
NOTES:
- this currently returns mask with all except ocean masked
- NEEDS UPDATE
"""
from funcs4GEOSC import get_land_map
if debug:
print 'ocean_mask called for: ', res
# Create a mask from land/water/ice indices
m = np.ma.masked_not_equal(get_land_map(res=res), 0)
if debug:
print mask, mask.shape,
return m.mask
