anchor.diff:          do, p=1, nRotorPoints
           
               ! initialize the TI_area_ratio to 0.0 for each turbine
  -            ti_area_ratio = 0.0_dp
  +            TI_area_ratio = 0.0_dp
       
               ! initialize deficit summation term to zero
               deficit_sum = 0.0_dp
  
positive.code:                         ! vertical spread at far wake onset
                          call sigmaz_func(kz_local, deltax0_dp, rotorDiameter(turb), sigmaz_dp)
  
                          print *, "inputs in parent: ", deltay, deltaz, Ct(turb), yaw(turb), sigmay_dp, sigmaz_dp, &
                                           & rotorDiameter(turb), x, discontinuity_point, sigmay_dp, sigmaz_dp, &
                                           & wake_model_version, kz_local, x0, &
                                           & opt_exp_fac
  
                          ! velocity deficit in the nearwake (linear model)
                          call deltav_near_wake_lin_func(deltay, deltaz, &
  
negative.code:     s = 1000.0_dp
      
      g = (log(exp(s*x) + exp(s*y)))/s
      !print *, "g1 = ", g
      
      !g = (x*exp(s*x)+y*exp(s*y))/(exp(s*x)+exp(s*y))
      !print *, "g2 = ", g
      
  end subroutine smooth_max
  
  