Changeset 2174 for palm/trunk/SOURCE/time_integration.f90

Timestamp:

Mar 13, 2017 8:18:57 AM (8 years ago)

Author:

maronga

Message:

nesting extended for cloud physics quantities, bugfixes in nesting, change of default most_method

File:

• palm/trunk/SOURCE/time_integration.f90 (modified) (5 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Added support for nesting with cloud microphysics
 ! 
 ! Former revisions:
@@ -250 +250 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  diss, dzu, e, e_p, nr_p, prho, pt, pt_p, pt_init, q_init, q,    &
-               ql, ql_c, ql_v, ql_vp, qr_p, q_p, ref_state, rho_ocean, s, s_p, sa_p, &
-               tend, u, u_p, v, vpt, v_p, w, w_p
+        ONLY:  diss, dzu, e, e_p, nr, nr_p, prho, pt, pt_p, pt_init, q_init, q,&
+               ql, ql_c, ql_v, ql_vp, qr, qr_p, q_p, ref_state, rho_ocean,     &
+               s, s_p, sa_p, tend, u, u_p, v, vpt, v_p, w, w_p
 
     USE calc_mean_profile_mod,                                                 &
@@ -603 +603 @@
              CALL pmci_datatrans( nesting_mode )
 
-             IF ( TRIM( nesting_mode ) == 'two-way' .OR.                               &
+             IF ( TRIM( nesting_mode ) == 'two-way' .OR.                       &
                   nesting_mode == 'vertical' )  THEN
 !
@@ -611 +611 @@
                 CALL exchange_horiz( v, nbgp )
                 CALL exchange_horiz( w, nbgp )
-                IF ( .NOT. neutral            )  CALL exchange_horiz( pt, nbgp )
-                IF ( humidity                 )  CALL exchange_horiz( q, nbgp  )
-                IF ( passive_scalar           )  CALL exchange_horiz( s, nbgp  )
-                IF ( .NOT. constant_diffusion )  CALL exchange_horiz( e, nbgp  )
+                IF ( .NOT. neutral )  CALL exchange_horiz( pt, nbgp )
+
+                IF ( humidity )  THEN
+
+                   CALL exchange_horiz( q, nbgp )
+
+                   IF ( cloud_physics  .AND.  microphysics_seifert )  THEN
+!                        CALL exchange_horiz( qc, nbgp )
+                       CALL exchange_horiz( qr, nbgp )
+!                        CALL exchange_horiz( nc, nbgp )
+                       CALL exchange_horiz( nr, nbgp )
+                   ENDIF
+
+                ENDIF
+
+                IF ( passive_scalar )  CALL exchange_horiz( s, nbgp )
+                IF ( .NOT. constant_diffusion )  CALL exchange_horiz( e, nbgp )
              ENDIF
 !
@@ -686 +699 @@
              CALL calc_liquid_water_content
           ENDIF
-! 
+!
 !--       If required, compute virtual potential temperature 
           IF ( humidity )  THEN