anchor.diff:         !
      end select
      !
  -   ! estimate memory 
  -   ! memory_limit-memory_now
  +   if (Nspecies>0) then 
  +      write(out,"(10x,'Pressure = ',e18.7)") pressure
  +      write(out,"(10x,'Voigt parameters:  gamma       n         T0            P0')")
  +      do i =1,Nspecies
  +        write(out,"(21x,a,4f12.4)") trim(species(i)%name),species(i)%gamma,species(i)%n,species(i)%t0,species(i)%p0
  +      enddo 
  +   endif
      !
      if (hitran.and.partfunc<0.0) then
        write(out,"('For HITRAN partition function must be defined using PF or QSTAT keywords')")
  
positive.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
            if (Nspecies>0) then 
               write(out,"(10x,'Pressure = ',e18.7)") pressure
               write(out,"(10x,'Voigt parameters:  gamma       n         T0            P0')")
               do i =1,Nspecies
                 write(out,"(21x,a,4f12.4)") trim(species(i)%name),species(i)%gamma,species(i)%n,species(i)%t0,species(i)%p0
               enddo 
            endif
            !
         endif
         !
  
negative.code:            !
          elseif (hitran) then 
             !
             ! HITRAN format
             !
             read(tunit,"(i3,f12.6,e10.3,e10.3,10x,f10.4,12x,a55,23x,2f7.1)",end=20) imol,tranfreq,abscoef,acoef,energyi,ch_q,gf,gi
             !
             energyf = tranfreq + energyi
             !
  