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

Timestamp:

Mar 27, 2018 3:52:42 PM (7 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/turbulence_closure_mod.f90 (modified) (8 diffs)

palm/trunk/SOURCE/turbulence_closure_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Further todo's
+! 
+! 2936 2018-03-27 14:49:27Z suehring
 ! - defined l_grid only within this module
 ! - Moved l_wall definition from modules.f90
@@ -69 +72 @@
 !>       add OpenMP directives whereever possible
 !>       remove debug output variables (dummy1, dummy2, dummy3)
+!> @todo Move initialization of wall-mixing length from init_grid
+!> @todo Check for random disturbances
 !> @note <Enter notes on the module>
 !> @bug  TKE-e closure still crashes due to too small dt
@@ -284 +289 @@
 
     USE control_parameters,                                                    &
-        ONLY:  message_string, neutral, turbulent_inflow, turbulent_outflow
+        ONLY:  message_string, nest_domain, neutral, turbulent_inflow,         &
+               turbulent_outflow
 
     IMPLICIT NONE
@@ -302 +308 @@
              rans_tke_e = .TRUE.
 
-             IF ( INDEX( initializing_actions, 'set_1d-model_profiles' )       &
-                  == 0 )  THEN
+             IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) == 0  &
+                  .AND.  .NOT.  nest_domain )  THEN
                 message_string = 'Initializing without 1D model while ' //     &
                                  'using TKE-e closure&' //                     &
@@ -826 +832 @@
         ONLY:  use_sgs_for_particles, wang_kernel
 
+    USE pmc_interface,                                                         &
+        ONLY:  nested_run
+
     IMPLICIT NONE
 
@@ -847 +856 @@
     ALLOCATE( e_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
 #endif
-
+!
+!-- Allocate arrays required for dissipation. 
+!-- Please note, if it is a nested run, arrays need to be allocated even if
+!-- they do not necessarily need to be transferred, which is attributed to 
+!-- the design of the model coupler which allocates memory for each variable. 
     IF ( rans_tke_e  .OR.  use_sgs_for_particles  .OR.  wang_kernel  .OR.      &
-         collision_turbulence )  THEN
+         collision_turbulence  .OR.  nested_run )  THEN
 #if defined( __nopointer )
        ALLOCATE( diss(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
@@ -858 +871 @@
 #else
        ALLOCATE( diss_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
-       IF ( rans_tke_e )  THEN
+       IF ( rans_tke_e  .OR.  nested_run )  THEN
           ALLOCATE( diss_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
           ALLOCATE( diss_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
@@ -871 +884 @@
 
     IF ( rans_tke_e  .OR.  use_sgs_for_particles  .OR.     &
-         wang_kernel  .OR.  collision_turbulence )  THEN
+         wang_kernel  .OR.  collision_turbulence  .OR.  nested_run )  THEN
        diss => diss_1
-       IF ( rans_tke_e )  THEN
+       IF ( rans_tke_e  .OR.  nested_run )  THEN
        diss_p => diss_2; tdiss_m => diss_3
        ENDIF