Remove workaround for borehole field initiation with two far-apart boreholes.

        for k in range(nsegments):

            la = gt.load_aggregation.ClaessonJaved(sim_params.dt, sim_params.tmax)

            LoadAgg.append(la)

        # TODO: JD: the workaround below was not reuired anymore?

        # The field contains only one borehole

        # boreField = [gt.boreholes.Borehole(custom_pipe.b.H, custom_pipe.b.D, custom_pipe.b.r_b, x=0., y=0.)]

        # The field contains only one borehole, but needs one extra at very large distance to be correct, ie. gfunc plateaus at 6.7)

        boreField = [

            gt.boreholes.Borehole(

                custom_pipe.b.H,

                custom_pipe.b.D,

                custom_pipe.b.r_b,

                x=0.0,

                y=0.0,

                tilt=np.radians(0.1),

            ),

            gt.boreholes.Borehole(

                custom_pipe.b.H,

                custom_pipe.b.D,

                custom_pipe.b.r_b,

                x=1000.0,

                y=0.0,

                tilt=np.radians(0.1),

            ),

        ]

        boreField = [gt.boreholes.Borehole(custom_pipe.b.H, custom_pipe.b.D, custom_pipe.b.r_b, x=0., y=0.)]

        # Get time values needed for g-function evaluation

        time_req = LoadAgg[0].get_times_for_simulation()