anchor.diff:        end if
       end do headers
   
  -    if(ntwin.eq.1) then
  +    if(ntwin == 1) then
         data1: do
           read(inUnit,end=530,err=531, iostat=iostat) &
                ia,ifipa,b,c,d,e,f,g,h,p
  
positive.code:                ia,ifipa,b,c,d,e,f,g,h,p
          end if
        end do data1
      else if (ntwin.eq.2) then
        data2: do
          read(inUnit,end=530,err=531, iostat=iostat) &
               ia,ifipa,b,c,d,e,f,g,h,p, ilapa,b,c1,d1,e1,f1,g1,h1,p1
  
negative.code: 
      integer nPairs
      integer iPair
      character(len=100) ofnameBuffer !TODO: More inteligent scaling (or at least a check) would be nice
                                      ! However the actual names are hard coded so whatever.
  
      nPairs = numSTFpairs(ifname)
  
  