 Timestamp:
 Mar 12, 2020 9:35:23 AM (10 months ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/wind_turbine_model_mod.f90
r4457 r4459 26 26 !  27 27 ! $Id$ 28 ! avoid division by zero in tip loss correction factor calculation 29 ! 30 ! 4457 20200311 14:20:43Z raasch 28 31 ! use statement for exchange horiz added 29 32 ! … … 2410 2413 ! coefficient and not to the drag coefficient in our case 2411 2414 ! 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 2416 2423 turb_cl(rseg) = tl_factor * turb_cl(rseg) 2417 2424
Note: See TracChangeset
for help on using the changeset viewer.