acflxlvl.f

Go to the documentation of this file.

      subroutine acflxlvl(bufx,sumgxx)
c     compute actinic flux at all the lavels in the atmosphere
c     include the common statements.....................................
      implicit real*8 (a-h,o-z)
      include 'afrt_rt2.cmn'
      real*8 bufx(4,50,46)
c
c***********************************************************************
c
      anophi=dfloat(nophi)
      sumg = 0.0d0
      do j = 1,nmum1
         psumg = 0.0d0
         do k = 1,2
            psumg=bufx(k,j,1)+bufx(k,j,jpart)+psumg
            do i=2,nophi/2
               psumg=psumg+2.0d0*bufx(k,j,i)
            enddo
         enddo
         dmuxx=dabs(dcmu(j))
         prodxx1=(dmuxx*2.0*pi)/anophi
         prodxx=psumg*prodxx1
         sumg=sumg+prodxx
c        write(*,600)j,sumg,prodxx,prodxx1,psumg,dmuxx,anophi
600      format(i3,1x,1p6e11.3)
      enddo
      sumgxx=sumg
c     write(*,*)'sumgxx',sumgxx
c
      return
      end
c***********************************************************************