Changeset 2851 for palm/trunk/SOURCE/check_parameters.f90

Timestamp:

Mar 5, 2018 2:39:31 PM (7 years ago)

Author:

maronga

Message:

bugfix: calculation of output time levels for restart runs (parallel NetCDF)

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (2 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix: calculation of output time levels in case of restart runs (parallel
+! NetCDF)
+! 
+! 2836 2018-02-26 13:40:05Z Giersch
 ! dt_dopr_listing is set to the default value zero 
 ! 
@@ -3673 +3677 @@
        CALL check_dt_do( dt_do2d_yz, 'dt_do2d_yz' )
 
-       ntdim_3d(0)    = INT( ( end_time - skip_time_do3d ) / dt_do3d )
+       ntdim_3d(0)    = CEILING( ( end_time - skip_time_do3d ) / dt_do3d )
        IF ( do3d_at_begin ) ntdim_3d(0) = ntdim_3d(0) + 1
-       ntdim_3d(1)    = INT( ( end_time - skip_time_data_output_av )           &
+       ntdim_3d(1)    = CEILING( ( end_time - skip_time_data_output_av )           &
                              / dt_data_output_av )
-       ntdim_2d_xy(0) = INT( ( end_time - skip_time_do2d_xy ) / dt_do2d_xy )
-       ntdim_2d_xz(0) = INT( ( end_time - skip_time_do2d_xz ) / dt_do2d_xz )
-       ntdim_2d_yz(0) = INT( ( end_time - skip_time_do2d_yz ) / dt_do2d_yz )
+       ntdim_2d_xy(0) = CEILING( ( end_time - skip_time_do2d_xy ) / dt_do2d_xy )
+       ntdim_2d_xz(0) = CEILING( ( end_time - skip_time_do2d_xz ) / dt_do2d_xz )
+       ntdim_2d_yz(0) = CEILING( ( end_time - skip_time_do2d_yz ) / dt_do2d_yz )
        IF ( do2d_at_begin )  THEN
           ntdim_2d_xy(0) = ntdim_2d_xy(0) + 1