Changeset 2698 for palm/trunk/SOURCE/init_3d_model.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/init_3d_model.f90 (modified) (4 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Bugfix in get_topography_top_index
 ! 
 ! Former revisions:
@@ -466 +466 @@
     USE surface_mod,                                                           &
         ONLY :  init_surface_arrays, init_surfaces, surf_def_h, surf_lsm_h,    &
-                surf_usm_h, get_topography_top_index
+                surf_usm_h, get_topography_top_index_ji
    
     USE transpose_indices
@@ -1580 +1580 @@
           DO  i = nxlg, nxrg
              DO  j = nysg, nyng
-                nz_u_shift = get_topography_top_index( j, i, 'u' )
-                nz_v_shift = get_topography_top_index( j, i, 'v' )
-                nz_w_shift = get_topography_top_index( j, i, 'w' )
-                nz_s_shift = get_topography_top_index( j, i, 's' )
+                nz_u_shift = get_topography_top_index_ji( j, i, 'u' )
+                nz_v_shift = get_topography_top_index_ji( j, i, 'v' )
+                nz_w_shift = get_topography_top_index_ji( j, i, 'w' )
+                nz_s_shift = get_topography_top_index_ji( j, i, 's' )
 
                 u(nz_u_shift:nzt+1,j,i)  = u(0:nzt+1-nz_u_shift,j,i)                
@@ -1625 +1625 @@
           IF ( complex_terrain )  THEN
              IF ( nxlg <= 0 .AND. nxrg >= 0 .AND. nysg <= 0 .AND. nyng >= 0 )  THEN
-                nz_u_shift_l = get_topography_top_index( 0, 0, 'u' )
-                nz_v_shift_l = get_topography_top_index( 0, 0, 'v' )
-                nz_w_shift_l = get_topography_top_index( 0, 0, 'w' )
-                nz_s_shift_l = get_topography_top_index( 0, 0, 's' )
+                nz_u_shift_l = get_topography_top_index_ji( 0, 0, 'u' )
+                nz_v_shift_l = get_topography_top_index_ji( 0, 0, 'v' )
+                nz_w_shift_l = get_topography_top_index_ji( 0, 0, 'w' )
+                nz_s_shift_l = get_topography_top_index_ji( 0, 0, 's' )
              ELSE
                 nz_u_shift_l = 0