Changeset 3274 for palm/trunk/SOURCE/header.f90

Timestamp:

Sep 24, 2018 3:42:55 PM (6 years ago)

Author:

knoop

Message:

Modularization of all bulk cloud physics code components

File:

• palm/trunk/SOURCE/header.f90 (modified) (9 diffs)

palm/trunk/SOURCE/header.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Modularization of all bulk cloud physics code components
+! 
+! 3241 2018-09-12 15:02:00Z raasch
 ! unused variables removed
 ! 
@@ -378 +381 @@
     USE arrays_3d,                                                             &
         ONLY:  pt_init, q_init, s_init, sa_init, ug, vg, w_subs, zu, zw
-        
+
+    USE basic_constants_and_equations_mod,                                     &
+        ONLY:  g, kappa, l_v, r_d
+
+    USE bulk_cloud_model_mod,                                                  &
+        ONLY:  bulk_cloud_model, bcm_header
+
     USE control_parameters
-        
-    USE cloud_parameters,                                                      &
-        ONLY:  cp, l_v, r_d
 
     USE cpulog,                                                                &
@@ -413 +419 @@
     USE lsf_nudging_mod,                                                       &
         ONLY:  lsf_nudging_header
-
-    USE microphysics_mod,                                                      &
-        ONLY:  cloud_water_sedimentation, collision_turbulence,                &
-               c_sedimentation, limiter_sedimentation, nc_const,               &
-               ventilation_effect
 
     USE model_1d_mod,                                                          &
@@ -742 +743 @@
     IF ( neutral )  WRITE ( io, 131 )  pt_surface
     IF ( humidity )  THEN
-       IF ( .NOT. cloud_physics )  THEN
+       IF ( .NOT. bulk_cloud_model )  THEN
           WRITE ( io, 129 )
        ELSE
@@ -1902 +1903 @@
 
 !
-!-- Cloud physics parameters
-    IF ( cloud_physics )  THEN
-       WRITE ( io, 415 )
-       WRITE ( io, 416 ) surface_pressure, r_d, rho_surface, cp, l_v
-       IF ( microphysics_seifert )  THEN
-          WRITE ( io, 510 ) 1.0E-6_wp * nc_const
-          WRITE ( io, 511 ) c_sedimentation
-       ENDIF
-    ENDIF
+!-- Output of microphysics information
+    IF ( bulk_cloud_model )  CALL bcm_header( io )
 
 !
 !-- Cloud physcis parameters / quantities / numerical methods
     WRITE ( io, 430 )
-    IF ( humidity .AND. .NOT. cloud_physics .AND. .NOT. cloud_droplets)  THEN
+    IF ( humidity .AND. .NOT. bulk_cloud_model .AND. .NOT. cloud_droplets)  THEN
        WRITE ( io, 431 )
-    ELSEIF ( humidity  .AND.  cloud_physics )  THEN
-       WRITE ( io, 432 )
-       IF ( cloud_top_radiation )  WRITE ( io, 132 )
-       IF ( microphysics_kessler )  THEN
-          WRITE ( io, 133 )
-       ELSEIF ( microphysics_seifert )  THEN
-          IF ( cloud_water_sedimentation )  WRITE ( io, 506 )
-          WRITE ( io, 505 )
-          IF ( collision_turbulence )  WRITE ( io, 507 )
-          IF ( ventilation_effect )  WRITE ( io, 508 )
-          IF ( limiter_sedimentation )  WRITE ( io, 509 )
-       ENDIF
-    ELSEIF ( humidity  .AND.  cloud_droplets )  THEN
+    ENDIF
+    IF ( cloud_top_radiation )  WRITE ( io, 132 )
+    IF ( humidity  .AND.  cloud_droplets )  THEN
        WRITE ( io, 433 )
        IF ( curvature_solution_effects )  WRITE ( io, 434 )
@@ -2105 +2089 @@
 132 FORMAT ('     Parameterization of long-wave radiation processes via'/ &
             '     effective emissivity scheme')
-133 FORMAT ('     Precipitation parameterization via Kessler-Scheme')
 134 FORMAT (' --> Additional prognostic equation for a passive scalar')
 135 FORMAT (' --> Solve perturbation pressure via ',A,' method (', &
@@ -2327 +2310 @@
 413 FORMAT ('       Reference density in buoyancy terms: ',F8.3,' kg/m**3')
 414 FORMAT ('       Reference temperature in buoyancy terms: ',F8.4,' K')
-415 FORMAT (/' Cloud physics parameters:'/ &
-             ' ------------------------'/)
-416 FORMAT ('    Surface pressure   :   p_0   = ',F7.2,' hPa'/      &
-            '    Gas constant       :   R     = ',F5.1,' J/(kg K)'/ &
-            '    Density of air     :   rho_0 =',F6.3,' kg/m**3'/  &
-            '    Specific heat cap. :   c_p   = ',F6.1,' J/(kg K)'/ &
-            '    Vapourization heat :   L_v   =',E9.2,' J/kg')
 420 FORMAT (/'    Characteristic levels of the initial temperature profile:'// &
             '       Height:        ',A,'  m'/ &
@@ -2376 +2352 @@
               ' ----------------------------------'/)
 431 FORMAT ('    Humidity is considered, bu no condensation')
-432 FORMAT ('    Bulk scheme with liquid water potential temperature and'/ &
-            '    total water content is used.'/ &
-            '    Condensation is parameterized via 0% - or 100% scheme.')
 433 FORMAT ('    Cloud droplets treated explicitly using the Lagrangian part', &
                  'icle model')
@@ -2473 +2446 @@
 503 FORMAT (' --> Momentum advection via Wicker-Skamarock-Scheme 5th order')
 504 FORMAT (' --> Scalar advection via Wicker-Skamarock-Scheme 5th order')
-505 FORMAT ('    Precipitation parameterization via Seifert-Beheng-Scheme')
-506 FORMAT ('    Cloud water sedimentation parameterization via Stokes law')
-507 FORMAT ('    Turbulence effects on precipitation process')
-508 FORMAT ('    Ventilation effects on evaporation of rain drops')
-509 FORMAT ('    Slope limiter used for sedimentation process')
-510 FORMAT ('    Droplet density    :   N_c   = ',F6.1,' 1/cm**3')
-511 FORMAT ('    Sedimentation Courant number:                  '/&
-            '                               C_s   =',F4.1,'        ')
 512 FORMAT (/' Date:                 ',A8,6X,'Run:       ',A34/      &
             ' Time:                 ',A8,6X,'Run-No.:   ',I2.2/     &