anchor.diff:     FUNCTION Mean ( Ary, AryLen )
   
   
  +   !NOTE: We should make AryLen an optional argument and use SIZE( Ary ) if it is not present.
  +
         ! This routine calculates the mean value of an array.
   
   
  
positive.code: 
        ! Global numeric-related variables.
  
     REAL(DbKi)                   :: D2R_D                                        ! Factor to convert degrees to radians in double precision
     REAL(DbKi)                   :: Inf_D                                        ! IEEE value for NaN (not-a-number) in double precision
     REAL(DbKi)                   :: Inv2Pi_D                                     ! 0.5/Pi (1/(2*Pi)) in double precision
     REAL(DbKi)                   :: NaN_D                                        ! IEEE value for Inf (infinity) in double precision
     REAL(DbKi)                   :: Pi_D                                         ! Ratio of a circle's circumference to its diameter in double precision
     REAL(DbKi)                   :: PiBy2_D                                      ! Pi/2 in double precision
     REAL(DbKi)                   :: R2D_D                                        ! Factor to convert radians to degrees in double precision
     REAL(DbKi)                   :: RPM2RPS_D                                    ! Factor to convert revolutions per minute to radians per second in double precision
     REAL(DbKi)                   :: RPS2RPM_D                                    ! Factor to convert radians per second to revolutions per minute in double precision
     REAL(DbKi)                   :: TwoByPi_D                                    ! 2/Pi in double precision
     REAL(DbKi)                   :: TwoPi_D                                      ! 2*Pi in double precision
  
  
     REAL(ReKi)                   :: D2R                                          ! Factor to convert degrees to radians
     REAL(ReKi)                   :: Inf                                          ! IEEE value for NaN (not-a-number)
     REAL(ReKi)                   :: Inv2Pi                                       ! 0.5/Pi = 1 / (2*pi)
     REAL(ReKi)                   :: NaN                                          ! IEEE value for Inf (infinity)
     REAL(ReKi)                   :: Pi                                           ! Ratio of a circle's circumference to its diameter
     REAL(ReKi)                   :: PiBy2                                        ! Pi/2
     REAL(ReKi)                   :: R2D                                          ! Factor to convert radians to degrees
     REAL(ReKi)                   :: RPM2RPS                                      ! Factor to convert revolutions per minute to radians per second
     REAL(ReKi)                   :: RPS2RPM                                      ! Factor to convert radians per second to revolutions per minute
     REAL(ReKi)                   :: TwoByPi                                      ! 2/Pi
     REAL(ReKi)                   :: TwoPi                                        ! 2*Pi
  
     INTEGER, PARAMETER           :: BITS_IN_ADDR  = C_INTPTR_T*8                 ! The number of bits in an address (32-bit or 64-bit).
  
  
  !=======================================================================
  
negative.code: !
  !    This file is part of the NWTC Subroutine Library.
  !
  !    NWTC_Num is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as
  !    published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.
  !
  !    This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty
  