subroutine xso_convective use modglobal, only : ie, rk3step, dt, dxi, nsv use modfields, only :sv0, svm, uouttot real rk3coef integer n rk3coef = dt/(4.-dble(rk3step)) do n = 1, nsv sv0(ie + 1, :, :, n) = sv0(ie + 1, :, :, n) - (sv0(ie + 1, :, :, n) - sv0(ie, :, :, n))*dxi*rk3coef*uouttot svm(ie + 1, :, :, n) = svm(ie + 1, :, :, n) - (svm(ie + 1, :, :, n) - svm(ie, :, :, n))*dxi*rk3coef*uouttot end do end subroutine xso_convective