NASA Ocean Color

ocssw V2022
       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***********************************************************************
const
subroutine const(NGAUSS, NMAX, MMAX, P, X, W, AN, ANN, S, SS, NP, EPS)
Definition: ampld.lp.f:924
#define pi
Definition: vincenty.c:23
float ei(float x)
iniclz
subroutine iniclz
Definition: iniclz.f:2