Changeset 3575 for palm/trunk


Ignore:
Timestamp:
Nov 28, 2018 3:06:32 PM (6 years ago)
Author:
kanani
Message:

Bugfix for output time levels in parallel NetCDF output with spinup (check_parameters)

File:
1 edited

Legend:

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

    r3569 r3575  
    2525! -----------------
    2626! $Id$
     27! Bugfix for output time levels in parallel NetCDF output with spinup
     28!
     29! 3569 2018-11-27 17:03:40Z kanani
    2730! Sort chemistry routine calls into list.
    2831! dom_dwd_user, Schrempf:
     
    34363439!--    Set needed time levels (ntdim) to
    34373440!--    saved time levels + to be saved time levels.
    3438        ntdim_3d(0) = do3d_time_count(0) + CEILING(                             &
    3439                      ( end_time - MAX( skip_time_do3d,                         &
    3440                                        simulated_time_at_begin )               &
     3441       ntdim_3d(0) = do3d_time_count(0) + CEILING(                                    &
     3442                     ( end_time - MAX(                                                &
     3443                         MERGE(skip_time_do3d, skip_time_do3d + spinup_time,          &
     3444                               data_output_during_spinup ),                           &
     3445                         simulated_time_at_begin )                                    &
    34413446                     ) / dt_do3d )
    34423447       IF ( do3d_at_begin ) ntdim_3d(0) = ntdim_3d(0) + 1
    34433448
    3444        ntdim_3d(1) = do3d_time_count(1) + CEILING(                             &
    3445                      ( end_time - MAX( skip_time_data_output_av,               &
    3446                                        simulated_time_at_begin )               &
     3449       ntdim_3d(1) = do3d_time_count(1) + CEILING(                                    &
     3450                     ( end_time - MAX(                                                &
     3451                         MERGE(skip_time_data_output_av, skip_time_data_output_av     &
     3452                              + spinup_time, data_output_during_spinup ),             &
     3453                         simulated_time_at_begin )                                    &
    34473454                     ) / dt_data_output_av )
    34483455
    3449        ntdim_2d_xy(0) = do2d_xy_time_count(0) + CEILING(                       &
    3450                         ( end_time - MAX( skip_time_do2d_xy,                   &
    3451                                           simulated_time_at_begin )            &
     3456       ntdim_2d_xy(0) = do2d_xy_time_count(0) + CEILING(                              &
     3457                        ( end_time - MAX(                                             &
     3458                           MERGE(skip_time_do2d_xy, skip_time_do2d_xy + spinup_time,  &
     3459                               data_output_during_spinup ),                           &
     3460                                          simulated_time_at_begin )                   &
    34523461                        ) / dt_do2d_xy )
    34533462
    3454        ntdim_2d_xz(0) = do2d_xz_time_count(0) + CEILING(                       &
    3455                         ( end_time - MAX( skip_time_do2d_xz,                   &
    3456                                           simulated_time_at_begin )            &
     3463       ntdim_2d_xz(0) = do2d_xz_time_count(0) + CEILING(                              &
     3464                        ( end_time - MAX(                                             &
     3465                         MERGE(skip_time_do2d_xz, skip_time_do2d_xz + spinup_time,    &
     3466                               data_output_during_spinup ),                           &
     3467                         simulated_time_at_begin )                                    &
    34573468                        ) / dt_do2d_xz )
    34583469
    3459        ntdim_2d_yz(0) = do2d_yz_time_count(0) + CEILING(                       &
    3460                         ( end_time - MAX( skip_time_do2d_yz,                   &
    3461                                           simulated_time_at_begin )            &
     3470       ntdim_2d_yz(0) = do2d_yz_time_count(0) + CEILING(                              &
     3471                        ( end_time - MAX(                                             &
     3472                         MERGE(skip_time_do2d_yz, skip_time_do2d_yz + spinup_time,    &
     3473                               data_output_during_spinup ),                           &
     3474                         simulated_time_at_begin )                                    &
    34623475                        ) / dt_do2d_yz )
    34633476
Note: See TracChangeset for help on using the changeset viewer.