Changeset 4625 for palm


Ignore:
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:
1 edited

Legend:

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

    r4586 r4625  
    2525! -----------------
    2626! $Id$
     27! for turbulence_closure = '1.5-order-dai': copy value of distance_to_wall to
     28! ghostpoints at dirichlet/radiation boundaries
     29!
     30! 4586 2020-07-01 16:16:43Z gronemeier
    2731! Renamed rif to richardson_number (Richardson Flux Number to Gradient Richardson Number)
    2832!
     
    19191923        CALL exchange_horiz( distance_to_wall, nbgp )
    19201924        !$ACC ENTER DATA COPYIN(distance_to_wall(nzb:nzt+1,nysg:nyng,nxlg:nxrg))
     1925
     1926!
     1927!--     Copy distance_to_wall to ghostpoints at inflow and outflow in case of
     1928!--     non-cyclic boundary conditions at inflow and outflow. This can not be done inside
     1929!--     exchange_horiz() because it is not compatible with prescribing mean inflow profiles.
     1930        IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
     1931           DO  i = nbgp, 1, -1
     1932             distance_to_wall(:,:,nxl-i) = distance_to_wall(:,:,nxl)
     1933           ENDDO
     1934        ENDIF
     1935        IF ( bc_dirichlet_r  .OR.  bc_radiation_r  )  THEN
     1936           DO  i = 1, nbgp
     1937              distance_to_wall(:,:,nxr+i) = distance_to_wall(:,:,nxr)
     1938           ENDDO
     1939        ENDIF
     1940        IF ( bc_dirichlet_s  .OR.  bc_radiation_s  )  THEN
     1941           DO  i = nbgp, 1, -1
     1942             distance_to_wall(:,nys-i,:) = distance_to_wall(:,nys,:)
     1943           ENDDO
     1944        ENDIF
     1945        IF ( bc_dirichlet_n  .OR.  bc_radiation_n  )  THEN
     1946           DO  i = 1, nbgp
     1947             distance_to_wall(:,nyn+i,:) = distance_to_wall(:,nyn,:)
     1948           ENDDO
     1949        ENDIF
    19211950     ENDIF
    19221951
Note: See TracChangeset for help on using the changeset viewer.