anchor.diff:  !     The mode is chosen by using either rrtmg_sw.nomcica.f90 (to not use
   !     McICA) or rrtmg_sw.f90 (to use McICA) to interface with a GCM.
   !
  -!    1) Standard, single forward model calculation (imca = 0); this is 
  +!    1) Standard, single forward model calculation (imca = 0); this is
   !       valid only for clear sky or fully overcast clouds
  -!    2) Monte Carlo Independent Column Approximation (McICA, Pincus et al., 
  +!    2) Monte Carlo Independent Column Approximation (McICA, Pincus et al.,
   !       JC, 2003) method is applied to the forward model calculation (imca = 1)
   !       This method is valid for clear sky or partial cloud conditions.
   !
  
positive.code:       real(kind=rb) :: zdpgcp                 ! flux to heating conversion ratio
  
  ! Atmosphere
        real(kind=rb) :: pavel(nlay+1)          ! layer pressures (mb) 
        real(kind=rb) :: tavel(nlay+1)          ! layer temperatures (K)
        real(kind=rb) :: pz(0:nlay+1)           ! level (interface) pressures (hPa, mb)
        real(kind=rb) :: tz(0:nlay+1)           ! level (interface) temperatures (K)
  
negative.code: ! Atmosphere/clouds/aerosol - spcvrt,spcvmc
        real(kind=rb) :: ztauc(nlay+1,nbndsw)     ! cloud optical depth
        real(kind=rb) :: ztaucorig(nlay+1,nbndsw) ! unscaled cloud optical depth
        real(kind=rb) :: zasyc(nlay+1,nbndsw)     ! cloud asymmetry parameter
                                                  !  (first moment of phase function)
        real(kind=rb) :: zomgc(nlay+1,nbndsw)     ! cloud single scattering albedo
        real(kind=rb) :: ztaua(nlay+1,nbndsw)     ! total aerosol optical depth
        real(kind=rb) :: zasya(nlay+1,nbndsw)     ! total aerosol asymmetry parameter
        real(kind=rb) :: zomga(nlay+1,nbndsw)     ! total aerosol single scattering albedo
  
        real(kind=rb) :: zcldfmc(nlay+1,ngptsw)   ! cloud fraction [mcica]
        real(kind=rb) :: ztaucmc(nlay+1,ngptsw)   ! cloud optical depth [mcica]
        real(kind=rb) :: ztaormc(nlay+1,ngptsw)   ! unscaled cloud optical depth [mcica]
        real(kind=rb) :: zasycmc(nlay+1,ngptsw)   ! cloud asymmetry parameter [mcica]
        real(kind=rb) :: zomgcmc(nlay+1,ngptsw)   ! cloud single scattering albedo [mcica]
  
        real(kind=rb) :: zbbfu(nlay+2)          ! temporary upward shortwave flux (w/m2)
  