glint_corr.f

Go to the documentation of this file.

      subroutine glint_corr
c     removes the direct beam from the buffer fio
c***********************************************************************
      implicit real*8 (a-h,o-z)
      include 'afrt_rt2.cmn'
c***********************************************************************
c
      do i=nx,nmum1
c       write(6,*)'Before'
c       call radnce(pi,conv,cosmu,the,fio,i,jpart,jphi)
        do j=1,jpart
            do k=1,4
               abcz=dexp(-tautot/dabs(cosmu(i)))
               fio(k,i,j)=fio(k,i,j)-fglint(k,i,j)*abcz
            enddo
        enddo
c       write(6,*)'After'
c       call radnce(pi,conv,cosmu,the,fio,i,jpart,jphi)
      enddo
c
      return
      end
c***********************************************************************