anchor.diff:            end if
           end if
         end do
  -      write(lout,"(a,i0,a,e22.15)") "ELENS> New theta at r2 for elens #",j, &
  +      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,e22.15)") "ELENS>   ...considering also geom. norm. fact.: ", &
  +      if ( elens_type(j)>=2 ) write(lout,"(a,1pe22.15)") "ELENS>   ...considering also geom. norm. fact.: ", &
              elens_geo_norm(j)
       end if
     end if
  
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 do
      ! - report geometrical factor
      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
  