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

Timestamp:

Dec 3, 2018 12:38:40 PM (6 years ago)

Author:

suehring

Message:

Determine number of coupled arrays dynamically

File:

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

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Number of coupled arrays is determined dynamically (instead of a fixed value
+! of 32)
+! 
+! 3524 2018-11-14 13:36:44Z raasch
 ! declaration statements rearranged to avoid compile time errors
 ! 
@@ -804 +808 @@
     CALL pmci_setup_coordinates
 !
+!-- Determine the number of coupled arrays
+    CALL pmci_num_arrays
+!
 !-- Initialize the child (must be called before pmc_setup_parent)
     CALL pmci_setup_child
@@ -3660 +3667 @@
 #endif
  END SUBROUTINE pmci_setup_coordinates
+ 
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> In this subroutine the number of coupled arrays is determined. 
+!------------------------------------------------------------------------------! 
+  SUBROUTINE pmci_num_arrays  
+               
+#if defined( __parallel ) 
+    USE pmc_general,                                                           &
+        ONLY:  pmc_max_array
+
+    IMPLICIT NONE
+!
+!-- The number of coupled arrays depends on the model settings. At least
+!-- 5 arrays need to be coupled (u, v, w, e, diss).  Please note, actually
+!-- e and diss (TKE and dissipation rate) are only required if RANS-RANS
+!-- nesting is applied, but memory is allocated nevertheless. This is because
+!-- the information whether they are needed or not is retrieved at a later
+!-- point in time. In case e and diss are not needed, they are also not 
+!-- exchanged between parent and child. 
+    pmc_max_array = 5
+!
+!-- pt
+    IF ( .NOT. neutral )  pmc_max_array = pmc_max_array + 1
+    
+    IF ( humidity )  THEN
+!
+!--    q
+       pmc_max_array = pmc_max_array + 1
+!
+!--    qc, nc
+       IF ( bulk_cloud_model  .AND.  microphysics_morrison )                   &
+          pmc_max_array = pmc_max_array + 2
+!
+!--    qr, nr
+       IF ( bulk_cloud_model  .AND.  microphysics_seifert )                    &
+          pmc_max_array = pmc_max_array + 2
+    ENDIF
+!
+!-- s
+    IF ( passive_scalar )  pmc_max_array = pmc_max_array + 1
+!
+!-- nr_part, part_adr
+    IF ( particle_advection )  pmc_max_array = pmc_max_array + 2
+!
+!-- Chemistry, depends on number of species
+    IF ( air_chemistry  .AND.  nest_chemistry )                                &
+       pmc_max_array = pmc_max_array + nspec
+
+
+#endif
+ END SUBROUTINE pmci_num_arrays
 
 
@@ -4081 +4141 @@
        ENDIF
 
-       IF ( air_chemistry  .AND.  nest_chemistry)  THEN
+       IF ( air_chemistry  .AND.  nest_chemistry )  THEN
           DO  n = 1, nspec
              CALL pmci_interp_tril_all ( chem_species(n)%conc,                 &