Changeset 1007 for palm/trunk/SOURCE/sum_up_3d_data.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/sum_up_3d_data.f90 (modified) (2 diffs)

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-!
+! Bugfix in calculation of ql_vp
 !
 ! Former revisions:
@@ -388 +388 @@
 
           CASE ( 'ql_vp' )
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
-                   DO  k = nzb, nzt+1
-                      ql_vp_av(k,j,i) = ql_vp_av(k,j,i) + ql_vp(k,j,i)
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   DO  k = nzb, nzt+1
+                      psi = prt_start_index(k,j,i)
+                      DO  n = psi, psi+prt_count(k,j,i)-1
+                         ql_vp_av(k,j,i) = ql_vp_av(k,j,i) + &
+                                           particles(n)%weight_factor / &
+                                           prt_count(k,j,i)
+                      ENDDO
                    ENDDO
                 ENDDO