anchor.diff:  
           ntwin = ilapa-ifipa+1
           if (.not. (ntwin == 1 .or. ntwin == 2)) then
  -          write(myOutUnit,'(a)' ) "Error when reading header; ifipa=",ifipa, "ilapa=",ilapa,"ntwin=",ntwin
  +          write(textOutUnit,'(a)' ) "Error when reading header; ifipa=",ifipa, "ilapa=",ilapa,"ntwin=",ntwin
             stop 1
           end if
   
  
positive.code:     close(inUnit)
  
      if(itopa <= 0) then
        write(myOutUnit,'(a)')    "ERR in numSTFpairs: Corrupt itopa in first header of file '" // ifname // "'"
        write(myOutUnit,'(a,i0,a)') "itopa = ", itopa, " <= 0"
        stop 1
      end if
  
      if (modulo(itopa,2) /= 0) then
        write(myOutUnit,'(a)')    "ERR in numSTFpairs: Number of particles not divisible by 2."
        write(myOutUnit,'(a,i6)') "ERR in numSTFpairs: itopa = ", itopa
        stop 1
      end if
  
  
negative.code:     end if
  
      !Should never reach this point
      write(textOutUnit,'(a)') "Control reached an unreachable point?!?"
      stop 1
  
      !Error handling
  410 continue
      write(textOutUnit,'(a,i6)') "ERR while opening input file; iostat=",iostat
      stop 1
  
  411 continue
      write(textOutUnit,'(a,i6)') "ERR while opening output file; iostat=",iostat
      stop 1
  
  511 continue
      write(textOutUnit,'(a,i6)') "END while reading header; iostat=",iostat
      stop 1
  
  520 continue
      write(textOutUnit,'(a,i6)') "ERR while reading header; iostat=",iostat
      stop 1
  
  530 continue
  