anchor.diff:    logical :: accepted_work = .false.
     logical :: completed_work = .true.
     logical :: all_done = .false.
  -  
  -   type(VoigtKampffCollection),save :: fast_voigt
  +  !
  +  type(VoigtKampffCollection),save :: fast_voigt
     
     
     !
  
positive.code:              !
               select case (trim(cutoff_model))
                  !
               case ("EXP")
                  !
                  continue
                  !
               case ("EXP-WEAK")
                  !
                  abscoef_ref =cmcoef*acoef*gtot(ilevelf)*exp(-beta_ref*energyi)*(1.0_rk-exp(-beta_ref*tranfreq))/&
                               (tranfreq**2*partfunc_ref)
                  !
                  cutoff = S_crit*exp(-tranfreq/nu_crit)
                  !
                  if ((abscoef>=cutoff.or.abscoef_ref>=cutoff).or.energyi<=ener_cut) then
                    cycle loop_swap
                  endif
                  !
  
negative.code:                 ilevelf = ilevelf_ram(iswap)
                  ileveli = ileveli_ram(iswap)
                  !
                  !
                  select case (trim(cutoff_model))
                     !
                  case ("EXP","EXP-WEAK","EXP-STRONG")
                    cutoff = S_crit*exp(-tranfreq/nu_crit)
                  case ("GRID")
                    inu  = min(max(1,nint(tranfreq/ ( (freqr-freql)/cutoffs%npoints ) )),cutoffs%npoints)
                    cutoff = cutoffs%intensity(inu)*cutoffs%factor
                  end select 
                  !
                  ileveli = ileveli_ram(iswap)
                  energyi = energies(ileveli)
                  !
                  if (abscoef<cutoff.and.energyi>cutoffs%enercut) then
                      write(wunit,my_fmt) ilevelf,ileveli,acoef !,tranfreq
                  else
                      write(sunit,my_fmt) ilevelf,ileveli,acoef,tranfreq
  