subroutine yso_convective use modglobal, only : je, rk3step, dt, dyi, nsv use modfields, only :sv0, svm, v0, vouttot real rk3coef integer n rk3coef = dt/(4.-dble(rk3step)) do n = 1, nsv sv0(:, je + 1, :, n) = sv0(:, je + 1, :, n) - (sv0(:, je + 1, :, n) - sv0(:, je, :, n))*dyi*rk3coef*vouttot svm(:, je + 1, :, n) = svm(:, je + 1, :, n) - (svm(:, je + 1, :, n) - svm(:, je, :, n))*dyi*rk3coef*vouttot end do end subroutine yso_convective