Changeset 724 for palm/tags/release-3.8/SOURCE/init_pegrid.f90

Timestamp:

Apr 11, 2011 6:42:09 AM (13 years ago)

Author:

raasch

Message:

bugfixes of r722 merged to release 3.8

Location:

palm/tags/release-3.8

Files:

• . (modified) (1 prop)
• SOURCE (modified) (1 prop)
• SOURCE/init_pegrid.f90 (modified) (8 diffs)

palm/tags/release-3.8/SOURCE/init_pegrid.f90

@@ -12 +12 @@
 ! -----------------
 ! $Id$
+!
+! 722 2011-04-11 06:21:09Z raasch
+! Bugfix: bc_lr/ns_cyc/dirrad/raddir replaced by bc_lr/ns, because variables
+!         are not yet set here; grid_level set to 0
 !
 ! 709 2011-03-30 09:31:40Z raasch
@@ -193 +197 @@
 !
 !-- If necessary, set horizontal boundary conditions to non-cyclic
-    IF ( .NOT. bc_lr_cyc )  cyclic(1) = .FALSE.
-    IF ( .NOT. bc_ns_cyc )  cyclic(2) = .FALSE.
+    IF ( bc_lr /= 'cyclic' )  cyclic(1) = .FALSE.
+    IF ( bc_ns /= 'cyclic' )  cyclic(2) = .FALSE.
 
 !
@@ -977 +981 @@
     ENDIF
 
-    grid_level = maximum_grid_level
+!
+!-- Default level 0 tells exchange_horiz that all ghost planes have to be
+!-- exchanged. grid_level is adjusted in poismg, where only one ghost plane
+!-- is required.
+    grid_level = 0
 
 #if defined( __parallel )
@@ -1067 +1075 @@
 !-- horizontal boundary conditions.
     IF ( pleft == MPI_PROC_NULL )  THEN
-       IF ( bc_lr_dirrad )  THEN
+       IF ( bc_lr == 'dirichlet/radiation' )  THEN
           inflow_l  = .TRUE.
-       ELSEIF ( bc_lr_raddir )  THEN
+       ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
           outflow_l = .TRUE.
        ENDIF
@@ -1075 +1083 @@
 
     IF ( pright == MPI_PROC_NULL )  THEN
-       IF ( bc_lr_dirrad )  THEN
+       IF ( bc_lr == 'dirichlet/radiation' )  THEN
           outflow_r = .TRUE.
-       ELSEIF ( bc_lr_raddir )  THEN
+       ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
           inflow_r  = .TRUE.
        ENDIF
@@ -1083 +1091 @@
 
     IF ( psouth == MPI_PROC_NULL )  THEN
-       IF ( bc_ns_dirrad )  THEN
+       IF ( bc_ns == 'dirichlet/radiation' )  THEN
           outflow_s = .TRUE.
-       ELSEIF ( bc_ns_raddir )  THEN
+       ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
           inflow_s  = .TRUE.
        ENDIF
@@ -1091 +1099 @@
 
     IF ( pnorth == MPI_PROC_NULL )  THEN
-       IF ( bc_ns_dirrad )  THEN
+       IF ( bc_ns == 'dirichlet/radiation' )  THEN
           inflow_n  = .TRUE.
-       ELSEIF ( bc_ns_raddir )  THEN
+       ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
           outflow_n = .TRUE.
        ENDIF
@@ -1123 +1131 @@
 
 #else
-    IF ( bc_lr_dirrad )  THEN
+    IF ( bc_lr == 'dirichlet/radiation' )  THEN
        inflow_l  = .TRUE.
        outflow_r = .TRUE.
-    ELSEIF ( bc_lr_raddir )  THEN
+    ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
        outflow_l = .TRUE.
        inflow_r  = .TRUE.
     ENDIF
 
-    IF ( bc_ns_dirrad )  THEN
+    IF ( bc_ns == 'dirichlet/radiation' )  THEN
        inflow_n  = .TRUE.
        outflow_s = .TRUE.
-    ELSEIF ( bc_ns_raddir )  THEN
+    ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
        outflow_n = .TRUE.
        inflow_s  = .TRUE.