flxbtm.f

Go to the documentation of this file.

      subroutine flxbtm                                                 
c
c     subroutine flxbtm computes the downward diffused flux at the      
c     base of the atmosphere for rg=0                                   
c***********************************************************************
c.....include common blocks                                             
      implicit real*8 (a-h,o-z)
      include 'afrt_rt2.cmn'
c***********************************************************************
c     compute the downward diffused flux                                
      sum = 0.0                                                         
      do 5003 j = 1, nx-1
      psum = 0.0                                                        
      do 4008 k = 1,2                                                   
      psum=fio(k,j,1)+fio(k,j,jpart)+psum                               
      do 4008 i=2,nophi/2
 4008 psum=psum+2.0d0*fio(k,j,i)                                           
 5003 sum=sum+psum*dcmusq(j)*6.2831854d0/float(nophi)
      sumdwn=sum                                                        
      sumdpi = sumdwn*pi                                                
      factr=eo(1)+eo(2)                                                 
c     write(6,100)factr,sumdwn,sumdpi
100   format('factr,sumdwn,sumdpi',1p3e12.4)
      return                                                            
      end                                                               
c***********************************************************************