@@ -1523,8 +1523,8 @@ contains
15231523 do i = momxb, E_idx
15241524 rhs_vf(i)%sf(j, k, l) = &
15251525 rhs_vf(i)%sf(j, k, l) + 1._wp / dx(j)* &
1526- (flux_src_n (i)%sf(j - 1 , k, l) &
1527- - flux_src_n (i)%sf(j, k, l))
1526+ (flux_src_n_in (i)%sf(j - 1 , k, l) &
1527+ - flux_src_n_in (i)%sf(j, k, l))
15281528 end do
15291529 end do
15301530 end do
@@ -1540,15 +1540,15 @@ contains
15401540 do i = chemxb, chemxe
15411541 rhs_vf(i)%sf(j, k, l) = &
15421542 rhs_vf(i)%sf(j, k, l) + 1._wp / dx(j)* &
1543- (flux_src_n (i)%sf(j - 1 , k, l) &
1544- - flux_src_n (i)%sf(j, k, l))
1543+ (flux_src_n_in (i)%sf(j - 1 , k, l) &
1544+ - flux_src_n_in (i)%sf(j, k, l))
15451545 end do
15461546
15471547 if (.not. viscous) then
15481548 rhs_vf(E_idx)%sf(j, k, l) = &
15491549 rhs_vf(E_idx)%sf(j, k, l) + 1._wp / dx(j)* &
1550- (flux_src_n (E_idx)%sf(j - 1 , k, l) &
1551- - flux_src_n (E_idx)%sf(j, k, l))
1550+ (flux_src_n_in (E_idx)%sf(j - 1 , k, l) &
1551+ - flux_src_n_in (E_idx)%sf(j, k, l))
15521552 end if
15531553 end do
15541554 end do
@@ -1631,8 +1631,8 @@ contains
16311631 do i = momxb, E_idx
16321632 rhs_vf(i)%sf(j, k, l) = &
16331633 rhs_vf(i)%sf(j, k, l) + 1._wp / dy(k)* &
1634- (flux_src_n (i)%sf(j, k - 1 , l) &
1635- - flux_src_n (i)%sf(j, k, l))
1634+ (flux_src_n_in (i)%sf(j, k - 1 , l) &
1635+ - flux_src_n_in (i)%sf(j, k, l))
16361636 end do
16371637 end do
16381638 end do
@@ -1648,14 +1648,14 @@ contains
16481648 do i = chemxb, chemxe
16491649 rhs_vf(i)%sf(j, k, l) = &
16501650 rhs_vf(i)%sf(j, k, l) + 1._wp / dy(k)* &
1651- (flux_src_n (i)%sf(j, k - 1 , l) &
1652- - flux_src_n (i)%sf(j, k, l))
1651+ (flux_src_n_in (i)%sf(j, k - 1 , l) &
1652+ - flux_src_n_in (i)%sf(j, k, l))
16531653 end do
16541654 if (.not. viscous) then
16551655 rhs_vf(E_idx)%sf(j, k, l) = &
16561656 rhs_vf(E_idx)%sf(j, k, l) + 1._wp / dy(k)* &
1657- (flux_src_n (E_idx)%sf(j, k - 1 , l) &
1658- - flux_src_n (E_idx)%sf(j, k, l))
1657+ (flux_src_n_in (E_idx)%sf(j, k - 1 , l) &
1658+ - flux_src_n_in (E_idx)%sf(j, k, l))
16591659 end if
16601660 end do
16611661 end do
@@ -1742,8 +1742,8 @@ contains
17421742 do i = momxb, E_idx
17431743 rhs_vf(i)%sf(j, k, l) = &
17441744 rhs_vf(i)%sf(j, k, l) + 1._wp / dz(l)* &
1745- (flux_src_n (i)%sf(j, k, l - 1 ) &
1746- - flux_src_n (i)%sf(j, k, l))
1745+ (flux_src_n_in (i)%sf(j, k, l - 1 ) &
1746+ - flux_src_n_in (i)%sf(j, k, l))
17471747 end do
17481748 end do
17491749 end do
@@ -1759,14 +1759,14 @@ contains
17591759 do i = chemxb, chemxe
17601760 rhs_vf(i)%sf(j, k, l) = &
17611761 rhs_vf(i)%sf(j, k, l) + 1._wp / dz(l)* &
1762- (flux_src_n (i)%sf(j, k, l - 1 ) &
1763- - flux_src_n (i)%sf(j, k, l))
1762+ (flux_src_n_in (i)%sf(j, k, l - 1 ) &
1763+ - flux_src_n_in (i)%sf(j, k, l))
17641764 end do
17651765 if (.not. viscous) then
17661766 rhs_vf(E_idx)%sf(j, k, l) = &
17671767 rhs_vf(E_idx)%sf(j, k, l) + 1._wp / dz(l)* &
1768- (flux_src_n (E_idx)%sf(j, k, l - 1 ) &
1769- - flux_src_n (E_idx)%sf(j, k, l))
1768+ (flux_src_n_in (E_idx)%sf(j, k, l - 1 ) &
1769+ - flux_src_n_in (E_idx)%sf(j, k, l))
17701770 end if
17711771 end do
17721772 end do
0 commit comments