Changeset 3932 for palm/trunk/SOURCE/pmc_interface_mod.f90

Timestamp:

Apr 24, 2019 5:31:34 PM (4 years ago)

Author:

suehring

Message:

Add missing if statements for call of pmc_set_dataarray_name for TKE and dissipation; minor bugfix for nesting of chemical species

File:

• palm/trunk/SOURCE/pmc_interface_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Add missing if statements for call of pmc_set_dataarray_name for TKE and 
+! dissipation.
 ! 
 ! Former revisions:
@@ -1341 +1342 @@
 !--    design of model coupler, however, data array names must be already 
 !--    available at this point.
-       CALL pmc_set_dataarray_name( 'coarse', 'e'  ,'fine', 'e',  ierr )
+       IF ( (        rans_mode_parent  .AND.         rans_mode )  .OR.         &
+            (  .NOT. rans_mode_parent  .AND.  .NOT.  rans_mode  .AND.          &
+               .NOT. constant_diffusion ) )  THEN
+          CALL pmc_set_dataarray_name( 'coarse', 'e'  ,'fine', 'e',  ierr )
+       ENDIF
 !
 !--    Nesting of dissipation rate only if both parent and child are in RANS
-!--    mode and TKE-epsilo closure is applied. Please so also comment for TKE
+!--    mode and TKE-epsilo closure is applied. Please see also comment for TKE
 !--    above.
-       CALL pmc_set_dataarray_name( 'coarse', 'diss'  ,'fine', 'diss',  ierr )
+       IF ( rans_mode_parent  .AND.  rans_mode  .AND.  rans_tke_e )  THEN
+          CALL pmc_set_dataarray_name( 'coarse', 'diss'  ,'fine', 'diss',  ierr )
+       ENDIF
 
        IF ( .NOT. neutral )  THEN
@@ -1511 +1518 @@
 
        DO  WHILE ( pmc_c_getnextarray( myname ) )
+!
 !--       Note that cg%nz is not the original nz of parent, but the highest
 !--       parent-grid level needed for nesting.