Changeset 4168 for palm/trunk/SOURCE/lagrangian_particle_model_mod.f90
- Timestamp:
- Aug 16, 2019 1:50:17 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/lagrangian_particle_model_mod.f90
r4145 r4168 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Replace function get_topography_top_index by topo_top_ind 28 ! 29 ! 4145 2019-08-06 09:55:22Z schwenkel 27 30 ! Some reformatting 28 31 ! … … 227 230 USE indices, & 228 231 ONLY: nx, nxl, nxlg, nxrg, nxr, ny, nyn, nys, nyng, nysg, nz, nzb, & 229 nzb_max, nzt, wall_flags_0,nbgp, ngp_2dh_outer 232 nzb_max, nzt,nbgp, ngp_2dh_outer, & 233 topo_top_ind, & 234 wall_flags_0 230 235 231 236 USE kinds … … 259 264 260 265 USE surface_mod, & 261 ONLY: get_topography_top_index_ji, surf_def_h, surf_lsm_h, surf_usm_h,& 262 bc_h 266 ONLY: bc_h, & 267 surf_def_h, & 268 surf_lsm_h, & 269 surf_usm_h 263 270 264 271 #if defined( __parallel ) && !defined( __mpifh ) … … 1402 1409 !-- Determine surface level. Therefore, check for 1403 1410 !-- upward-facing wall on w-grid. 1404 k_surf = get_topography_top_index_ji( jp, ip, 'w')1411 k_surf = topo_top_ind(jp,ip,3) 1405 1412 IF ( seed_follows_topography ) THEN 1406 1413 ! … … 3469 3476 !-- above topography (Prandtl-layer height) 3470 3477 !-- Determine vertical index of topography top 3471 k_wall = get_topography_top_index_ji( jp, ip, 's')3478 k_wall = topo_top_ind(jp,ip,0) 3472 3479 3473 3480 IF ( constant_flux_layer .AND. zv(n) - zw(k_wall) < z_p ) THEN … … 3558 3565 ! 3559 3566 !-- Determine vertical index of topography top 3560 k_wall = get_topography_top_index_ji( jp,ip, 's')3567 k_wall = topo_top_ind(jp,ip,0) 3561 3568 3562 3569 IF ( constant_flux_layer .AND. zv(n) - zw(k_wall) < z_p ) THEN
Note: See TracChangeset
for help on using the changeset viewer.