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