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

palm/trunk/SOURCE/check_parameters.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: check data output intervals to be /= 0.0 in case of parallel NetCDF4
 ! 
 ! Former revisions:
@@ -3835 +3835 @@
     IF ( netcdf_data_format > 4 )  THEN
 
+!
+!--    Check if any of the follwoing data output interval is 0.0s, which is
+!--    not allowed for parallel output.
+       CALL check_dt_do( dt_do3d,    'dt_do3d'    )
+       CALL check_dt_do( dt_do2d_xy, 'dt_do2d_xy' )
+       CALL check_dt_do( dt_do2d_xz, 'dt_do2d_xz' )
+       CALL check_dt_do( dt_do2d_yz, 'dt_do2d_yz' )
+
        ntdim_3d(0)    = INT( ( end_time - skip_time_do3d ) / dt_do3d )
        IF ( do3d_at_begin ) ntdim_3d(0) = ntdim_3d(0) + 1
@@ -3850 +3858 @@
        ntdim_2d_xz(1) = ntdim_3d(1)
        ntdim_2d_yz(1) = ntdim_3d(1)
-
+             write(9,*) "ntdim_3d(av)=",ntdim_3d(0)
+             CALL local_flush(9)
     ENDIF
 
@@ -4333 +4342 @@
     CALL user_check_parameters
 
+ CONTAINS
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Check the length of data output intervals. In case of parallel NetCDF output
+!> the time levels of the output files need to be fixed. Therefore setting the
+!> output interval to 0.0s (usually used to output each timestep) is not 
+!> possible as long as a non-fixed timestep is used.
+!------------------------------------------------------------------------------!
+
+    SUBROUTINE check_dt_do( dt_do, dt_do_name )
+
+       IMPLICIT NONE
+
+       CHARACTER (LEN=*), INTENT (IN) :: dt_do_name !< parin variable name
+
+       REAL(wp), INTENT (INOUT)       :: dt_do      !< data output interval
+
+       IF ( dt_do == 0.0_wp )  THEN
+          IF ( dt_fixed )  THEN
+             WRITE( message_string, '(A,F9.4,A)' )  'Output at every '  //  &
+                    'timestep is desired (' // dt_do_name // ' = 0.0).&'//  &
+                    'Setting the output interval to the fixed timestep '//  &
+                    'dt = ', dt, 's.'
+             CALL message( 'check_parameters', 'PA0060', 0, 0, 0, 6, 0 )
+             dt_do = dt
+          ELSE
+             message_string = dt_do_name // ' = 0.0 while using a ' //      &
+                              'variable timestep and parallel netCDF4 ' //  &
+                              'is not allowed.'
+             CALL message( 'check_parameters', 'PA0081', 1, 2, 0, 6, 0 )
+          ENDIF
+       ENDIF
+
+    END SUBROUTINE check_dt_do
 
  END SUBROUTINE check_parameters