anchor.diff:        character(len=25) ch1
         dimension qwc(3),clo(3),clop(3),di0(2),dip0(2)
         dimension ta(6,6)
  -
  +      
  +      interface
  +         integer function dtostr(x,results)
  +           double precision, intent(in) :: x
  +           character*(*), intent(out) :: results
  +         end function dtostr
  +      end interface
  +      
         integer :: stat
         
         INTEGER :: cmdarg_i, cmdarg_length, cmdarg_status
  
positive.code: !  Read a fort.90 and print with correct binary/decimal conversion
  !-----------------------------------------------------------------------
        implicit none
        integer errno,dtostr
        character(len=80) sixtit,commen
        integer icode
        double precision dpscor,sigcor
  
negative.code: !----------------------------------------------------------------------
        end
        integer function dtostr(x,results)
  ! Uses the dtoa_c.c version of dtoa via the dtoaf.c interface copied from roundctl
        use, intrinsic :: iso_fortran_env, only : stdout=>output_unit
        implicit none
        double precision, intent(in) :: x
        character*(*), intent(out) :: results
        integer dtoaf 
        integer ilen,mode,ndigits,decpoint,mysign
        integer i,l,d,e
  