anchor.diff:      deltavr = deltav0m
   
       ! linearized gaussian magnitude term for near wake
  -    deltav = (((deltav0m - deltavr)/x0) * x + deltavr) *       &
  +    deltav = (((deltav0m - deltavr)/discontinuity_point) * x + deltavr) *       &
           exp(-0.5_dp * (deltay / sigmay) ** 2) *                                 &
           exp(-0.5_dp * (deltaz / sigmaz) ** 2)
                   
  
positive.code:     ! adjusted wind speed for wind shear
      adjusted_wind_speed = u_ref*((point_z-z_0)/(z_ref-z_0))**shear_exp
      
  end subroutine wind_shear_func
negative.code: ! calculates the velocity difference between hub velocity and free stream for a given wake
  ! for use in the near wake region only
  subroutine deltav_near_wake_lin_func(deltay, deltaz, Ct, yaw,  &
                                   & sigmay, sigmaz, rotor_diameter_ust, x, &
                                   & discontinuity_point, sigmay0, sigmaz0, deltav) 
                         
      implicit none
  
  