Beispiel #1
0
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
Beispiel #2
0
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
Beispiel #3
0
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
Beispiel #4
0
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