Changeset 2773 for palm/trunk/SOURCE/check_parameters.f90

Timestamp:

Jan 30, 2018 2:12:54 PM (6 years ago)

Author:

suehring

Message:

Nesting for chemical species implemented; Bugfix passive scalar boundary condition after anterpolation; Timeseries output of surface temperature; Enable initialization of 3D topography (was commented out so far)

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (3 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Check for consistent initialization in nesting mode added.
+! 
+! 2766 2018-01-22 17:17:47Z kanani
 ! Removed preprocessor directive __chem
 ! 
@@ -666 +669 @@
     INTEGER(iwp) ::  lsp                             !< running index for chem spcs.
 
+    LOGICAL     ::  check_nest                       !< flag used to check initialization in case of nesting
     LOGICAL     ::  found                            !<
 
@@ -1340 +1344 @@
        CALL message( 'check_parameters', 'PA0033', 1, 2, 0, 6, 0 )
     ENDIF
+!
+!-- In case of nested run assure that all domains are initialized the same
+!-- way, i.e. if at least at one domain is initialized with soil and 
+!-- atmospheric data provided by COSMO, all domains must be initialized the 
+!-- same way, to assure that soil and atmospheric quantities are 
+!-- consistent. 
+    IF ( nested_run )  THEN
+       check_nest = .TRUE.
+#if defined( __parallel )
+       CALL MPI_ALLREDUCE( TRIM( initializing_actions ) == 'inifor',           &
+                           check_nest, 1, MPI_LOGICAL,                         &
+                           MPI_LAND, MPI_COMM_WORLD, ierr )
+
+       IF ( TRIM( initializing_actions ) == 'inifor'  .AND.                    &
+            .NOT.  check_nest )  THEN
+          message_string = 'In case of nesting, if at least in one ' //        &
+                           'domain initializing_actions = inifor, '  //        &
+                           'all domains need to be initialized that way.'
+          CALL message( 'netcdf_data_input_mod', 'PA0430', 3, 2, 0, 6, 0 )
+       ENDIF
+#endif
+    ENDIF
 
     IF ( cloud_physics  .AND.  .NOT.  humidity )  THEN