Changeset 2769

Timestamp:

Jan 25, 2018 9:22:24 AM (7 years ago)

Author:

raasch

Message:

bugfix for calculating number of required output time levels in case of output at the beginning of a restart run

File:

• palm/trunk/SOURCE/netcdf_interface_mod.f90 (modified) (10 diffs)

palm/trunk/SOURCE/netcdf_interface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! bugfix for calculating number of required output time levels in case of output
+! at the beginning of a restart run
+! 
+! 2766 2018-01-22 17:17:47Z kanani
 ! Removed preprocessor directive __chem
 ! 
@@ -494 +498 @@
                data_output_pr, domask, dopr_n,                                 &
                dopr_time_count, dopts_time_count, dots_time_count,             &
-               do2d, do2d_xz_time_count, do3d,                                 &
+               do2d, do2d_at_begin, do2d_xz_time_count, do3d, do3d_at_begin,   &
                do2d_yz_time_count, dt_data_output_av, dt_do2d_xy, dt_do2d_xz,  &
                dt_do2d_yz, dt_do3d, mask_size, do2d_xy_time_count,             &
@@ -1810 +1814 @@
                                                 simulated_time_at_begin )  &
                               ) / dt_do3d )
+                IF ( do3d_at_begin )  ntdim_3d(0) = ntdim_3d(0) + 1
              ELSE
                 ntdim_3d(1) = do3d_time_count(1) + CEILING(                &
@@ -1838 +1843 @@
                                              simulated_time_at_begin )  &
                            ) / dt_do3d )
+                   IF ( do3d_at_begin )  ntdim_3d(0) = ntdim_3d(0) + 1
                 ELSE
                    ntdim_3d(1) = CEILING(                               &
@@ -2576 +2582 @@
                                                    simulated_time_at_begin )  &
                                  ) / dt_do2d_xy )
+                IF ( do2d_at_begin )  ntdim_2d_xy(0) = ntdim_2d_xy(0) + 1
              ELSE
                 ntdim_2d_xy(1) = do2d_xy_time_count(1) + CEILING(             &
@@ -2604 +2611 @@
                                              simulated_time_at_begin )  &
                            ) / dt_do2d_xy )
+                   IF ( do2d_at_begin )  ntdim_2d_xy(0) = ntdim_2d_xy(0) + 1
                 ELSE
                    ntdim_2d_xy(1) = CEILING(                            &
@@ -3226 +3234 @@
                                                    simulated_time_at_begin )  &
                                  ) / dt_do2d_xz )
+                IF ( do2d_at_begin )  ntdim_2d_xz(0) = ntdim_2d_xz(0) + 1
              ELSE
                 ntdim_2d_xz(1) = do2d_xz_time_count(1) + CEILING(             &
@@ -3254 +3263 @@
                                              simulated_time_at_begin )  &
                            ) / dt_do2d_xz )
+                   IF ( do2d_at_begin )  ntdim_2d_xz(0) = ntdim_2d_xz(0) + 1
                 ELSE
                    ntdim_2d_xz(1) = CEILING(                            &
@@ -3879 +3889 @@
                                                    simulated_time_at_begin )  &
                                  ) / dt_do2d_yz )
+                IF ( do2d_at_begin )  ntdim_2d_yz(0) = ntdim_2d_yz(0) + 1
              ELSE
                 ntdim_2d_yz(1) = do2d_yz_time_count(1) + CEILING(             &
@@ -3907 +3918 @@
                                              simulated_time_at_begin )  &
                            ) / dt_do2d_yz )
+                   IF ( do2d_at_begin )  ntdim_2d_yz(0) = ntdim_2d_yz(0) + 1
                 ELSE
                    ntdim_2d_yz(1) = CEILING(                            &