Changeset 1308 for palm/trunk/SOURCE/data_output_3d.f90

Timestamp:

Mar 13, 2014 2:58:42 PM (10 years ago)

Author:

fricke

Message:

Adjustments for parallel NetCDF output for lccrayh/lccrayb (Cray XC30 systems)

File:

• palm/trunk/SOURCE/data_output_3d.f90 (modified) (3 diffs)

palm/trunk/SOURCE/data_output_3d.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! Check, if the limit of the time dimension is exceeded for parallel output
+! To increase the performance for parallel output, the following is done:
+! - Update of time axis is only done by PE0
 !
 ! Former revisions:
@@ -144 +146 @@
 !-- Return, if nothing to output
     IF ( do3d_no(av) == 0 )  RETURN
+!
+!-- For parallel netcdf output the time axis must be limited. Return, if this
+!-- limit is exceeded. This could be the case, if the simulated time exceeds 
+!-- the given end time by the length of the given output interval.
+    IF ( netcdf_data_format > 4 )  THEN
+       IF ( do3d_time_count(av) + 1 > ntdim_3d(av) )  THEN
+          WRITE ( message_string, * ) 'Output of 3d data is not given at t=',  &
+                                      simulated_time, '&because the maximum ', & 
+                                      'number of output time levels is ',      &
+                                      'exceeded.'
+          CALL message( 'data_output_3d', 'PA0387', 0, 1, 0, 6, 0 )          
+          RETURN
+       ENDIF
+    ENDIF
 
     CALL cpu_log (log_point(14),'data_output_3d','start')
@@ -174 +190 @@
 !
 !-- Update the netCDF time axis
+!-- In case of parallel output, this is only done by PE0 to increase the
+!-- performance.
 #if defined( __netcdf )
-    IF ( myid == 0  .OR.  netcdf_data_format > 4 )  THEN
-       do3d_time_count(av) = do3d_time_count(av) + 1
+    do3d_time_count(av) = do3d_time_count(av) + 1
+    IF ( myid == 0 )  THEN
        IF ( netcdf_output )  THEN
           nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_time_3d(av), &