anchor.diff:  !**FUNCTION NEWDATE : CONVERTS DATES BETWEEN TWO OF THE FOLLOWING
   !FORMATS: PRINTABLE DATE, CMC DATE-TIME STAMP, TRUE DATE
   !
  -      FUNCTION NEWDATE(DAT1,DAT2,DAT3,MODE)
  +      INTEGER FUNCTION NEWDATE(DAT1,DAT2,DAT3,MODE)
         IMPLICIT NONE
  -      INTEGER NEWDATE,DAT1,DAT2(*),DAT3,MODE
  +      INTEGER DAT1,DAT2(*),DAT3,MODE
   !
   !AUTHOR   - E. BEAUCHESNE  -  JUN 96
   !
  
positive.code:       if (y2 > y1 .or. (y1 == y2 .and. m2 > m1)) inc=-1
        do annee = y2,y1,inc
          if (bissextile(annee)) then
            ier = newdate(tdate28f,annee*10000+0228,limite,print2true)
            ier = newdate(tdate29f,annee*10000+0229,0,print2true)
            if (tdate29f <= tdate28f) print *,'Error tdate29f < tdate28f'
            if (inc > 0) then
              if ((tdate2 <= tdate28f) .and. (tdate1L >= tdate29f)) then
                ndays = ndays+inc
                tdate1L = tdate1L+addit*inc
  
negative.code:                ndays = Calendar_Adjust(tdate1,tdate2,'E',adding)
                 tdate1 = tdate1 + (ndays*24)
               endif
               result=newdate(tdate1,idate,runnum,-6) ; idate1=idate(1)
             else
               result=newdate(tdate1,idate,runnum,-1) ; idate1=idate(1)
             endif
             if (result.ne.0)  then
                print *,'after if adding,if rounding',tdate1
  