Changeset 1992 for palm/trunk/SOURCE/check_parameters.f90
- Timestamp:
- Aug 12, 2016 3:14:59 PM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/check_parameters.f90
r1985 r1992 19 19 ! Current revisions: 20 20 ! ----------------- 21 ! 21 ! New checks for top_scalarflux 22 ! Bugfixes concerning data output of profiles in case of passive_scalar 22 23 ! 23 24 ! Former revisions: … … 1706 1707 message_string = 'boundary condition: bc_sa_t = "' // & 1707 1708 TRIM( bc_sa_t ) // '" is not allowed with ' // & 1708 ' constant_top_salinityflux = .TRUE.'1709 'top_salinityflux /= 0.0' 1709 1710 CALL message( 'check_parameters', 'PA0070', 1, 2, 0, 6, 0 ) 1710 1711 ENDIF … … 1749 1750 CALL check_bc_scalars( 's', bc_s_b, ibc_s_b, 'PA0073', 'PA0074', & 1750 1751 constant_scalarflux, s_surface_initial_change ) 1752 1753 IF ( top_scalarflux == 9999999.9_wp ) constant_top_scalarflux = .FALSE. 1754 ! 1755 !-- A fixed scalar concentration at the top implies Dirichlet boundary 1756 !-- condition for scalar. Hence, in this case specification of a constant 1757 !-- scalar flux is forbidden. 1758 IF ( ( ibc_s_t == 0 .OR. ibc_s_t == 2 ) .AND. constant_top_scalarflux & 1759 .AND. top_scalarflux /= 0.0_wp ) THEN 1760 message_string = 'boundary condition: bc_s_t = "' // & 1761 TRIM( bc_sa_t ) // '" is not allowed with ' // & 1762 'top_scalarflux /= 0.0' 1763 CALL message( 'check_parameters', 'PA0441', 1, 2, 0, 6, 0 ) 1764 ENDIF 1751 1765 ENDIF 1752 1766 ! … … 2210 2224 CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 ) 2211 2225 ELSE 2212 dopr_index(i) = 412226 dopr_index(i) = 117 2213 2227 dopr_unit(i) = 'kg/m3' 2214 hom(:,2, 41,:) = SPREAD( zu, 2, statistic_regions+1 )2228 hom(:,2,117,:) = SPREAD( zu, 2, statistic_regions+1 ) 2215 2229 IF ( data_output_pr(i)(1:1) == '#' ) THEN 2216 dopr_initial_index(i) = 262217 hom(:,2, 26,:)= SPREAD( zu, 2, statistic_regions+1 )2218 hom(nzb,2, 26,:)= 0.0_wp ! because zu(nzb) is negative2230 dopr_initial_index(i) = 117 2231 hom(:,2,117,:) = SPREAD( zu, 2, statistic_regions+1 ) 2232 hom(nzb,2,117,:) = 0.0_wp ! because zu(nzb) is negative 2219 2233 data_output_pr(i) = data_output_pr(i)(2:) 2220 2234 ENDIF … … 2331 2345 CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 ) 2332 2346 ELSE 2333 dopr_index(i) = 482347 dopr_index(i) = 119 2334 2348 dopr_unit(i) = 'kg/m3 m/s' 2335 hom(:,2, 48,:) = SPREAD( zw, 2, statistic_regions+1 )2349 hom(:,2,119,:) = SPREAD( zw, 2, statistic_regions+1 ) 2336 2350 ENDIF 2337 2351 … … 2343 2357 CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 ) 2344 2358 ELSE 2345 dopr_index(i) = 492359 dopr_index(i) = 116 2346 2360 dopr_unit(i) = 'kg/m3 m/s' 2347 hom(:,2, 49,:) = SPREAD( zw, 2, statistic_regions+1 )2361 hom(:,2,116,:) = SPREAD( zw, 2, statistic_regions+1 ) 2348 2362 ENDIF 2349 2363 … … 2355 2369 CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 ) 2356 2370 ELSE 2357 dopr_index(i) = 502371 dopr_index(i) = 120 2358 2372 dopr_unit(i) = 'kg/m3 m/s' 2359 hom(:,2, 50,:) = SPREAD( zw, 2, statistic_regions+1 )2373 hom(:,2,120,:) = SPREAD( zw, 2, statistic_regions+1 ) 2360 2374 ENDIF 2361 2375 … … 3916 3930 ELSEIF ( bc_t == 'neumann' ) THEN 3917 3931 ibc_t = 1 3932 ELSEIF ( bc_t == 'initial_gradient' ) THEN 3933 ibc_t = 2 3918 3934 ELSEIF ( bc_t == 'nested' ) THEN 3919 3935 ibc_t = 3
Note: See TracChangeset
for help on using the changeset viewer.