Ignore:
Timestamp:
May 5, 2008 2:09:05 PM (13 years ago)
Author:
raasch
Message:

bugfixes for turbulent inflow in init_pegrid, inflow_turbulence, and init_3d_model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/inflow_turbulence.f90

    r151 r163  
    9191    i = recycling_plane
    9292
    93     IF ( i >= nxl  .AND.  i <= nxr )  THEN
     93    IF ( myidx == id_recycling )  THEN
    9494
    9595       DO  j = nys-1, nyn+1
     
    110110!-- For parallel runs, send the disturbances to the respective inflow PE
    111111#if defined( __parallel )
    112     IF ( i >= nxl  .AND.  i <= nxr  .AND.  myid /= id_inflow )  THEN
     112    IF ( myidx == id_recycling  .AND.  myidx /= id_inflow )  THEN
    113113
    114 !       print*, '*** sending id = ', myid, '   send to:', id_inflow
    115114       CALL MPI_SEND( inflow_dist(nzb,nys-1,1), ngp_ifd, MPI_REAL, &
    116115                      id_inflow, 1, comm1dx, ierr )
    117116
    118     ELSEIF ( ( i < nxl  .OR.  i > nxr )  .AND.  myid == id_inflow )  THEN
     117    ELSEIF ( myidx /= id_recycling  .AND.  myidx == id_inflow )  THEN
    119118
    120 !       print*, '*** receiving id = ', myid
     119       inflow_dist = 0.0
    121120       CALL MPI_RECV( inflow_dist(nzb,nys-1,1), ngp_ifd, MPI_REAL, &
    122                       MPI_ANY_SOURCE, 1, comm1dx, status, ierr )
     121                      id_recycling, 1, comm1dx, status, ierr )
    123122
    124123    ENDIF
Note: See TracChangeset for help on using the changeset viewer.