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

Timestamp:

Jul 27, 2007 9:09:17 AM (16 years ago)

Author:

raasch

Message:

preliminary version for coupled runs

File:

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

palm/trunk/SOURCE/palm.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Get coupling mode from environment variable
 !
 ! Former revisions:
@@ -43 +43 @@
     USE arrays_3d
     USE constants
+    USE control_parameters
     USE cpulog
     USE dvrp_variables
@@ -53 +54 @@
     USE spectrum
     USE statistics
-    USE control_parameters
 
     IMPLICIT NONE
@@ -63 +63 @@
     INTEGER           ::  i, run_description_header_i(80)
 
-    version = 'PALM 3.3'
+    version = 'PALM 3.3a'
 
 #if defined( __parallel )
@@ -75 +75 @@
 #endif
 
+#if defined( __mpi2 )
+!
+!-- Get information about the coupling mode from the environment variable
+!-- which has been set by the mpiexec command
+    CALL local_getenv( 'coupling_mode', 13, coupling_mode, i )
+    IF ( i == 0 )  coupling_mode = 'uncoupled'
+    IF ( coupling_mode == 'ocean_to_atmosphere' )  coupling_char = '_O'
+#endif
+
 !
 !-- Initialize measuring of the CPU-time remaining to the run
@@ -104 +113 @@
 !
 !-- Open a file for debug output
-    OPEN( 9, FILE='DEBUG'//myid_char, FORM='FORMATTED' )
+    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 )
+#endif
 
 !
@@ -167 +184 @@
     CALL cpu_log( log_point(4), 'last actions', 'stop' )
 
+#if defined( __mpi2 )
+!
+!-- Test exchange via intercommunicator
+    IF ( coupling_mode == 'atmosphere_to_ocean' )  THEN
+       i = 12345 + myid
+       CALL MPI_SEND( i, 1, MPI_INTEGER, myid, 11, comm_inter, ierr )
+    ELSEIF ( coupling_mode == 'ocean_to_atmosphere' )  THEN
+       CALL MPI_RECV( i, 1, MPI_INTEGER, myid, 11, comm_inter, status, ierr )
+       PRINT*, '### myid: ', myid, '   received from atmosphere:  i = ', i
+    ENDIF
+#endif
+
 !
 !-- Take final CPU-time for CPU-time analysis