anchor.diff:      REAL(ReKi), intent(in)       ::   spacing                       ! the distance from the downstream turbine to the upstream turbine
       REAL(ReKi), intent(in)       ::   angle                         ! the angle between the investigated turbine and the line connecting the upwind turbine and wind origin
       
  -    REAL       ::   shifted_velocity                   ! the output
  +    REAL(ReKi)       ::   shifted_velocity                   ! the output
                                                     !   the velocity at a certain point
       
  -    REAL       ::   distance                      ! the distance from the point to the meandered wake center
  -    REAL       ::   y0                            ! wake center position on y axis
  -    REAL       ::   z0                            ! wake center position on z axis
  -    REAL       ::   unit                          ! single unit length  R/ppR
  -    REAL       ::   scale_factor
  +    REAL(ReKi)       ::   distance                      ! the distance from the point to the meandered wake center
  +    REAL(ReKi)       ::   y0                            ! wake center position on y axis
  +    REAL(ReKi)       ::   z0                            ! wake center position on z axis
  +    REAL(ReKi)       ::   unit                          ! single unit length  R/ppR
  +    REAL(ReKi)       ::   scale_factor
       INTEGER    ::   p1
       INTEGER    ::   p2
       INTEGER    ::   time_position                 ! to define which plane's wake center is used
  -    REAL       ::   Yshifted
  -    REAL       ::   Zshifted
  +    REAL(ReKi)       ::   Yshifted
  +    REAL(ReKi)       ::   Zshifted
       
       
       !ALLOCATE  (Uwake(NINT( m%DWDD%ppR*Rdomain ))) ! the axis symmetrical velocity
  
positive.code: 
      REAL(DbKi)   ::  timestep                     ! upper limit = usable time + grid width /  mean wind speed              %%% will change wrt wind speed  
                                                    ! = second / 0.05    timestep >= 1
      REAL         ::  y_0                          ! wake center point
      REAL         ::  z_0
      INTEGER      ::  wake_radius                  ! b(:) in cal_mixl
      REAL         ::  filter_velocity (3)          ! only v,w components
      
      INTEGER      ::  number_counter               ! counter : how many points are in the circle
      INTEGER      ::  radius_length                ! wake radius (meters)
      INTEGER      ::  y_axis                       
      INTEGER      ::  z_axis
      REAL         ::  temp_filter_velocity (3)     ! interpolation function, has u,v,w three components
      REAL         ::  temp_wind_velocity (3)
      
      INTEGER( IntKi )                   :: ErrStat           ! Error status of the operation
      CHARACTER                          :: ErrMsg            ! Error message if ErrStat /= ErrID_None
  
negative.code:     INTEGER            ::  counter1
      INTEGER            ::  counter2
      INTEGER            ::  initial_timestep
      REAL(ReKi)               ::  y_axis_turbine
      REAL(ReKi)               ::  z_axis_turbine
      REAL(ReKi)               ::  distance                  ! the distance between one point to the meandered wake center
      REAL(ReKi)               ::  TI_downstream_node        ! the TI at a specfic point in the inbestigated cross plane
      REAL(ReKi)               ::  TI_node_temp
      REAL(ReKi)               ::  TI_node
      REAL(ReKi)               ::  TI_accumulation
      REAL(ReKi)               ::  TI_apprant_accumulation
      REAl(ReKi)               ::  TI_average                ! THE AVERAGE TI OF THE CROSS PLANE
      REAL(ReKi)               ::  TI_apprant                ! The TI due to the meadering
      REAL(ReKi)               ::  HubHt 
      REAL(ReKi)               ::  wake_center_y
      REAL(ReKi)               ::  wake_center_z
      REAL(ReKi)               ::  Rscale
      REAL(ReKi)               ::  y_coor
      REAL(ReKi)               ::  z_coor  
      REAL(ReKi)               ::  zero_spacing
      REAL(ReKi)               ::  temp1,temp2,temp3
      REAL(ReKi)               ::       c_uw
      
     !-------------------------------------------------------------------------------------------------
     ! calculate the TI at each node at the downstream turbine plane from the wake deficit calculation
  