Hello, There are some minor issues here.

[Timemakesmegrow]

Can this work achieve mature code for a 4f system, such as determining the positions of each component?
Can this work solely use CZT in segmented propagation?

[Timemakesmegrow]

When there is a bug in the original code for generating field images, resulting in two colorbars in the phase that cause inconsistent image sizes.
This is the revised code.
def draw_field(self,
logarithm: float = 0.,
normalize: str = 'maximum',
title: str = "",
cut_value: float = 0.,
colormap_kind: str = '',
percentage_intensity: float = 0.,
**kwargs):

    amplitude, intensity, phase = field_parameters(self.u, has_amplitude_sign=True)
    intensity = reduce_matrix_size(self.reduce_matrix, self.x, self.y, intensity)
    phase     = reduce_matrix_size(self.reduce_matrix, self.x, self.y, phase)

    if percentage_intensity is None:
        percentage_intensity = percentage_intensity_config

    phase[intensity < percentage_intensity * intensity.max()] = 0

    intensity = normalize_draw(intensity, logarithm, normalize, cut_value)
    phase_deg = phase / degrees

    extension = [self.x[0], self.x[-1], self.y[0], self.y[-1]]

    fig = plt.figure()
    fig.suptitle(title)

    # ----- intensity subplot -----
    ax1 = fig.add_subplot(1, 2, 1)
    h1 = ax1.imshow(intensity, interpolation='bilinear',
                    aspect='auto', origin='lower', extent=extension)
    h1.set_cmap(self.CONF_DRAWING["color_intensity"])
    ax1.set_xlabel(r"$x  (\mu m)$")
    ax1.set_ylabel(r"$y  (\mu m)$")
    ax1.set_title("intensity")
    ax1.axis('scaled')
    ax1.axis(extension)
    fig.colorbar(h1, ax=ax1, orientation='horizontal', shrink=0.66, pad=0.12)

    # ----- phase subplot -----
    ax2 = fig.add_subplot(1, 2, 2)
    h2 = ax2.imshow(phase_deg, interpolation='bilinear',
                    aspect='auto', origin='lower', extent=extension)
    h2.set_cmap(self.CONF_DRAWING["color_phase"])
    h2.set_clim(-180, 180)
    ax2.set_xlabel(r"$x  (\mu m)$")
    ax2.set_ylabel(r"$y  (\mu m)$")
    ax2.set_title("phase")
    ax2.axis('scaled')
    ax2.axis(extension)
    fig.colorbar(h2, ax=ax2, orientation='horizontal', shrink=0.66, pad=0.12)

    plt.subplots_adjust(0.05, 0.18, 0.98, 0.90, 0.25, 0.25)
    return (h1, h2)