There was a bug in the AVEINT subroutine, which could result in wrong averaging of
angular variables sometimes. If you used previous distribution please correct 
aveint_new.F or aveint_period.F, correct the last line in the following part of
the code:

C
C AL 4/2/2024: For NAVINT.GT.1, XNOWP at the end of the NSTEPS(ITYPE) period
C is the average value corresponding to the entire NSTEPS(ITYPE) period; 
C for NAVINT.EQ.1 it is the current value.
C
        XNOWP = (XAVEAVE(1)+XNOWP)/NAVINT
        IF (NSTEPS(1).LE.1) THEN


to

        IF (NSTEPS(ITYPE).LE.1) THEN
