focn_below.f

Go to the documentation of this file.

      subroutine focn_below
c
c     compute the water leaving radiance due to uniform illumination 
c    from below aflat ocean surface
c
c**********************************************************************
c
c     include the common blocks
c
      implicit real*8 (a-h,o-z)
      include 'common_all.cmn'
c
c**********************************************************************
c     compute the upwelling radiance just below the 
c     ocean surface due to uniform illumination (f=1)
c
      radu=1.0d0/pi
c     compute the radiance just above the ocean surface
      do i=1,jjjj
         ii=nmum1-i+1
         do j=1,jpart
            radxi(ii,j)=((1.0d0-rfwat(i))*radu)/(xrw**2+xiw*2)
            radxi(i,j)=0.0d0
         enddo
      enddo
c
c     initialize the buffers
      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,nolyrp
         write(55,rec=il)fio
      enddo
c
      if(iref.eq.1 .and. itrans.eq.1)then
        do i=jjj,nmum1
           do j=1,jpart
              do k=1,2
                 fio(k,i,j)=0.5d0*radxi(i,j)
                 fio(k+2,i,j)=0.0d0
              enddo
           enddo
        enddo
        write(55,rec=nolyrp)fio
      endif
c
      return
      end
c**********************************************************************