Changeset 4742 for palm/trunk/SOURCE/init_masks.f90

Timestamp:

Oct 14, 2020 3:11:02 PM (4 years ago)

Author:

schwenkel

Message:

Implement snow and graupel (bulk microphysics)

File:

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

palm/trunk/SOURCE/init_masks.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! Implement snow and graupel (bulk microphysics)
+!
+! 4648 2020-08-25 07:52:08Z raasch
 ! file re-formatted to follow the PALM coding standard
 !
@@ -66 +69 @@
     USE bulk_cloud_model_mod,                                                                      &
         ONLY: bulk_cloud_model, microphysics_ice_phase, microphysics_morrison,                     &
-              microphysics_seifert
+              microphysics_seifert, snow, graupel
 
 
@@ -281 +284 @@
                 unit = '1/m3'
 
+             CASE ( 'ng' )
+                IF ( .NOT. bulk_cloud_model )  THEN
+                   WRITE ( message_string, * ) 'output of "', TRIM( var ),                         &
+                        '" requires bulk_cloud_model = .TRUE.'
+                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
+                 ELSEIF ( .NOT. microphysics_ice_phase  .OR.  .NOT.  graupel )  THEN
+                   message_string = 'output of "' // TRIM( var ) // '" ' //                        &
+                         'requires  microphysics_ice_phase = .TRUE.'
+                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
+                ENDIF
+                unit = '1/m3'
+
+
              CASE ( 'nr' )
                 IF ( .NOT. bulk_cloud_model )  THEN
@@ -293 +309 @@
                 unit = '1/m3'
 
+             CASE ( 'ns' )
+                IF ( .NOT. bulk_cloud_model )  THEN
+                   WRITE ( message_string, * ) 'output of "', TRIM( var ),                         &
+                        '" requires bulk_cloud_model = .TRUE.'
+                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
+                 ELSEIF ( .NOT. microphysics_seifert  .OR.  .NOT.  snow)  THEN
+                   message_string = 'output of "' // TRIM( var ) // '"' //                         &
+                         'requires cloud_scheme = seifert_beheng'
+                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
+                ENDIF
+                unit = '1/m3'
+
+
              CASE ( 'pc', 'pr' )
                 IF ( .NOT. particle_advection )  THEN
@@ -347 +376 @@
                 unit = 'kg/kg'
 
+             CASE ( 'qg' )
+                IF ( .NOT. bulk_cloud_model )  THEN
+                   message_string = 'output of "' // TRIM( var ) // '" ' //                        &
+                            'requires bulk_cloud_model = .TRUE.'
+                   CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
+                ELSEIF ( .NOT. microphysics_ice_phase  .OR.  .NOT.  graupel )  THEN
+                   message_string = 'output of "' // TRIM( var ) // '" ' //                        &
+                            'requires microphysics_ice_phase = .TRUE.'
+                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
+                ENDIF
+                unit = 'kg/kg'
+
              CASE ( 'qi' )
                 IF ( .NOT. bulk_cloud_model )  THEN
@@ -370 +411 @@
                 ENDIF
                 unit = 'kg/kg'
+
+             CASE ( 'qs' )
+                IF ( .NOT. bulk_cloud_model )  THEN
+                   message_string = 'output of "' // TRIM( var ) // '" ' //                        &
+                            'requires bulk_cloud_model = .TRUE.'
+                   CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
+                ELSEIF ( .NOT. microphysics_seifert  .OR.  snow )  THEN
+                   message_string = 'output of "' // TRIM( var ) // '" ' //                        &
+                            'requires cloud_scheme = seifert_beheng'
+                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
+                ENDIF
+                unit = 'kg/kg'
+
 
              CASE ( 'rho_sea_water' )