Changeset 2215 for palm

Timestamp:

Apr 28, 2017 12:49:49 PM (8 years ago)

Author:

suehring

Message:

Normalization of spectra output

File:

• palm/trunk/SOURCE/spectra_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/spectra_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Normalization of spectra output adjusted
 !
 ! Former revisions:
@@ -559 +559 @@
 !
 !--    Subtract horizontal mean from the array, for which spectra have to be
-!--    calculated
+!--    calculated. Moreover, calculate variance of the respective quantitiy, 
+!--    later used for normalizing spectra output. 
        var_d_l(:) = 0.0_wp
        DO  i = nxl, nxr
@@ -698 +699 @@
 !
 !--       Normalize spectra by variance
-          sum_spec_dum = SUM( sums_spectra(:,n) )
+          sum_spec_dum = SUM( sums_spectra(1:nx/2,n) )
+
           IF ( sum_spec_dum /= 0.0_wp )  THEN
-             sums_spectra(:,n) = sums_spectra(:,n) * var_d(k) / sum_spec_dum
+             sums_spectra(1:nx/2,n) = sums_spectra(1:nx/2,n) *                 &
+                                      var_d(k) / sum_spec_dum
           ENDIF
           n = n + 1
@@ -839 +842 @@
 !
 !--       Normalize spectra by variance
-          sum_spec_dum = SUM( sums_spectra(:,n) )
-          IF ( SUM(sums_spectra(:,n)) /= 0.0_wp )  THEN
-             sums_spectra(:,n) = sums_spectra(:,n) *                           &
-                                 var_d(k) / SUM(sums_spectra(:,n))
+          sum_spec_dum = SUM( sums_spectra(1:ny/2,n) )
+          IF ( sum_spec_dum /= 0.0_wp )  THEN
+             sums_spectra(1:ny/2,n) = sums_spectra(1:ny/2,n) *                 &
+                                      var_d(k) / sum_spec_dum
           ENDIF
           n = n + 1