anchor.diff:      end if
       elens_lFull(ielens(iElem))=elens_r1(ielens(iElem))==zero
       elens_lZeroThick(ielens(iElem))=elens_r1(ielens(iElem))==elens_r2(ielens(iElem))
  -    if ( elens_lFull(ielens(iElem)).and.elens_lZeroThick(ielens(iElem)) ) then
  -      write(lerr,"(a)") "ELENS> ERROR corner case of elens being a wire: R1=R2=0."
  -      write(lerr,"(a)") "ELENS>       theta_R2 looses meaning. Try using a wire."
  -      iErr = .true.
  -      return
  -    end if
       
       if(st_debug) then
         call sixin_echoVal("name",    bez(iElem),                   "ELENS",iLine)
  
positive.code:           ! rr>R2: formula: r2/r
            frr=elens_r2(ielens(ix))/rr
          else
            ! rr=R2: formula: 1
            frr=one
          end if
          
        endif
  
negative.code:       end do
        write(lout,"(a,i0,a,1pe22.15)") "ELENS> New theta at r2 for elens #",j, &
             " named "//trim(bez(jj))//": ",elens_theta_r2(j)
        ! take into account geometrical normalisation factor only in case lens is not a wire
        if ( elens_type(j)>=2 .and. (.not.elens_lFull(j) .or. .not.elens_lZeroThick(j)) ) then
          write(lout,"(a,1pe22.15)") "ELENS>   ...considering also geom. norm. fact.: ", elens_geo_norm(j)
        end if
      end if
    end if
  
  