Ignore:
Timestamp:
Sep 19, 2012 2:30:36 PM (9 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r979 r1007  
    44! Current revisions:
    55! -----------------
    6 !
     6! Bugfix in calculation of ql_vp
    77!
    88! Former revisions:
     
    388388
    389389          CASE ( 'ql_vp' )
    390              DO  i = nxlg, nxrg
    391                 DO  j = nysg, nyng
    392                    DO  k = nzb, nzt+1
    393                       ql_vp_av(k,j,i) = ql_vp_av(k,j,i) + ql_vp(k,j,i)
     390             DO  i = nxl, nxr
     391                DO  j = nys, nyn
     392                   DO  k = nzb, nzt+1
     393                      psi = prt_start_index(k,j,i)
     394                      DO  n = psi, psi+prt_count(k,j,i)-1
     395                         ql_vp_av(k,j,i) = ql_vp_av(k,j,i) + &
     396                                           particles(n)%weight_factor / &
     397                                           prt_count(k,j,i)
     398                      ENDDO
    394399                   ENDDO
    395400                ENDDO
Note: See TracChangeset for help on using the changeset viewer.