Changeset 3884 for palm/trunk/SOURCE


Ignore:
Timestamp:
Apr 10, 2019 1:31:55 PM (2 years ago)
Author:
Giersch
Message:

id_recycling is only calculated in case of tubulent inflow

File:
1 edited

Legend:

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

    r3761 r3884  
    2525! -----------------
    2626! $Id$
     27! id_recycling is only calculated in case of tubulent inflow
     28!
     29! 3761 2019-02-25 15:31:42Z raasch
    2730! unused variable removed
    2831!
     
    263266               message_string, mg_switch_to_pe0_level, momentum_advec,         &
    264267               psolver, outflow_source_plane, recycling_width, scalar_advec,   &
    265                subdomain_size, turbulent_outflow, y_shift
     268               subdomain_size, turbulent_inflow, turbulent_outflow, y_shift
    266269
    267270    USE grid_variables,                                                        &
     
    13521355!-- Broadcast the id of the recycling plane
    13531356!-- WARNING: needs to be adjusted in case of inflows other than from left side!
    1354     IF ( NINT( recycling_width / dx ) >= nxl  .AND. &
    1355          NINT( recycling_width / dx ) <= nxr )  THEN
    1356        id_recycling_l = myidx
    1357     ELSE
    1358        id_recycling_l = 0
    1359     ENDIF
    1360     IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
    1361     CALL MPI_ALLREDUCE( id_recycling_l, id_recycling, 1, MPI_INTEGER, MPI_SUM, &
    1362                         comm1dx, ierr )
     1357    IF ( turbulent_inflow ) THEN
     1358   
     1359       IF ( NINT( recycling_width / dx ) >= nxl  .AND. &
     1360            NINT( recycling_width / dx ) <= nxr )  THEN
     1361          id_recycling_l = myidx
     1362       ELSE
     1363          id_recycling_l = 0
     1364       ENDIF
     1365       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
     1366       CALL MPI_ALLREDUCE( id_recycling_l, id_recycling, 1, MPI_INTEGER, MPI_SUM, &
     1367                           comm1dx, ierr )
     1368                           
     1369    ENDIF
    13631370
    13641371!
Note: See TracChangeset for help on using the changeset viewer.