Distributed Aperture

A element aperture can be distributed on every step of integrations. This requires an NODE_LAYOUT.

Again the main program contains

call  put_an_s_aperture(0.00001d0,0.00002d0,0.8d0)
call Track_and_test_s_aperture

Here 0.8d0 indicates that a single aperture will be put 80 % down the magnet. The aperture is tiny: $(x,y)=(0.00001,0.00002)$. The subroutine Track_and_test_s_aperture is called. The following questions are asked.

write(6,*) "aperture => True or false "
read(5,*) aperture_flag
write(6,*) " Check_madx_aperture (magnet aperture) => True or false "
read(5,*) check_madx_aperture
write(6,*) " s_aperture_check ( Extended magnet aperture) => True or false "
read(5,*) s_aperture_check
write(6,*) " Give value of the absolute_aperture "
read(5,*) absolute_aperture

WRITE(6,*) " GIVE A 'BIG' INITIAL VALUE (ONE NUMBER) FOR X(1:6)= R (example 0.1) " 
READ(5,*) R
X0(1:6)=R
x=x0

s_aperture_check powers the distributed apertures.

N.B. The flag check_madx_aperture is necessary because distributed apertures are made of NST+1 regular element apertures.