Changeset 2698 for palm/trunk/SOURCE/turbulence_closure_mod.f90

Timestamp:

Dec 14, 2017 6:46:24 PM (6 years ago)

Author:

suehring

Message:

Particle reflections at downward-facing walls; revision of particle speed interpolations at walls; bugfixes in get_topography_index and in date constants

File:

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

palm/trunk/SOURCE/turbulence_closure_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix in get_topography_top_index
 ! 
 ! Former revisions:
@@ -870 +870 @@
 
     USE surface_mod,                                                           &
-        ONLY:  get_topography_top_index
+        ONLY:  get_topography_top_index_ji
 
     IMPLICIT NONE
@@ -972 +972 @@
           DO  i = nxlg, nxrg
              DO  j = nysg, nyng
-                nz_s_shift = get_topography_top_index( j, i, 's' )
+                nz_s_shift = get_topography_top_index_ji( j, i, 's' )
 
                 e(nz_s_shift:nzt+1,j,i)  =  e(0:nzt+1-nz_s_shift,j,i)
@@ -991 +991 @@
           IF ( complex_terrain )  THEN
              IF ( nxlg <= 0 .AND. nxrg >= 0 .AND. nysg <= 0 .AND. nyng >= 0 )  THEN
-                nz_s_shift_l = get_topography_top_index( 0, 0, 's' )
+                nz_s_shift_l = get_topography_top_index_ji( 0, 0, 's' )
              ELSE
                 nz_s_shift_l = 0