anchor.diff:           ! A NEGATIVE EXPONENTIAL DlSTRIBUTION WlTH DENSITY PROPORTIONAL
            ! TO EXP(-random_exponential), USING INVERSION.
   
  -         REAL  :: fn_val
  -         REAL  :: Lambda
  +         REAL(ReKi)  :: fn_val
  +         REAL(ReKi)  :: Lambda
   
            !     Local variable
  -         REAL  :: r
  +         REAL(ReKi)  :: r
   
            DO
              CALL RANDOM_NUMBER(r)
  
positive.code: 
           IMPLICIT NONE
  
           REAL :: mean
           REAL :: var
           REAL :: fn_val
  
           !Local variables
           REAL :: k
           REAL :: Lambda
           REAL :: u
  
           k = wbpar (mean, var)
           lambda = mean / NWTC_gamma(1+1/k)
  
negative.code:          ! A NEGATIVE EXPONENTIAL DlSTRIBUTION WlTH DENSITY PROPORTIONAL
           ! TO EXP(-random_exponential), USING INVERSION.
  
           REAL(ReKi)  :: fn_val
           REAL(ReKi)  :: Sigma
  
           !     Local variable
           REAL(ReKi)  :: r
  
           DO
             CALL RANDOM_NUMBER(r)
  