Changeset 163 for palm/trunk/SOURCE/inflow_turbulence.f90

Timestamp:

May 5, 2008 2:09:05 PM (16 years ago)

Author:

raasch

Message:

bugfixes for turbulent inflow in init_pegrid, inflow_turbulence, and init_3d_model

File:

• palm/trunk/SOURCE/inflow_turbulence.f90 (modified) (2 diffs)

palm/trunk/SOURCE/inflow_turbulence.f90

@@ -91 +91 @@
     i = recycling_plane
 
-    IF ( i >= nxl  .AND.  i <= nxr )  THEN
+    IF ( myidx == id_recycling )  THEN
 
        DO  j = nys-1, nyn+1
@@ -110 +110 @@
 !-- For parallel runs, send the disturbances to the respective inflow PE
 #if defined( __parallel )
-    IF ( i >= nxl  .AND.  i <= nxr  .AND.  myid /= id_inflow )  THEN
+    IF ( myidx == id_recycling  .AND.  myidx /= id_inflow )  THEN
 
-!       print*, '*** sending id = ', myid, '   send to:', id_inflow
        CALL MPI_SEND( inflow_dist(nzb,nys-1,1), ngp_ifd, MPI_REAL, &
                       id_inflow, 1, comm1dx, ierr )
 
-    ELSEIF ( ( i < nxl  .OR.  i > nxr )  .AND.  myid == id_inflow )  THEN
+    ELSEIF ( myidx /= id_recycling  .AND.  myidx == id_inflow )  THEN
 
-!       print*, '*** receiving id = ', myid
+       inflow_dist = 0.0
        CALL MPI_RECV( inflow_dist(nzb,nys-1,1), ngp_ifd, MPI_REAL, &
-                      MPI_ANY_SOURCE, 1, comm1dx, status, ierr )
+                      id_recycling, 1, comm1dx, status, ierr )
 
     ENDIF