anchor.diff:                write(ans,9999) yr2,mo2,dy2,hr,mn,sc
             endif
         endif
  - 9997 format(F20.8)
  + 9997 format(1PE20.10)
    9998 format(I6,I4,1F6.1)
    9999 format(I6,I4,4F6.1)
         RETURN
  
positive.code:      2    year.eq.1582.and.month.eq.10.and.day.lt.15.0d0) then
            B = 0.0d0
        else
            A = floor(dble(yrp)/100)
            B = 2.0d0 - A + floor(A*0.25d0)
        endif
        if (yrp.lt.0) then
            C = floor(365.25d0*dble(yrp) - 0.75d0)
        else
            C = floor(365.25*dble(yrp))
        endif
        D = floor(30.6001d0*(dble(mop)+1.0d0))
        jd = B+C+D+day+1720994.5
        RETURN
        END
  
  
negative.code:               read(arg2,*) ndy
            endif
        endif
        if (opt.eq.1.or.opt.eq.3.and.long.eq.1) then
            dy1 = dy1 + hr1/2.4d1 + mn1/1.44d3 + sc1/8.64d4
            dy2 = dy2 + hr2/2.4d1 + mn2/1.44d3 + sc2/8.64d4
        endif
  