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

Timestamp:

Mar 26, 2018 7:13:22 PM (6 years ago)

Author:

suehring

Message:

Synchronize parent and child model after initialization and spinup phase; Check for consistent setting of spinup times in parent and child model; remove obsolete masking of tendency arrays during initialization

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Check spinup in case of nested runs, spinup time and timestep must be 
+! identical to assure synchronuous simulations. 
+! 
+! 2932 2018-03-26 09:39:22Z maronga
 ! renamed particles_par to particle_parameters
 ! 
@@ -703 +707 @@
     LOGICAL     ::  found                            !<
 
+    REAL(wp)    ::  dt_spinup_max                    !< maximum spinup timestep in nested domains
     REAL(wp)    ::  dum                              !<
     REAL(wp)    ::  gradient                         !<
     REAL(wp)    ::  remote = 0.0_wp                  !<
+    REAL(wp)    ::  spinup_time_max                  !< maximum spinup time in nested domains
     REAL(wp)    ::  time_to_be_simulated_from_reference_point  !<
 
@@ -1397 +1403 @@
 #endif
     ENDIF
+!
+!-- In case of spinup and nested run, spinup end time must be identical 
+!-- in order to have synchronously running simulations. 
+    IF ( nested_run  .AND.  spinup )  THEN
+#if defined( __parallel )
+       CALL MPI_ALLREDUCE( spinup_time, spinup_time_max, 1, MPI_REAL,          &
+                           MPI_MAX, MPI_COMM_WORLD, ierr )
+       CALL MPI_ALLREDUCE( dt_spinup,   dt_spinup_max,   1, MPI_REAL,          &
+                           MPI_MAX, MPI_COMM_WORLD, ierr )
+
+       IF ( spinup_time /= spinup_time_max  .OR.  dt_spinup /= dt_spinup_max ) &
+       THEN
+          message_string = 'In case of nesting, spinup_time and ' //           &
+                           'dt_spinup must be identical in all parent ' //     &
+                           'and child domains.'
+          CALL message( 'check_parameters', 'PA0489', 3, 2, 0, 6, 0 )
+       ENDIF
+#endif
+    ENDIF
 
     IF ( cloud_physics  .AND.  .NOT.  humidity )  THEN