Changeset 4528 for palm/trunk/SOURCE

Timestamp:

May 11, 2020 2:14:09 PM (5 years ago)

Author:

oliver.maas

Message:

added wind_turbine_parameters-namelist parameter smearing_kernel_size

File:

• palm/trunk/SOURCE/wind_turbine_model_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/wind_turbine_model_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! added namelist parameter smearing_kernel_size
+! 
+! 4497 2020-04-15 10:20:51Z raasch
 ! file re-formatted to follow the PALM coding standard
 !
@@ -231 +234 @@
     REAL(wp) ::  segment_width_radial       = 0.5_wp  !< width of the segments, the rotor area is divided into
                                                       !< (in radial direction, as factor of MIN(dx,dy,dz))
+    
+    REAL(wp) ::  smearing_kernel_size       = 2.0_wp  !< size of the smearing kernel as multiples of dx
+
     REAL(wp) ::  time_turbine_on            = 0.0_wp  !< time at which turbines are started
     REAL(wp) ::  tilt_angle                 = 0.0_wp  !< vertical tilt_angle of the rotor [degree] ( positive = backwards )
@@ -530 +536 @@
                 nacelle_radius, rotor_radius, segment_length_tangential, segment_width_radial,     &
                 region_2_slope, speed_control, tilt_angle, time_turbine_on,                        &
-                tower_cd, pitch_rate,                                                              &
+                smearing_kernel_size, tower_cd, pitch_rate,                                        &
                 yaw_control, yaw_speed, tip_loss_correction
 !                , nacelle_cd
@@ -1068 +1074 @@
 !-- The kernel is set equal to twice the grid spacing which has turned out to be a reasonable
 !-- value (see e.g. Troldborg et al. (2013), Wind Energy, DOI: 10.1002/we.1608):
-    eps_kernel = 2.0_wp * dx
+    eps_kernel = smearing_kernel_size * dx
 !
 !-- The zero point (eps_min) of the polynomial function must be the following if the integral of