Changeset 291 for palm/trunk/SOURCE/check_parameters.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/check_parameters.f90 (modified) (6 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! Coupling with independent precursor runs.
 ! Check particle_color, particle_dvrpsize, color_interval, dvrpsize_interval
 ! Bugfix: pressure included for profile output
@@ -115 +116 @@
          position, prec
     LOGICAL ::  found, ldum
-    REAL    ::  gradient, maxn, maxp, remote = 0.0
+    REAL    ::  gradient, maxn, maxp, remote = 0.0, &
+                simulation_time_since_reference
 
 !
@@ -188 +190 @@
        ENDIF
 
-       CALL MPI_SEND( end_time, 1, MPI_REAL, target_id, 14, comm_inter, ierr )
+       simulation_time_since_reference = end_time - coupling_start_time
+       CALL MPI_SEND( simulation_time_since_reference, 1, MPI_REAL, target_id, &
+                      14, comm_inter, ierr )
        CALL MPI_RECV( remote, 1, MPI_REAL, target_id, 14, comm_inter, &
                       status, ierr )
-       IF ( end_time /= remote )  THEN
+       IF ( simulation_time_since_reference /= remote )  THEN
           WRITE( message_string, * ) 'coupling mode "', TRIM( coupling_mode ), &
-                 '": end_time = ', end_time, '& is not equal to ',             &
-                 'end_time_remote = ', remote
+                 '": simulation_time_since_reference = ',                      &
+                 simulation_time_since_reference, '& is not equal to ',        &
+                 'simulation_time_since_reference_remote = ', remote
           CALL message( 'check_parameters', 'PA0008', 1, 2, 0, 6, 0 )
        ENDIF
@@ -383 +388 @@
 
 !
+!-- Check whether an (uncoupled) atmospheric run has been declared as an
+!-- ocean run (this setting is done via mrun-option -y)
+    IF ( TRIM( coupling_mode ) == 'uncoupled'  .AND.  &
+         TRIM( coupling_char ) == '_O' )  THEN
+       message_string = 'ocean = .F. does not allow coupling_char = "' // &
+                        TRIM( coupling_char ) // '" set by mrun-option "-y"'
+       CALL message( 'check_parameters', 'PA0317', 1, 2, 0, 6, 0 )
+    ENDIF
+
+!
 !-- Check whether there are any illegal values
 !-- Pressure solver:
@@ -1027 +1042 @@
 
 !
+!-- Store reference time for coupled runs and change the coupling flag,
+!-- if ...
+    IF ( simulated_time == 0.0 )  THEN
+       IF ( coupling_start_time == 0.0 )  THEN
+          time_since_reference_point = 0.0
+       ELSEIF ( time_since_reference_point < 0.0 )  THEN
+          run_coupled = .FALSE.
+       ENDIF
+    ENDIF
+
+!
 !-- Set wind speed in the Galilei-transformed system
     IF ( galilei_transformation )  THEN
@@ -2813 +2839 @@
 !-- internal parameter for steering restart events)
     IF ( restart_time /= 9999999.9 )  THEN
-       IF ( restart_time > simulated_time )  time_restart = restart_time
+       IF ( restart_time > time_since_reference_point )  THEN
+          time_restart = restart_time
+       ENDIF
     ELSE
 !