anchor.diff:                          ! 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
  +                       !  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, &
  
positive.code:         do, p=1, nRotorPoints
          
              ! initialize the TI_area_ratio to 0.0 for each turbine
              ti_area_ratio = 0.0_dp
      
              ! initialize deficit summation term to zero
              deficit_sum = 0.0_dp
  
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
  
  