Changeset 1762 for palm/trunk/SOURCE/init_pegrid.f90

Timestamp:

Feb 25, 2016 12:31:13 PM (8 years ago)

Author:

hellstea

Message:

Introduction of nested domain system

File:

• palm/trunk/SOURCE/init_pegrid.f90 (modified) (11 diffs)

palm/trunk/SOURCE/init_pegrid.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! Introduction of nested domain feature
 ! 
 ! Former revisions:
@@ -129 +129 @@
                io_blocks, io_group, maximum_grid_level,                        &
                maximum_parallel_io_streams, message_string,                    &
-               mg_switch_to_pe0_level, momentum_advec, neutral, psolver,       &
+               mg_switch_to_pe0_level, momentum_advec, nest_bound_l,           &
+               nest_bound_n, nest_bound_r, nest_bound_s, neutral, psolver,     &
                outflow_l, outflow_n, outflow_r, outflow_s, recycling_width,    &
                scalar_advec, subdomain_size 
@@ -147 +148 @@
     USE pegrid
  
+#if defined( PMC_ACTIVE )
+    USE pmc_interface,                                                         &
+        ONLY:  cpl_npex,cpl_npey
+#endif
+
     USE transpose_indices,                                                     &
         ONLY:  nxl_y, nxl_yd, nxl_z, nxr_y, nxr_yd, nxr_z, nyn_x, nyn_z, nys_x,&
@@ -206 +212 @@
     CALL location_message( 'creating virtual PE grids + MPI derived data types', &
                            .FALSE. )
+#if defined( PMC_ACTIVE )
+!
+!-- In case of nested-domain runs, the processor grid is explicitly given
+!-- by the user in the nestpar-NAMELIST
+    pdims(1) = cpl_npex
+    pdims(2) = cpl_npey
+#else
 !
 !-- Determine the processor topology or check it, if prescribed by the user
@@ -243 +256 @@
 
     ENDIF
+#endif
 
 !
@@ -1058 +1072 @@
 #if defined( __parallel ) && ! defined ( __check )
 !
-!-- Setting of flags for inflow/outflow conditions in case of non-cyclic
+!-- Setting of flags for inflow/outflow/nesting conditions in case of non-cyclic
 !-- horizontal boundary conditions.
     IF ( pleft == MPI_PROC_NULL )  THEN
@@ -1065 +1079 @@
        ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
           outflow_l = .TRUE.
+#if defined( PMC_ACTIVE )
+       ELSEIF ( bc_lr == 'nested' )  THEN
+          nest_bound_l = .TRUE.
+#endif
        ENDIF
     ENDIF
@@ -1073 +1091 @@
        ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
           inflow_r  = .TRUE.
+#if defined( PMC_ACTIVE )
+       ELSEIF ( bc_lr == 'nested' )  THEN
+          nest_bound_r = .TRUE.
+#endif
        ENDIF
     ENDIF
@@ -1081 +1103 @@
        ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
           inflow_s  = .TRUE.
+#if defined( PMC_ACTIVE )
+       ELSEIF ( bc_ns == 'nested' )  THEN
+          nest_bound_s = .TRUE.
+#endif
        ENDIF
     ENDIF
@@ -1089 +1115 @@
        ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
           outflow_n = .TRUE.
+#if defined( PMC_ACTIVE )
+       ELSEIF ( bc_ns == 'nested' )  THEN
+          nest_bound_n = .TRUE.
+#endif
        ENDIF
     ENDIF
@@ -1139 +1169 @@
 !-- At the inflow or outflow, u or v, respectively, have to be calculated for
 !-- one more grid point.
-    IF ( inflow_l .OR. outflow_l )  THEN
+    IF ( inflow_l .OR. outflow_l .OR. nest_bound_l )  THEN
        nxlu = nxl + 1
     ELSE
        nxlu = nxl
     ENDIF
-    IF ( inflow_s .OR. outflow_s )  THEN
+    IF ( inflow_s .OR. outflow_s .OR. nest_bound_s )  THEN
        nysv = nys + 1
     ELSE