Changeset 2292 for palm/trunk/SOURCE/swap_timelevel.f90

Timestamp:

Jun 20, 2017 9:51:42 AM (7 years ago)

Author:

schwenkel

Message:

implementation of new bulk microphysics scheme

File:

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

palm/trunk/SOURCE/swap_timelevel.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implementation of new microphysic scheme: cloud_scheme = 'morrison' 
+! includes two more prognostic equations for cloud drop concentration (nc)  
+! and cloud water content (qc). 
+! 
+! 2233 2017-05-30 18:08:54Z suehring
 !
 ! 2232 2017-05-30 17:47:52Z suehring
@@ -103 +108 @@
 #if defined( __nopointer )
     USE arrays_3d,                                                             &
-        ONLY:  e, e_p, nr, nr_p, pt, pt_p, q, q_p, qr, qr_p, s, s_p, sa, sa_p, &
-               u, u_p, v, v_p, w, w_p
+        ONLY:  e, e_p, nc, nc_p, nr, nr_p, pt, pt_p, q, q_p, qc, qc_p qr, qr_p,&
+               s, s_p, sa, sa_p, u, u_p, v, v_p, w, w_p
 #else
     USE arrays_3d,                                                             &
-        ONLY:  e, e_1, e_2, e_p, nr, nr_1, nr_2, nr_p, pt, pt_1, pt_2, pt_p, q,&
-               q_1, q_2, q_p, qr, qr_1, qr_2, qr_p, s, s_1, s_2, s_p, sa, sa_1,&
-               sa_2, sa_p, u, u_1, u_2, u_p, v, v_1, v_2, v_p, w, w_1, w_2, w_p
+        ONLY:  e, e_1, e_2, e_p, nc, nc_1, nc_2, nc_p, nr, nr_1, nr_2, nr_p,   &
+               pt, pt_1, pt_2, pt_p, q, q_1, q_2, q_p, qc, qc_1, qc_2, qc_p,   &
+               qr, qr_1, qr_2, qr_p, s, s_1, s_2, s_p, sa, sa_1, sa_2, sa_p,   &
+               u, u_1, u_2, u_p, v, v_1, v_2, v_p, w, w_1, w_2, w_p
 
 #endif
@@ -121 +127 @@
     USE control_parameters,                                                    &
         ONLY:  cloud_physics, constant_diffusion, humidity, land_surface,      &
-               microphysics_seifert, neutral, ocean, passive_scalar,           &
-               timestep_count, urban_surface
+               microphysics_morrison, microphysics_seifert, neutral, ocean,    &
+               passive_scalar, timestep_count, urban_surface
 
     USE indices,                                                               &
@@ -174 +180 @@
 
     IF ( humidity )  THEN
-       q = q_p             
+       q = q_p
+       IF ( cloud_physics  .AND.  microphysics_morrison )  THEN
+          qc = qc_p
+          nc = nc_p
+       ENDIF             
        IF ( cloud_physics  .AND.  microphysics_seifert )  THEN
           qr = qr_p
@@ -214 +224 @@
           IF ( humidity )  THEN
              q => q_1;    q_p => q_2
+             IF ( cloud_physics  .AND.  microphysics_morrison )  THEN
+                qc => qc_1;    qc_p => qc_2
+                nc => nc_1;    nc_p => nc_2
+             ENDIF
              IF ( cloud_physics  .AND.  microphysics_seifert )  THEN
                 qr => qr_1;    qr_p => qr_2
@@ -241 +255 @@
           IF ( humidity )  THEN
              q => q_2;    q_p => q_1
+             IF ( cloud_physics  .AND.  microphysics_morrison )  THEN
+                qc => qc_2;    qc_p => qc_1
+                nc => nc_2;    nc_p => nc_1
+             ENDIF
              IF ( cloud_physics  .AND.  microphysics_seifert )  THEN
                 qr => qr_2;    qr_p => qr_1