anchor.diff:          wake_overlap = 0.0_dp
       end if
       
  -    !print *, "wake overlap in func: ", wake_overlap
  +    if (wake_verlap > 1.0_dp) then
  +        print *, "wake overlap in func: ", wake_overlap
  +        STOP 1
  +    end if
                                
   end subroutine overlap_area_func
   
  
positive.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
  
negative.code:         !    TI_dst = TI_tmp
          !end if
  !         print *, "before: ", TI_dst, TI_tmp
          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
  