Changeset 3987 for palm/trunk/SOURCE/init_3d_model.f90
- Timestamp:
- May 22, 2019 9:52:13 AM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_3d_model.f90
r3939 r3987 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Convert most location messages to debug messages to reduce output in 28 ! job logfile to a minimum 29 ! 30 ! 27 31 ! unused variable removed 28 32 ! … … 737 741 INTEGER(iwp) :: nz_s_shift_l !< topography-top index on scalar-grid, used to vertically shift initial profiles 738 742 739 CALL location_message( 'init_3d_model', 'start' ) 740 CALL location_message( 'allocating arrays', 'start' ) 743 744 CALL location_message( 'model initialization', 'start' ) 745 746 IF ( debug_output ) CALL debug_message( 'allocating arrays', 'start' ) 741 747 ! 742 748 !-- Allocate arrays … … 1087 1093 intermediate_timestep_count = 0 ! needed when simulated_time = 0.0 1088 1094 1089 CALL location_message( 'allocating arrays', 'finished' )1095 IF ( debug_output ) CALL debug_message( 'allocating arrays', 'end' ) 1090 1096 1091 1097 ! … … 1110 1116 !-- Initialization with provided input data derived from larger-scale model 1111 1117 IF ( INDEX( initializing_actions, 'inifor' ) /= 0 ) THEN 1112 CALL location_message( 'initializing with INIFOR', 'start' )1118 IF ( debug_output ) CALL debug_message( 'initializing with INIFOR', 'start' ) 1113 1119 ! 1114 1120 !-- Read initial 1D profiles or 3D data from NetCDF file, depending … … 1313 1319 CALL init_surfaces 1314 1320 1315 CALL location_message( 'initializing with INIFOR', 'finished' )1321 IF ( debug_output ) CALL debug_message( 'initializing with INIFOR', 'end' ) 1316 1322 ! 1317 1323 !-- Initialization via computed 1D-model profiles 1318 1324 ELSEIF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN 1319 1325 1320 CALL location_message( 'initializing with 1D model profiles', 'start' )1326 IF ( debug_output ) CALL debug_message( 'initializing with 1D model profiles', 'start' ) 1321 1327 ! 1322 1328 !-- Use solutions of the 1D model as initial profiles, … … 1380 1386 CALL init_surfaces 1381 1387 1382 CALL location_message( 'initializing with 1D model profiles', 'finished' )1388 IF ( debug_output ) CALL debug_message( 'initializing with 1D model profiles', 'end' ) 1383 1389 1384 1390 ELSEIF ( INDEX(initializing_actions, 'set_constant_profiles') /= 0 ) & 1385 1391 THEN 1386 1392 1387 CALL location_message( 'initializing with constant profiles', 'start' )1393 IF ( debug_output ) CALL debug_message( 'initializing with constant profiles', 'start' ) 1388 1394 1389 1395 ! … … 1448 1454 CALL init_surfaces 1449 1455 1450 CALL location_message( 'initializing with constant profiles', 'finished' )1456 IF ( debug_output ) CALL debug_message( 'initializing with constant profiles', 'end' ) 1451 1457 1452 1458 ELSEIF ( INDEX(initializing_actions, 'by_user') /= 0 ) & 1453 1459 THEN 1454 1460 1455 CALL location_message( 'initializing by user', 'start' )1461 IF ( debug_output ) CALL debug_message( 'initializing by user', 'start' ) 1456 1462 ! 1457 1463 !-- Pre-initialize surface variables, i.e. setting start- and end-indices … … 1463 1469 CALL user_init_3d_model 1464 1470 1465 CALL location_message( 'initializing by user', 'finished' )1466 1467 ENDIF 1468 1469 CALL location_message( 'initializing statistics, boundary conditions, etc.', 'start' )1471 IF ( debug_output ) CALL debug_message( 'initializing by user', 'end' ) 1472 1473 ENDIF 1474 1475 IF ( debug_output ) CALL debug_message( 'initializing statistics, boundary conditions, etc.', 'start' ) 1470 1476 1471 1477 ! … … 1601 1607 ENDIF 1602 1608 1603 CALL location_message( 'initializing statistics, boundary conditions, etc.', 'finished' )1609 IF ( debug_output ) CALL debug_message( 'initializing statistics, boundary conditions, etc.', 'end' ) 1604 1610 1605 1611 ELSEIF ( TRIM( initializing_actions ) == 'read_restart_data' .OR. & … … 1607 1613 THEN 1608 1614 1609 CALL location_message( 'initializing in case of restart / cyclic_fill', 'start' )1615 IF ( debug_output ) CALL debug_message( 'initializing in case of restart / cyclic_fill', 'start' ) 1610 1616 ! 1611 1617 !-- Initialize surface elements and its attributes, e.g. heat- and … … 1867 1873 IF ( passive_scalar ) ts_m = 0.0_wp 1868 1874 1869 CALL location_message( 'initializing in case of restart / cyclic_fill', 'finished' )1875 IF ( debug_output ) CALL debug_message( 'initializing in case of restart / cyclic_fill', 'end' ) 1870 1876 1871 1877 ELSE … … 2228 2234 TRIM( initializing_actions ) /= 'cyclic_fill' ) THEN 2229 2235 2236 IF ( debug_output ) CALL debug_message( 'creating disturbances + applying pressure solver', 'start' ) 2230 2237 ! 2231 2238 !-- Needed for both disturb_field and pres … … 2235 2242 !$ACC COPY(v(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) 2236 2243 2237 CALL location_message( 'creating initial disturbances', 'start' )2238 2244 CALL disturb_field( 'u', tend, u ) 2239 2245 CALL disturb_field( 'v', tend, v ) 2240 CALL location_message( 'creating initial disturbances', 'finished' )2241 2246 2242 2247 !$ACC DATA & … … 2255 2260 !$ACC COPYIN(bc_h(1)%k(1:bc_h(1)%ns)) 2256 2261 2257 CALL location_message( 'applying pressure solver', 'start' )2258 2262 n_sor = nsor_ini 2259 2263 CALL pres 2260 2264 n_sor = nsor 2261 CALL location_message( 'applying pressure solver', 'finished' )2262 2265 2263 2266 !$ACC END DATA 2264 2267 !$ACC END DATA 2268 2269 IF ( debug_output ) CALL debug_message( 'creating disturbances + applying pressure solver', 'end' ) 2265 2270 2266 2271 ENDIF … … 2466 2471 2467 2472 2468 CALL location_message( ' init_3d_model', 'finished' )2473 CALL location_message( 'model initialization', 'finished' ) 2469 2474 2470 2475 END SUBROUTINE init_3d_model
Note: See TracChangeset
for help on using the changeset viewer.