def pressure(x, y, z, tt):
    sum_pres = 0
    t_ind = find_time(tt)
    nx = len(dx)
    ny = len(dy)
    nz = len(dz)
    nt = len(time)
    print("dx =  ", dx)
    if t_ind >= 2:
        dt = np.array(make_dt(i) for i in range(min(nt, t_ind)))  # Is it necessary to substract 1?

        for k in range(min(nt, t_ind)):  # Is it necessary to substract 1?
            for i in range(nx):
                x_frac = dx[k] * i
                for j in range(nz):
                    z_frac = dz[k] * j
                    if time[k] >= exposure_time(x_frac, z_frac):
                        sum_pres = sum_pres + 2 * 2 * (D * mu / k) * leakoff * (
                                1 / math.sqrt(time(k) - exposure_time(x_frac, z_frac) + dt[k] / 4)) * pressure_point(
                            x - x_frac, y, z - z_frac, D, tt - time[k]) * dx[k] * dz[k] * dt[k]
                    else: print ("noooo")

        if tt < nt:
            sum_pres = sum_pres + 2 * coeff * leakoff * (
                    1 / math.sqrt(time[len(time) - 1] - exposure_time(x_frac, x_frac) + 0.01)) * pressure_point(
                x - x_frac, y,
                z - z_frac, D,
                tt - time[t_ind - 1]) * dx[k] * dz[k] * (tt - time[t_ind - 1])

    return sum_pres + p0