Changeset 2365 for palm/trunk/SOURCE/init_coupling.f90

Timestamp:

Aug 21, 2017 2:59:59 PM (7 years ago)

Author:

kanani

Message:

Vertical nesting implemented (SadiqHuq?)

File:

• palm/trunk/SOURCE/init_coupling.f90 (modified) (6 diffs)

palm/trunk/SOURCE/init_coupling.f90

@@ -25 +25 @@
 ! ------------------
 ! $Id$
+! Vertical nesting implemented (SadiqHuq)
+! 
+! 2298 2017-06-29 09:28:18Z raasch
 ! MPI2 coupling removed
 ! 
@@ -67 +70 @@
     USE pegrid
 
+    USE vertical_nesting_mod
+
     IMPLICIT NONE
 
@@ -96 +101 @@
        IF ( TRIM( coupling_mode ) == 'coupled_run' )  THEN
           i = 1
+       ELSEIF ( TRIM( coupling_mode ) == 'vnested_twi' )  THEN
+          i = 9
        ELSE
           i = 0
@@ -111 +118 @@
     i = bc_data(0)
 
-    IF ( i == 0 ) THEN
+    IF ( i == 0 )  THEN
        coupling_mode = 'uncoupled'
 !
@@ -121 +128 @@
           CLOSE ( 90 )
        ENDIF
+    ELSEIF ( i == 9 )  THEN
+
+!
+!--    Set a flag to identify runs with vertical nesting
+       vnested = .TRUE.
+       
+       comm_inter = MPI_COMM_WORLD
+       
+!
+!--    Split the total available PE's into two groups
+!--    numprocs for Coarse and Fine Grid are specified via mrun argument -N 
+       IF ( myid < bc_data(1) )  THEN
+          inter_color     = 0
+          numprocs        = bc_data(1)
+          coupling_mode   = 'vnested_crse'
+       ELSE
+          inter_color     = 1
+          numprocs        = bc_data(2)
+          coupling_mode   = 'vnested_fine'
+       ENDIF
+       
+       CALL MPI_COMM_SPLIT( MPI_COMM_WORLD, inter_color, 0, comm_palm, ierr )
+       comm2d = comm_palm
+       
+       OPEN( 90, FILE='VNESTING_PORT_OPENED', FORM='FORMATTED' )
+       WRITE ( 90, '(''TRUE'')' )
+       CLOSE ( 90 )
+      
     ELSE
        comm_inter = MPI_COMM_WORLD
@@ -154 +189 @@
     ENDIF
 
+    IF (  TRIM( coupling_mode ) == 'vnested_fine' )  THEN
+!
+!-- Set file extension for vertical nesting
+       coupling_char = '_N'
+    ENDIF
+
  END SUBROUTINE init_coupling