Ignore:
Timestamp:
Dec 28, 2018 5:58:49 PM (3 years ago)
Author:
kanani
Message:

Bugfix: replace simulated_time by time_since_reference_point where required

File:
1 edited

Legend:

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

    r3637 r3646  
    2525! -----------------
    2626! $Id$
     27! Bugfix: use time_since_reference_point instead of simulated_time (relevant
     28! when using wall/soil spinup)
     29!
     30! 3637 2018-12-20 01:51:36Z knoop
    2731! Implementation of the PALM module interface
    2832!
     
    298302               ibc_uv_b, io_blocks, io_group, message_string,                  &
    299303               ntdim_2d_xy, ntdim_2d_xz, ntdim_2d_yz,                          &
    300                psolver, section, simulated_time,                               &
     304               psolver, section,                                               &
    301305               time_since_reference_point
    302306
     
    510514            ntdim_2d_xy(av) )  THEN
    511515          WRITE ( message_string, * ) 'Output of xy cross-sections is not ',   &
    512                           'given at t=', simulated_time, '&because the',       &
     516                          'given at t=', time_since_reference_point, 's because the',       &
    513517                          ' maximum number of output time levels is exceeded.'
    514518          CALL message( 'data_output_2d', 'PA0384', 0, 1, 0, 6, 0 )
     
    519523            ntdim_2d_xz(av) )  THEN
    520524          WRITE ( message_string, * ) 'Output of xz cross-sections is not ',   &
    521                           'given at t=', simulated_time, '&because the',       &
     525                          'given at t=', time_since_reference_point, 's because the',       &
    522526                          ' maximum number of output time levels is exceeded.'
    523527          CALL message( 'data_output_2d', 'PA0385', 0, 1, 0, 6, 0 )
     
    528532            ntdim_2d_yz(av) )  THEN
    529533          WRITE ( message_string, * ) 'Output of yz cross-sections is not ',   &
    530                           'given at t=', simulated_time, '&because the',       &
     534                          'given at t=', time_since_reference_point, 's because the',       &
    531535                          ' maximum number of output time levels is exceeded.'
    532536          CALL message( 'data_output_2d', 'PA0386', 0, 1, 0, 6, 0 )
     
    697701             CASE ( 'pc_xy', 'pc_xz', 'pc_yz' )  ! particle concentration
    698702                IF ( av == 0 )  THEN
    699                    IF ( simulated_time >= particle_advection_start )  THEN
     703                   IF ( time_since_reference_point >= particle_advection_start )  THEN
    700704                      tend = prt_count
    701705!                      CALL exchange_horiz( tend, nbgp )
     
    722726             CASE ( 'pr_xy', 'pr_xz', 'pr_yz' )  ! mean particle radius (effective radius)
    723727                IF ( av == 0 )  THEN
    724                    IF ( simulated_time >= particle_advection_start )  THEN
     728                   IF ( time_since_reference_point >= particle_advection_start )  THEN
    725729                      DO  i = nxl, nxr
    726730                         DO  j = nys, nyn
     
    844848             CASE ( 'ql_vp_xy', 'ql_vp_xz', 'ql_vp_yz' )
    845849                IF ( av == 0 )  THEN
    846                    IF ( simulated_time >= particle_advection_start )  THEN
     850                   IF ( time_since_reference_point >= particle_advection_start )  THEN
    847851                      DO  i = nxl, nxr
    848852                         DO  j = nys, nyn
     
    14221426!--                In case of parallel output, this is only done by PE0
    14231427!--                to increase the performance.
    1424                    IF ( simulated_time /= do2d_xy_last_time(av) )  THEN
     1428                   IF ( time_since_reference_point /= do2d_xy_last_time(av) )  THEN
    14251429                      do2d_xy_time_count(av) = do2d_xy_time_count(av) + 1
    1426                       do2d_xy_last_time(av)  = simulated_time
     1430                      do2d_xy_last_time(av)  = time_since_reference_point
    14271431                      IF ( myid == 0 )  THEN
    14281432                         IF ( .NOT. data_output_2d_on_each_pe  &
     
    16111615!--                In case of parallel output, this is only done by PE0
    16121616!--                to increase the performance.
    1613                    IF ( simulated_time /= do2d_xz_last_time(av) )  THEN
     1617                   IF ( time_since_reference_point /= do2d_xz_last_time(av) )  THEN
    16141618                      do2d_xz_time_count(av) = do2d_xz_time_count(av) + 1
    1615                       do2d_xz_last_time(av)  = simulated_time
     1619                      do2d_xz_last_time(av)  = time_since_reference_point
    16161620                      IF ( myid == 0 )  THEN
    16171621                         IF ( .NOT. data_output_2d_on_each_pe  &
     
    18291833!--                In case of parallel output, this is only done by PE0
    18301834!--                to increase the performance.
    1831                    IF ( simulated_time /= do2d_yz_last_time(av) )  THEN
     1835                   IF ( time_since_reference_point /= do2d_yz_last_time(av) )  THEN
    18321836                      do2d_yz_time_count(av) = do2d_yz_time_count(av) + 1
    1833                       do2d_yz_last_time(av)  = simulated_time
     1837                      do2d_yz_last_time(av)  = time_since_reference_point
    18341838                      IF ( myid == 0 )  THEN
    18351839                         IF ( .NOT. data_output_2d_on_each_pe  &
Note: See TracChangeset for help on using the changeset viewer.