anchor.diff:          !
         endif 
         !
  -      indecies(itemp) = i
  +      indices(itemp) = i
         !
         j = 2*int(jrot(i))+1
         energy = energies(i)
  
positive.code:           !
            if (mod(npoints,2)==0) npoints = npoints + 1
            !
          case ("ABSORPTION","EMISSION")
            !
            specttype = trim(w)
            !
  
negative.code:    use  input
     !
     integer(ik) :: info,ipoint,nlevels,i,itemp,enunit,tunit,sunit,j,ilog,ib,ie,j0,ilevelf,ileveli,indexi,indexf,iline,maxitems,kitem,l
     real(rk)    :: beta,ln2,dtemp,dfreq,temp0,beta0,intband,dpwcoef,x0,tranfreq,tranfreq_i,abscoef,dfreq_,xp,xm,de,lor,b,lor2,dfreq_2,halfwidth0
     real(rk)    :: cmcoef,emcoef,energy,energyf,energyi,elow,jf,ji,acoef,j0rk,alpha
     character(len=cl) :: ioname
     character(wl) :: string_tmp
     !
     real(rk),allocatable :: freq(:),intens(:),jrot(:),pf(:,:),energies(:),Asum(:)
     integer(ik),allocatable :: gtot(:),indices(:)
     character(len=20),allocatable :: quantum_numbers(:,:)
     !
     character(len=9) :: npoints_fmt  !text variable containing formats for reads/writes
  