Changeset 4459 for palm


Ignore:
Timestamp:
Mar 12, 2020 9:35:23 AM (15 months ago)
Author:
oliver.maas
Message:

avoid division by zero in tip loss correction factor calculation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/wind_turbine_model_mod.f90

    r4457 r4459  
    2626! -----------------
    2727! $Id$
     28! avoid division by zero in tip loss correction factor calculation
     29!
     30! 4457 2020-03-11 14:20:43Z raasch
    2831! use statement for exchange horiz added
    2932!
     
    24102413!--                  coefficient and not to the drag coefficient in our case
    24112414!--                 
    2412                      tl_factor = ( 2.0 / pi ) *                                &
    2413                           ACOS( EXP( -1.0 * ( 3.0 * ( rotor_radius(inot) - cur_r ) /     &
    2414                           ( 2.0 * cur_r * abs( sin( phi_rel(rseg) ) ) ) ) ) )
    2415                          
     2415                     IF ( phi_rel(rseg) == 0.0_wp ) THEN
     2416                        tl_factor = 1.0_wp
     2417                     ELSE
     2418                        tl_factor = ( 2.0 / pi ) *                                &
     2419                             ACOS( EXP( -1.0 * ( 3.0 * ( rotor_radius(inot) - cur_r ) /     &
     2420                             ( 2.0 * cur_r * abs( sin( phi_rel(rseg) ) ) ) ) ) )
     2421                     ENDIF
     2422                     
    24162423                     turb_cl(rseg)  = tl_factor * turb_cl(rseg)                                 
    24172424                                 
Note: See TracChangeset for help on using the changeset viewer.