Changeset 4168 for palm/trunk/SOURCE/virtual_measurement_mod.f90
- Timestamp:
- Aug 16, 2019 1:50:17 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/virtual_measurement_mod.f90
r3988 r4168 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Replace function get_topography_top_index by topo_top_ind 28 ! 29 ! 3988 2019-05-22 11:32:37Z kanani 27 30 ! Add variables to enable steering of output interval for virtual measurements 28 31 ! … … 115 118 116 119 USE indices, & 117 ONLY: nzb, nzt, nxl, nxr, nys, nyn, nx, ny, wall_flags_0120 ONLY: nzb, nzt, nxl, nxr, nys, nyn, nx, ny, topo_top_ind, wall_flags_0 118 121 119 122 USE kinds … … 435 438 netcdf_data_input_get_dimension_length, & 436 439 netcdf_data_input_att, netcdf_data_input_var 437 438 USE surface_mod, &439 ONLY: get_topography_top_index_ji440 440 441 441 IMPLICIT NONE … … 763 763 !-- Determine vertical index which correspond to the observation 764 764 !-- height. 765 ksurf = get_topography_top_index_ji( js, is, 's')765 ksurf = topo_top_ind(js,is,0) 766 766 ks = MINLOC( ABS( zu - zw(ksurf) - z_ag(t,n) ), DIM = 1 ) - 1 767 767 ! … … 822 822 vmea(l)%j(ns) = j 823 823 vmea(l)%k(ns) = k 824 vmea(l)%z_ag(ns) = zu(k) - & 825 zw(get_topography_top_index_ji( j, i, 's' )) 824 vmea(l)%z_ag(ns) = zu(k) - zw(topo_top_ind(j,i,0)) 826 825 ENDIF 827 826 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.