Changeset 108 for palm/trunk/SOURCE/local_stop.f90

Timestamp:

Aug 24, 2007 3:10:38 PM (17 years ago)

Author:

letzel

Message:

• Improved coupler: evaporation - salinity-flux coupling for humidity = .T.,

avoid MPI hangs when coupled runs terminate, add DOC/app/chapter_3.8;

• Optional calculation of km and kh from initial TKE e_init;
• Default initialization of km,kh = 0.00001 for ocean = .T.;
• Allow data_output_pr= q, wq, w"q", w*q* for humidity = .T.;
• Bugfix: Rayleigh damping for ocean fixed.

File:

• palm/trunk/SOURCE/local_stop.f90 (modified) (2 diffs)

palm/trunk/SOURCE/local_stop.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! modifications to terminate coupled runs
 ! 
 !
@@ -24 +25 @@
 
     USE pegrid
+    USE control_parameters
 
 #if defined( __parallel )
-    CALL MPI_FINALIZE( ierr )
-#endif
+    IF ( coupling_mode == 'uncoupled' )  THEN
+       CALL MPI_FINALIZE( ierr )
+    ELSE
+
+       SELECT CASE ( terminate_coupled_remote )
+
+          CASE ( 0 )
+             IF ( myid == 0 )  THEN
+                PRINT*, '+++ local_stop:'
+                PRINT*, '    local model "', TRIM( coupling_mode ), &
+                     '" stops now'
+             ENDIF
+!
+!--          Inform the remote model of the termination and its reason, provided
+!--          the remote model has not already been informed of another 
+!--          termination reason (terminate_coupled > 0) before.
+             IF ( terminate_coupled == 0 )  THEN
+                terminate_coupled = 1
+                CALL MPI_SENDRECV( &
+                     terminate_coupled,        1, MPI_INTEGER, myid,  0, &
+                     terminate_coupled_remote, 1, MPI_INTEGER, myid,  0, &
+                     comm_inter, status, ierr )
+             ENDIF
+             CALL MPI_FINALIZE( ierr )
+
+          CASE ( 1 )
+             IF ( myid == 0 )  THEN
+                PRINT*, '+++ local_stop:'
+                PRINT*, '    remote model "', TRIM( coupling_mode_remote ), &
+                     '" stopped'
+             ENDIF
+             CALL MPI_FINALIZE( ierr )
+
+          CASE ( 2 )
+             IF ( myid == 0 )  THEN
+                PRINT*, '+++ local_stop:'
+                PRINT*, '    remote model "', TRIM( coupling_mode_remote ), &
+                     '" terminated'
+                PRINT*, '    with stop_dt = .T.'
+             ENDIF
+             stop_dt = .TRUE.
+
+          CASE ( 3 )
+             IF ( myid == 0 )  THEN
+                PRINT*, '+++ local_stop:'
+                PRINT*, '    remote model "', TRIM( coupling_mode_remote ), &
+                     '" terminated'
+                PRINT*, '    with terminate_run = .T. (CPU-time limit)'
+             ENDIF
+             terminate_run = .TRUE.
+
+          CASE ( 4 )
+             IF ( myid == 0 )  THEN
+                PRINT*, '+++ local_stop:'
+                PRINT*, '    remote model "', TRIM( coupling_mode_remote ), &
+                     '" terminated'
+                PRINT*, '    with terminate_run = .T. (restart)'
+             ENDIF
+             terminate_run = .TRUE.
+             time_restart = time_restart + dt_restart
+
+          CASE ( 5 )
+             IF ( myid == 0 )  THEN
+                PRINT*, '+++ local_stop:'
+                PRINT*, '    remote model "', TRIM( coupling_mode_remote ), &
+                     '" terminated'
+                PRINT*, '    with terminate_run = .T. (single restart)'
+             ENDIF
+             terminate_run = .TRUE.
+             time_restart = 9999999.9
+
+       END SELECT
+
+    ENDIF
+
+#else
 
     STOP
 
+#endif
+
  END SUBROUTINE local_stop