Changeset 1115 for palm/trunk/SOURCE/prognostic_equations.f90

Timestamp:

Mar 26, 2013 6:16:16 PM (11 years ago)

Author:

hoffmann

Message:

optimization of two-moments cloud physics

File:

• palm/trunk/SOURCE/prognostic_equations.f90 (modified) (12 diffs)

palm/trunk/SOURCE/prognostic_equations.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! optimized cloud physics: calculation of microphysical tendencies transfered 
+! to microphysics.f90; qr and nr are only calculated if precipitation is required
 !
 ! Former revisions:
@@ -229 +230 @@
     IF ( ( ws_scheme_mom .OR. ws_scheme_sca )  .AND.  &
          intermediate_timestep_count == 1 )  CALL ws_statistics
-
 !
 !-- Loop over all prognostic equations
@@ -285 +285 @@
 
              CALL user_actions( i, j, 'u-tendency' )
-
 !
 !--          Prognostic equation for u-velocity component
@@ -341 +340 @@
 
              CALL user_actions( i, j, 'v-tendency' )
-
 !
 !--          Prognostic equation for v-velocity component
@@ -422 +420 @@
              ENDIF
           ENDIF
-
-!
-!--       If required, calculate tendencies for total water content, rain water 
-!--       content, rain drop concentration and liquid temperature
-          IF ( cloud_physics  .AND.  icloud_scheme == 0 )  THEN
-
-             tend_q(:,j,i)  = 0.0
-             tend_qr(:,j,i) = 0.0
-             tend_nr(:,j,i) = 0.0
-             tend_pt(:,j,i) = 0.0
-!
-!--          Droplet size distribution (dsd) properties are needed for the
-!--          computation of selfcollection, breakup, evaporation and
-!--          sedimentation of rain
-             IF ( precipitation )  THEN
-                CALL dsd_properties( i,j )
-                CALL autoconversion( i,j )
-                CALL accretion( i,j )
-                CALL selfcollection_breakup( i,j )
-                CALL evaporation_rain( i,j )
-                CALL sedimentation_rain( i,j )
-             ENDIF
-
-             IF ( drizzle )  CALL sedimentation_cloud( i,j )
-
-          ENDIF
-
+!
+!--       If required, calculate tendencies for total water content, liquid water 
+!--       potential temperature, rain water content and rain drop concentration
+          IF ( cloud_physics  .AND.  icloud_scheme == 0 )  CALL microphysics_control( i, j )
 !
 !--       If required, compute prognostic equation for potential temperature
@@ -650 +625 @@
 !--          If required, calculate prognostic equations for rain water content 
 !--          and rain drop concentration
-             IF ( cloud_physics  .AND.  icloud_scheme == 0 )  THEN
+             IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.              &
+                  precipitation )  THEN
 !
 !--             Calculate prognostic equation for rain water content
@@ -667 +643 @@
                 ENDIF
                 CALL diffusion_s( i, j, qr, qrsws, qrswst, wall_qrflux )
-
 !
 !--             Using microphysical tendencies (autoconversion, accretion, 
@@ -673 +648 @@
                 tend(:,j,i) = tend(:,j,i) + tend_qr(:,j,i)
 
-!
-!--             If required, compute influence of large-scale subsidence/ascent
-                IF ( large_scale_subsidence )  THEN
-                   CALL subsidence( i, j, tend, qr, qr_init )
-                ENDIF
-
-!              CALL user_actions( i, j, 'qr-tendency' )
-
+                CALL user_actions( i, j, 'qr-tendency' )
 !
 !--             Prognostic equation for rain water content
                 DO  k = nzb_s_inner(j,i)+1, nzt
-                   qr_p(k,j,i) = qr(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) +    &
-                                                     tsc(3) * tqr_m(k,j,i) )     &
-                                         - tsc(5) * rdf_sc(k) *                  &
-                                        ( qr(k,j,i) - qr_init(k) )
-                   IF ( qr_p(k,j,i) < 0.0 )  qr_p(k,j,i) = 0.1 * qr(k,j,i)
+                   qr_p(k,j,i) = qr(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) +  &
+                                                       tsc(3) * tqr_m(k,j,i) ) &
+                                           - tsc(5) * rdf_sc(k) * qr(k,j,i)
+                   IF ( qr_p(k,j,i) < 0.0 )  qr_p(k,j,i) = 0.0
                 ENDDO
 !
@@ -711 +678 @@
                 IF ( timestep_scheme(1:5) == 'runge' )  THEN
                    IF ( ws_scheme_sca )  THEN
-                      CALL advec_s_ws( i, j, nr, 'nr', flux_s_nr,       & 
-                                       diss_s_nr, flux_l_nr, diss_l_nr, &
-                                       i_omp_start, tn )
+                      CALL advec_s_ws( i, j, nr, 'nr', flux_s_nr,    & 
+                                    diss_s_nr, flux_l_nr, diss_l_nr, &
+                                    i_omp_start, tn )
                    ELSE 
                       CALL advec_s_pw( i, j, nr )
@@ -721 +688 @@
                 ENDIF
                 CALL diffusion_s( i, j, nr, nrsws, nrswst, wall_nrflux )
-
+!
 !--             Using microphysical tendencies (autoconversion, accretion, 
 !--             selfcollection, breakup, evaporation; 
@@ -727 +694 @@
                 tend(:,j,i) = tend(:,j,i) + tend_nr(:,j,i)
 
-!
-!--             If required, compute influence of large-scale subsidence/ascent
-                IF ( large_scale_subsidence )  THEN
-                   CALL subsidence( i, j, tend, nr, nr_init )
-                ENDIF
-
-!                CALL user_actions( i, j, 'nr-tendency' )
-
+                CALL user_actions( i, j, 'nr-tendency' )
 !
 !--             Prognostic equation for rain drop concentration
                 DO  k = nzb_s_inner(j,i)+1, nzt
-                   nr_p(k,j,i) = nr(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) + &
-                                                  tsc(3) * tnr_m(k,j,i) )     &
-                                      - tsc(5) * rdf_sc(k) *                  &
-                                        ( nr(k,j,i) - nr_init(k) )
-                   IF ( nr_p(k,j,i) < 0.0 )  nr_p(k,j,i) = 0.1 * nr(k,j,i)
+                   nr_p(k,j,i) = nr(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) +  &
+                                                       tsc(3) * tnr_m(k,j,i) ) &
+                                           - tsc(5) * rdf_sc(k) * nr(k,j,i)
+                   IF ( nr_p(k,j,i) < 0.0 )  nr_p(k,j,i) = 0.0
                 ENDDO
 !
@@ -762 +721 @@
              ENDIF
 
-          ENDIF
-
+         ENDIF
 !
 !--       If required, compute prognostic equation for turbulent kinetic