- Timestamp:
- Oct 17, 2019 11:29:38 AM (5 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 6 edited
-
boundary_conds.f90 (modified) (23 diffs)
-
bulk_cloud_model_mod.f90 (modified) (6 diffs)
-
chemistry_model_mod.f90 (modified) (4 diffs)
-
module_interface.f90 (modified) (7 diffs)
-
salsa_mod.f90 (modified) (4 diffs)
-
time_integration.f90 (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/boundary_conds.f90
r4182 r4268 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Removing bulk cloud variables to respective module 28 ! 29 ! 4182 2019-08-22 15:20:23Z scharf 27 30 ! Corrected "Former revisions" section 28 31 ! … … 64 67 USE arrays_3d, & 65 68 ONLY: c_u, c_u_m, c_u_m_l, c_v, c_v_m, c_v_m_l, c_w, c_w_m, c_w_m_l, & 66 dzu, nc_p, nr_p, pt, pt_init, pt_p, q,&67 q_p, qc_p, qr_p, s, s_p, sa, sa_p, u, u_init, u_m_l, u_m_n,&69 dzu, pt, pt_init, pt_p, q, & 70 q_p, s, s_p, sa, sa_p, u, u_init, u_m_l, u_m_n, & 68 71 u_m_r, u_m_s, u_p, v, v_init, v_m_l, v_m_n, v_m_r, v_m_s, v_p, & 69 72 w, w_p, w_m_l, w_m_n, w_m_r, w_m_s 70 71 USE bulk_cloud_model_mod, &72 ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert73 73 74 74 USE chemistry_model_mod, & … … 102 102 103 103 USE salsa_mod, & 104 ONLY: salsa_boundary_conds 104 ONLY: salsa_boundary_conds 105 105 106 106 USE surface_mod, & … … 136 136 !$ACC PRESENT(bc_h, w_p) 137 137 DO m = 1, bc_h(l)%ns 138 i = bc_h(l)%i(m) 138 i = bc_h(l)%i(m) 139 139 j = bc_h(l)%j(m) 140 140 k = bc_h(l)%k(m) … … 174 174 !$OMP PARALLEL DO PRIVATE( i, j, k ) 175 175 DO m = 1, bc_h(l)%ns 176 i = bc_h(l)%i(m) 176 i = bc_h(l)%i(m) 177 177 j = bc_h(l)%j(m) 178 178 k = bc_h(l)%k(m) … … 223 223 !$OMP PARALLEL DO PRIVATE( i, j, k ) 224 224 DO m = 1, bc_h(l)%ns 225 i = bc_h(l)%i(m) 225 i = bc_h(l)%i(m) 226 226 j = bc_h(l)%j(m) 227 227 k = bc_h(l)%k(m) … … 253 253 !$OMP PARALLEL DO PRIVATE( i, j, k ) 254 254 DO m = 1, bc_h(l)%ns 255 i = bc_h(l)%i(m) 255 i = bc_h(l)%i(m) 256 256 j = bc_h(l)%j(m) 257 257 k = bc_h(l)%k(m) … … 259 259 ENDDO 260 260 ENDDO 261 261 262 262 ELSE 263 263 264 264 DO l = 0, 1 265 265 !$OMP PARALLEL DO PRIVATE( i, j, k ) 266 266 DO m = 1, bc_h(l)%ns 267 i = bc_h(l)%i(m) 267 i = bc_h(l)%i(m) 268 268 j = bc_h(l)%j(m) 269 269 k = bc_h(l)%k(m) … … 279 279 q_p(nzt+1,:,:) = q_p(nzt,:,:) + bc_q_t_val * dzu(nzt+1) 280 280 ENDIF 281 282 IF ( bulk_cloud_model .AND. microphysics_morrison ) THEN283 !284 !-- Surface conditions cloud water (Dirichlet)285 !-- Run loop over all non-natural and natural walls. Note, in wall-datatype286 !-- the k coordinate belongs to the atmospheric grid point, therefore, set287 !-- qr_p and nr_p at upward (k-1) and downward-facing (k+1) walls288 DO l = 0, 1289 !$OMP PARALLEL DO PRIVATE( i, j, k )290 DO m = 1, bc_h(l)%ns291 i = bc_h(l)%i(m)292 j = bc_h(l)%j(m)293 k = bc_h(l)%k(m)294 qc_p(k+bc_h(l)%koff,j,i) = 0.0_wp295 nc_p(k+bc_h(l)%koff,j,i) = 0.0_wp296 ENDDO297 ENDDO298 !299 !-- Top boundary condition for cloud water (Dirichlet)300 qc_p(nzt+1,:,:) = 0.0_wp301 nc_p(nzt+1,:,:) = 0.0_wp302 303 ENDIF304 305 IF ( bulk_cloud_model .AND. microphysics_seifert ) THEN306 !307 !-- Surface conditions rain water (Dirichlet)308 !-- Run loop over all non-natural and natural walls. Note, in wall-datatype309 !-- the k coordinate belongs to the atmospheric grid point, therefore, set310 !-- qr_p and nr_p at upward (k-1) and downward-facing (k+1) walls311 DO l = 0, 1312 !$OMP PARALLEL DO PRIVATE( i, j, k )313 DO m = 1, bc_h(l)%ns314 i = bc_h(l)%i(m)315 j = bc_h(l)%j(m)316 k = bc_h(l)%k(m)317 qr_p(k+bc_h(l)%koff,j,i) = 0.0_wp318 nr_p(k+bc_h(l)%koff,j,i) = 0.0_wp319 ENDDO320 ENDDO321 !322 !-- Top boundary condition for rain water (Dirichlet)323 qr_p(nzt+1,:,:) = 0.0_wp324 nr_p(nzt+1,:,:) = 0.0_wp325 326 ENDIF327 281 ENDIF 328 282 ! … … 336 290 !-- s_p at k-1 337 291 IF ( ibc_s_b == 0 ) THEN 338 292 339 293 DO l = 0, 1 340 294 !$OMP PARALLEL DO PRIVATE( i, j, k ) 341 295 DO m = 1, bc_h(l)%ns 342 i = bc_h(l)%i(m) 296 i = bc_h(l)%i(m) 343 297 j = bc_h(l)%j(m) 344 298 k = bc_h(l)%k(m) … … 346 300 ENDDO 347 301 ENDDO 348 302 349 303 ELSE 350 304 351 305 DO l = 0, 1 352 306 !$OMP PARALLEL DO PRIVATE( i, j, k ) 353 307 DO m = 1, bc_h(l)%ns 354 i = bc_h(l)%i(m) 308 i = bc_h(l)%i(m) 355 309 j = bc_h(l)%j(m) 356 310 k = bc_h(l)%k(m) … … 411 365 IF ( humidity ) THEN 412 366 q_p(:,nys-1,:) = q_p(:,nys,:) 413 IF ( bulk_cloud_model .AND. microphysics_morrison ) THEN414 qc_p(:,nys-1,:) = qc_p(:,nys,:)415 nc_p(:,nys-1,:) = nc_p(:,nys,:)416 ENDIF417 IF ( bulk_cloud_model .AND. microphysics_seifert ) THEN418 qr_p(:,nys-1,:) = qr_p(:,nys,:)419 nr_p(:,nys-1,:) = nr_p(:,nys,:)420 ENDIF421 367 ENDIF 422 368 IF ( passive_scalar ) s_p(:,nys-1,:) = s_p(:,nys,:) … … 425 371 IF ( humidity ) THEN 426 372 q_p(:,nyn+1,:) = q_p(:,nyn,:) 427 IF ( bulk_cloud_model .AND. microphysics_morrison ) THEN428 qc_p(:,nyn+1,:) = qc_p(:,nyn,:)429 nc_p(:,nyn+1,:) = nc_p(:,nyn,:)430 ENDIF431 IF ( bulk_cloud_model .AND. microphysics_seifert ) THEN432 qr_p(:,nyn+1,:) = qr_p(:,nyn,:)433 nr_p(:,nyn+1,:) = nr_p(:,nyn,:)434 ENDIF435 373 ENDIF 436 374 IF ( passive_scalar ) s_p(:,nyn+1,:) = s_p(:,nyn,:) … … 439 377 IF ( humidity ) THEN 440 378 q_p(:,:,nxl-1) = q_p(:,:,nxl) 441 IF ( bulk_cloud_model .AND. microphysics_morrison ) THEN442 qc_p(:,:,nxl-1) = qc_p(:,:,nxl)443 nc_p(:,:,nxl-1) = nc_p(:,:,nxl)444 ENDIF445 IF ( bulk_cloud_model .AND. microphysics_seifert ) THEN446 qr_p(:,:,nxl-1) = qr_p(:,:,nxl)447 nr_p(:,:,nxl-1) = nr_p(:,:,nxl)448 ENDIF449 379 ENDIF 450 380 IF ( passive_scalar ) s_p(:,:,nxl-1) = s_p(:,:,nxl) … … 453 383 IF ( humidity ) THEN 454 384 q_p(:,:,nxr+1) = q_p(:,:,nxr) 455 IF ( bulk_cloud_model .AND. microphysics_morrison ) THEN456 qc_p(:,:,nxr+1) = qc_p(:,:,nxr)457 nc_p(:,:,nxr+1) = nc_p(:,:,nxr)458 ENDIF459 IF ( bulk_cloud_model .AND. microphysics_seifert ) THEN460 qr_p(:,:,nxr+1) = qr_p(:,:,nxr)461 nr_p(:,:,nxr+1) = nr_p(:,:,nxr)462 ENDIF463 385 ENDIF 464 386 IF ( passive_scalar ) s_p(:,:,nxr+1) = s_p(:,:,nxr) … … 479 401 u_p(:,-1,:) = u(:,0,:) 480 402 v_p(:,0,:) = v(:,1,:) 481 w_p(:,-1,:) = w(:,0,:) 403 w_p(:,-1,:) = w(:,0,:) 482 404 ELSEIF ( .NOT. use_cmax ) THEN 483 405 … … 528 450 IF ( denom /= 0.0_wp ) THEN 529 451 c_w(k,i) = -c_max * ( w(k,0,i) - w_m_s(k,0,i) ) / ( denom * tsc(2) ) 452 IF ( c_w(k,i) < 0.0_wp ) THEN 453 c_w(k,i) = 0.0_wp 454 ELSEIF ( c_w(k,i) > c_max ) THEN 455 c_w(k,i) = c_max 456 ENDIF 457 ELSE 458 c_w(k,i) = c_max 459 ENDIF 460 461 c_u_m_l(k) = c_u_m_l(k) + c_u(k,i) 462 c_v_m_l(k) = c_v_m_l(k) + c_v(k,i) 463 c_w_m_l(k) = c_w_m_l(k) + c_w(k,i) 464 465 ENDDO 466 ENDDO 467 468 #if defined( __parallel ) 469 IF ( collective_wait ) CALL MPI_BARRIER( comm1dx, ierr ) 470 CALL MPI_ALLREDUCE( c_u_m_l(nzb+1), c_u_m(nzb+1), nzt-nzb, MPI_REAL, & 471 MPI_SUM, comm1dx, ierr ) 472 IF ( collective_wait ) CALL MPI_BARRIER( comm1dx, ierr ) 473 CALL MPI_ALLREDUCE( c_v_m_l(nzb+1), c_v_m(nzb+1), nzt-nzb, MPI_REAL, & 474 MPI_SUM, comm1dx, ierr ) 475 IF ( collective_wait ) CALL MPI_BARRIER( comm1dx, ierr ) 476 CALL MPI_ALLREDUCE( c_w_m_l(nzb+1), c_w_m(nzb+1), nzt-nzb, MPI_REAL, & 477 MPI_SUM, comm1dx, ierr ) 478 #else 479 c_u_m = c_u_m_l 480 c_v_m = c_v_m_l 481 c_w_m = c_w_m_l 482 #endif 483 484 c_u_m = c_u_m / (nx+1) 485 c_v_m = c_v_m / (nx+1) 486 c_w_m = c_w_m / (nx+1) 487 488 ! 489 !-- Save old timelevels for the next timestep 490 IF ( intermediate_timestep_count == 1 ) THEN 491 u_m_s(:,:,:) = u(:,0:1,:) 492 v_m_s(:,:,:) = v(:,1:2,:) 493 w_m_s(:,:,:) = w(:,0:1,:) 494 ENDIF 495 496 ! 497 !-- Calculate the new velocities 498 DO k = nzb+1, nzt+1 499 DO i = nxlg, nxrg 500 u_p(k,-1,i) = u(k,-1,i) - dt_3d * tsc(2) * c_u_m(k) * & 501 ( u(k,-1,i) - u(k,0,i) ) * ddy 502 503 v_p(k,0,i) = v(k,0,i) - dt_3d * tsc(2) * c_v_m(k) * & 504 ( v(k,0,i) - v(k,1,i) ) * ddy 505 506 w_p(k,-1,i) = w(k,-1,i) - dt_3d * tsc(2) * c_w_m(k) * & 507 ( w(k,-1,i) - w(k,0,i) ) * ddy 508 ENDDO 509 ENDDO 510 511 ! 512 !-- Bottom boundary at the outflow 513 IF ( ibc_uv_b == 0 ) THEN 514 u_p(nzb,-1,:) = 0.0_wp 515 v_p(nzb,0,:) = 0.0_wp 516 ELSE 517 u_p(nzb,-1,:) = u_p(nzb+1,-1,:) 518 v_p(nzb,0,:) = v_p(nzb+1,0,:) 519 ENDIF 520 w_p(nzb,-1,:) = 0.0_wp 521 522 ! 523 !-- Top boundary at the outflow 524 IF ( ibc_uv_t == 0 ) THEN 525 u_p(nzt+1,-1,:) = u_init(nzt+1) 526 v_p(nzt+1,0,:) = v_init(nzt+1) 527 ELSE 528 u_p(nzt+1,-1,:) = u_p(nzt,-1,:) 529 v_p(nzt+1,0,:) = v_p(nzt,0,:) 530 ENDIF 531 w_p(nzt:nzt+1,-1,:) = 0.0_wp 532 533 ENDIF 534 535 ENDIF 536 537 IF ( bc_radiation_n ) THEN 538 539 IF ( use_cmax ) THEN 540 u_p(:,ny+1,:) = u(:,ny,:) 541 v_p(:,ny+1,:) = v(:,ny,:) 542 w_p(:,ny+1,:) = w(:,ny,:) 543 ELSEIF ( .NOT. use_cmax ) THEN 544 545 c_max = dy / dt_3d 546 547 c_u_m_l = 0.0_wp 548 c_v_m_l = 0.0_wp 549 c_w_m_l = 0.0_wp 550 551 c_u_m = 0.0_wp 552 c_v_m = 0.0_wp 553 c_w_m = 0.0_wp 554 555 ! 556 !-- Calculate the phase speeds for u, v, and w, first local and then 557 !-- average along the outflow boundary. 558 DO k = nzb+1, nzt+1 559 DO i = nxl, nxr 560 561 denom = u_m_n(k,ny,i) - u_m_n(k,ny-1,i) 562 563 IF ( denom /= 0.0_wp ) THEN 564 c_u(k,i) = -c_max * ( u(k,ny,i) - u_m_n(k,ny,i) ) / ( denom * tsc(2) ) 565 IF ( c_u(k,i) < 0.0_wp ) THEN 566 c_u(k,i) = 0.0_wp 567 ELSEIF ( c_u(k,i) > c_max ) THEN 568 c_u(k,i) = c_max 569 ENDIF 570 ELSE 571 c_u(k,i) = c_max 572 ENDIF 573 574 denom = v_m_n(k,ny,i) - v_m_n(k,ny-1,i) 575 576 IF ( denom /= 0.0_wp ) THEN 577 c_v(k,i) = -c_max * ( v(k,ny,i) - v_m_n(k,ny,i) ) / ( denom * tsc(2) ) 578 IF ( c_v(k,i) < 0.0_wp ) THEN 579 c_v(k,i) = 0.0_wp 580 ELSEIF ( c_v(k,i) > c_max ) THEN 581 c_v(k,i) = c_max 582 ENDIF 583 ELSE 584 c_v(k,i) = c_max 585 ENDIF 586 587 denom = w_m_n(k,ny,i) - w_m_n(k,ny-1,i) 588 589 IF ( denom /= 0.0_wp ) THEN 590 c_w(k,i) = -c_max * ( w(k,ny,i) - w_m_n(k,ny,i) ) / ( denom * tsc(2) ) 530 591 IF ( c_w(k,i) < 0.0_wp ) THEN 531 592 c_w(k,i) = 0.0_wp … … 567 628 !-- Save old timelevels for the next timestep 568 629 IF ( intermediate_timestep_count == 1 ) THEN 569 u_m_s(:,:,:) = u(:,0:1,:)570 v_m_s(:,:,:) = v(:,1:2,:)571 w_m_s(:,:,:) = w(:,0:1,:)572 ENDIF573 574 !575 !-- Calculate the new velocities576 DO k = nzb+1, nzt+1577 DO i = nxlg, nxrg578 u_p(k,-1,i) = u(k,-1,i) - dt_3d * tsc(2) * c_u_m(k) * &579 ( u(k,-1,i) - u(k,0,i) ) * ddy580 581 v_p(k,0,i) = v(k,0,i) - dt_3d * tsc(2) * c_v_m(k) * &582 ( v(k,0,i) - v(k,1,i) ) * ddy583 584 w_p(k,-1,i) = w(k,-1,i) - dt_3d * tsc(2) * c_w_m(k) * &585 ( w(k,-1,i) - w(k,0,i) ) * ddy586 ENDDO587 ENDDO588 589 !590 !-- Bottom boundary at the outflow591 IF ( ibc_uv_b == 0 ) THEN592 u_p(nzb,-1,:) = 0.0_wp593 v_p(nzb,0,:) = 0.0_wp594 ELSE595 u_p(nzb,-1,:) = u_p(nzb+1,-1,:)596 v_p(nzb,0,:) = v_p(nzb+1,0,:)597 ENDIF598 w_p(nzb,-1,:) = 0.0_wp599 600 !601 !-- Top boundary at the outflow602 IF ( ibc_uv_t == 0 ) THEN603 u_p(nzt+1,-1,:) = u_init(nzt+1)604 v_p(nzt+1,0,:) = v_init(nzt+1)605 ELSE606 u_p(nzt+1,-1,:) = u_p(nzt,-1,:)607 v_p(nzt+1,0,:) = v_p(nzt,0,:)608 ENDIF609 w_p(nzt:nzt+1,-1,:) = 0.0_wp610 611 ENDIF612 613 ENDIF614 615 IF ( bc_radiation_n ) THEN616 617 IF ( use_cmax ) THEN618 u_p(:,ny+1,:) = u(:,ny,:)619 v_p(:,ny+1,:) = v(:,ny,:)620 w_p(:,ny+1,:) = w(:,ny,:)621 ELSEIF ( .NOT. use_cmax ) THEN622 623 c_max = dy / dt_3d624 625 c_u_m_l = 0.0_wp626 c_v_m_l = 0.0_wp627 c_w_m_l = 0.0_wp628 629 c_u_m = 0.0_wp630 c_v_m = 0.0_wp631 c_w_m = 0.0_wp632 633 !634 !-- Calculate the phase speeds for u, v, and w, first local and then635 !-- average along the outflow boundary.636 DO k = nzb+1, nzt+1637 DO i = nxl, nxr638 639 denom = u_m_n(k,ny,i) - u_m_n(k,ny-1,i)640 641 IF ( denom /= 0.0_wp ) THEN642 c_u(k,i) = -c_max * ( u(k,ny,i) - u_m_n(k,ny,i) ) / ( denom * tsc(2) )643 IF ( c_u(k,i) < 0.0_wp ) THEN644 c_u(k,i) = 0.0_wp645 ELSEIF ( c_u(k,i) > c_max ) THEN646 c_u(k,i) = c_max647 ENDIF648 ELSE649 c_u(k,i) = c_max650 ENDIF651 652 denom = v_m_n(k,ny,i) - v_m_n(k,ny-1,i)653 654 IF ( denom /= 0.0_wp ) THEN655 c_v(k,i) = -c_max * ( v(k,ny,i) - v_m_n(k,ny,i) ) / ( denom * tsc(2) )656 IF ( c_v(k,i) < 0.0_wp ) THEN657 c_v(k,i) = 0.0_wp658 ELSEIF ( c_v(k,i) > c_max ) THEN659 c_v(k,i) = c_max660 ENDIF661 ELSE662 c_v(k,i) = c_max663 ENDIF664 665 denom = w_m_n(k,ny,i) - w_m_n(k,ny-1,i)666 667 IF ( denom /= 0.0_wp ) THEN668 c_w(k,i) = -c_max * ( w(k,ny,i) - w_m_n(k,ny,i) ) / ( denom * tsc(2) )669 IF ( c_w(k,i) < 0.0_wp ) THEN670 c_w(k,i) = 0.0_wp671 ELSEIF ( c_w(k,i) > c_max ) THEN672 c_w(k,i) = c_max673 ENDIF674 ELSE675 c_w(k,i) = c_max676 ENDIF677 678 c_u_m_l(k) = c_u_m_l(k) + c_u(k,i)679 c_v_m_l(k) = c_v_m_l(k) + c_v(k,i)680 c_w_m_l(k) = c_w_m_l(k) + c_w(k,i)681 682 ENDDO683 ENDDO684 685 #if defined( __parallel )686 IF ( collective_wait ) CALL MPI_BARRIER( comm1dx, ierr )687 CALL MPI_ALLREDUCE( c_u_m_l(nzb+1), c_u_m(nzb+1), nzt-nzb, MPI_REAL, &688 MPI_SUM, comm1dx, ierr )689 IF ( collective_wait ) CALL MPI_BARRIER( comm1dx, ierr )690 CALL MPI_ALLREDUCE( c_v_m_l(nzb+1), c_v_m(nzb+1), nzt-nzb, MPI_REAL, &691 MPI_SUM, comm1dx, ierr )692 IF ( collective_wait ) CALL MPI_BARRIER( comm1dx, ierr )693 CALL MPI_ALLREDUCE( c_w_m_l(nzb+1), c_w_m(nzb+1), nzt-nzb, MPI_REAL, &694 MPI_SUM, comm1dx, ierr )695 #else696 c_u_m = c_u_m_l697 c_v_m = c_v_m_l698 c_w_m = c_w_m_l699 #endif700 701 c_u_m = c_u_m / (nx+1)702 c_v_m = c_v_m / (nx+1)703 c_w_m = c_w_m / (nx+1)704 705 !706 !-- Save old timelevels for the next timestep707 IF ( intermediate_timestep_count == 1 ) THEN708 630 u_m_n(:,:,:) = u(:,ny-1:ny,:) 709 631 v_m_n(:,:,:) = v(:,ny-1:ny,:) … … 731 653 u_p(nzb,ny+1,:) = 0.0_wp 732 654 v_p(nzb,ny+1,:) = 0.0_wp 733 ELSE 655 ELSE 734 656 u_p(nzb,ny+1,:) = u_p(nzb+1,ny+1,:) 735 657 v_p(nzb,ny+1,:) = v_p(nzb+1,ny+1,:) … … 757 679 u_p(:,:,0) = u(:,:,1) 758 680 v_p(:,:,-1) = v(:,:,0) 759 w_p(:,:,-1) = w(:,:,0) 681 w_p(:,:,-1) = w(:,:,0) 760 682 ELSEIF ( .NOT. use_cmax ) THEN 761 683 … … 870 792 u_p(nzb,:,0) = 0.0_wp 871 793 v_p(nzb,:,-1) = 0.0_wp 872 ELSE 794 ELSE 873 795 u_p(nzb,:,0) = u_p(nzb+1,:,0) 874 796 v_p(nzb,:,-1) = v_p(nzb+1,:,-1) … … 896 818 u_p(:,:,nx+1) = u(:,:,nx) 897 819 v_p(:,:,nx+1) = v(:,:,nx) 898 w_p(:,:,nx+1) = w(:,:,nx) 820 w_p(:,:,nx+1) = w(:,:,nx) 899 821 ELSEIF ( .NOT. use_cmax ) THEN 900 822 … … 1009 931 u_p(nzb,:,nx+1) = 0.0_wp 1010 932 v_p(nzb,:,nx+1) = 0.0_wp 1011 ELSE 933 ELSE 1012 934 u_p(nzb,:,nx+1) = u_p(nzb+1,:,nx+1) 1013 935 v_p(nzb,:,nx+1) = v_p(nzb+1,:,nx+1) -
palm/trunk/SOURCE/bulk_cloud_model_mod.f90
r4182 r4268 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Introducing bcm_boundary_conditions 28 ! 29 ! 4182 2019-08-22 15:20:23Z scharf 27 30 ! Corrected "Former revisions" section 28 31 ! … … 114 117 bc_radiation_n, & 115 118 bc_radiation_r, & 116 bc_radiation_s, & 119 bc_radiation_s, & 117 120 debug_output, & 118 121 dt_3d, dt_do2d_xy, intermediate_timestep_count, & … … 247 250 bcm_exchange_horiz, & 248 251 bcm_prognostic_equations, & 252 bcm_boundary_conditions, & 249 253 bcm_3d_data_averaging, & 250 254 bcm_data_output_2d, & … … 313 317 INTERFACE bcm_exchange_horiz 314 318 MODULE PROCEDURE bcm_exchange_horiz 315 END INTERFACE bcm_exchange_horiz 319 END INTERFACE bcm_exchange_horiz 316 320 317 321 INTERFACE bcm_prognostic_equations … … 319 323 MODULE PROCEDURE bcm_prognostic_equations_ij 320 324 END INTERFACE bcm_prognostic_equations 325 326 INTERFACE bcm_boundary_conditions 327 MODULE PROCEDURE bcm_boundary_conditions 328 END INTERFACE bcm_boundary_conditions 321 329 322 330 INTERFACE bcm_swap_timelevel … … 2014 2022 2015 2023 !------------------------------------------------------------------------------! 2024 ! Description: Boundary conditions of the bulk cloud module variables 2025 !------------------------------------------------------------------------------! 2026 SUBROUTINE bcm_boundary_conditions 2027 2028 IMPLICIT NONE 2029 2030 INTEGER(iwp) :: i !< 2031 INTEGER(iwp) :: j !< 2032 INTEGER(iwp) :: k !< 2033 INTEGER(iwp) :: m !< 2034 INTEGER(iwp) :: l !< 2035 2036 IF ( microphysics_morrison ) THEN 2037 ! 2038 !-- Surface conditions cloud water (Dirichlet) 2039 !-- Run loop over all non-natural and natural walls. Note, in wall-datatype 2040 !-- the k coordinate belongs to the atmospheric grid point, therefore, set 2041 !-- qr_p and nr_p at upward (k-1) and downward-facing (k+1) walls 2042 DO l = 0, 1 2043 !$OMP PARALLEL DO PRIVATE( i, j, k ) 2044 DO m = 1, bc_h(l)%ns 2045 i = bc_h(l)%i(m) 2046 j = bc_h(l)%j(m) 2047 k = bc_h(l)%k(m) 2048 qc_p(k+bc_h(l)%koff,j,i) = 0.0_wp 2049 nc_p(k+bc_h(l)%koff,j,i) = 0.0_wp 2050 ENDDO 2051 ENDDO 2052 ! 2053 !-- Top boundary condition for cloud water (Dirichlet) 2054 qc_p(nzt+1,:,:) = 0.0_wp 2055 nc_p(nzt+1,:,:) = 0.0_wp 2056 2057 ENDIF 2058 2059 IF ( microphysics_seifert ) THEN 2060 ! 2061 !-- Surface conditions rain water (Dirichlet) 2062 !-- Run loop over all non-natural and natural walls. Note, in wall-datatype 2063 !-- the k coordinate belongs to the atmospheric grid point, therefore, set 2064 !-- qr_p and nr_p at upward (k-1) and downward-facing (k+1) walls 2065 DO l = 0, 1 2066 !$OMP PARALLEL DO PRIVATE( i, j, k ) 2067 DO m = 1, bc_h(l)%ns 2068 i = bc_h(l)%i(m) 2069 j = bc_h(l)%j(m) 2070 k = bc_h(l)%k(m) 2071 qr_p(k+bc_h(l)%koff,j,i) = 0.0_wp 2072 nr_p(k+bc_h(l)%koff,j,i) = 0.0_wp 2073 ENDDO 2074 ENDDO 2075 ! 2076 !-- Top boundary condition for rain water (Dirichlet) 2077 qr_p(nzt+1,:,:) = 0.0_wp 2078 nr_p(nzt+1,:,:) = 0.0_wp 2079 2080 ENDIF 2081 2082 ! 2083 !-- Lateral boundary conditions for scalar quantities at the outflow. 2084 !-- Lateral oundary conditions for TKE and dissipation are set 2085 !-- in tcm_boundary_conds. 2086 IF ( bc_radiation_s ) THEN 2087 IF ( microphysics_morrison ) THEN 2088 qc_p(:,nys-1,:) = qc_p(:,nys,:) 2089 nc_p(:,nys-1,:) = nc_p(:,nys,:) 2090 ENDIF 2091 IF ( microphysics_seifert ) THEN 2092 qr_p(:,nys-1,:) = qr_p(:,nys,:) 2093 nr_p(:,nys-1,:) = nr_p(:,nys,:) 2094 ENDIF 2095 ELSEIF ( bc_radiation_n ) THEN 2096 IF ( microphysics_morrison ) THEN 2097 qc_p(:,nyn+1,:) = qc_p(:,nyn,:) 2098 nc_p(:,nyn+1,:) = nc_p(:,nyn,:) 2099 ENDIF 2100 IF ( microphysics_seifert ) THEN 2101 qr_p(:,nyn+1,:) = qr_p(:,nyn,:) 2102 nr_p(:,nyn+1,:) = nr_p(:,nyn,:) 2103 ENDIF 2104 ELSEIF ( bc_radiation_l ) THEN 2105 IF ( microphysics_morrison ) THEN 2106 qc_p(:,:,nxl-1) = qc_p(:,:,nxl) 2107 nc_p(:,:,nxl-1) = nc_p(:,:,nxl) 2108 ENDIF 2109 IF ( microphysics_seifert ) THEN 2110 qr_p(:,:,nxl-1) = qr_p(:,:,nxl) 2111 nr_p(:,:,nxl-1) = nr_p(:,:,nxl) 2112 ENDIF 2113 ELSEIF ( bc_radiation_r ) THEN 2114 IF ( microphysics_morrison ) THEN 2115 qc_p(:,:,nxr+1) = qc_p(:,:,nxr) 2116 nc_p(:,:,nxr+1) = nc_p(:,:,nxr) 2117 ENDIF 2118 IF ( microphysics_seifert ) THEN 2119 qr_p(:,:,nxr+1) = qr_p(:,:,nxr) 2120 nr_p(:,:,nxr+1) = nr_p(:,:,nxr) 2121 ENDIF 2122 ENDIF 2123 2124 END SUBROUTINE bcm_boundary_conditions 2125 2126 !------------------------------------------------------------------------------! 2016 2127 ! 2017 2128 ! Description: -
palm/trunk/SOURCE/chemistry_model_mod.f90
r4230 r4268 27 27 ! ----------------- 28 28 ! $Id$ 29 ! Moving module specific boundary conditions from time_integration to module 30 ! 31 ! 4230 2019-09-11 13:58:14Z suehring 29 32 ! Bugfix, initialize mean profiles also in restart runs. Also initialize 30 33 ! array used for Runge-Kutta tendecies in restart runs. … … 380 383 END INTERFACE chem_boundary_conds 381 384 385 INTERFACE chem_boundary_conditions 386 MODULE PROCEDURE chem_boundary_conditions 387 END INTERFACE chem_boundary_conditions 388 382 389 INTERFACE chem_check_data_output 383 390 MODULE PROCEDURE chem_check_data_output … … 537 544 538 545 539 PUBLIC chem_3d_data_averaging, chem_boundary_conds, 546 PUBLIC chem_3d_data_averaging, chem_boundary_conds, chem_boundary_conditions, & 540 547 chem_boundary_conds_decycle, chem_check_data_output, & 541 548 chem_check_data_output_pr, chem_check_parameters, & … … 799 806 END SUBROUTINE chem_boundary_conds 800 807 808 !------------------------------------------------------------------------------! 809 ! Description: 810 ! ------------ 811 !> Subroutine for boundary conditions 812 !------------------------------------------------------------------------------! 813 SUBROUTINE chem_boundary_conditions 814 815 IMPLICIT NONE 816 817 INTEGER(iwp) :: lsp !< 818 INTEGER(iwp) :: lsp_usr !< 819 820 ! 821 !-- Boundary conditions for prognostic quantitites of other modules: 822 !-- Here, only decycling is carried out 823 824 DO lsp = 1, nvar 825 lsp_usr = 1 826 DO WHILE ( TRIM( cs_name( lsp_usr ) ) /= 'novalue' ) 827 IF ( TRIM( chem_species(lsp)%name ) == TRIM( cs_name(lsp_usr) ) ) THEN 828 CALL chem_boundary_conds( chem_species(lsp)%conc_p, & 829 chem_species(lsp)%conc_pr_init ) 830 ENDIF 831 lsp_usr = lsp_usr + 1 832 ENDDO 833 ENDDO 834 835 836 END SUBROUTINE chem_boundary_conditions 801 837 802 838 !------------------------------------------------------------------------------! -
palm/trunk/SOURCE/module_interface.f90
r4182 r4268 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Introduction of module_interface_boundary_conditions 28 ! 29 ! 4182 2019-08-22 15:20:23Z scharf 27 30 ! Corrected "Former revisions" section 28 31 ! … … 245 248 bcm_exchange_horiz, & 246 249 bcm_prognostic_equations, & 250 bcm_boundary_conditions, & 247 251 bcm_swap_timelevel, & 248 252 bcm_3d_data_averaging, & … … 266 270 chem_non_advective_processes, & 267 271 chem_prognostic_equations, & 272 chem_boundary_conditions, & 268 273 chem_swap_timelevel, & 269 274 chem_3d_data_averaging, & … … 418 423 salsa_exchange_horiz_bounds, & 419 424 salsa_prognostic_equations, & 425 salsa_boundary_conditions, & 420 426 salsa_swap_timelevel, & 421 427 salsa_3d_data_averaging, & … … 520 526 module_interface_exchange_horiz, & 521 527 module_interface_prognostic_equations, & 528 module_interface_boundary_conditions, & 522 529 module_interface_swap_timelevel, & 523 530 module_interface_3d_data_averaging, & … … 598 605 MODULE PROCEDURE module_interface_swap_timelevel 599 606 END INTERFACE module_interface_swap_timelevel 607 608 INTERFACE module_interface_boundary_conditions 609 MODULE PROCEDURE module_interface_boundary_conditions 610 END INTERFACE module_interface_boundary_conditions 600 611 601 612 INTERFACE module_interface_3d_data_averaging … … 1273 1284 END SUBROUTINE module_interface_prognostic_equations_ij 1274 1285 1286 !------------------------------------------------------------------------------! 1287 ! Description: 1288 ! ------------ 1289 !> Compute module-specific boundary conditions 1290 !------------------------------------------------------------------------------! 1291 SUBROUTINE module_interface_boundary_conditions 1292 1293 1294 IF ( debug_output_timestep ) CALL debug_message( 'module-specific boundary_conditions', 'start' ) 1295 1296 IF ( bulk_cloud_model ) CALL bcm_boundary_conditions 1297 IF ( air_chemistry ) CALL chem_boundary_conditions 1298 IF ( salsa ) CALL salsa_boundary_conditions 1299 1300 IF ( debug_output_timestep ) CALL debug_message( 'module-specific boundary_conditions', 'end' ) 1301 1302 1303 END SUBROUTINE module_interface_boundary_conditions 1275 1304 1276 1305 !------------------------------------------------------------------------------! -
palm/trunk/SOURCE/salsa_mod.f90
r4256 r4268 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Moving module specific boundary conditions from time_integration to module 29 ! 30 ! 4256 2019-10-07 10:08:52Z monakurppa 28 31 ! Document previous changes: use global variables nx, ny and nz in salsa_header 29 32 ! … … 773 776 END INTERFACE salsa_boundary_conds 774 777 778 INTERFACE salsa_boundary_conditions 779 MODULE PROCEDURE salsa_boundary_conditions 780 END INTERFACE salsa_boundary_conditions 781 775 782 INTERFACE salsa_check_data_output 776 783 MODULE PROCEDURE salsa_check_data_output … … 866 873 ! 867 874 !-- Public functions: 868 PUBLIC salsa_boundary_conds, salsa_check_data_output, salsa_check_parameters, 875 PUBLIC salsa_boundary_conds, salsa_check_data_output, salsa_check_parameters, salsa_boundary_conditions, & 869 876 salsa_3d_data_averaging, salsa_data_output_2d, salsa_data_output_3d, & 870 877 salsa_data_output_mask, salsa_define_netcdf_grid, salsa_diagnostics, salsa_driver, & … … 7954 7961 END SUBROUTINE salsa_tendency 7955 7962 7963 7964 !------------------------------------------------------------------------------! 7965 ! Description: 7966 ! ------------ 7967 !> Boundary conditions for prognostic variables in SALSA from module interface 7968 !------------------------------------------------------------------------------! 7969 SUBROUTINE salsa_boundary_conditions 7970 7971 IMPLICIT NONE 7972 7973 INTEGER(iwp) :: ib !< index for aerosol size bins 7974 INTEGER(iwp) :: ic !< index for aerosol mass bins 7975 INTEGER(iwp) :: icc !< additional index for aerosol mass bins 7976 INTEGER(iwp) :: ig !< index for salsa gases 7977 7978 ! 7979 !-- Boundary conditions for prognostic quantitites of other modules: 7980 !-- Here, only decycling is carried out 7981 IF ( time_since_reference_point >= skip_time_do_salsa ) THEN 7982 7983 DO ib = 1, nbins_aerosol 7984 CALL salsa_boundary_conds( aerosol_number(ib)%conc_p, aerosol_number(ib)%init ) 7985 DO ic = 1, ncomponents_mass 7986 icc = ( ic - 1 ) * nbins_aerosol + ib 7987 CALL salsa_boundary_conds( aerosol_mass(icc)%conc_p, aerosol_mass(icc)%init ) 7988 ENDDO 7989 ENDDO 7990 IF ( .NOT. salsa_gases_from_chem ) THEN 7991 DO ig = 1, ngases_salsa 7992 CALL salsa_boundary_conds( salsa_gas(ig)%conc_p, salsa_gas(ig)%init ) 7993 ENDDO 7994 ENDIF 7995 7996 ENDIF 7997 7998 END SUBROUTINE salsa_boundary_conditions 7999 7956 8000 !------------------------------------------------------------------------------! 7957 8001 ! Description: -
palm/trunk/SOURCE/time_integration.f90
r4227 r4268 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Removing module specific boundary conditions an put them into their modules 28 ! 29 ! 4227 2019-09-10 18:04:34Z gronemeier 27 30 ! implement new palm_date_time_mod 28 31 ! … … 254 257 255 258 USE module_interface, & 256 ONLY: module_interface_actions, module_interface_swap_timelevel 259 ONLY: module_interface_actions, module_interface_swap_timelevel, & 260 module_interface_boundary_conditions 257 261 258 262 USE multi_agent_system_mod, & … … 583 587 !-- Execute all other module actions routunes 584 588 CALL module_interface_actions( 'before_timestep' ) 585 589 586 590 ! 587 591 !-- Start of intermediate step loop … … 750 754 !-- velocities at the outflow in case of a non-cyclic lateral wall) 751 755 CALL boundary_conds 752 753 ! 754 !-- Boundary conditions for prognostic quantitites of other modules: 755 !-- Here, only decycling is carried out 756 IF ( air_chemistry ) THEN 757 758 DO lsp = 1, nvar 759 lsp_usr = 1 760 DO WHILE ( TRIM( cs_name( lsp_usr ) ) /= 'novalue' ) 761 IF ( TRIM( chem_species(lsp)%name ) == TRIM( cs_name(lsp_usr) ) ) THEN 762 CALL chem_boundary_conds( chem_species(lsp)%conc_p, & 763 chem_species(lsp)%conc_pr_init ) 764 ENDIF 765 lsp_usr = lsp_usr + 1 766 ENDDO 767 ENDDO 768 769 ENDIF 770 771 IF ( salsa .AND. time_since_reference_point >= skip_time_do_salsa ) THEN 772 773 DO ib = 1, nbins_aerosol 774 CALL salsa_boundary_conds( aerosol_number(ib)%conc_p, aerosol_number(ib)%init ) 775 DO ic = 1, ncomponents_mass 776 icc = ( ic - 1 ) * nbins_aerosol + ib 777 CALL salsa_boundary_conds( aerosol_mass(icc)%conc_p, aerosol_mass(icc)%init ) 778 ENDDO 779 ENDDO 780 IF ( .NOT. salsa_gases_from_chem ) THEN 781 DO ig = 1, ngases_salsa 782 CALL salsa_boundary_conds( salsa_gas(ig)%conc_p, salsa_gas(ig)%init ) 783 ENDDO 784 ENDIF 785 786 ENDIF 787 756 ! 757 !-- Boundary conditions for module-specific variables 758 CALL module_interface_boundary_conditions 788 759 ! 789 760 !-- Incrementing timestep counter … … 1009 980 1010 981 ELSE 1011 ! 982 ! 1012 983 !-- Mass (volume) flux correction to ensure global mass conservation for child domains. 1013 984 IF ( child_domain ) THEN … … 1018 989 ENDIF 1019 990 ENDIF 1020 991 1021 992 CALL pres 1022 993 … … 1084 1055 CALL lsm_energy_balance( .FALSE., 3 ) 1085 1056 CALL lsm_soil_model( .FALSE., 3, .TRUE. ) 1086 1057 1087 1058 ! 1088 1059 !-- At the end, set boundary conditons for potential temperature … … 1091 1062 CALL lsm_boundary_condition 1092 1063 1093 1064 1094 1065 CALL cpu_log( log_point(54), 'land_surface', 'stop' ) 1095 1066 ENDIF … … 1099 1070 IF (urban_surface) THEN 1100 1071 CALL cpu_log( log_point(74), 'urban_surface', 'start' ) 1101 1072 1102 1073 CALL usm_surface_energy_balance( .FALSE. ) 1103 1074 IF ( usm_material_model ) THEN
Note: See TracChangeset
for help on using the changeset viewer.
