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

Timestamp:

Feb 25, 2016 12:31:13 PM (9 years ago)

Author:

hellstea

Message:

Introduction of nested domain system

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! Previous version restored.
+! Introduction of nested domain feature
 ! 
 ! Former revisions:
@@ -1755 +1755 @@
 !-- Lateral boundary conditions
     IF ( bc_lr /= 'cyclic'  .AND.  bc_lr /= 'dirichlet/radiation'  .AND. &
-         bc_lr /= 'radiation/dirichlet' )  THEN
+         bc_lr /= 'radiation/dirichlet'  .AND.  bc_lr /= 'nested' )  THEN
        message_string = 'unknown boundary condition: bc_lr = "' // &
                         TRIM( bc_lr ) // '"'
@@ -1761 +1761 @@
     ENDIF
     IF ( bc_ns /= 'cyclic'  .AND.  bc_ns /= 'dirichlet/radiation'  .AND. &
-         bc_ns /= 'radiation/dirichlet' )  THEN
+         bc_ns /= 'radiation/dirichlet'  .AND.  bc_ns /= 'nested' )  THEN
        message_string = 'unknown boundary condition: bc_ns = "' // &
                         TRIM( bc_ns ) // '"'
@@ -1843 +1843 @@
     IF ( bc_p_t == 'dirichlet' )  THEN
        ibc_p_t = 0
-    ELSEIF ( bc_p_t == 'neumann' )  THEN
+!-- TO_DO: later set bc_p_t to neumann before, in case of nested domain
+    ELSEIF ( bc_p_t == 'neumann' .OR. bc_p_t == 'nested' )  THEN
        ibc_p_t = 1
     ELSE
@@ -1873 +1874 @@
     ELSEIF ( bc_pt_t == 'initial_gradient' )  THEN
        ibc_pt_t = 2
+    ELSEIF ( bc_pt_t == 'nested' )  THEN
+       ibc_pt_t = 3
     ELSE
        message_string = 'unknown boundary condition: bc_pt_t = "' // &
@@ -2015 +2018 @@
        ELSEIF ( bc_q_t == 'neumann' )  THEN
           ibc_q_t = 1
+       ELSEIF ( bc_q_t == 'nested' )  THEN
+          ibc_q_t = 3
        ELSE
           message_string = 'unknown boundary condition: bc_' // TRIM( sq ) // &
@@ -2098 +2103 @@
        ELSEIF ( bc_uv_t == 'neumann' )  THEN
           ibc_uv_t = 1
+       ELSEIF ( bc_uv_t == 'nested' )  THEN
+          ibc_uv_t = 3
        ELSE
           message_string = 'unknown boundary condition: bc_uv_t = "' // &
@@ -3861 +3868 @@
        ntdim_2d_xz(1) = ntdim_3d(1)
        ntdim_2d_yz(1) = ntdim_3d(1)
-             write(9,*) "ntdim_3d(av)=",ntdim_3d(0)
-             CALL local_flush(9)
+
     ENDIF