anchor.diff:        integer, intent(in) :: i,n
         real(fp), intent(in) :: h0, h1
         real(fp) :: t,p
  -      real(fp), parameter :: k = 1.5
  +      real(fp), parameter :: k = 1.0
         t = real(i - 1, fp) / (n - 1)
         p = 1 - exp(-(h1 - h0) / (k * par % hscale))
         h = h0 - k * (par % hscale) * log(1 - t * p)
  
positive.code:         end do calc_hray
  
          Q(:) = (h_ray - hterr) / (adist * radearth * t)
          attn(j) = th(minval(Q) / real(0.01,fp) + 0.5)
        end do along_heights
  
      end block trace_ray
  
    contains
  
negative.code:       do i = 1,size(sky,1)
          call arr2geo(gt, real(i,fp), real(j,fp), lat, lng)
          call onepoint(par, lat, lng, src, maph, gth, sky(i,j,:), hobs(i,j))
          write (*,'(2F10.4,F8.1,4F8.3)') lat, lng, hobs(i,j), ity2mag(sky(i,j,:))
        end do
      end do
      !$omp end parallel do
  