anchor.diff:        end if
       end do
       ! - report geometrical factor
  -    write(lout,"(a,i0,a,e22.15)") "ELENS> Geom. norm. fact. for elens #",j, &
  +    write(lout,"(a,i0,a,1pe22.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
  
positive.code:     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, &
             ". Max is ",elens_radial_dim
        goto 30
      end if
  
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
  