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

Timestamp:

Mar 22, 2007 9:54:05 AM (17 years ago)

Author:

raasch

Message:

preliminary update for changes concerning non-cyclic boundary conditions

File:

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

palm/trunk/SOURCE/init_pegrid.f90

@@ -5 +5 @@
 ! -----------------
 ! uxrp, vynp eliminated,
-! dirichlet/neumann changed to dirichlet/radiation, etc.
+! dirichlet/neumann changed to dirichlet/radiation, etc.,
+! poisfft_init is only called if fft-solver is switched on
 !
 ! Former revisions:
@@ -493 +494 @@
     nzta = nz
     nnz  = nz
-
-!
-!-- For non-cyclic boundaries extend array u (v) by one gridpoint
-!    IF ( bc_lr /= 'cyclic' )  uxrp = 1
-!    IF ( bc_ns /= 'cyclic' )  vynp = 1
 
 !
@@ -757 +753 @@
           inflow_r  = .TRUE.
        ENDIF
-!       uxrp      = 1
     ENDIF
 
@@ -774 +769 @@
           outflow_n = .TRUE.
        ENDIF
-!       vynp      = 1
-    ENDIF
-
-!
-!-- Additional MPI derived data type for the exchange of ghost points along x
-!-- needed in case of non-cyclic boundary conditions along y on the northmost
-!-- processors (for the exchange of the enlarged v array)
-    IF ( bc_ns /= 'cyclic'  .AND.  pnorth == MPI_PROC_NULL )  THEN
-       ngp_yz_p = ( nzt - nzb + 2 ) * ( nyn + vynp - nys + 3 )
-       CALL MPI_TYPE_VECTOR( nxr-nxl+3, nzt-nzb+2, ngp_yz_p, &
-                             MPI_REAL, type_xz_p, ierr )
-       CALL MPI_TYPE_COMMIT( type_xz_p, ierr )
-    ENDIF
+    ENDIF
+
 #else
     IF ( bc_lr == 'dirichlet/radiation' )  THEN
        inflow_l  = .TRUE.
        outflow_r = .TRUE.
-!       uxrp      = 1
     ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
        outflow_l = .TRUE.
@@ -803 +786 @@
        outflow_n = .TRUE.
        inflow_s  = .TRUE.
-!       vynp      = 1
     ENDIF
 #endif
@@ -809 +791 @@
     IF ( psolver == 'poisfft_hybrid' )  THEN
        CALL poisfft_hybrid_ini
-    ELSE
+    ELSEIF ( psolver == 'poisfft' )  THEN
        CALL poisfft_init
     ENDIF