anchor.diff:    REAL(dp), DIMENSION(nbdirs) :: result2d
     INTEGER :: nd
     INTEGER :: nbdirs
  -!   PRINT*, turbine_y, turbine_z, rotor_diameter, wake_center_y, &
  -! & wake_center_z, wake_diameter, wake_overlap
  +!     print *, turbine_y, turbine_z, rotor_diameter, &
  +!                             wake_center_y, wake_center_z, wake_diameter, &
  +!                             wake_overlap
   ! distance between wake center and rotor center
     IF (wake_center_z .GT. turbine_z + tol .OR. wake_center_z .LT. &
   &     turbine_z - tol) THEN
  
positive.code:       END IF
      END DO
      ovdyd = SQRT(arg1)
    ELSE IF (wake_center_y .GT. turbine_y) THEN
      DO nd=1,nbdirs
  ! potential source of gradient issues, abs() did not cause a problem in FLORIS
        ovdydd(nd) = -turbine_yd(nd)
      END DO
      ovdyd = wake_center_y - turbine_y
    ELSE IF (turbine_y .GT. wake_center_y) THEN
      DO nd=1,nbdirs
        ovdydd(nd) = turbine_yd(nd)
      END DO
  
negative.code:   REAL(dp), DIMENSION(nbdirs), INTENT(OUT) :: wake_overlapd
  ! local
    REAL(dp), PARAMETER :: pi=3.141592653589793_dp, tol=0.000001_dp
    REAL(dp) :: ovdyd, ovr, ovrr, ovl, ovz, ovz2
    REAL(dp), DIMENSION(nbdirs) :: ovdydd, ovrd, ovrrd, ovld, ovzd, &
  & ovz2d
  ! load intrinsic functions
    INTRINSIC ACOS, SQRT
    REAL(dp) :: arg1
  