Ignore:
Timestamp:
Jan 25, 2018 9:22:24 AM (3 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/netcdf_interface_mod.f90

    r2766 r2769  
    2525! -----------------
    2626! $Id$
     27! bugfix for calculating number of required output time levels in case of output
     28! at the beginning of a restart run
     29!
     30! 2766 2018-01-22 17:17:47Z kanani
    2731! Removed preprocessor directive __chem
    2832!
     
    494498               data_output_pr, domask, dopr_n,                                 &
    495499               dopr_time_count, dopts_time_count, dots_time_count,             &
    496                do2d, do2d_xz_time_count, do3d,                                 &
     500               do2d, do2d_at_begin, do2d_xz_time_count, do3d, do3d_at_begin,   &
    497501               do2d_yz_time_count, dt_data_output_av, dt_do2d_xy, dt_do2d_xz,  &
    498502               dt_do2d_yz, dt_do3d, mask_size, do2d_xy_time_count,             &
     
    18101814                                                simulated_time_at_begin )  &
    18111815                              ) / dt_do3d )
     1816                IF ( do3d_at_begin )  ntdim_3d(0) = ntdim_3d(0) + 1
    18121817             ELSE
    18131818                ntdim_3d(1) = do3d_time_count(1) + CEILING(                &
     
    18381843                                             simulated_time_at_begin )  &
    18391844                           ) / dt_do3d )
     1845                   IF ( do3d_at_begin )  ntdim_3d(0) = ntdim_3d(0) + 1
    18401846                ELSE
    18411847                   ntdim_3d(1) = CEILING(                               &
     
    25762582                                                   simulated_time_at_begin )  &
    25772583                                 ) / dt_do2d_xy )
     2584                IF ( do2d_at_begin )  ntdim_2d_xy(0) = ntdim_2d_xy(0) + 1
    25782585             ELSE
    25792586                ntdim_2d_xy(1) = do2d_xy_time_count(1) + CEILING(             &
     
    26042611                                             simulated_time_at_begin )  &
    26052612                           ) / dt_do2d_xy )
     2613                   IF ( do2d_at_begin )  ntdim_2d_xy(0) = ntdim_2d_xy(0) + 1
    26062614                ELSE
    26072615                   ntdim_2d_xy(1) = CEILING(                            &
     
    32263234                                                   simulated_time_at_begin )  &
    32273235                                 ) / dt_do2d_xz )
     3236                IF ( do2d_at_begin )  ntdim_2d_xz(0) = ntdim_2d_xz(0) + 1
    32283237             ELSE
    32293238                ntdim_2d_xz(1) = do2d_xz_time_count(1) + CEILING(             &
     
    32543263                                             simulated_time_at_begin )  &
    32553264                           ) / dt_do2d_xz )
     3265                   IF ( do2d_at_begin )  ntdim_2d_xz(0) = ntdim_2d_xz(0) + 1
    32563266                ELSE
    32573267                   ntdim_2d_xz(1) = CEILING(                            &
     
    38793889                                                   simulated_time_at_begin )  &
    38803890                                 ) / dt_do2d_yz )
     3891                IF ( do2d_at_begin )  ntdim_2d_yz(0) = ntdim_2d_yz(0) + 1
    38813892             ELSE
    38823893                ntdim_2d_yz(1) = do2d_yz_time_count(1) + CEILING(             &
     
    39073918                                             simulated_time_at_begin )  &
    39083919                           ) / dt_do2d_yz )
     3920                   IF ( do2d_at_begin )  ntdim_2d_yz(0) = ntdim_2d_yz(0) + 1
    39093921                ELSE
    39103922                   ntdim_2d_yz(1) = CEILING(                            &
Note: See TracChangeset for help on using the changeset viewer.