Changeset 2130 for palm/trunk/SOURCE


Ignore:
Timestamp:
Jan 24, 2017 4:25:39 PM (7 years ago)
Author:
raasch
Message:

bugfix: in case of nested runs the stop condition in case of too small timesteps is communicated to all parent/child processes

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/Makefile

    r2119 r2130  
    2020# Current revisions:
    2121# ------------------
    22 #
     22# dependency for timestep updated
    2323#
    2424# Former revisions:
     
    584584        urban_surface_mod.o virtual_flight_mod.o wind_turbine_model_mod.o
    585585time_to_string.o: mod_kinds.o
    586 timestep.o: modules.o cpulog_mod.o mod_kinds.o microphysics_mod.o
     586timestep.o: modules.o cpulog_mod.o mod_kinds.o microphysics_mod.o \
     587 pmc_interface_mod.o
    587588timestep_scheme_steering.o: modules.o mod_kinds.o
    588589transpose.o: modules.o cpulog_mod.o mod_kinds.o
  • palm/trunk/SOURCE/timestep.f90

    r2119 r2130  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! bugfix: in case of nested runs the stop condition in case of too small
     23! timesteps is communicated to all parent/child processes,
     24! some formatting done
    2325!
    2426! Former revisions:
     
    123125    USE pegrid
    124126
     127    USE pmc_interface,                                                         &
     128        ONLY:  nested_run
     129
    125130    USE statistics,                                                            &
    126131        ONLY:  flow_statistics_called, hom, u_max, u_max_ijk, v_max, v_max_ijk,&
     
    132137    INTEGER(iwp) ::  j !<
    133138    INTEGER(iwp) ::  k !<
     139
     140    LOGICAL ::  stop_dt_local !< local switch for controlling the time stepping
    134141
    135142    REAL(wp) ::  div               !<
     
    193200             ENDDO
    194201          ENDDO
    195           uv_gtrans_l(1) = u_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb), KIND=wp )
    196           uv_gtrans_l(2) = v_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb), KIND=wp )
     202          uv_gtrans_l(1) = u_gtrans_l /                                        &
     203                           REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb), KIND=wp )
     204          uv_gtrans_l(2) = v_gtrans_l /                                        &
     205                           REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb), KIND=wp )
    197206#if defined( __parallel )
    198207          IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
    199           CALL MPI_ALLREDUCE( uv_gtrans_l, uv_gtrans, 2, MPI_REAL, MPI_SUM, &
     208          CALL MPI_ALLREDUCE( uv_gtrans_l, uv_gtrans, 2, MPI_REAL, MPI_SUM,    &
    200209                              comm2d, ierr )
    201210          u_gtrans = uv_gtrans(1) / REAL( numprocs, KIND=wp )
     
    229238          DO  j = nys, nyn
    230239             DO  k = nzb+1, nzt
    231                 dt_u_l = MIN( dt_u_l, ( dx     / ( ABS( u(k,j,i) - u_gtrans ) + 1.0E-10_wp ) ) )
    232                 dt_v_l = MIN( dt_v_l, ( dy     / ( ABS( v(k,j,i) - v_gtrans ) + 1.0E-10_wp ) ) )
    233                 dt_w_l = MIN( dt_w_l, ( dzu(k) / ( ABS( w(k,j,i) )            + 1.0E-10_wp ) ) )
     240                dt_u_l = MIN( dt_u_l, ( dx     /                               &
     241                                 ( ABS( u(k,j,i) - u_gtrans ) + 1.0E-10_wp ) ) )
     242                dt_v_l = MIN( dt_v_l, ( dy     /                               &
     243                                 ( ABS( v(k,j,i) - v_gtrans ) + 1.0E-10_wp ) ) )
     244                dt_w_l = MIN( dt_w_l, ( dzu(k) /                               &
     245                                 ( ABS( w(k,j,i) )            + 1.0E-10_wp ) ) )
    234246             ENDDO
    235247          ENDDO
     
    254266!--    Compute time step according to the diffusion criterion.
    255267!--    First calculate minimum grid spacing which only depends on index k
    256 !--    Note: also at k=nzb+1 a full grid length is being assumed, although
    257 !--          in the Prandtl-layer friction term only dz/2 is used.
    258 !--          Experience from the old model seems to justify this.
    259268       dt_diff_l = 999999.0_wp
    260269
     
    268277          DO  j = nys, nyn
    269278             DO  k = nzb+1, nzt
    270                 dt_diff_l = MIN( dt_diff_l, dxyz2_min(k) / &
    271                                        ( MAX( kh(k,j,i), km(k,j,i) ) + 1E-20_wp ) )
     279                dt_diff_l = MIN( dt_diff_l, dxyz2_min(k) /                     &
     280                                    ( MAX( kh(k,j,i), km(k,j,i) ) + 1E-20_wp ) )
    272281             ENDDO
    273282          ENDDO
     
    276285#if defined( __parallel )
    277286       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
    278        CALL MPI_ALLREDUCE( dt_diff_l, dt_diff, 1, MPI_REAL, MPI_MIN, comm2d, &
     287       CALL MPI_ALLREDUCE( dt_diff_l, dt_diff, 1, MPI_REAL, MPI_MIN, comm2d,   &
    279288                           ierr )
    280289#else
     
    286295!--    step minus a reduction (cfl_factor) to be on the safe side.
    287296!--    The time step must not exceed the maximum allowed value.
    288        dt_3d = cfl_factor * MIN( dt_diff, dt_u, dt_v, dt_w,   &
    289                                  dt_precipitation )
     297       dt_3d = cfl_factor * MIN( dt_diff, dt_u, dt_v, dt_w, dt_precipitation )
    290298       dt_3d = MIN( dt_3d, dt_max )
    291299
     
    303311          stop_dt = .TRUE.
    304312
    305           WRITE( message_string, * ) 'Time step has reached minimum limit.',  &
    306                '&dt              = ', dt_3d, ' s  Simulation is terminated.', &
    307                '&old_dt          = ', old_dt, ' s',                           &
    308                '&dt_u            = ', dt_u, ' s',                             &
    309                '&dt_v            = ', dt_v, ' s',                             &
    310                '&dt_w            = ', dt_w, ' s',                             &
    311                '&dt_diff         = ', dt_diff, ' s',                          &
    312                '&u_max           = ', u_max, ' m/s   k=', u_max_ijk(1),       &
    313                '  j=', u_max_ijk(2), '  i=', u_max_ijk(3),                    &
    314                '&v_max           = ', v_max, ' m/s   k=', v_max_ijk(1),       &
    315                '  j=', v_max_ijk(2), '  i=', v_max_ijk(3),                    &
    316                '&w_max           = ', w_max, ' m/s   k=', w_max_ijk(1),       &
     313          WRITE( message_string, * ) 'Time step has reached minimum limit.',   &
     314               '&dt              = ', dt_3d, ' s  Simulation is terminated.',  &
     315               '&old_dt          = ', old_dt, ' s',                            &
     316               '&dt_u            = ', dt_u, ' s',                              &
     317               '&dt_v            = ', dt_v, ' s',                              &
     318               '&dt_w            = ', dt_w, ' s',                              &
     319               '&dt_diff         = ', dt_diff, ' s',                           &
     320               '&u_max           = ', u_max, ' m/s   k=', u_max_ijk(1),        &
     321               '  j=', u_max_ijk(2), '  i=', u_max_ijk(3),                     &
     322               '&v_max           = ', v_max, ' m/s   k=', v_max_ijk(1),        &
     323               '  j=', v_max_ijk(2), '  i=', v_max_ijk(3),                     &
     324               '&w_max           = ', w_max, ' m/s   k=', w_max_ijk(1),        &
    317325               '  j=', w_max_ijk(2), '  i=', w_max_ijk(3)
    318326          CALL message( 'timestep', 'PA0312', 0, 1, 0, 6, 0 )
     
    324332          IF ( coupling_mode /= 'uncoupled' .AND. terminate_coupled == 0 )  THEN
    325333             terminate_coupled = 2
    326              IF ( myid == 0 ) THEN
     334             IF ( myid == 0 )  THEN
    327335                CALL MPI_SENDRECV( &
    328                      terminate_coupled,        1, MPI_INTEGER, target_id,  0, &
    329                      terminate_coupled_remote, 1, MPI_INTEGER, target_id,  0, &
     336                     terminate_coupled,        1, MPI_INTEGER, target_id,  0,  &
     337                     terminate_coupled_remote, 1, MPI_INTEGER, target_id,  0,  &
    330338                     comm_inter, status, ierr )
    331339             ENDIF
    332              CALL MPI_BCAST( terminate_coupled_remote, 1, MPI_INTEGER, 0, comm2d, ierr)
     340             CALL MPI_BCAST( terminate_coupled_remote, 1, MPI_INTEGER, 0,      &
     341                             comm2d, ierr)
    333342          ENDIF
    334343#endif
     344       ENDIF
     345
     346!
     347!--    In case of nested runs all parent/child processes have to terminate if
     348!--    one process has set the stop flag, i.e. they need to set the stop flag
     349!--    too.
     350       IF ( nested_run )  THEN
     351          stop_dt_local = stop_dt
     352          CALL MPI_ALLREDUCE( stop_dt_local, stop_dt, 1, MPI_LOGICAL, MPI_LOR, &
     353                              MPI_COMM_WORLD, ierr )
    335354       ENDIF
    336355
Note: See TracChangeset for help on using the changeset viewer.