Changeset 206 for palm/trunk/SOURCE/palm.f90

Timestamp:

Oct 13, 2008 2:59:11 PM (13 years ago)

Author:

raasch

Message:

ocean-atmosphere coupling realized with MPI-1, adjustments in mrun, mbuild, subjob for lcxt4

File:

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

palm/trunk/SOURCE/palm.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Initialization of coupled runs modified for MPI-1 and moved to external
+! subroutine init_coupling
 !
 ! Former revisions:
@@ -77 +78 @@
     CALL MPI_INIT( ierr )
     CALL MPI_COMM_SIZE( MPI_COMM_WORLD, numprocs, ierr )
+    CALL MPI_COMM_RANK( MPI_COMM_WORLD, myid, ierr )
     comm_palm = MPI_COMM_WORLD
     comm2d    = MPI_COMM_WORLD
-#endif
-
-#if defined( __mpi2 )
-!
-!-- Get information about the coupling mode from the environment variable
-!-- which has been set by the mpiexec command.
-!-- This method is currently not used because the mpiexec command is not
-!-- available on some machines
-!    CALL local_getenv( 'coupling_mode', 13, coupling_mode, i )
-!    IF ( i == 0 )  coupling_mode = 'uncoupled'
-!    IF ( coupling_mode == 'ocean_to_atmosphere' )  coupling_char = '_O'
-
-!
-!-- Get information about the coupling mode from standard input (PE0 only) and
-!-- distribute it to the other PEs
-    CALL MPI_COMM_RANK( MPI_COMM_WORLD, myid, ierr )
-    IF ( myid == 0 )  THEN
-       READ (*,*,ERR=10,END=10)  coupling_mode
-10     IF ( TRIM( coupling_mode ) == 'atmosphere_to_ocean' )  THEN
-          i = 1
-       ELSEIF ( TRIM( coupling_mode ) ==  'ocean_to_atmosphere' )  THEN
-          i = 2
-       ELSE
-          i = 0
-       ENDIF
-    ENDIF
-    CALL MPI_BCAST( i, 1, MPI_INTEGER, 0, MPI_COMM_WORLD, ierr )
-    IF ( i == 0 )  THEN
-       coupling_mode = 'uncoupled'
-    ELSEIF ( i == 1 )  THEN
-       coupling_mode = 'atmosphere_to_ocean'
-    ELSEIF ( i == 2 )  THEN
-       coupling_mode = 'ocean_to_atmosphere'
-    ENDIF
-    IF ( coupling_mode == 'ocean_to_atmosphere' )  coupling_char = '_O'
+
+!
+!-- Initialize PE topology in case of coupled runs
+    CALL init_coupling
 #endif
 
@@ -124 +95 @@
     CALL cpu_log( log_point(1), 'total', 'start' )
     CALL cpu_log( log_point(2), 'initialisation', 'start' )
+
+!
+!-- Open a file for debug output
+    WRITE (myid_char,'(''_'',I4.4)')  myid
+    OPEN( 9, FILE='DEBUG'//TRIM( coupling_char )//myid_char, FORM='FORMATTED' )
 
 !
@@ -132 +108 @@
 #if defined( __parallel )
     CALL MPI_COMM_RANK( comm_palm, myid, ierr )
-#endif
-
-!
-!-- Open a file for debug output
-    WRITE (myid_char,'(''_'',I4.4)')  myid
-    OPEN( 9, FILE='DEBUG'//TRIM( coupling_char )//myid_char, FORM='FORMATTED' )
-
-#if defined( __mpi2 )
 !
 !-- TEST OUTPUT (TO BE REMOVED)
     WRITE(9,*) '*** coupling_mode = "', TRIM( coupling_mode ), '"'
     CALL LOCAL_FLUSH( 9 )
-    print*, '*** PE', myid, '  ', TRIM( coupling_mode )
+    PRINT*, '*** PE', myid, ' Global target PE:', target_id, &
+            TRIM( coupling_mode )
 #endif
 
@@ -220 +189 @@
 #if defined( __mpi2 )
 !
-!-- Test exchange via intercommunicator
+!-- Test exchange via intercommunicator in case of a MPI-2 coupling
     IF ( coupling_mode == 'atmosphere_to_ocean' )  THEN
        i = 12345 + myid
@@ -240 +209 @@
 
  END PROGRAM palm
-