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

Timestamp:

Apr 17, 2020 4:14:16 PM (4 years ago)

Author:

schwenkel

Message:

Implementation of ice microphysics

File:

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

palm/trunk/SOURCE/init_masks.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implementation of ice microphysics
+! 
+! 4444 2020-03-05 15:59:50Z raasch
 ! bugfix: cpp-directives for serial mode added
 ! 
@@ -58 +61 @@
 
     USE bulk_cloud_model_mod,                                                  &
-        ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert
+        ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert,   &
+              microphysics_ice_extension
 
     USE control_parameters,                                                    &
@@ -268 +272 @@
                 unit = '1/m3'
 
+             CASE ( 'ni' )
+                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_extension ) THEN
+                   message_string = 'output of "' // TRIM( var ) // '" ' //    &
+                         'requires  microphysics_ice_extension = .TRUE.'
+                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
+                ENDIF
+                unit = '1/m3'
+
              CASE ( 'nr' )
                 IF ( .NOT. bulk_cloud_model )  THEN
@@ -331 +347 @@
                         '" requires bulk_cloud_model = .TRUE.'
                    CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
+                ENDIF
+                unit = 'kg/kg'
+
+             CASE ( 'qi' )
+                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_extension ) THEN
+                   message_string = 'output of "' // TRIM( var ) // '" ' //    &
+                            'requires microphysics_ice_extension = .TRUE.'
+                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
                 ENDIF
                 unit = 'kg/kg'