Changeset 1933 for palm/trunk/SOURCE/parin.f90

Timestamp:

Jun 13, 2016 7:12:51 AM (8 years ago)

Author:

hellstea

Message:

last commit documented

File:

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

palm/trunk/SOURCE/parin.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-!
+! 
 ! 
 ! Former revisions:
@@ -25 +25 @@
 ! $Id$
 !
+! 1917 2016-05-27 14:28:12Z witha
+! Initial version of purely vertical nesting introduced.
+! 
 ! 1914 2016-05-26 14:44:07Z witha
 ! Added call to wind turbine model for reading of &wind_turbine_par
@@ -245 +248 @@
 
     USE pmc_interface,                                                         &
-        ONLY:  nested_run
+        ONLY:  nested_run, nesting_mode
 
     USE profil_parameter,                                                      &
@@ -429 +432 @@
 
 !
-!--       In case of nested runs, explicitly set nesting boundary conditions
-!--       except for the root domain. This will overwrite the user settings.
-          IF ( nest_domain )  THEN
-             bc_lr   = 'nested'
-             bc_ns   = 'nested'
-             bc_uv_t = 'nested'
-             bc_pt_t = 'nested'
-             bc_q_t  = 'nested'
-             bc_p_t  = 'neumann'
+!--       In case of nested runs, explicitly set nesting boundary conditions. 
+!--       This will overwrite the user settings. bc_lr and bc_ns always need
+!--       to be cyclic for vertical nesting.
+          IF ( nesting_mode == 'vertical' )  THEN
+             IF (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' )  THEN
+                WRITE ( message_string, *)  'bc_lr and bc_ns were set to ,',   &
+                                      'cyclic for vertical nesting'
+                CALL message( 'parin', 'PA0428', 0, 0, 0, 6, 0 )
+                bc_lr   = 'cyclic'
+                bc_ns   = 'cyclic'
+             ENDIF
+             IF ( nest_domain )  THEN
+                bc_uv_t = 'nested'
+                bc_pt_t = 'nested'
+                bc_q_t  = 'nested'
+                bc_p_t  = 'neumann'
+             ENDIF
+          ELSE
+          
+!
+!--       For other nesting modes only set boundary conditions for 
+!--       nested domains.
+             IF ( nest_domain )  THEN
+                bc_lr   = 'nested'
+                bc_ns   = 'nested'
+                bc_uv_t = 'nested'
+                bc_pt_t = 'nested'
+                bc_q_t  = 'nested'
+                bc_p_t  = 'neumann'
+             ENDIF
           ENDIF
+            
 !
 !--       Check validity of lateral boundary conditions. This has to be done