- Timestamp:
- Apr 9, 2018 3:14:01 PM (7 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 3 edited
-
Makefile (modified) (2 diffs)
-
init_grid.f90 (modified) (4 diffs)
-
netcdf_data_input_mod.f90 (modified) (12 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/Makefile
r2938 r2955 25 25 # ----------------- 26 26 # $Id$ 27 # Add log-points to measure CPU time of NetCDF data input. 28 # 29 # 2938 2018-03-27 15:52:42Z suehring 27 30 # No initialization of child domains via dynamic input file, except for soil 28 31 # moisture and temperature … … 1169 1172 mod_kinds.o 1170 1173 netcdf_data_input_mod.o: \ 1174 cpulog_mod.o \ 1171 1175 mod_kinds.o \ 1172 1176 modules.o -
palm/trunk/SOURCE/init_grid.f90
r2927 r2955 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Improve topography filter routine and add ghost-point exchange for building 28 ! ID and building type. 29 ! 30 ! 2927 2018-03-23 15:13:00Z suehring 27 31 ! Bugfix, setting boundary conditions for topography index array. 28 32 ! … … 1186 1190 INTEGER(iwp) :: num_wall !< number of surrounding vertical walls for a single grid point 1187 1191 1188 INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE :: topo_tmp !< temporary 3D-topography used to fill holes 1189 INTEGER(iwp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: topo_3d !< 3D-topography array merging buildings and orography 1192 INTEGER(iwp), DIMENSION(nysg:nyng,nxlg:nxrg) :: var_exchange_int !< dummy array for exchanging ghost-points 1193 INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE :: topo_tmp !< temporary 3D-topography used to fill holes 1194 INTEGER(iwp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: topo_3d !< 3D-topography array merging buildings and orography 1190 1195 ! 1191 1196 !-- Before checking for holes, set lateral boundary conditions for … … 1201 1206 num_hole = 0 1202 1207 CALL exchange_horiz_int( topo_3d, nys, nyn, nxl, nxr, nzt, nbgp ) 1208 ! 1209 !-- Exchange also building ID and type. Note, building_type is an one-byte 1210 !-- variable. 1211 IF ( building_id_f%from_file ) & 1212 CALL exchange_horiz_2d_int( building_id_f%var, nys, nyn, nxl, nxr, nbgp ) 1213 IF ( building_type_f%from_file ) THEN 1214 var_exchange_int = INT( building_type_f%var, KIND = 4 ) 1215 CALL exchange_horiz_2d_int( var_exchange_int, nys, nyn, nxl, nxr, nbgp ) 1216 building_type_f%var = INT( var_exchange_int, KIND = 1 ) 1217 ENDIF 1203 1218 1204 1219 topo_tmp = topo_3d … … 1343 1358 IF ( nxl == 0 ) topo_3d(:,:,-1) = topo_3d(:,:,0) 1344 1359 IF ( nxr == nx ) topo_3d(:,:,nx+1) = topo_3d(:,:,nx) 1360 ENDIF 1361 ! 1362 !-- Exchange building ID and type. Note, building_type is an one-byte variable. 1363 IF ( building_id_f%from_file ) & 1364 CALL exchange_horiz_2d_int( building_id_f%var, nys, nyn, nxl, nxr, nbgp ) 1365 IF ( building_type_f%from_file ) THEN 1366 var_exchange_int = INT( building_type_f%var, KIND = 4 ) 1367 CALL exchange_horiz_2d_int( var_exchange_int, nys, nyn, nxl, nxr, nbgp ) 1368 building_type_f%var = INT( var_exchange_int, KIND = 1 ) 1345 1369 ENDIF 1346 1370 -
palm/trunk/SOURCE/netcdf_data_input_mod.f90
r2953 r2955 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Extend checks for consistent setting of buildings, its ID and type. 28 ! Add log-points to measure CPU time of NetCDF data input. 29 ! 30 ! 2953 2018-04-09 11:26:02Z suehring 27 31 ! Bugfix in checks for initialization data 28 32 ! … … 108 112 ONLY: coupling_char, io_blocks, io_group 109 113 114 USE cpulog, & 115 ONLY: cpu_log, log_point_s 116 110 117 USE kinds 111 118 … … 595 602 !-- If not static input file is available, skip this routine 596 603 IF ( .NOT. input_pids_static ) RETURN 604 ! 605 !-- Measure CPU time 606 CALL cpu_log( log_point_s(82), 'NetCDF input', 'start' ) 597 607 ! 598 608 !-- Read plant canopy variables. … … 1215 1225 #endif 1216 1226 ENDDO 1227 ! 1228 !-- End of CPU measurement 1229 CALL cpu_log( log_point_s(82), 'NetCDF input', 'stop' ) 1217 1230 ! 1218 1231 !-- Exchange 1 ghost points for surface variables. Please note, ghost point … … 1689 1702 1690 1703 REAL(wp) :: dum !< dummy variable to skip columns while reading topography file 1691 1704 ! 1705 !-- CPU measurement 1706 CALL cpu_log( log_point_s(83), 'NetCDF/ASCII input topo', 'start' ) 1692 1707 1693 1708 DO ii = 0, io_blocks-1 … … 1890 1905 ENDDO 1891 1906 ! 1907 !-- End of CPU measurement 1908 CALL cpu_log( log_point_s(83), 'NetCDF/ASCII input topo', 'stop' ) 1909 ! 1892 1910 !-- Check for minimum requirement to setup building topography. If buildings 1893 1911 !-- are provided, also an ID and a type are required. … … 2002 2020 IF ( bc_lr_cyc .AND. nxl == 0 ) nxlu = 1 2003 2021 IF ( bc_ns_cyc .AND. nys == 0 ) nysv = 1 2022 2023 ! 2024 !-- CPU measurement 2025 CALL cpu_log( log_point_s(85), 'NetCDF input init', 'start' ) 2004 2026 2005 2027 DO ii = 0, io_blocks-1 … … 2348 2370 #endif 2349 2371 ENDDO 2372 ! 2373 !-- End of CPU measurement 2374 CALL cpu_log( log_point_s(85), 'NetCDF input init', 'stop' ) 2350 2375 ! 2351 2376 !-- Finally, check if the input data has any fill values. Please note, … … 2471 2496 !-- Skip input if no forcing from larger-scale models is applied. 2472 2497 IF ( .NOT. forcing ) RETURN 2498 2499 ! 2500 !-- CPU measurement 2501 CALL cpu_log( log_point_s(86), 'NetCDF input forcing', 'start' ) 2473 2502 2474 2503 DO ii = 0, io_blocks-1 … … 2837 2866 2838 2867 ! 2868 !-- End of CPU measurement 2869 CALL cpu_log( log_point_s(86), 'NetCDF input forcing', 'stop' ) 2870 2871 ! 2839 2872 !-- Finally, after data input set control flag indicating that vertical 2840 2873 !-- inter- and/or extrapolation is required. … … 3250 3283 ENDIF 3251 3284 ! 3252 !-- Check if building_type is set at each building 3285 !-- Check if building_type is set at each building and vice versa. 3253 3286 IF ( building_type_f%from_file .AND. buildings_f%from_file ) THEN 3254 3287 IF ( buildings_f%lod == 1 ) THEN 3255 3288 IF ( buildings_f%var_2d(j,i) /= buildings_f%fill1 .AND. & 3256 3289 building_type_f%var(j,i) == building_type_f%fill ) THEN 3257 WRITE( message_string, * ) 'Each building requires ' // & 3258 'a type in case the ' // & 3259 'urban-surface model is ' // & 3260 'applied. i, j = ', i, j 3290 3291 WRITE( message_string, * ) 'Each location where a ' // & 3292 'building is set requires a type ' // & 3293 '( and vice versa ) in case the ' // & 3294 'urban-surface model is applied. ' // & 3295 'i, j = ', i, j 3261 3296 CALL message( 'netcdf_data_input_mod', 'NDI035', & 3262 3297 2, 2, 0, 6, 0 ) … … 3264 3299 ENDIF 3265 3300 IF ( buildings_f%lod == 2 ) THEN 3301 IF ( ANY( buildings_f%var_3d(:,j,i) == 1 ) .AND. & 3302 building_type_f%var(j,i) == building_type_f%fill ) THEN 3303 WRITE( message_string, * ) 'Each location where a ' // & 3304 'building is set requires a type ' // & 3305 '( and vice versa ) in case the ' // & 3306 'urban-surface model is applied. ' // & 3307 'i, j = ', i, j 3308 CALL message( 'netcdf_data_input_mod', 'NDI035', & 3309 2, 2, 0, 6, 0 ) 3310 ENDIF 3311 ENDIF 3312 ENDIF 3313 ! 3314 !-- Check if at each location where a building is present also an ID 3315 !-- is set and vice versa. 3316 IF ( buildings_f%from_file ) THEN 3317 IF ( buildings_f%lod == 1 ) THEN 3318 IF ( buildings_f%var_2d(j,i) /= buildings_f%fill1 .AND. & 3319 building_id_f%var(j,i) == building_id_f%fill ) THEN 3320 WRITE( message_string, * ) 'Each location where a ' // & 3321 'building is set requires an ID ' // & 3322 '( and vice versa ). i, j = ', i, j 3323 CALL message( 'netcdf_data_input_mod', 'NDI036', & 3324 2, 2, 0, 6, 0 ) 3325 ENDIF 3326 ELSEIF ( buildings_f%lod == 2 ) THEN 3327 IF ( ANY( buildings_f%var_3d(:,j,i) == 1 ) .AND. & 3328 building_id_f%var(j,i) == building_id_f%fill ) THEN 3329 WRITE( message_string, * ) 'Each location where a ' // & 3330 'building is set requires an ID ' // & 3331 '( and vice versa ). i, j = ', i, j 3332 CALL message( 'netcdf_data_input_mod', 'NDI036', & 3333 2, 2, 0, 6, 0 ) 3334 ENDIF 3335 ENDIF 3336 ENDIF 3337 ! 3338 !-- Check if at each location where a building ID or a -type is set 3339 !-- also a bulding is defined. 3340 IF ( buildings_f%from_file ) THEN 3341 IF ( buildings_f%lod == 1 ) THEN 3342 IF ( buildings_f%var_2d(j,i) /= buildings_f%fill1 .AND. & 3343 building_id_f%var(j,i) == building_id_f%fill ) THEN 3344 WRITE( message_string, * ) 'Each building grid point '// & 3345 'requires an ID.', i, j 3346 CALL message( 'netcdf_data_input_mod', 'NDI036', & 3347 2, 2, 0, 6, 0 ) 3348 ENDIF 3349 ELSEIF ( buildings_f%lod == 2 ) THEN 3266 3350 IF ( ANY( buildings_f%var_3d(:,j,i) == 1 ) & 3267 .AND. building_type_f%var(j,i) == building_type_f%fill ) THEN 3268 WRITE( message_string, * ) 'Each building requires ' // & 3269 'a type in case the ' // & 3270 'urban-surface model is ' // & 3271 'applied. i, j = ', i, j 3351 .AND. building_id_f%var(j,i) == building_id_f%fill ) THEN 3352 WRITE( message_string, * ) 'Each building grid point '// & 3353 'requires an ID.', i, j 3272 3354 CALL message( 'netcdf_data_input_mod', 'NDI036', & 3273 3355 2, 2, 0, 6, 0 )
Note: See TracChangeset
for help on using the changeset viewer.
