Changeset 2292 for palm/trunk/SOURCE/check_parameters.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/check_parameters.f90 (modified) (5 diffs)

palm/trunk/SOURCE/check_parameters.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). 
+! 
+! 2274 2017-06-09 13:27:48Z Giersch
 ! Changed error messages 
 ! 
@@ -983 +988 @@
        microphysics_seifert    = .TRUE.
        microphysics_kessler    = .FALSE.
+       microphysics_morrison  = .FALSE.
        precipitation           = .TRUE.
     ELSEIF ( cloud_scheme == 'kessler' )  THEN
@@ -988 +994 @@
        microphysics_seifert    = .FALSE.
        microphysics_kessler    = .TRUE.
+       microphysics_morrison   = .FALSE.
+       precipitation           = .TRUE.
+    ELSEIF ( cloud_scheme == 'morrison' )  THEN
+       microphysics_sat_adjust = .FALSE.
+       microphysics_seifert    = .TRUE.
+       microphysics_kessler    = .FALSE.
+       microphysics_morrison   = .TRUE.
        precipitation           = .TRUE.
     ELSE
@@ -2766 +2779 @@
              hom(:,2,122,:) = SPREAD( zw, 2, statistic_regions+1 )
 
+          CASE ( 'nc' )
+             IF (  .NOT.  cloud_physics )  THEN
+                message_string = 'data_output_pr = ' //                        &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for cloud_physics = .FALSE.'
+                CALL message( 'check_parameters', 'PA0094', 1, 2, 0, 6, 0 )
+             ELSEIF ( .NOT.  microphysics_morrison )  THEN
+                message_string = 'data_output_pr = ' //                        &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for cloud_scheme /= morrison'
+                CALL message( 'check_parameters', 'PA0358', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 89
+                dopr_unit(i)  = '1/m3'
+                hom(:,2,89,:)  = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
           CASE ( 'nr' )
              IF (  .NOT.  cloud_physics )  THEN
@@ -3061 +3091 @@
              ENDIF
              unit = 'K'
+
+          CASE ( 'nc' )
+             IF (  .NOT.  cloud_physics )  THEN
+                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+                         'res cloud_physics = .TRUE.'
+                CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
+             ELSEIF ( .NOT.  microphysics_morrison )  THEN
+                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+                         'res = microphysics morrison '
+                CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
+             ENDIF
+             unit = '1/m3'
 
           CASE ( 'nr' )