Changeset 2938 for palm/trunk/SOURCE/boundary_conds.f90

Timestamp:

Mar 27, 2018 3:52:42 PM (6 years ago)

Author:

suehring

Message:

Nesting in RANS-LES and RANS-RANS mode enabled; synthetic turbulence generator at all lateral boundaries in nesting or non-cyclic forcing mode; revised Inifor initialization in nesting mode

File:

• palm/trunk/SOURCE/boundary_conds.f90 (modified) (5 diffs)

palm/trunk/SOURCE/boundary_conds.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Set boundary condition for TKE and TKE dissipation rate in case of nesting
+! and if parent model operates in RANS mode but child model in LES mode.
+! mode
+! 
+! 2793 2018-02-07 10:54:33Z suehring
 ! Removed preprocessor directive __chem
 ! 
@@ -193 +198 @@
                inflow_s, intermediate_timestep_count,                          &
                microphysics_morrison, microphysics_seifert, nest_domain,       &
-               nest_bound_l, nest_bound_s, nudging, ocean, outflow_l,          &
-               outflow_n, outflow_r, outflow_s, passive_scalar, rans_tke_e,    &
-               tsc, use_cmax
+               nest_bound_l, nest_bound_n, nest_bound_r, nest_bound_s, nudging,&
+               ocean, outflow_l, outflow_n, outflow_r, outflow_s,              &
+               passive_scalar, rans_mode, rans_tke_e, tsc, use_cmax
 
     USE grid_variables,                                                        &
@@ -209 +214 @@
 
     USE pmc_interface,                                                         &
-        ONLY : nesting_mode
+        ONLY : nesting_mode, rans_mode_parent
 
     USE surface_mod,                                                           &
@@ -321 +326 @@
 
 !
-!-- Boundary conditions for TKE
-!-- Generally Neumann conditions with de/dz=0 are assumed
+!-- Boundary conditions for TKE.
+!-- Generally Neumann conditions with de/dz=0 are assumed.
     IF ( .NOT. constant_diffusion )  THEN
 
@@ -343 +348 @@
        IF ( .NOT. nest_domain )  THEN
           e_p(nzt+1,:,:) = e_p(nzt,:,:)
-       ENDIF
-    ENDIF
-
-!
-!-- Boundary conditions for TKE dissipation rate
+!
+!--    Nesting case: if parent operates in RANS mode and child in LES mode,
+!--    no TKE is transfered. This case, set Neumann conditions at lateral and 
+!--    top child boundaries. 
+!--    If not ( both either in RANS or in LES mode ), TKE boundary condition
+!--    is treated in the nesting. 
+       ELSE
+
+          IF ( rans_mode_parent  .AND.  .NOT. rans_mode )  THEN
+
+
+
+             e_p(nzt+1,:,:) = e_p(nzt,:,:)
+             IF ( nest_bound_l )  e_p(:,:,nxl-1) = e_p(:,:,nxl)
+             IF ( nest_bound_r )  e_p(:,:,nxr+1) = e_p(:,:,nxr)
+             IF ( nest_bound_s )  e_p(:,nys-1,:) = e_p(:,nys,:)
+             IF ( nest_bound_n )  e_p(:,nyn+1,:) = e_p(:,nyn,:)
+
+          ENDIF
+       ENDIF
+    ENDIF
+
+!
+!-- Boundary conditions for TKE dissipation rate. 
     IF ( rans_tke_e .AND. .NOT. nest_domain )  THEN
        diss_p(nzt+1,:,:) = diss_p(nzt,:,:)