anchor.diff:        call chr_fromReal(ejv(j),  roundBuf(177:200),17,3,rErr)
         if(withIons) then
           call chr_fromReal(nucm(j), roundBuf(202:225),17,3,rErr)
  -        write(fileUnit, "(i8,1x,i8,2(1x,l4),a225,2(1x,i4))") &
  -        partID(j),parentID(j),llostp(j),isPrim,roundBuf(1:225),naa(j),nzz(j)
  +        write(fileUnit, "(i8,1x,i8,2(1x,l4),a225,3(1x,i4))") &
  +        partID(j),parentID(j),llostp(j),isPrim,roundBuf(1:225),naa(j),nzz(j),nqq(j)
         else
           write(fileUnit, "(i8,1x,i8,2(1x,l4),a200)") &
             partID(j),parentID(j),llostp(j),isPrim,roundBuf(1:200)
  
positive.code:     end do
    else
      call f_open(unit=fileUnit,file=fileName,formatted=.false.,mode="w",status="replace",access="stream")
      write(fileUnit) napx,napxo,npart,iDummy ! 4x32bit
      do j=1,npart
        ! These have to be set explicitly as ifort converts logical to integer differently than gfortran and nagfor
        if(partID(j) <= napxo) then
  
negative.code:       else
          iLost = 0
        end if
        write(fileUnit) partID(j),parentID(j),iLost,iPrim      ! 4x32 bit
        write(fileUnit) xv1(j),xv2(j),yv1(j),yv2(j)            ! 4x64 bit
        write(fileUnit) sigmv(j),dpsv(j),ejfv(j),ejv(j)        ! 4x64 bit
        if(withIons) then
          write(fileUnit) nucm(j),naa(j),nzz(j),nqq(j),0_int16 ! 64 bit + 4x16 bit
        end if
      end do
    end if
  