grefl.f

Go to the documentation of this file.

      subroutine grefl
c
c     compute the stokes parameters for the upwelling radiation just
c     above the ground/ocean when the atmosphere is illuminated from below
c
c*************************************************************************
c.....includes the common blocks
      implicit real*8 (a-h,o-z)
      include 'afrt_rt2.cmn'
c
c**************************************************************************
c
c     write(6,*)'welcome to subroutine grfl'
c     initialize the buffers
c
      do i=1,nmum1
         do j=1,jpart
            do k=1,4
               fio(k,i,j)=0.0d0
            enddo
         enddo
      enddo
c
c     initialize unit 55
c
      do il=1,nolyr+1
         write(55,rec=il)fio
      enddo
c
      if(iref.eq.0)then
        do i=nx,nmum1
           do j=1,jpart
              do k=1,2
                 fio(k,i,j)=0.5d0/pi
                 fio(k+2,i,j)=0.0d0
              enddo
           enddo
        enddo
        write(55,rec=nolyr+1)fio
      endif
c
      if(iref.eq.1 .or. iref.eq.2 .or. iref.eq.3)then
        read(53,rec=nolyr+1)fio
c       write(55,rec=nolyr+1)fio
      endif
c
c
c     write(6,*)'leaving subroutine grefl'
      return
      end
c***********************************************************************