Changeset 137 for palm/trunk/SOURCE

Timestamp:

Nov 28, 2007 8:50:10 AM (17 years ago)

Author:

letzel

Message:

Bugfix: summation of sums_l_l in diffusivities.

Location:

palm/trunk/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (2 diffs)
• diffusivities.f90 (modified) (2 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -51 +51 @@
 ------
 
+Bugfix: summation of sums_l_l in diffusivities.
+
 Several bugfixes in the ocean part: Initial density rho is calculated
 (init_ocean). Error in initializing u_init and v_init removed
@@ -67 +69 @@
 Bugfix: absolute value of f must be used when calculating the Blackadar mixing length (init_1d_model)
 
-advec_particles, check_parameters, combine_plot_fields, diffusion_s, init_ocean, init_1d_model, poisfft, production_e
+advec_particles, check_parameters, combine_plot_fields, diffusion_s, diffusivities, init_ocean, init_1d_model, poisfft, production_e

palm/trunk/SOURCE/diffusivities.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! Bugfix for summation of sums_l_l for flow_statistics
 ! Vertical scalar profiles now based on nzb_s_inner and ngp_2dh_s_inner.
 ! 
@@ -140 +141 @@
 !
 !--       Summation for averaged profile (cf. flow_statistics)
+!--       (the IF statement still requires a performance check on NEC machines)
           DO  sr = 0, statistic_regions
-             IF ( rmask(j,i,sr) /= 0.0 )  THEN
+             IF ( rmask(j,i,sr) /= 0.0 .AND.  &
+                  i >= nxl .AND. i <= nxr .AND. j >= nys .AND. j <= nyn )  THEN
                 DO  k = nzb_s_inner(j,i)+1, nzt
                    sums_l_l(k,sr,tn) = sums_l_l(k,sr,tn) + l(k)