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

Timestamp:

Nov 8, 2016 3:00:55 PM (7 years ago)

Author:

gronemeier

Message:

Implement turbulent outflow condition

File:

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

palm/trunk/SOURCE/init_pegrid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Implement turbulent outflow condition
 ! 
 ! Former revisions:
@@ -169 +169 @@
                nest_bound_n, nest_bound_r, nest_bound_s, nest_domain, neutral, &
                psolver, outflow_l, outflow_n, outflow_r, outflow_s,            &
-               recycling_width, scalar_advec, subdomain_size 
+               outflow_source_plane, recycling_width, scalar_advec,            &
+               subdomain_size, turbulent_outflow
 
     USE grid_variables,                                                        &
@@ -199 +200 @@
     INTEGER(iwp) ::  i                        !<
     INTEGER(iwp) ::  id_inflow_l              !<
+    INTEGER(iwp) ::  id_outflow_l             !< local value of id_outflow
+    INTEGER(iwp) ::  id_outflow_source_l      !< local value of id_outflow_source
     INTEGER(iwp) ::  id_recycling_l           !<
     INTEGER(iwp) ::  ind(5)                   !<
@@ -1170 +1173 @@
     ENDIF
 
-       
 !
 !-- Broadcast the id of the inflow PE
@@ -1194 +1196 @@
     CALL MPI_ALLREDUCE( id_recycling_l, id_recycling, 1, MPI_INTEGER, MPI_SUM, &
                         comm1dx, ierr )
+
+!
+!-- Broadcast the id of the outflow PE and outflow-source plane
+    IF ( turbulent_outflow )  THEN
+
+       IF ( outflow_r )  THEN
+          id_outflow_l = myidx
+       ELSE
+          id_outflow_l = 0
+       ENDIF
+       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
+       CALL MPI_ALLREDUCE( id_outflow_l, id_outflow, 1, MPI_INTEGER, MPI_SUM, &
+                           comm1dx, ierr )
+
+       IF ( NINT( outflow_source_plane / dx ) >= nxl  .AND. &
+            NINT( outflow_source_plane / dx ) <= nxr )  THEN
+          id_outflow_source_l = myidx
+       ELSE
+          id_outflow_source_l = 0
+       ENDIF
+       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
+       CALL MPI_ALLREDUCE( id_outflow_source_l, id_outflow_source, 1, &
+                           MPI_INTEGER, MPI_SUM, comm1dx, ierr )
+
+    ENDIF
 
     CALL location_message( 'finished', .TRUE. )