Time Step Integration Sequence

+ ATMOS_DYN
  + ATMOS_DYN_tinteg_large
    | {var_dyn}0 = {var_dyn}
    - RK loop
      | {var_dyn} = {var_dyn}0
      - ATMOS_DYN_tstep_large
        | nstep = ceiling( dt_large/dt_small )
    | calculation pres0 from {var_dyn}0, {var_tracer}
    | calculation damping for var_tracer from {var_tracer}
    + nstep loop
    | | calculation damping for var_dyn from {var_dyn}
    | | calculation numerical diffusion for var_dyn from {var_dyn}
    | | calculation divergence for var_dyn from {var_dyn}
         | | {var_dyn}00 = {var_dyn}0 at the first nstep loop, else {var_dyn}00 = {var_dyn}
    | - ATMOS_DYN_tinteg_short
    |   - RK loop
    |     - ATMOS_DYN_tstep_short
    |       - {var_dyn} = {var_dyn}00 + f({var_dyn})
    | calculation momentum flux for tracer advection (nstep loop average)
    | calculation numerical diffusion for var_tracer from {var_tracer}
    - ATMOS_DYN_tinteg_tracer
      - RK loop
        - ATMOS_DYN_tstep_tracer
          - {var_tracer} = {var_tracer}0 + g({var_tracer}, mflux)
