anchor.diff:          !    TI_dst = TI_tmp
           !end if
   !         print *, "before: ", TI_dst, TI_tmp
  -        call smooth_max(TI_dst, TI_tmp, TI_dst)
  +        TI_dst_in = TI_dst
  +        call smooth_max(TI_dst_in, TI_tmp, TI_dst)
   !         print *, "after:: ", TI_dst, TI_tmp
   
       ! Niayifar and Porte Agel 2015, 2016 using max on area TI ratio
  
positive.code:     
      ! local
      real(dp) :: axial_induction_ust, beta, epsilon, sigma, wake_diameter, wake_overlap
      real(dp) :: TI_added, TI_tmp, rotor_area_dst, ti_area_ratio_tmp
      real(dp), parameter :: pi = 3.141592653589793_dp
      
      ! out  
  
negative.code: 
  ! adjust wind speed for wind shear
  subroutine wind_shear_func(point_z, u_ref, z_ref, z_0, shear_exp, adjusted_wind_speed)
  
      implicit none
      
      ! define precision to be the standard for a double precision ! on local system
  