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

Timestamp:

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

Author:

gronemeier

Message:

Implement turbulent outflow condition

File:

• palm/trunk/SOURCE/modules.f90 (modified) (4 diffs)

palm/trunk/SOURCE/modules.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Implement turbulent outflow condition
 ! 
 ! Former revisions:
@@ -783 +783 @@
                 terminate_run = .FALSE., transpose_compute_overlap = .FALSE., &
                 turbulent_inflow = .FALSE.
+    LOGICAL ::  turbulent_outflow = .FALSE. !< flag for turbulent outflow condition
     LOGICAL ::  urban_surface = .FALSE.  !< flag for urban surface model
     LOGICAL ::  use_cmax = .TRUE., use_initial_profile_as_reference = .FALSE., &
@@ -832 +833 @@
                  molecular_viscosity = 1.461E-5_wp, &
                  old_dt = 1.0E-10_wp, omega = 7.29212E-5_wp, omega_sor = 1.8_wp, &
+                 outflow_source_plane = -9999999.9_wp, & !< x-position of outflow-source plane (turbulent outflow method)
                  particle_maximum_age = 9999999.9_wp, &
                  phi = 55.0_wp, prandtl_number = 1.0_wp, &
@@ -1195 +1197 @@
     CHARACTER(LEN=2) ::  send_receive = 'al'
     CHARACTER(LEN=7) ::  myid_char = ''
-    INTEGER(iwp)          ::  acc_rank, comm1dx, comm1dy, comm2d, comm_inter,       &
-                              comm_palm, id_inflow = 0, id_recycling = 0, ierr,     &
-                              myid = 0, myidx = 0, myidy = 0, ndim = 2, ngp_a,      &
-                              ngp_o, ngp_xy, ngp_y, npex = -1, npey = -1,           &
-                              numprocs = 1, numprocs_previous_run = -1,             &
-                              num_acc_per_node = 0, pleft, pnorth, pright, psouth,  &
-                              req_count = 0, sendrecvcount_xy, sendrecvcount_yz,    &
-                              sendrecvcount_zx, sendrecvcount_zyd,                  &
-                              sendrecvcount_yxd, target_id, tasks_per_node = -9999, &
-                              threads_per_task = 1, type_x, type_xy,    &
-                              type_y
+    INTEGER(iwp)     ::  acc_rank, comm1dx, comm1dy, comm2d, comm_inter,       &
+                         comm_palm, id_inflow = 0,                             &
+                         id_outflow = 0,        & !< myidx of procs at outflow (turbulent outflow method)
+                         id_outflow_source = 0, & !< myidx of procs including ouflow source plane (turbulent outflow method)
+                         id_recycling = 0, ierr,                               &
+                         myid = 0, myidx = 0, myidy = 0, ndim = 2, ngp_a,      &
+                         ngp_o, ngp_xy, ngp_y, npex = -1, npey = -1,           &
+                         numprocs = 1, numprocs_previous_run = -1,             &
+                         num_acc_per_node = 0, pleft, pnorth, pright, psouth,  &
+                         req_count = 0, sendrecvcount_xy, sendrecvcount_yz,    &
+                         sendrecvcount_zx, sendrecvcount_zyd,                  &
+                         sendrecvcount_yxd, target_id, tasks_per_node = -9999, &
+                         threads_per_task = 1, type_x, type_xy,    &
+                         type_y
 
     INTEGER(iwp)          ::  pdims(2) = 1, req(100)