Changeset 2747 for palm/trunk
- Timestamp:
- Jan 15, 2018 12:44:17 PM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_grid.f90
r2718 r2747 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Bugfix, topography height is rounded to the nearest discrete grid level 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 806 809 807 810 USE arrays_3d, & 808 ONLY: z w811 ONLY: zu, zw 809 812 810 813 USE control_parameters, & … … 1018 1021 !-- Hence, in this case, de-flag the grid point and give it 1019 1022 !-- urban type instead. 1020 IF ( z w(k) - ocean_offset <= terrain_height_f%var(j,i) ) THEN1023 IF ( zu(k) - ocean_offset <= terrain_height_f%var(j,i) ) THEN 1021 1024 topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 ) 1022 1025 topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 1 ) … … 1033 1036 building_id_f%var(j,i) ), DIM = 1 ) 1034 1037 1035 IF ( z w(k) - ocean_offset <= &1038 IF ( zu(k) - ocean_offset <= & 1036 1039 oro_max(nr) + buildings_f%var_2d(j,i) ) THEN 1037 1040 topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 ) … … 1081 1084 DO j = nys, nyn 1082 1085 DO k = nzb, nzt 1083 IF ( z w(k) - ocean_offset <= buildings_f%var_2d(j,i) ) THEN1086 IF ( zu(k) - ocean_offset <= buildings_f%var_2d(j,i) ) THEN 1084 1087 topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 ) 1085 1088 topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 ) !indicates terrain
Note: See TracChangeset
for help on using the changeset viewer.