Changeset 801 for palm/trunk/SOURCE/modules.f90

Timestamp:

Jan 10, 2012 5:30:36 PM (12 years ago)

Author:

suehring

Message:

Bugfix concerning OpenMP parallelization. Calculation of turbulent fluxes in advec_ws is now thread-safe.

File:

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

palm/trunk/SOURCE/modules.f90

@@ -7 +7 @@
 ! particle arrays (particles, parrticles_temp) implemented as pointers,
 ! +particles_1, particles_2, sort_count
+!
+! Dimesion of sums_wsus_ws_l, ! sums_wsvs_ws_l, sums_us2_ws_l, sums_vs2_ws_l,
+! sums_ws2_ws_l, sums_wspts_ws_l, sums_wsqs_ws_l, sums_wssas_ws_l increased.
+! for thread-safe summation in advec_ws.
 !
 ! Former revisions:
@@ -1439 +1443 @@
     REAL ::     u_max, v_max, w_max
     REAL, DIMENSION(:), ALLOCATABLE       ::  sums_divnew_l, sums_divold_l,   &
+                                              weight_substep, weight_pres
+    REAL, DIMENSION(:,:), ALLOCATABLE     ::  sums, sums_wsts_bc_l, ts_value, &
                                               sums_wsus_ws_l, sums_wsvs_ws_l, &
                                               sums_us2_ws_l, sums_vs2_ws_l,   &
                                               sums_ws2_ws_l,                  &
                                               sums_wspts_ws_l,                &
-                                              sums_wssas_ws_l,sums_wsqs_ws_l, &                 
-                                              weight_substep, weight_pres
-    REAL, DIMENSION(:,:), ALLOCATABLE     ::  sums, sums_wsts_bc_l, ts_value
+                                              sums_wssas_ws_l,sums_wsqs_ws_l
     REAL, DIMENSION(:,:,:), ALLOCATABLE   ::  hom_sum, rmask, spectrum_x, &
                                               spectrum_y, sums_l, sums_l_l, &