- Timestamp:
- Feb 5, 2019 5:21:16 PM (6 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 3 edited
-
boundary_conds.f90 (modified) (3 diffs)
-
module_interface.f90 (modified) (3 diffs)
-
virtual_measurement_mod.f90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/boundary_conds.f90
r3655 r3717 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Bugfix, do not set boundary conditions for potential temperature in neutral 28 ! runs. 29 ! 30 ! 3655 2019-01-07 16:51:22Z knoop 27 31 ! OpenACC port for SPEC 28 32 ! … … 235 239 humidity, ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, ibc_s_b, & 236 240 ibc_s_t, ibc_uv_b, ibc_uv_t, intermediate_timestep_count, & 237 nesting_offline, n udging, ocean_mode, passive_scalar, rans_mode,&238 rans_ tke_e, tsc, salsa, use_cmax241 nesting_offline, neutral, nudging, ocean_mode, passive_scalar, & 242 rans_mode, rans_tke_e, tsc, salsa, use_cmax 239 243 240 244 USE grid_variables, & … … 328 332 !-- the sea surface temperature of the coupled ocean model. 329 333 !-- Dirichlet 330 IF ( ibc_pt_b == 0 ) THEN 331 DO l = 0, 1 332 ! 333 !-- Set kb, for upward-facing surfaces value at topography top (k-1) is set, 334 !-- for downward-facing surfaces at topography bottom (k+1). 335 kb = MERGE( -1, 1, l == 0 ) 336 !$OMP PARALLEL DO PRIVATE( i, j, k ) 337 DO m = 1, bc_h(l)%ns 338 i = bc_h(l)%i(m) 339 j = bc_h(l)%j(m) 340 k = bc_h(l)%k(m) 341 pt_p(k+kb,j,i) = pt(k+kb,j,i) 342 ENDDO 343 ENDDO 344 ! 345 !-- Neumann, zero-gradient 346 ELSEIF ( ibc_pt_b == 1 ) THEN 347 DO l = 0, 1 348 ! 349 !-- Set kb, for upward-facing surfaces value at topography top (k-1) is set, 350 !-- for downward-facing surfaces at topography bottom (k+1). 351 kb = MERGE( -1, 1, l == 0 ) 352 !$OMP PARALLEL DO PRIVATE( i, j, k ) 353 !$ACC PARALLEL LOOP PRIVATE(i, j, k) & 354 !$ACC PRESENT(bc_h, pt_p) 355 DO m = 1, bc_h(l)%ns 356 i = bc_h(l)%i(m) 357 j = bc_h(l)%j(m) 358 k = bc_h(l)%k(m) 359 pt_p(k+kb,j,i) = pt_p(k,j,i) 360 ENDDO 361 ENDDO 362 ENDIF 363 364 ! 365 !-- Temperature at top boundary 366 IF ( ibc_pt_t == 0 ) THEN 367 pt_p(nzt+1,:,:) = pt(nzt+1,:,:) 368 ! 369 !-- In case of nudging adjust top boundary to pt which is 370 !-- read in from NUDGING-DATA 371 IF ( nudging ) THEN 372 pt_p(nzt+1,:,:) = pt_init(nzt+1) 373 ENDIF 374 ELSEIF ( ibc_pt_t == 1 ) THEN 375 pt_p(nzt+1,:,:) = pt_p(nzt,:,:) 376 ELSEIF ( ibc_pt_t == 2 ) THEN 377 !$ACC KERNELS PRESENT(pt_p, dzu) 378 pt_p(nzt+1,:,:) = pt_p(nzt,:,:) + bc_pt_t_val * dzu(nzt+1) 379 !$ACC END KERNELS 334 IF ( .NOT. neutral ) THEN 335 IF ( ibc_pt_b == 0 ) THEN 336 DO l = 0, 1 337 ! 338 !-- Set kb, for upward-facing surfaces value at topography top (k-1) 339 !-- is set, for downward-facing surfaces at topography bottom (k+1). 340 kb = MERGE( -1, 1, l == 0 ) 341 !$OMP PARALLEL DO PRIVATE( i, j, k ) 342 DO m = 1, bc_h(l)%ns 343 i = bc_h(l)%i(m) 344 j = bc_h(l)%j(m) 345 k = bc_h(l)%k(m) 346 pt_p(k+kb,j,i) = pt(k+kb,j,i) 347 ENDDO 348 ENDDO 349 ! 350 !-- Neumann, zero-gradient 351 ELSEIF ( ibc_pt_b == 1 ) THEN 352 DO l = 0, 1 353 ! 354 !-- Set kb, for upward-facing surfaces value at topography top (k-1) 355 !-- is set, for downward-facing surfaces at topography bottom (k+1). 356 kb = MERGE( -1, 1, l == 0 ) 357 !$OMP PARALLEL DO PRIVATE( i, j, k ) 358 !$ACC PARALLEL LOOP PRIVATE(i, j, k) & 359 !$ACC PRESENT(bc_h, pt_p) 360 DO m = 1, bc_h(l)%ns 361 i = bc_h(l)%i(m) 362 j = bc_h(l)%j(m) 363 k = bc_h(l)%k(m) 364 pt_p(k+kb,j,i) = pt_p(k,j,i) 365 ENDDO 366 ENDDO 367 ENDIF 368 369 ! 370 !-- Temperature at top boundary 371 IF ( ibc_pt_t == 0 ) THEN 372 pt_p(nzt+1,:,:) = pt(nzt+1,:,:) 373 ! 374 !-- In case of nudging adjust top boundary to pt which is 375 !-- read in from NUDGING-DATA 376 IF ( nudging ) THEN 377 pt_p(nzt+1,:,:) = pt_init(nzt+1) 378 ENDIF 379 ELSEIF ( ibc_pt_t == 1 ) THEN 380 pt_p(nzt+1,:,:) = pt_p(nzt,:,:) 381 ELSEIF ( ibc_pt_t == 2 ) THEN 382 !$ACC KERNELS PRESENT(pt_p, dzu) 383 pt_p(nzt+1,:,:) = pt_p(nzt,:,:) + bc_pt_t_val * dzu(nzt+1) 384 !$ACC END KERNELS 385 ENDIF 380 386 ENDIF 381 387 -
palm/trunk/SOURCE/module_interface.f90
r3711 r3717 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Add check_parameters routine for virtual measurements 28 ! 29 ! 3711 2019-01-31 13:44:26Z knoop 27 30 ! Introduced module_interface_init_checks for post-init checks 28 31 ! … … 297 300 298 301 USE virtual_measurement_mod, & 299 ONLY: vm_init, & 302 ONLY: vm_check_parameters, & 303 vm_init, & 300 304 vm_last_actions, & 301 305 vm_parin … … 492 496 IF ( syn_turb_gen ) CALL stg_check_parameters 493 497 IF ( urban_surface ) CALL usm_check_parameters 498 IF ( virtual_measurement ) CALL vm_check_parameters 494 499 IF ( wind_turbine ) CALL wtm_check_parameters 495 500 IF ( user_module_enabled ) CALL user_check_parameters -
palm/trunk/SOURCE/virtual_measurement_mod.f90
r3706 r3717 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Additional check + error numbers adjusted 28 ! 29 ! 3706 2019-01-29 20:02:26Z suehring 27 30 ! unused variables removed 28 31 ! … … 313 316 314 317 USE netcdf_data_input_mod, & 315 ONLY: input_pids_static 318 ONLY: input_pids_static, input_pids_vm 316 319 317 320 IMPLICIT NONE 318 321 322 ! 323 !-- Virtual measurements require a setup file. 324 IF ( virtual_measurement .AND. .NOT. input_pids_vm ) THEN 325 message_string = 'If virtual measurements are taken, a setup input ' // & 326 'file for the site locations is mandatory.' 327 CALL message( 'vm_check_parameters', 'PA0533', 1, 2, 0, 6, 0 ) 328 ENDIF 319 329 ! 320 330 !-- In case virtual measurements are taken, a static input file is required. … … 325 335 message_string = 'If virtual measurements are taken, a static input ' //& 326 336 'file is mandatory.' 327 CALL message( 'vm_check_parameters', 'PA0 000', 1, 2, 0, 6, 0 )337 CALL message( 'vm_check_parameters', 'PA0534', 1, 2, 0, 6, 0 ) 328 338 ENDIF 329 339 … … 522 532 TRIM( vmea(l)%feature_type ) // & 523 533 ' is not allowed.' 524 CALL message( 'vm_init', 'PA0 000', 1, 2, 0, 6, 0 )534 CALL message( 'vm_init', 'PA0535', 1, 2, 0, 6, 0 ) 525 535 ENDIF 526 536 !
Note: See TracChangeset
for help on using the changeset viewer.
