Changeset 3981 for palm/trunk/SOURCE


Ignore:
Timestamp:
May 15, 2019 2:43:01 PM (2 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:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/lpm_init.f90

    r3655 r3981  
    2525! -----------------
    2626! $Id$
     27! Bugfix in particle nesting, initialize TKE gradient with zero
     28!
     29! 3655 2019-01-07 16:51:22Z knoop
    2730! Changed number concentration unit to SI standard
    2831!
     
    431434                 de_dy(nzb:nzt+1,nysg:nyng,nxlg:nxrg), &
    432435                 de_dz(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     436       de_dx = 0.0_wp
     437       de_dy = 0.0_wp
     438       de_dz = 0.0_wp
    433439
    434440       sgs_wf_part = 1.0_wp / 3.0_wp
  • palm/trunk/SOURCE/lpm_init_sgs_tke.f90

    r3655 r3981  
    2525! -----------------
    2626! $Id$
     27! Bugfix in particle nesting, set boundary conditions at non-cyclic boundaries
     28!
     29! 3655 2019-01-07 16:51:22Z knoop
    2730! unused variables removed
    2831!
     
    8184        ONLY:  de_dx, de_dy, de_dz, diss, e, u, v, w, zu
    8285
     86    USE control_parameters,                                                    &
     87        ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r, bc_dirichlet_s
     88       
    8389    USE grid_variables,                                                        &
    8490        ONLY:  ddx, ddy
     
    206212    CALL exchange_horiz( de_dz, nbgp )
    207213    CALL exchange_horiz( diss, nbgp  )
    208 
    209 
     214!
     215!-- Set boundary conditions at non-periodic boundaries. Note, at non-period
     216!-- boundaries zero-gradient boundary conditions are set, so that the gradient
     217!-- normal to the boundary is zero.
     218    IF ( bc_dirichlet_l )  THEN
     219       de_dx(:,:,-1) = 0.0_wp
     220       de_dy(:,:,-1) = de_dy(:,:,0)
     221       de_dz(:,:,-1) = de_dz(:,:,0)
     222    ENDIF
     223    IF ( bc_dirichlet_r )  THEN
     224       de_dx(:,:,nxr+1) = 0.0_wp
     225       de_dy(:,:,nxr+1) = de_dy(:,:,nxr)
     226       de_dz(:,:,nxr+1) = de_dz(:,:,nxr)
     227    ENDIF
     228    IF ( bc_dirichlet_n )  THEN
     229       de_dx(:,nyn+1,:) = de_dx(:,nyn,:)
     230       de_dy(:,nyn+1,:) = 0.0_wp
     231       de_dz(:,nyn+1,:) = de_dz(:,nyn,:)
     232    ENDIF
     233    IF ( bc_dirichlet_s )  THEN
     234       de_dx(:,nys-1,:) = de_dx(:,nys,:)
     235       de_dy(:,nys-1,:) = 0.0_wp
     236       de_dz(:,nys-1,:) = de_dz(:,nys,:)
     237    ENDIF 
    210238!
    211239!-- Calculate the horizontally averaged profiles of SGS TKE and resolved
Note: See TracChangeset for help on using the changeset viewer.