anchor.diff:          endif
           dp1=dp00
           dp0=dp00
  -        ddp1 = zero ! -Wmaybe-uninitialized
  -        do 250 ib=1,1
  +
  +! ================================================================================================================================ !
  +
             dp10=dp1
  -!-----------------------------------------------------------------------
             if(idp /= 1 .or. iation /= 1) iclo6=0
             if(iclo6 == 1 .or. iclo6 == 2) then
  -            if(ib == 1) then
  -              if(iclo6r == 0) then
  -                clo6(1)  = clo(1)
  -                clop6(1) = clop(1)
  -                clo6(2)  = clo(2)
  -                clop6(2) = clop(2)
  -                clo6(3)  = zero
  -                clop6(3) = zero
  -              else
  -                write(lout,"(a)") "MAINCR> Reading closed orbit guess from fort.33"
  -                call readFort33
  -              end if
  -              call clorb(zero)
  -              call betalf(zero,qw)
  -              call phasad(zero,qwc)
  -              sigm(1) = clo6(3)
  -              dps(1)  = clop6(3)
  -              call qmodda(3,qwc)
  -              if(ilin >= 2) then
  -                nlinoo = nlin
  -                nlin   = nlino
  -                ilinc  = 1
  -                call mydaini(2,2,6,3,6,1)
  -                nlin   = nlinoo
  -              end if
  -              dp1 = dp10+clop6(3)
  +            if(iclo6r == 0) then
  +              clo6(1)  = clo(1)
  +              clop6(1) = clop(1)
  +              clo6(2)  = clo(2)
  +              clop6(2) = clop(2)
  +              clo6(3)  = zero
  +              clop6(3) = zero
  +            else
  +              write(lout,"(a)") "MAINCR> Reading closed orbit guess from fort.33"
  +              call readFort33
               end if
  +            call clorb(zero)
  +            call betalf(zero,qw)
  +            call phasad(zero,qwc)
  +            sigm(1) = clo6(3)
  +            dps(1)  = clop6(3)
  +            call qmodda(3,qwc)
  +            if(ilin >= 2) then
  +              nlinoo = nlin
  +              nlin   = nlino
  +              ilinc  = 1
  +              call mydaini(2,2,6,3,6,1)
  +              nlin   = nlinoo
  +            end if
  +            dp1 = dp10+clop6(3)
  +
               if(iqmod6 == 1) then
                 do ncrr=1,iu
                   ix=ic(ncrr)
  
positive.code:               enddo
              endif
  
              do 190 ib1=1,napx
                ib3=ib1+(ib-1)*napx
  !--beam-beam element
                clo6v(1,ib3)=clo6(1)
                clo6v(2,ib3)=clo6(2)
                clo6v(3,ib3)=clo6(3)
                clop6v(1,ib3)=clop6(1)
                clop6v(2,ib3)=clop6(2)
                clop6v(3,ib3)=clop6(3)
                di0xs(ib3)=di0(1)
                di0zs(ib3)=di0(2)
                dip0xs(ib3)=dip0(1)
                dip0zs(ib3)=dip0(2)
                qwcs(ib3,1)=qwc(1)
                qwcs(ib3,2)=qwc(2)
                qwcs(ib3,3)=qwc(3)
  
                do i2=1,6
                  do j2=1,6
                    tas(ib3,i2,j2)=tasm(i2,j2)
                  end do
                end do
  
    190       continue
            else
              if(idp.eq.1.and.iation.eq.1) then
                ncorruo=ncorru
                ncorru=1
                call clorb(zero)
  #ifdef DEBUG
  !     call dumpbin('aclorb',1,1)
  !     call abend('after  clorb                                      ')
  #endif
                call betalf(zero,qw)
                call phasad(zero,qwc)
  #ifdef DEBUG
  !     call dumpbin('abetphas',1,1)
  !     call abend('after  abetphas                                   ')
  #endif
  !--beam-beam element
                if(nbeam.ge.1) then
                nd=3
                nd2=6
  #include "include/beamcou.f90"
                endif
                ncorru=ncorruo
                iqmodc=3
                call mydaini(2,2,6,3,6,1)
  #ifdef DEBUG
  !     call dumpbin('bmydaini',999,9999)
  !     call abend('before mydaini                                    ')
  #endif
                do i=1,2
                  qwc(i)=real(int(qwc(i)),fPrec)+wxys(i)
                enddo
                if(ilin.ge.2) then
  #ifdef DEBUG
  !     call dumpbin('bmydaini',999,9999)
  !     call abend('before mydaini                                    ')
  #endif
                  nlinoo=nlin
                  nlin=nlino
                  ilinc=1
                  call mydaini(2,2,6,3,6,1)
  #ifdef DEBUG
  !     call dumpbin('amydaini',999,9999)
  !     call abend('after  mydaini                                    ')
  #endif
                  nlin=nlinoo
                endif
              else
  
negative.code: !     call abend('after bib1                                        ')
  #endif
            do 220 ib1=1,napx
              iar=ib1
  
              do ib2=1,6
                do ib3=1,6
  