Changeset 1007 for palm/trunk/SOURCE/data_output_3d.f90

Timestamp:

Sep 19, 2012 2:30:36 PM (12 years ago)

Author:

franke

Message:

Bugfixes

---

Missing calculation of mean particle weighting factor for output added. (data_output_2d, data_output_3d, data_output_mask, sum_up_3d_data)
Calculation of mean particle radius for output now considers the weighting factor. (data_output_mask)
Calculation of sugrid-scale buoyancy flux for humidity and cloud droplets corrected. (flow_statistics)
Factor in calculation of enhancement factor for collision efficencies corrected. (lpm_collision_kernels)
Calculation of buoyancy production now considers the liquid water mixing ratio in case of cloud droplets. (production_e)

Changes

---

Calculation of buoyancy flux for humidity in case of WS-scheme is now using turbulent fluxes of WS-scheme. (flow_statistics)
Calculation of the collision kernels now in SI units. (lpm_collision_kernels)

File:

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

palm/trunk/SOURCE/data_output_3d.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: missing calculation of ql_vp added
 !
 ! Former revisions:
@@ -295 +295 @@
           CASE ( 'ql_vp' )
              IF ( av == 0 )  THEN
-                to_be_resorted => ql_vp
-             ELSE
+                DO  i = nxl, nxr
+                   DO  j = nys, nyn
+                      DO  k = nzb, nz_do3d
+                         psi = prt_start_index(k,j,i)
+                         DO  n = psi, psi+prt_count(k,j,i)-1
+                            tend(k,j,i) = tend(k,j,i) + &
+                                          particles(n)%weight_factor / &
+                                          prt_count(k,j,i)
+                         ENDDO
+                      ENDDO
+                   ENDDO
+                ENDDO
+                CALL exchange_horiz( tend, nbgp )
+                DO  i = nxlg, nxrg
+                   DO  j = nysg, nyng
+                      DO  k = nzb, nz_do3d
+                         local_pf(i,j,k) = tend(k,j,i)
+                      ENDDO
+                   ENDDO
+                ENDDO
+                resorted = .TRUE.
+             ELSE
+                CALL exchange_horiz( ql_vp_av, nbgp )
                 to_be_resorted => ql_vp_av
              ENDIF