Changeset 3832 for palm/trunk/SOURCE/wind_turbine_model_mod.f90

Timestamp:

Mar 28, 2019 1:16:58 PM (5 years ago)

Author:

raasch

Message:

some routines instrumented with openmp directives, loop reordering for performance optimization

File:

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

palm/trunk/SOURCE/wind_turbine_model_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! instrumented with openmp directives
+! 
+! 3725 2019-02-07 10:11:02Z raasch
 ! unused variables removed
 ! 
@@ -340 +343 @@
     REAL(wp) ::  w_int_l                     !< 
     REAL(wp) ::  w_int_u                     !<
+    !$OMP THREADPRIVATE (cur_r, phi_rotor, pre_factor, torque_seg, u_int_l, u_int_u, u_rot, &
+    !$OMP&               v_int_l, v_int_u, w_int_l, w_int_u)
 !
 !-  Tendencies from the nacelle and tower thrust
@@ -346 +351 @@
     REAL(wp) ::  tend_nac_y = 0.0_wp  !< 
     REAL(wp) ::  tend_tow_y = 0.0_wp  !<
+    !$OMP THREADPRIVATE (tend_nac_x, tend_tow_x, tend_nac_y, tend_tow_y)
 
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  alpha_attack !< 
@@ -1731 +1737 @@
 !
 !--          Loop over rings of each turbine:
+
+             !$OMP PARALLEL PRIVATE (ring, rseg, thrust_seg, torque_seg_y, torque_seg_z, sin_rot,  &
+             !$OMP&                  cos_rot, re, rbx, rby, rbz, ii, jj, kk, aa, bb, cc, dd, gg)
+             !$OMP DO
              DO ring = 1, nrings(inot)
 
@@ -1919 +1929 @@
                 ENDDO
              ENDDO
+             !$OMP END PARALLEL
 
           ENDDO
@@ -1960 +1971 @@
 !
 !--          Loop over rings of each turbine:
+             !$OMP PARALLEL PRIVATE (ring, rseg, sin_rot, cos_rot, re, rea, ren, rote, rota, rotn, &
+             !$OMP&                  vtheta, phi_rel, lct, rad_d, alpha_attack, vrel,              &
+             !$OMP&                  chord, iialpha, iir, turb_cl, tl_factor, thrust_seg,          &
+             !$OMP&                  torque_seg_y, turb_cd, torque_seg_z, thrust_ring,             &
+             !$OMP&                  torque_ring_y, torque_ring_z)
+             !$OMP DO
              DO ring = 1, nrings(inot)
 !
@@ -2184 +2201 @@
 !
 !--                Add the segment thrust to the thrust of the whole rotor
+                   !$OMP CRITICAL
                    thrust_rotor(inot) = thrust_rotor(inot) +                   &
                                         thrust_seg(rseg)                   
@@ -2189 +2207 @@
 
                    torque_total(inot) = torque_total(inot) + (torque_seg * cur_r)
+                   !$OMP END CRITICAL
 
                 ENDDO   !-- end of loop over ring segments
@@ -2201 +2220 @@
 
              ENDDO   !-- end of loop over rings
+             !$OMP END PARALLEL
 
 
@@ -2297 +2317 @@
                                         eps_min ) / dy )
 
+             !$OMP PARALLEL PRIVATE (i, j, k, ring, rseg, flag, dist_u_3d, dist_v_3d, dist_w_3d)
+             !$OMP DO
              DO i = MAX( nxl, i_hub(inot) - i_smear(inot) ),                   &
                     MIN( nxr, i_hub(inot) + i_smear(inot) )
@@ -2377 +2399 @@
                 ENDDO           ! End of loop over j
              ENDDO              ! End of loop over i
+             !$OMP END PARALLEL
 
              CALL cpu_log( log_point_s(63), 'wtm_smearing', 'stop' )