anchor.diff:           iout = 99
            do iplon = 1, ncol
               do iband = 1,16
  +               ! CLIMLAB
                  istart = iband
                  iend   = iband
   
  
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:                            totuclfl, totdclfl, fnetc, htrc, &
                             idrv, dplankbnd_dt, dtotuflux_dt, dtotuclfl_dt )
  
  ! CLIMLAB write to spectrally-resolved output
                 olr_sr(iplon,iband) = totuflux(nlayers)
  
  !  Transfer up and down fluxes and heating rate to output arrays.
  