Changeset 4079 for palm/trunk/SOURCE/nesting_offl_mod.f90
 Timestamp:
 Jul 9, 2019 6:04:41 PM (4 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/nesting_offl_mod.f90
r4022 r4079 25 25 !  26 26 ! $Id$ 27 !  Set boundary condition for w at nzt+1 at the lateral boundaries, even 28 ! though these won't enter the numerical solution. However, due to the mass 29 ! conservation these values might some up to very large values which will 30 ! occur in the runcontrol file 31 !  Bugfix in offline nesting of chemical species 32 !  Do not set Neumann conditions for TKE and passive scalar 33 ! 34 ! 4022 20190612 11:52:39Z suehring 27 35 ! Detection of boundarylayer depth in stable boundary layer on basis of 28 36 ! boundary data improved … … 353 361 BTEST( wall_flags_0(k,j,1), 3 ) ) 354 362 ENDDO 363 w(nzt,j,1) = w(nzt1,j,1) 355 364 ENDDO 356 365 … … 427 436 BTEST( wall_flags_0(k,j,nxr+1), 3 ) ) 428 437 ENDDO 438 w(nzt,j,nxr+1) = w(nzt1,j,nxr+1) 429 439 ENDDO 430 440 … … 504 514 BTEST( wall_flags_0(k,1,i), 3 ) ) 505 515 ENDDO 516 w(nzt,1,i) = w(nzt1,1,i) 506 517 ENDDO 507 518 … … 580 591 BTEST( wall_flags_0(k,nyn+1,i), 3 ) ) 581 592 ENDDO 593 w(nzt,nyn+1,i) = w(nzt1,nyn+1,i) 582 594 ENDDO 583 595 … … 710 722 DO j = nys, nyn 711 723 chem_species(n)%conc(nzt+1,j,i) = interpolate_in_time( & 712 nest_offl%chem_ north(0,j,i,n), &713 nest_offl%chem_ north(1,j,i,n), &724 nest_offl%chem_top(0,j,i,n), & 725 nest_offl%chem_top(1,j,i,n), & 714 726 fac_dt ) 715 727 ENDDO … … 727 739 IF ( bc_dirichlet_n ) diss(:,nyn+1,:) = diss(:,nyn,:) 728 740 ENDIF 729 IF ( .NOT. constant_diffusion ) THEN730 IF ( bc_dirichlet_l ) e(:,:,nxl1) = e(:,:,nxl)731 IF ( bc_dirichlet_r ) e(:,:,nxr+1) = e(:,:,nxr)732 IF ( bc_dirichlet_s ) e(:,nys1,:) = e(:,nys,:)733 IF ( bc_dirichlet_n ) e(:,nyn+1,:) = e(:,nyn,:)734 e(nzt+1,:,:) = e(nzt,:,:)735 ENDIF736 IF ( passive_scalar ) THEN737 IF ( bc_dirichlet_l ) s(:,:,nxl1) = s(:,:,nxl)738 IF ( bc_dirichlet_r ) s(:,:,nxr+1) = s(:,:,nxr)739 IF ( bc_dirichlet_s ) s(:,nys1,:) = s(:,nys,:)740 IF ( bc_dirichlet_n ) s(:,nyn+1,:) = s(:,nyn,:)741 ENDIF741 ! IF ( .NOT. constant_diffusion ) THEN 742 ! IF ( bc_dirichlet_l ) e(:,:,nxl1) = e(:,:,nxl) 743 ! IF ( bc_dirichlet_r ) e(:,:,nxr+1) = e(:,:,nxr) 744 ! IF ( bc_dirichlet_s ) e(:,nys1,:) = e(:,nys,:) 745 ! IF ( bc_dirichlet_n ) e(:,nyn+1,:) = e(:,nyn,:) 746 ! e(nzt+1,:,:) = e(nzt,:,:) 747 ! ENDIF 748 ! IF ( passive_scalar ) THEN 749 ! IF ( bc_dirichlet_l ) s(:,:,nxl1) = s(:,:,nxl) 750 ! IF ( bc_dirichlet_r ) s(:,:,nxr+1) = s(:,:,nxr) 751 ! IF ( bc_dirichlet_s ) s(:,nys1,:) = s(:,nys,:) 752 ! IF ( bc_dirichlet_n ) s(:,nyn+1,:) = s(:,nyn,:) 753 ! ENDIF 742 754 743 755 CALL exchange_horiz( u, nbgp ) … … 755 767 ENDDO 756 768 ENDIF 757 769 ! 770 ! Set top boundary condition at all horizontal grid points, also at the 771 ! lateral boundary grid points. 772 w(nzt+1,:,:) = w(nzt,:,:) 758 773 ! 759 774 ! In case of Rayleigh damping, where the profiles u_init, v_init
Note: See TracChangeset
for help on using the changeset viewer.