anchor.diff:  
            !Calculate Weibull distribution parameters due to mean value and variance of the data
            IMPLICIT NONE
  -
  -         REAL :: mean
  -         REAL :: var
  -         REAL :: k 
  -         REAL :: k1, F1, dFdk
  -         REAL :: error 
  +                                                                                                                                                                                       
  +         REAL(ReKi) :: mean
  +         REAL(ReKi) :: var
  +         REAL(ReKi) :: k 
  +         REAL(ReKi) :: k1, F1, dFdk
  +         REAL(ReKi) :: error 
   
            INTEGER :: I
            
  
positive.code:          ! A NEGATIVE EXPONENTIAL DlSTRIBUTION WlTH DENSITY PROPORTIONAL
           ! TO EXP(-random_exponential), USING INVERSION.
  
           REAL  :: fn_val
           REAL  :: Lambda
  
           !     Local variable
           REAL  :: r
  
           DO
             CALL RANDOM_NUMBER(r)
  
negative.code: 
           IMPLICIT NONE
  
           REAL(ReKi) :: mean
           REAL(ReKi) :: var
           REAL(ReKi) :: fn_val
  
           !Local variables                                                                                                                                                                     
           REAL(ReKi) :: k
           REAL(ReKi) :: Lambda
           REAL(ReKi) :: u
  
           k = wbpar (mean, var)
           lambda = mean / NWTC_gamma(1+1/k)
  