Changeset 3981 for palm/trunk

Timestamp:

May 15, 2019 2:43:01 PM (6 years ago)

Author:

suehring

Message:

Bugfix in particlel nesting, TKE-gradients at ghost points at non-cyclic boundaries were not initialized

Location:

palm/trunk/SOURCE

Files:

• lpm_init.f90 (modified) (2 diffs)
• lpm_init_sgs_tke.f90 (modified) (3 diffs)

palm/trunk/SOURCE/lpm_init.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in particle nesting, initialize TKE gradient with zero
+! 
+! 3655 2019-01-07 16:51:22Z knoop
 ! Changed number concentration unit to SI standard
 ! 
@@ -431 +434 @@
                  de_dy(nzb:nzt+1,nysg:nyng,nxlg:nxrg), &
                  de_dz(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+       de_dx = 0.0_wp
+       de_dy = 0.0_wp
+       de_dz = 0.0_wp
 
        sgs_wf_part = 1.0_wp / 3.0_wp

palm/trunk/SOURCE/lpm_init_sgs_tke.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in particle nesting, set boundary conditions at non-cyclic boundaries
+! 
+! 3655 2019-01-07 16:51:22Z knoop
 ! unused variables removed
 ! 
@@ -81 +84 @@
         ONLY:  de_dx, de_dy, de_dz, diss, e, u, v, w, zu
 
+    USE control_parameters,                                                    &
+        ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r, bc_dirichlet_s
+        
     USE grid_variables,                                                        &
         ONLY:  ddx, ddy
@@ -206 +212 @@
     CALL exchange_horiz( de_dz, nbgp )
     CALL exchange_horiz( diss, nbgp  )
-
-
+!
+!-- Set boundary conditions at non-periodic boundaries. Note, at non-period
+!-- boundaries zero-gradient boundary conditions are set, so that the gradient
+!-- normal to the boundary is zero.
+    IF ( bc_dirichlet_l )  THEN
+       de_dx(:,:,-1) = 0.0_wp
+       de_dy(:,:,-1) = de_dy(:,:,0) 
+       de_dz(:,:,-1) = de_dz(:,:,0)
+    ENDIF
+    IF ( bc_dirichlet_r )  THEN
+       de_dx(:,:,nxr+1) = 0.0_wp
+       de_dy(:,:,nxr+1) = de_dy(:,:,nxr) 
+       de_dz(:,:,nxr+1) = de_dz(:,:,nxr)
+    ENDIF
+    IF ( bc_dirichlet_n )  THEN
+       de_dx(:,nyn+1,:) = de_dx(:,nyn,:)
+       de_dy(:,nyn+1,:) = 0.0_wp 
+       de_dz(:,nyn+1,:) = de_dz(:,nyn,:)
+    ENDIF
+    IF ( bc_dirichlet_s )  THEN
+       de_dx(:,nys-1,:) = de_dx(:,nys,:)
+       de_dy(:,nys-1,:) = 0.0_wp 
+       de_dz(:,nys-1,:) = de_dz(:,nys,:)
+    ENDIF  
 !
 !-- Calculate the horizontally averaged profiles of SGS TKE and resolved