anchor.diff:                                                        !   (lw aerosols/scattering not yet available)
   
   ! ----- Output -----
  -      real(kind=rb), intent(out) :: olr_sr(:,:)       ! Spectrally-decomposed OLR (W/m2)
  +      real(kind=rb), intent(out) :: olr_sr(:,:)       ! CLIMLAB Spectrally-decomposed OLR (W/m2)
                                                         !    Dimensions: (ncol,nbndlw)
         real(kind=rb), intent(out) :: uflx(:,:)         ! Total sky longwave upward flux (W/m2)
                                                         !    Dimensions: (ncol,nlay+1)
  
positive.code: ! to be used.  Clear sky calculation is done simultaneously.
  ! For McICA, RTRNMC is called for clear and cloudy calculations.
  
  !  Do a loop over bands by setting istart=iend = 1,2,3,4,5,...16
  
                 call rtrnmc(nlayers, istart, iend, iout, pz, semiss, ncbands, &
                             cldfmc, taucmc, planklay, planklev, plankbnd, &
  
negative.code:     endif
  
  
  !  CLIMLAB SECOND COLUMN LOOP, ALWAYS CALLED
  !  Calculates total fluxes
  !  Reset flags for normal calculation of total fluxes
        istart = 1
  