anchor.diff:         !
          write(ioname, '(a)') 'Stick spectrum'
          call IOstart(trim(ioname),sunit)
  -       open(unit=sunit,file=trim(output)//".stick",action='write',status='replace')
  -       !
  +       !
  +       if (stick_hitran) then 
  +         open(unit=sunit,file=trim(output)//".par",action='write',status='replace')
  +       else
  +         open(unit=sunit,file=trim(output)//".stick",action='write',status='replace')
  +       endif
          if (verbose>=2) then
             !
             if (cutoff_model=="HITRAN") then
  
positive.code:        !
         if (stick_hitran) then
           !
           write(out,"(i3,f12.6,e10.3,e10.3,f5.4,f5.3,f10.4,f4.2,f8.6)",advance="no") &
                       iso,tranfreq,abscoef,acoef,&
                       species(1)%gamma,species(2)%gamma,&
                       energyi,species(1)%n,species(1)%delta
  
negative.code:             !
              hitran_do = .false.
              stick_hitran = .true.
              proftype = 'STICK'
              specttype = 'ABSORPTION'
              !
              iE = 0
              iS = 0
  