Changeset 73 for palm/trunk/SOURCE/init_pegrid.f90

Timestamp:

Mar 20, 2007 8:33:14 AM (17 years ago)

Author:

raasch

Message:

preliminary changes for radiation conditions

File:

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

palm/trunk/SOURCE/init_pegrid.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! uxrp, vynp eliminated,
+! dirichlet/neumann changed to dirichlet/radiation, etc.
 !
 ! Former revisions:
@@ -495 +496 @@
 !
 !-- For non-cyclic boundaries extend array u (v) by one gridpoint
-    IF ( bc_lr /= 'cyclic' )  uxrp = 1
-    IF ( bc_ns /= 'cyclic' )  vynp = 1
+!    IF ( bc_lr /= 'cyclic' )  uxrp = 1
+!    IF ( bc_ns /= 'cyclic' )  vynp = 1
 
 !
@@ -743 +744 @@
 !-- by one gridpoint on the left/rightmost (northest/southest) processor
     IF ( pleft == MPI_PROC_NULL )  THEN
-       IF ( bc_lr == 'dirichlet/neumann' )  THEN
+       IF ( bc_lr == 'dirichlet/radiation' )  THEN
           inflow_l  = .TRUE.
-       ELSEIF ( bc_lr == 'neumann/dirichlet' )  THEN
+       ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
           outflow_l = .TRUE.
        ENDIF
@@ -751 +752 @@
 
     IF ( pright == MPI_PROC_NULL )  THEN
-       IF ( bc_lr == 'dirichlet/neumann' )  THEN
+       IF ( bc_lr == 'dirichlet/radiation' )  THEN
           outflow_r = .TRUE.
-       ELSEIF ( bc_lr == 'neumann/dirichlet' )  THEN
+       ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
           inflow_r  = .TRUE.
        ENDIF
-       uxrp      = 1
+!       uxrp      = 1
     ENDIF
 
     IF ( psouth == MPI_PROC_NULL )  THEN
-       IF ( bc_ns == 'dirichlet/neumann' )  THEN
+       IF ( bc_ns == 'dirichlet/radiation' )  THEN
           outflow_s = .TRUE.
-       ELSEIF ( bc_ns == 'neumann/dirichlet' )  THEN
+       ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
           inflow_s  = .TRUE.
        ENDIF
@@ -768 +769 @@
 
     IF ( pnorth == MPI_PROC_NULL )  THEN
-       IF ( bc_ns == 'dirichlet/neumann' )  THEN
+       IF ( bc_ns == 'dirichlet/radiation' )  THEN
           inflow_n  = .TRUE.
-       ELSEIF ( bc_ns == 'neumann/dirichlet' )  THEN
+       ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
           outflow_n = .TRUE.
        ENDIF
-       vynp      = 1
+!       vynp      = 1
     ENDIF
 
@@ -787 +788 @@
     ENDIF
 #else
-    IF ( bc_lr == 'dirichlet/neumann' )  THEN
+    IF ( bc_lr == 'dirichlet/radiation' )  THEN
        inflow_l  = .TRUE.
        outflow_r = .TRUE.
-       uxrp      = 1
-    ELSEIF ( bc_lr == 'neumann/dirichlet' )  THEN
+!       uxrp      = 1
+    ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
        outflow_l = .TRUE.
        inflow_r  = .TRUE.
     ENDIF
 
-    IF ( bc_ns == 'dirichlet/neumann' )  THEN
+    IF ( bc_ns == 'dirichlet/radiation' )  THEN
        inflow_n  = .TRUE.
        outflow_s = .TRUE.
-    ELSEIF ( bc_ns == 'neumann/dirichlet' )  THEN
+    ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
        outflow_n = .TRUE.
        inflow_s  = .TRUE.
-       vynp      = 1
+!       vynp      = 1
     ENDIF
 #endif