def calculateDeltas(gridTRNX, gridTRNY, gridTRNZ):
    """This function calculates the deltas for each axis
    
    Parameters
    ----------
    gridTRNX : array(I)
        Array containing x-coordinates of grid
    gridTRNY : array(J)
        Array containing y-coordinates of grid
    gridTRNZ : array(K)
        Array containing z-coordinates of grid
    
    Returns
    -------
    float
        Delta for x-axis
    float
        Delta for y-axis
    float
        Delta for z-axis
    """
    
    dx = (gridTRNX.max()-gridTRNX.min())/(gridTRNX.shape[0]-1)
    dy = (gridTRNY.max()-gridTRNY.min())/(gridTRNY.shape[0]-1)
    dz = (gridTRNZ.max()-gridTRNZ.min())/(gridTRNZ.shape[0]-1)
    return dx, dy, dz