anchor.diff:        nfile=190
         n=0
         open(nfile,file=fname,form='UNFORMATTED',status='OLD')
  -      
  +
  +      if (ofoutput .eqv. .TRUE.) then
  +         open(191,file=ofname,action='WRITE',status='REPLACE')
  +         ounit = 191
  +      end if
  +
   100   read(nfile,end=511,err=520) sixtit,commen,cdate,ctime,       &
        &progrm,ifipa,ilapa,itopa,icode,numl,qwc(1),qwc(2),qwc(3), clo(1), &
        &clop(1),clo(2),clop(2),clo(3),clop(3), di0(1),dip0(1),di0(2),dip0 &
  
positive.code:          endif
        endif
        
        write (*,*) 'Read header, record=', n
        ntwin=1
        if(ilapa.ne.ifipa) ntwin=2
        write (*,*) 'Header ntwin ',ntwin,ifipa,ilapa
        ! ifipa=0
        n=n+1 !Increase record number
        write(*,*) sixtit
        write(*,*) commen
        write (*,*) progrm,itopa,icode,numl
        errno=dtostr(qwc(1),ch1)
        write(*,*) ch1
        errno=dtostr(qwc(2),ch1)
        write(*,*) ch1
        errno=dtostr(qwc(3),ch1)
        write(*,*) ch1
        errno=dtostr(clo(1),ch1)
        write(*,*) ch1
        errno=dtostr(clop(1),ch1)
        write(*,*) ch1
        errno=dtostr(clo(2),ch1)
        write(*,*) ch1
        errno=dtostr(clop(2),ch1)
        write(*,*) ch1
        errno=dtostr(clo(3),ch1)
        write(*,*) ch1
        errno=dtostr(clop(3),ch1)
        write(*,*) ch1
        errno=dtostr(di0(1),ch1)
        write(*,*) ch1
        errno=dtostr(dip0(1),ch1)
        write(*,*) ch1
        errno=dtostr(di0(2),ch1)
        write(*,*) ch1
        errno=dtostr(dip0(2),ch1)
        write(*,*) ch1
        do i=1,6
          do j=1,6
            errno=dtostr(ta(i,j),ch1)
            write (*,*) ch1
          enddo
        enddo
        errno=dtostr(dmmac,ch1)
        write(*,*) ch1
        errno=dtostr(dnms,ch1)
        write(*,*) ch1
        errno=dtostr(dizu0,ch1)
        write(*,*) ch1
        errno=dtostr(dnumlr,ch1)
        write(*,*) ch1
        errno=dtostr(sigcor,ch1)
        write(*,*) ch1
        errno=dtostr(dpscor,ch1)
        write(*,*) ch1
  
        if (STF) then
           if (ifipa .lt. itopa-1) goto 100 !Read more headers
  
negative.code:                write(*,*) " first and last pair inclusive. Expects first number < last number"
                 write(*,*) "The --fname flag can be used to specify the name of the file to read from (max 100 characters)."
                 write(*,*) " If nothing is specified, it defaults to 'fort.190'."
                 write(*,*) "The --ofname flag can be used to specify the name of the file to write to (max 100 characters)."
                 write(*,*) " If nothing is specified, it defaults to stdout."
                 stop 1
              end if
           end if
  