Changeset 4459 for palm

Timestamp:

Mar 12, 2020 9:35:23 AM (4 years ago)

Author:

oliver.maas

Message:

avoid division by zero in tip loss correction factor calculation

File:

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

palm/trunk/SOURCE/wind_turbine_model_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! avoid division by zero in tip loss correction factor calculation
+! 
+! 4457 2020-03-11 14:20:43Z raasch
 ! use statement for exchange horiz added
 ! 
@@ -2410 +2413 @@
 !--                  coefficient and not to the drag coefficient in our case 
 !--                  
-                     tl_factor = ( 2.0 / pi ) *                                &
-                          ACOS( EXP( -1.0 * ( 3.0 * ( rotor_radius(inot) - cur_r ) /     &
-                          ( 2.0 * cur_r * abs( sin( phi_rel(rseg) ) ) ) ) ) )
-                          
+                     IF ( phi_rel(rseg) == 0.0_wp ) THEN
+                        tl_factor = 1.0_wp
+                     ELSE
+                        tl_factor = ( 2.0 / pi ) *                                &
+                             ACOS( EXP( -1.0 * ( 3.0 * ( rotor_radius(inot) - cur_r ) /     &
+                             ( 2.0 * cur_r * abs( sin( phi_rel(rseg) ) ) ) ) ) )
+                     ENDIF
+                     
                      turb_cl(rseg)  = tl_factor * turb_cl(rseg)