anchor.diff:       &     ATOL,' RTOL= ',RTOL, ' SF= ',SF, 'DN_SAFE= ',DN_SAFE
         ENDIF
         DO IPT=1,NPT
  -
            IF(IVERBOSE.GT.1) THEN
   CMWF DEBUG
               WRITE(LU_O,*)'T_I(',IPT,')= ',T_I(IPT),
  
positive.code: C ==> NO NEED TO REDUCE TIMESTEP
  C
  C PYADH HAD               IF(DABS(D2).LE.DN_SAFE) THEN
                 IF(DABS(D2*DL).LE.ZEROTOL) THEN
                    IF (IDEBUG.GT.0) THEN
  C MWF DEBUG
                       WRITE(6,*)'D2*DL.LE.ZEROTOL GOTO 150'
  
negative.code:       ELSE
           TERM4 = (SE - (1.D0 + S + 0.5*S**2 + ONESIXTH*S**3))/(VX**4)
        ENDIF
        IF (SMAG.LT.(120.D0*ZEROTOL)**(0.2D0)) THEN
           TERM5 = DT**(3.D0)*(ONESIXTH + ONE24*S)
        ELSE
           TERM5 = (SE - (1.D0 + S + 0.5*S*S))/(VX**3)
  