anchor.diff:            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
          !
  
positive.code:        !
     end select
     !
     ! estimate memory 
     ! memory_limit-memory_now
     !
     if (hitran.and.partfunc<0.0) then
       write(out,"('For HITRAN partition function must be defined using PF or QSTAT keywords')")
  
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
             !
  