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

Timestamp:

Feb 5, 2016 1:06:51 PM (8 years ago)

Author:

gronemeier

Message:

Bugfix:calculation of time levels for parallel NetCDF output

File:

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

palm/trunk/SOURCE/data_output_3d.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! Bugfix: test if time axis limit exceeds moved to point after call of check_open
 ! 
 ! Former revisions:
@@ -176 +176 @@
 !-- 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')
@@ -209 +195 @@
        CALL check_open( 106+av*10 )
     ENDIF
+
+!
+!-- 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 )
+          CALL cpu_log( log_point(14), 'data_output_3d', 'stop' )
+          RETURN
+       ENDIF
+    ENDIF
+    WRITE(9,*) '___hier4'
+    CALL local_flush(9)
 
 !