Changeset 3592

Timestamp:

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

Author:

suehring

Message:

Determine number of coupled arrays dynamically

Location:

palm/trunk/SOURCE

Files:

• pmc_general_mod.f90 (modified) (2 diffs)
• pmc_interface_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/pmc_general_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Determine number of coupled arrays dynamically.
+! 
+! 3241 2018-09-12 15:02:00Z raasch
 ! unused variable removed
 ! 
@@ -99 +102 @@
     SAVE
 
+    INTEGER(iwp), PUBLIC            :: pmc_max_array          !< max # of arrays which can be coupled - will be determined dynamically in pmc_interface
+    
     INTEGER(iwp), PARAMETER, PUBLIC :: da_desclen       =  8  !<
     INTEGER(iwp), PARAMETER, PUBLIC :: da_namelen       = 16  !<
     INTEGER(iwp), PARAMETER, PUBLIC :: pmc_da_name_err  = 10  !<
-    INTEGER(iwp), PARAMETER, PUBLIC :: pmc_max_array    = 32  !< max # of arrays which can be coupled
     INTEGER(iwp), PARAMETER, PUBLIC :: pmc_max_models   = 64  !<
     INTEGER(iwp), PARAMETER, PUBLIC :: pmc_status_ok    =  0  !<

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,                 &