Changeset 4486
 Timestamp:
 Apr 2, 2020 8:45:12 PM (13 months ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/turbulence_closure_mod.f90
r4481 r4486 25 25 !  26 26 ! $Id$ 27 ! Bugfix: include topography in calculation of distance_to_wall (1.5orderdai 28 ! closure) 29 ! 30 ! 4481 20200331 18:55:54Z maronga 27 31 !  added new LES closure after Dai et al. (2020), which provides much better grid 28 32 ! convergence in stable boundary layer runs. The implementation is experimental … … 1481 1485 DO k = nzb, nzt+1 1482 1486 delta(k,:,:) = gridsize_geometric_mean(k) 1483 1484 !1485 ! If Dai et al. (2020) closure is used, the distance to the wall must be stored1486 IF ( les_dai ) THEN1487 distance_to_wall(k,:,:) = zu(k)1488 ENDIF1489 1487 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 1490 1501 1491 1502 IF ( wall_adjustment ) THEN … … 1870 1881 !$ACC ENTER DATA COPYIN(delta(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) 1871 1882 1872 IF ( les_dai ) THEN1873 CALL exchange_horiz( distance_to_wall, nbgp )1874 !$ACC ENTER DATA COPYIN(distance_to_wall(nzb:nzt+1,nysg:nyng,nxlg:nxrg))1875 ENDIF1883 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 1876 1887 1877 1888
Note: See TracChangeset
for help on using the changeset viewer.