anchor.diff:             !
              if (stick_hitran) then 
                !
  -             write(out,"(i2,i1,f12.6,e10.3,e10.3,f5.4,f5.3,f10.4,f4.2,f8.6)",advance="no") &
  -                         hitran%molec_id,hitran%molec_iso_id,tranfreq,abscoef,acoef,&
  +             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
                !
                !
                nchars_ = 0 
                !
  -             do kitem = 2,maxitems 
  +             do kitem = 3,maxitems 
                  !
                  l = len(trim(quantum_numbers(kitem,ilevelf)))
                  !
  
positive.code:           !
            hitran_do = .true.
            !
            call read_line(eof) ; if (eof) exit
            call readu(w)
            !
            do while (trim(w)/="".and.trim(w)/="END")
              !
              select case(w)
              !
              case('WRITE')
              !
              hitran_do = .false.
              stick_hitran = .true.
              !
              case('ID')
                !
                call readi(hitran%molec_id)
                !
                if (hitran%molec_id<0) call report ("Illegal molec_id "//trim(w),.true.)
                !
              case ("ISO","ISO_ID")
                !
                call readi(hitran%molec_iso_id)
                !
              case default
                !
                call report ("Unrecognized unit name "//trim(w),.true.)
                !
              end select
              !
              call read_line(eof) ; if (eof) exit
              call readu(w)
              !
            enddo          !
          case ("SELECT","FILTER")
            !
            filter  = .true.
  
negative.code:              !
               nchars_ = 0 
               !
               do kitem = 3,maxitems 
                 !
                 !l = len(trim(quantum_numbers(kitem,ileveli)))
                 !
  