anchor.diff:      
       INTEGER           ::    num_element
       INTEGER           ::    i
  -    REAL              ::    r_t ( num_element )              ! The distance from the node to the hub
  -    REAL              ::    turbine_mean_velocity            ! turbine mean velocity
  -    REAL              ::    node_radius    ( num_element )
  -    REAL              ::    element_length ( num_element )
  -    REAL              ::    TI_normalization
  -    REAL              ::    FAST_Time
  +    REAL(ReKi)              ::    r_t ( num_element )              ! The distance from the node to the hub
  +    REAL(ReKi)              ::    turbine_mean_velocity            ! turbine mean velocity
  +    REAL(ReKi)              ::    node_radius    ( num_element )
  +    REAL(ReKi)              ::    element_length ( num_element )
  +    REAL(ReKi)              ::    TI_normalization
  +    REAL(ReKi)              ::    FAST_Time
       
       ! check if the meandering simulation time is valid
       IF (p%WakePosition_1 < FAST_Time/ ( (20*p%RotorR/p%p_p_r)/(0.32*p%Uambient) ) + 1 ) THEN
  
positive.code:       END DO
      END DO
      
      TI_apprant = ((TI_apprant_accumulation / REAL(counter2))**0.5)
    
      
     !-------------------------------------------------------------------------------------------------
  
negative.code:        k_wiener                 = 2*Shear_add_du_dz * delrad**2;
         DO I = 1,m%DWDD%n_r_vector
             One_div_du_dr_DWM(j-1,I) = du_dr_DWM(j-1,I) / (du_dr_DWM(j-1,I)**2 + k_wiener)
             visc_fac(I)              = max(1.0_ReKi, (du_dr_tot(j-1,I) * ABS(One_div_du_dr_DWM(j-1,I))))
             visc(j-1,I)              = visc(j-1,I) * visc_fac(I)
         END DO
         
  