Changeset 4502 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Apr 17, 2020 4:14:16 PM (4 years ago)

Author:

schwenkel

Message:

Implementation of ice microphysics

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (6 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implementation of ice microphysics
+! 
+! 4472 2020-03-24 12:21:00Z Giersch
 ! Calculations of the Kolmogorov lengt scale eta implemented
 !
@@ -84 +87 @@
     USE arrays_3d,                                                             &
         ONLY:  ddzu, ddzw, e, heatflux_output_conversion, hyp, km, kh,         &
-               momentumflux_output_conversion, nc, nr, p, prho, prr, pt, q,    &
-               qc, ql, qr, rho_air, rho_air_zw, rho_ocean, s,                  &
+               momentumflux_output_conversion, nc, ni, nr, p, prho, prr, pt, q,&
+               qc, qi, ql, qr, rho_air, rho_air_zw, rho_ocean, s,              &
                sa, u, ug, v, vg, vpt, w, w_subs, waterflux_output_conversion,  &
                zw, d_exner
@@ -93 +96 @@
 
     USE bulk_cloud_model_mod,                                                  &
-        ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert
+        ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert,   &
+              microphysics_ice_extension
 
     USE chem_modules,                                                          &
@@ -1414 +1418 @@
                                                                 flag
                          ENDIF
+                         IF ( microphysics_ice_extension )  THEN
+                            sums_l(k,124,tn) = sums_l(k,124,tn) + ni(k,j,i) *  &
+                                                                rmask(j,i,sr) *&
+                                                                flag
+                            sums_l(k,125,tn) = sums_l(k,125,tn) + qi(k,j,i) *  &
+                                                                rmask(j,i,sr) *&
+                                                                flag
+                         ENDIF
+
                          IF ( microphysics_seifert )  THEN
                             sums_l(k,73,tn) = sums_l(k,73,tn) + nr(k,j,i) *    &
@@ -1915 +1928 @@
              sums(k,116)           = sums(k,116)           / ngp_2dh_s_inner(k,sr)
              sums(k,118:pr_palm-2) = sums(k,118:pr_palm-2) / ngp_2dh_s_inner(k,sr)
-             sums(k,123)           = sums(k,123) * ngp_2dh_s_inner(k,sr)  / ngp_2dh(sr)
+             sums(k,123:125)       = sums(k,123:125) * ngp_2dh_s_inner(k,sr)  / ngp_2dh(sr)
           ENDIF
        ENDDO
@@ -2029 +2042 @@
        hom(:,1,72,sr) = hyp * 1E-2_wp  ! hyp in hPa
        hom(:,1,123,sr) = sums(:,123)   ! nc
+       hom(:,1,124,sr) = sums(:,124)   ! ni
+       hom(:,1,125,sr) = sums(:,125)   ! qi
        hom(:,1,73,sr) = sums(:,73)     ! nr
        hom(:,1,74,sr) = sums(:,74)     ! qr