Skip to content

cinei.utils

Calculate grid cell area based on latitude and resolution.

Parameters

lat : ndarray Latitude values in degrees res : float Resolution in degrees

Returns

ndarray Grid cell areas in square kilometers

Source code in cinei/utils.py 
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
def ll_area(lat, res):
    """
    Calculate grid cell area based on latitude and resolution.

    Parameters
    ----------
    lat : ndarray
        Latitude values in degrees
    res : float
        Resolution in degrees

    Returns
    -------
    ndarray
        Grid cell areas in square kilometers
    """
    Re = 6371.392  # Earth radius in km
    X = Re * np.cos(lat * (np.pi/180)) * (np.pi/180) * res
    Y = Re * (np.pi/180) * res
    return X * Y  # Bug fix: added return statement