anchor.diff:      ii=ii+1
       if(ii>elens_radial_dim) then
         iErr = 2
  -      write(lerr,"(a,i0,a,i0)") "ELENS> ERROR too many points in radial profile: ",ii, &
  +      write(lerr,"(2(a,i0))") "ELENS> ERROR too many points in radial profile: ",ii, &
              ". Max is ",elens_radial_dim
         goto 30
       end if
  
positive.code:       end if
      end do
      ! - report geometrical factor
      write(lout,"(a,i0,a,e22.15)") "ELENS> Geom. norm. fact. for elens #",j, &
           " named "//trim(bez(jj))//": ",elens_geo_norm(j)
      
      ! Compute elens theta at R2, if requested by user
  
negative.code:           end if
          end if
        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)
        if ( elens_type(j)>=2 ) write(lout,"(a,1pe22.15)") "ELENS>   ...considering also geom. norm. fact.: ", &
             elens_geo_norm(j)
      end if
    end if
  