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

Timestamp:

Oct 14, 2020 3:11:02 PM (4 years ago)

Author:

schwenkel

Message:

Implement snow and graupel (bulk microphysics)

File:

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

palm/trunk/SOURCE/flow_statistics.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! Implement snow and graupel (bulk microphysics)
+!
+! 4703 2020-09-28 09:21:45Z suehring
 ! Revise averaging of land-surface quantities
 !
@@ -103 +106 @@
     USE arrays_3d,                                                                                 &
         ONLY:  ddzu, ddzw, d_exner, e, heatflux_output_conversion, hyp, km, kh,                    &
-               momentumflux_output_conversion, nc, ni, nr, p, prho, prr, pt, q, qc, qi, ql, qr,    &
+               momentumflux_output_conversion,                                                     &
+               nc, ni, ng, nr, ns, p, prho, prr, pt, q, qc, qi, qg, ql, qr, qs,                    &
                rho_air, rho_air_zw, rho_ocean, s, sa, u, ug, v, vg, vpt, w, w_subs,                &
                waterflux_output_conversion, zw
@@ -111 +115 @@
 
     USE bulk_cloud_model_mod,                                                                      &
-        ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert, microphysics_ice_phase
+        ONLY: bulk_cloud_model, graupel, snow, microphysics_morrison, microphysics_seifert,        &
+              microphysics_ice_phase
 
     USE chem_modules,                                                                              &
@@ -1354 +1359 @@
                             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
+                            IF ( graupel  .AND.  snow )  THEN
+                               sums_l(k,126,tn) = sums_l(k,126,tn) + ng(k,j,i) * rmask(j,i,sr) * flag
+                               sums_l(k,127,tn) = sums_l(k,127,tn) + qg(k,j,i) * rmask(j,i,sr) * flag
+                               sums_l(k,128,tn) = sums_l(k,128,tn) + ns(k,j,i) * rmask(j,i,sr) * flag
+                               sums_l(k,129,tn) = sums_l(k,129,tn) + qs(k,j,i) * rmask(j,i,sr) * flag
+                            ENDIF
                          ENDIF
 
@@ -1821 +1832 @@
              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:125)       = sums(k,123:125) * ngp_2dh_s_inner(k,sr)  / ngp_2dh(sr)
+             sums(k,123:129)       = sums(k,123:129) * ngp_2dh_s_inner(k,sr)  / ngp_2dh(sr)
           ENDIF
        ENDDO
@@ -1931 +1942 @@
        hom(:,1,124,sr) = sums(:,124)   ! ni
        hom(:,1,125,sr) = sums(:,125)   ! qi
+       hom(:,1,126,sr) = sums(:,126)   ! ng
+       hom(:,1,127,sr) = sums(:,127)   ! qg
+       hom(:,1,128,sr) = sums(:,128)   ! ns
+       hom(:,1,129,sr) = sums(:,129)   ! qs
        hom(:,1,73,sr) = sums(:,73)     ! nr
        hom(:,1,74,sr) = sums(:,74)     ! qr