getrad.f

Go to the documentation of this file.

      subroutine getrad
c
c     include the common statements.....................................
      implicit real*8 (a-h,o-z)
      include 'afrt_rt2.cmn'
c
c************************************************************************
c
c     get the fluxes & radiances at the top and bottom of the atmosphere
c
      if(jpass.eq.1)then
           jread=53
      else
         jread=iwrt
      endif
c
c     get the radiances at the top and bottom of the atmosphere
c
      read(jread,rec=nolyr+1)fio
      call flxbtm
      fluxd(jpass,nolyr+1)=sumdwn
      do i=1,nx-1
         do j=1,jpart
            xint=fio(1,i,j)+fio(2,i,j)
            geoidn(jpass,i,j)=xint
         enddo
      enddo
c
      read(jread,rec=1)fio
      do i=nx,nmum1
         ij=i-(nx-1)
         do j=1,jpart
            xint=fio(1,i,j)+fio(2,i,j)
            geoiup(jpass,ij,j)=xint
         enddo
      enddo
      call flxtop
      fluxu(jpass,1)=sumup
c
      return
      end
c**********************************************************************