anchor.diff:           iout = 99
            do iplon = 1, ncol
               do iband = 1,16
  +               ! CLIMLAB
                  istart = iband
                  iend   = iband
   
  
positive.code:                            totuclfl, totdclfl, fnetc, htrc, &
                             idrv, dplankbnd_dt, dtotuflux_dt, dtotuclfl_dt )
  
                 olr_sr(iplon,iband) = totuflux(nlayers)
  
  !  Transfer up and down fluxes and heating rate to output arrays.
  
negative.code:     endif
  
  
  !  CLIMLAB SECOND COLUMN LOOP, ALWAYS CALLED
  !  Calculates total fluxes
  !  Reset flags for normal calculation of total fluxes
        istart = 1
  