iniclz.f

Go to the documentation of this file.

      subroutine iniclz                                                 
c                                                                       
c     the subroutine iniclz defines the various variable limits         
c     used in the second part of the radiative transfer program         
c                                                                       
c***********************************************************************
c.....include the common blocks                                         
      implicit real*8 (a-h,o-z)
      include 'afrt_rt2.cmn'
c***********************************************************************
c                                                                       
c     write(6,*)'entering iniclz'
      pi=dacos(-1.0d0)                                                  
      conv=pi/180.0d0                                                   
c                                                                       
      dtau=deltau                                                       
c                                                                       
      nmu=2*nx-1                                                        
      nmu2=2*nmu                                                        
      nmum1=nmu-1                                                       
      nmum12=2*nmum1                                                    
      jjjj=nmum1/2                                                      
      jjj=jjjj+1                                                        
      nangl=nmum1/2                                                     
      nangl1=nangl+1                                                    
c     write(6,*)'nx,dtheta',nx,dtheta
c                                                                       
      nophi=360.0d0/dphi
      anophi=nophi                                                      
      lphi=dphi+0.001d0
      ddphi=dphi*conv                                                   
      noph2=nophi/2                                                     
      jpart=noph2+1                                                     
      nomat=jpart*32                                                    
      nomatp=nomat+31                                                   
      nommu=nomat*nmum1                                                 
      jjpart=jpart+1                                                    
      jpart2=jpart/2                                                    
c                                                                       
      nopas=0                                                           
      jpass=1                                                           
      nopass=npass1                                                     
c                                                                       
      nolyrp=nolyr+1                                                    
      laym1=nolyr-1                                                     
c                                                                       
      do 949 i=1,nmum1                                                  
  949 c(i)=1.                                                           
      do 200 ij=1,nmum12                                                
      do 200 ik=1,4                                                     
      pp(ij,ik)=0.0d0                                                   
      qq(ij,ik)=0.0d0                                                   
200   continue                                                          
      do 300 i=1,2                                                      
      ei(i)=0.50d0                                                      
      ei(i+2)=0.0d0                                                     
 300  continue                                                          
c     write(6,*)'qs,qt',qs,qt
      const=wvlth**2/(4.0*pi**2*qs)                                     
      ratio=qt/qs                                                       
      qsqt=qs/qt
      conr=3.0d0/(8.0d0*pi)
c
      calb=0.0d0
      xrw=1.334d0
      xiw=0.0d0
c
c     write(6,*)'leaving iniclz'
      return                                                            
      end                                                               
c***********************************************************************