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

Timestamp:

Aug 21, 2017 2:59:59 PM (7 years ago)

Author:

kanani

Message:

Vertical nesting implemented (SadiqHuq?)

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Vertical grid nesting implemented: Check coupling mode. Generate file header
+! (SadiqHuq)
+! 
+! 2354 2017-08-17 10:49:36Z schwenkel
 ! Bugfix correlated to lsm_check_data_output_pr.
 ! If-statement for following checks is essential, otherwise units for lsm output 
@@ -561 +565 @@
     USE wind_turbine_model_mod,                                                &
         ONLY:  wtm_check_parameters, wind_turbine
+    USE vertical_nesting_mod,                                                  &
+        ONLY:  vnested, vnest_check_parameters
 
 
@@ -603 +609 @@
 !> @todo Check if any queries for other coupling modes (e.g. precursor_ocean) are missing
     IF ( coupling_mode /= 'uncoupled'            .AND.  &
+         coupling_mode /= 'vnested_crse'         .AND.  &
+         coupling_mode /= 'vnested_fine'         .AND.  &
          coupling_mode /= 'atmosphere_to_ocean'  .AND.  &
          coupling_mode /= 'ocean_to_atmosphere' )  THEN
@@ -611 +619 @@
 !
 !-- Check dt_coupling, restart_time, dt_restart, end_time, dx, dy, nx and ny
-    IF ( coupling_mode /= 'uncoupled')  THEN
+    IF ( coupling_mode /= 'uncoupled'  .AND.                                   &
+         coupling_mode(1:8) /= 'vnested_' )  THEN
 
        IF ( dt_coupling == 9999999.9_wp )  THEN
@@ -809 +818 @@
     IF ( coupling_mode == 'uncoupled' )  THEN
        coupling_string = ''
+    ELSEIF ( coupling_mode == 'vnested_crse' )  THEN
+       coupling_string = ' nested (coarse)'
+    ELSEIF ( coupling_mode == 'vnested_fine' )  THEN
+       coupling_string = ' nested (fine)'
     ELSEIF ( coupling_mode == 'atmosphere_to_ocean' )  THEN
        coupling_string = ' coupled (atmosphere)'
@@ -3858 +3871 @@
     ENDIF
 
+!-- Vertical nesting: check fine and coarse grid compatibility for data exchange
+    IF ( vnested )  CALL vnest_check_parameters
+
     CALL location_message( 'finished', .TRUE. )
 !