anchor.diff:     real(rk)    :: beta,ln2,ln22,dtemp,dfreq,temp0,beta0,intband,dpwcoef,x0,tranfreq,tranfreq_i,abscoef,dfreq_,xp,xm,de,lor,b,lor2,dfreq_2,halfwidth0,dnu_half,dxp
      real(rk)    :: cmcoef,emcoef,energy,energyf,energyi,elow,jf,ji,acoef,j0rk,bnorm,f,eta1,eta2,intens1,intens2,Va0,gammaV,a,wg,d
      real(rk)    :: sigma,alpha,gamma,y,x1,x2,voigt_,dx2,xi,L1,L2,acoef_
  -   integer(ik) :: Jmax,Jp,Jpp,Kp,Kpp,J_,ispecies
  +   integer(ik) :: Jmax,Jp,Jpp,Kp,Kpp,J_,ispecies,alloc_p
      real(rk)    :: gamma_,n_
      character(len=cl) :: ioname
      character(wl) :: string_tmp
  
positive.code:               !
                write(out,my_fmt,advance="no"), &
                tranfreq,abscoef,jrot(ilevelf),energyf, jrot(ileveli),energyi 
                !
                !write(out,"(a4)",advance="no"), " <- "
                !
  
negative.code:           write(out,"(10x,'Partition function = ',f17.4)") partfunc
            write(out,"(10x,'Spectrum type = ',a/)") trim(specttype)
            if (offset<0) write(out,"(10x,'Line truncated at ',f9.2/)") offset
            !
         endif
         !
  