anchor.diff:      
       ! wake combination method error
       else
  -        print *, "Invalid wake combination method. Must be one of [0,1,2,3]."
  +        !print *, "Invalid wake combination method. Must be one of [0,1,2,3]."
           stop 1
       end if                       
       
  
positive.code: 
      !print *, "here 0"
      
      do, loc=1, nSamples
           
          ! set combined velocity deficit to zero for loc
          deficit_sum = 0.0_dp
          
          do, turb=1, nTurbines ! at turbineX-locations
          
              ! set velocity deficit at loc due to turb to zero
              deltav = 0.0_dp
  
negative.code:     real(dp), dimension(nSamples), intent(out) :: wsArray
  
      intrinsic sin, cos, atan, max, sqrt, log
      
      !print *, "entering visualization code"
  
      ! bastankhah and porte agel 2016 define yaw to be positive clockwise, this is reversed
      yaw = - yawDeg*pi/180.0_dp
  