Changeset 4079 for palm/trunk/SOURCE/nesting_offl_mod.f90

Timestamp:

Jul 9, 2019 6:04:41 PM (5 years ago)

Author:

suehring

Message:

Implementation of a monotonic flux limiter for vertical advection term in Wicker-Skamarock scheme. The flux limiter is currently only applied for passive scalars (passive scalar, chemical species, aerosols) within the region up to the highest topography, in order to avoid the built-up of large concentrations within poorly resolved cavities in urban environments. To enable the limiter monotonic_limiter_z = .T. must be set. Note, the limiter is currently only implemented for the cache-optimized version of advec_ws. Further changes in offline nesting: Set boundary condition for w at nzt+1 at all lateral boundaries (even though these won't enter the numerical solution), in order to avoid high vertical velocities in the run-control file which might built-up due to the mass-conservation; bugfix in offline nesting for chemical species

File:

• palm/trunk/SOURCE/nesting_offl_mod.f90 (modified) (8 diffs)

palm/trunk/SOURCE/nesting_offl_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! - Set boundary condition for w at nzt+1 at the lateral boundaries, even 
+!   though these won't enter the numerical solution. However, due to the mass
+!   conservation these values might some up to very large values which will
+!   occur in the run-control file
+! - Bugfix in offline nesting of chemical species
+! - Do not set Neumann conditions for TKE and passive scalar
+! 
+! 4022 2019-06-12 11:52:39Z suehring
 ! Detection of boundary-layer depth in stable boundary layer on basis of 
 ! boundary data improved 
@@ -353 +361 @@
                                    BTEST( wall_flags_0(k,j,-1), 3 ) )
              ENDDO
+             w(nzt,j,-1) = w(nzt-1,j,-1)
           ENDDO
 
@@ -427 +436 @@
                                     BTEST( wall_flags_0(k,j,nxr+1), 3 ) )
              ENDDO
+             w(nzt,j,nxr+1) = w(nzt-1,j,nxr+1)
           ENDDO
 
@@ -504 +514 @@
                                   BTEST( wall_flags_0(k,-1,i), 3 ) )
              ENDDO
+             w(nzt,-1,i) = w(nzt-1,-1,i)
           ENDDO
 
@@ -580 +591 @@
                                       BTEST( wall_flags_0(k,nyn+1,i), 3 ) )
              ENDDO
+             w(nzt,nyn+1,i) = w(nzt-1,nyn+1,i)
           ENDDO
 
@@ -710 +722 @@
                    DO  j = nys, nyn
                       chem_species(n)%conc(nzt+1,j,i) = interpolate_in_time(   &
-                                              nest_offl%chem_north(0,j,i,n),   &
-                                              nest_offl%chem_north(1,j,i,n),   &
+                                              nest_offl%chem_top(0,j,i,n),   &
+                                              nest_offl%chem_top(1,j,i,n),   &
                                               fac_dt                       )
                    ENDDO
@@ -727 +739 @@
           IF (  bc_dirichlet_n )  diss(:,nyn+1,:) = diss(:,nyn,:)
        ENDIF
-       IF ( .NOT. constant_diffusion )  THEN
-          IF (  bc_dirichlet_l )  e(:,:,nxl-1) = e(:,:,nxl)
-          IF (  bc_dirichlet_r )  e(:,:,nxr+1) = e(:,:,nxr)
-          IF (  bc_dirichlet_s )  e(:,nys-1,:) = e(:,nys,:)
-          IF (  bc_dirichlet_n )  e(:,nyn+1,:) = e(:,nyn,:)
-          e(nzt+1,:,:) = e(nzt,:,:)
-       ENDIF
-       IF ( passive_scalar )  THEN
-          IF (  bc_dirichlet_l )  s(:,:,nxl-1) = s(:,:,nxl)
-          IF (  bc_dirichlet_r )  s(:,:,nxr+1) = s(:,:,nxr)
-          IF (  bc_dirichlet_s )  s(:,nys-1,:) = s(:,nys,:)
-          IF (  bc_dirichlet_n )  s(:,nyn+1,:) = s(:,nyn,:)
-       ENDIF
+!        IF ( .NOT. constant_diffusion )  THEN
+!           IF (  bc_dirichlet_l )  e(:,:,nxl-1) = e(:,:,nxl)
+!           IF (  bc_dirichlet_r )  e(:,:,nxr+1) = e(:,:,nxr)
+!           IF (  bc_dirichlet_s )  e(:,nys-1,:) = e(:,nys,:)
+!           IF (  bc_dirichlet_n )  e(:,nyn+1,:) = e(:,nyn,:)
+!           e(nzt+1,:,:) = e(nzt,:,:)
+!        ENDIF
+!        IF ( passive_scalar )  THEN
+!           IF (  bc_dirichlet_l )  s(:,:,nxl-1) = s(:,:,nxl)
+!           IF (  bc_dirichlet_r )  s(:,:,nxr+1) = s(:,:,nxr)
+!           IF (  bc_dirichlet_s )  s(:,nys-1,:) = s(:,nys,:)
+!           IF (  bc_dirichlet_n )  s(:,nyn+1,:) = s(:,nyn,:)
+!        ENDIF
 
        CALL exchange_horiz( u, nbgp )
@@ -755 +767 @@
           ENDDO
        ENDIF
-       
+!
+!--    Set top boundary condition at all horizontal grid points, also at the
+!--    lateral boundary grid points. 
+       w(nzt+1,:,:) = w(nzt,:,:)       
 !
 !--    In case of Rayleigh damping, where the profiles u_init, v_init