Changeset 4625

Timestamp:

Jul 24, 2020 12:13:16 PM (4 years ago)

Author:

oliver.maas

Message:

for turbulence_closure = 1.5-order-dai: copy value of distance_to_wall to ghostpoints at dirichlet/radiation boundaries

File:

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

palm/trunk/SOURCE/turbulence_closure_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! for turbulence_closure = '1.5-order-dai': copy value of distance_to_wall to
+! ghostpoints at dirichlet/radiation boundaries
+!
+! 4586 2020-07-01 16:16:43Z gronemeier
 ! Renamed rif to richardson_number (Richardson Flux Number to Gradient Richardson Number)
 !
@@ -1919 +1923 @@
         CALL exchange_horiz( distance_to_wall, nbgp )
         !$ACC ENTER DATA COPYIN(distance_to_wall(nzb:nzt+1,nysg:nyng,nxlg:nxrg))
+
+!
+!--     Copy distance_to_wall to ghostpoints at inflow and outflow in case of
+!--     non-cyclic boundary conditions at inflow and outflow. This can not be done inside
+!--     exchange_horiz() because it is not compatible with prescribing mean inflow profiles.
+        IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
+           DO  i = nbgp, 1, -1
+             distance_to_wall(:,:,nxl-i) = distance_to_wall(:,:,nxl)
+           ENDDO
+        ENDIF
+        IF ( bc_dirichlet_r  .OR.  bc_radiation_r  )  THEN
+           DO  i = 1, nbgp
+              distance_to_wall(:,:,nxr+i) = distance_to_wall(:,:,nxr)
+           ENDDO
+        ENDIF
+        IF ( bc_dirichlet_s  .OR.  bc_radiation_s  )  THEN
+           DO  i = nbgp, 1, -1
+             distance_to_wall(:,nys-i,:) = distance_to_wall(:,nys,:)
+           ENDDO
+        ENDIF
+        IF ( bc_dirichlet_n  .OR.  bc_radiation_n  )  THEN
+           DO  i = 1, nbgp
+             distance_to_wall(:,nyn+i,:) = distance_to_wall(:,nyn,:)
+           ENDDO
+        ENDIF
      ENDIF