anchor.diff:      end if
     end do
   
  +  ! Recompute from loaded arrays
  +  oidpsv(1:napxo) = one/(one + dpsv(1:napxo))
  +  rvv(1:napxo)    = (ejv(1:napxo)*e0f)/(e0*ejfv(1:napxo))
  +
  +  ! Aperture data
  +  aperv(1:napxo,1:2) = craperv(1:napxo,1:2)
  +
  +  ! Module data
  +  call hions_crstart
     call meta_crstart
     if(dynk_enabled) then
       call dynk_crstart
  
positive.code:   if(scatter_active) then
      call scatter_crstart
    end if
    call hions_crstart
    if(melens > 0) then
      call elens_crstart
    end if
  
negative.code: 
    if(dynk_enabled) then ! Store current settings of elements affected by DYNK
      if(st_debug) then
        write(93,"(a)") "CR_POINT> Filling DYNK sets"
        flush(93)
      end if
      do j=1,dynk_nSets_unique
  