Changeset 4486


Ignore:
Timestamp:
Apr 2, 2020 8:45:12 PM (18 months ago)
Author:
maronga
Message:

1.5-order-dai SGS schema modified for use with topography

File:
1 edited

Legend:

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

    r4481 r4486  
    2525! -----------------
    2626! $Id$
     27! Bugfix: include topography in calculation of distance_to_wall (1.5-order-dai
     28! closure)
     29!
     30! 4481 2020-03-31 18:55:54Z maronga
    2731! - added new LES closure after Dai et al. (2020), which provides much better grid
    2832!   convergence in stable boundary layer runs. The implementation is experimental
     
    14811485       DO  k = nzb, nzt+1
    14821486          delta(k,:,:) = gridsize_geometric_mean(k)
    1483 
    1484 !
    1485 !--       If Dai et al. (2020) closure is used, the distance to the wall must be stored
    1486           IF ( les_dai )  THEN
    1487              distance_to_wall(k,:,:) = zu(k)
    1488           ENDIF
    14891487       ENDDO
     1488
     1489!
     1490!--    If Dai et al. (2020) closure is used, the distance to the wall (distance to nearest upward facing surface)
     1491!--    must be stored
     1492       IF ( les_dai )  THEN
     1493          DO  i = nxl, nxr
     1494             DO  j = nys, nyn
     1495                DO  k = nzb, nzt+1
     1496                   distance_to_wall(k,j,i) = zu(k) - zw(topo_top_ind(j,i,0))
     1497                ENDDO
     1498             ENDDO
     1499          ENDDO
     1500       ENDIF
    14901501
    14911502       IF ( wall_adjustment )  THEN
     
    18701881    !$ACC ENTER DATA COPYIN(delta(nzb:nzt+1,nysg:nyng,nxlg:nxrg))
    18711882
    1872     IF ( les_dai )  THEN
    1873        CALL exchange_horiz( distance_to_wall, nbgp )
    1874        !$ACC ENTER DATA COPYIN(distance_to_wall(nzb:nzt+1,nysg:nyng,nxlg:nxrg))
    1875     ENDIF
     1883     IF ( les_dai )  THEN
     1884        CALL exchange_horiz( distance_to_wall, nbgp )
     1885        !$ACC ENTER DATA COPYIN(distance_to_wall(nzb:nzt+1,nysg:nyng,nxlg:nxrg))
     1886     ENDIF
    18761887
    18771888
Note: See TracChangeset for help on using the changeset viewer.