Changeset 1115 for palm/trunk/SOURCE/header.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/header.f90 (modified) (4 diffs)

palm/trunk/SOURCE/header.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! descriptions for Seifert-Beheng-cloud-physics-scheme added
 !
 ! Former revisions:
@@ -1237 +1237 @@
        WRITE ( io, 415 )
        WRITE ( io, 416 ) surface_pressure, r_d, rho_surface, cp, l_v
+       IF ( icloud_scheme == 0 )  THEN
+          WRITE ( io, 510 ) 1.0E-6 * nc_const
+          IF ( precipitation )  WRITE ( io, 511 ) c_sedimentation
+       ENDIF
     ENDIF
 
@@ -1458 +1462 @@
     ELSEIF ( humidity  .AND.  cloud_physics )  THEN
        WRITE ( io, 432 )
-       IF ( radiation )      WRITE ( io, 132 )
-       IF ( precipitation )  WRITE ( io, 133 )
+       IF ( radiation )  WRITE ( io, 132 )
+       IF ( icloud_scheme == 1 )  THEN
+          IF ( precipitation )  WRITE ( io, 133 )
+       ELSEIF ( icloud_scheme == 0 )  THEN
+          IF ( drizzle )  WRITE ( io, 506 )
+          IF ( precipitation )  THEN
+             WRITE ( io, 505 )
+             IF ( turbulence )  WRITE ( io, 507 )
+             IF ( ventilation_effect )  WRITE ( io, 508 )
+             IF ( limiter_sedimentation )  WRITE ( io, 509 )
+          ENDIF
+       ENDIF
     ELSEIF ( humidity  .AND.  cloud_droplets )  THEN
        WRITE ( io, 433 )
@@ -2020 +2034 @@
 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 ('    Drizzle 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   = ',F3.1,'        ')
 
  END SUBROUTINE header