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

Timestamp:

Apr 16, 2009 12:07:26 PM (15 years ago)

Author:

raasch

Message:

changes for coupling with independent precursor runs; z_i calculation with Sullivan criterion

File:

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

palm/trunk/SOURCE/init_coupling.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Coupling with independent precursor runs.
 !
 ! Former revisions:
@@ -27 +27 @@
 !-- Local variables
     INTEGER               ::  i, inter_color
-    INTEGER, DIMENSION(:) ::  bc_data(0:2) = 0
+    INTEGER, DIMENSION(:) ::  bc_data(0:3) = 0
 
 !
@@ -48 +48 @@
        READ (*,*,ERR=10,END=10)  coupling_mode, bc_data(1), bc_data(2)
 10     CONTINUE
-
 #if defined( __mpi2 )
        IF ( TRIM( coupling_mode ) == 'atmosphere_to_ocean' )  THEN
@@ -64 +63 @@
        ENDIF
 #endif
-    bc_data(0) = i
+       bc_data(0) = i
+
+!
+!--    Check if '_O' has to be used as file extension in an uncoupled ocean
+!--    run. This is required, if this run shall be continued as a coupled run.
+       IF ( TRIM( coupling_mode ) == 'precursor_ocean' )  bc_data(3) = 1
+
     ENDIF
 
-    CALL MPI_BCAST( bc_data(0), 3, MPI_INTEGER, 0, MPI_COMM_WORLD, ierr )
+    CALL MPI_BCAST( bc_data(0), 4, MPI_INTEGER, 0, MPI_COMM_WORLD, ierr )
     i = bc_data(0)
 
@@ -114 +119 @@
 #endif
 
-    IF ( coupling_mode == 'ocean_to_atmosphere' )  coupling_char = '_O'
+!
+!-- In case of a precursor ocean run (followed by a coupled run), or a
+!-- coupled atmosphere-ocean run, set the file extension for the ocean files
+    IF ( TRIM( coupling_mode ) == 'ocean_to_atmosphere' .OR. bc_data(3) == 1 ) &
+    THEN
+       coupling_char = '_O'
+    ENDIF
 
  END SUBROUTINE init_coupling