Changeset 691 for palm/trunk/SOURCE

Timestamp:

Mar 4, 2011 8:45:30 AM (13 years ago)

Author:

maronga

Message:

Bugfix for precursor atmosphere/ocean runs, re-adjustments for lcxt4

Location:

palm/trunk/SOURCE

Files:

• data_output_2d.f90 (modified) (4 diffs)
• data_output_3d.f90 (modified) (6 diffs)
• init_coupling.f90 (modified) (2 diffs)

palm/trunk/SOURCE/data_output_2d.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+! Replaced simulated_time by time_since_reference_point
 !
 ! Former revisions:
@@ -727 +728 @@
 #if defined( __netcdf )
                          IF ( netcdf_output  .AND.  myid == 0 )  THEN
-                            WRITE ( 21 )  simulated_time, &
+                            WRITE ( 21 )  time_since_reference_point, &
                                           do2d_xy_time_count(av), av
                          ENDIF
@@ -1023 +1024 @@
 #if defined( __netcdf )
                          IF ( netcdf_output  .AND.  myid == 0 )  THEN
-                            WRITE ( 22 )  simulated_time, &
+                            WRITE ( 22 )  time_since_reference_point, &
                                           do2d_xz_time_count(av), av
                          ENDIF
@@ -1324 +1325 @@
 #if defined( __netcdf )
                          IF ( netcdf_output  .AND.  myid == 0 )  THEN
-                            WRITE ( 23 )  simulated_time, &
+                            WRITE ( 23 )  time_since_reference_point, &
                                           do2d_yz_time_count(av), av
                          ENDIF

palm/trunk/SOURCE/data_output_3d.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+! Replaced simulated_time by time_since_reference_point
 !
 ! Former revisions:
@@ -242 +243 @@
                 to_be_resorted => q_av
              ENDIF
-             
+
           CASE ( 'ql' )
              IF ( av == 0 )  THEN
@@ -291 +292 @@
                 to_be_resorted => rho_av
              ENDIF
-             
+
           CASE ( 's' )
              IF ( av == 0 )  THEN
@@ -298 +299 @@
                 to_be_resorted => s_av
              ENDIF
-             
+
           CASE ( 'sa' )
              IF ( av == 0 )  THEN
@@ -305 +306 @@
                 to_be_resorted => sa_av
              ENDIF
-             
+
           CASE ( 'u' )
              IF ( av == 0 )  THEN
@@ -408 +409 @@
 !--             combine_plot_fields
                 IF ( myid == 0 )  THEN
-                   WRITE ( 30 )  simulated_time, do3d_time_count(av), av
+                   WRITE ( 30 )  time_since_reference_point,                   &
+                                 do3d_time_count(av), av
                 ENDIF
                 WRITE ( 30 )  nxlg, nxrg, nysg, nyng, nzb, nz_do3d

palm/trunk/SOURCE/init_coupling.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-!
+! Bugfix: combine_plot_fields did not work with data_output_2d_on_each_pe = .T.
+! for precursor ocean runs
 !
 ! Former revisions:
@@ -96 +97 @@
     IF ( i == 0 ) THEN
        coupling_mode = 'uncoupled'
+!
+!--    In case of a precursor ocean run, an additional flag file is created.
+!--    This is necessary for data_output_2d_on_each_pe = .T.
+       IF ( bc_data(3) == 1 )  THEN
+          OPEN( 90, FILE='PRECURSOR_OCEAN', FORM='FORMATTED' )
+          WRITE ( 90, '(''TRUE'')' )
+          CLOSE ( 90 )
+       ENDIF
     ELSE
        comm_inter = MPI_COMM_WORLD