Skip to content

Commit d0e10ac

Browse files
Dimitrios AdamDimitrios Adam
Dimitrios Adam
authored and
Dimitrios Adam
committed
Suggestions from AI
1 parent 100309e commit d0e10ac

File tree

2 files changed

+72
-88
lines changed

2 files changed

+72
-88
lines changed

‎src/pre_process/include/1dHardcodedIC.fpp‎

Lines changed: 0 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -29,7 +29,6 @@
2929
q_prim_vf(momxb)%sf(i, 0, 0) = 0.0_wp
3030
q_prim_vf(E_idx)%sf(i, 0, 0) = 1.01325_wp*(10.0_wp)**5
3131
q_prim_vf(advxb)%sf(i, 0, 0) = 1.0_wp
32-
q_prim_vf(advxb)%sf(i, 0, 0) = 1.0_wp
3332

3433
y1 = (0.195_wp - 0.142_wp)*profile_shape + 0.142_wp
3534
y2 = (0.0_wp - 0.1_wp)*profile_shape + 0.1_wp

‎src/simulation/m_rhs.fpp‎

Lines changed: 72 additions & 87 deletions
Original file line numberDiff line numberDiff line change
@@ -1514,41 +1514,36 @@ contains
15141514
end do
15151515
end if
15161516

1517-
if (surface_tension .or. viscous) then
1517+
if ((surface_tension .or. viscous) .or. chem_params%diffusion) then
15181518
$:GPU_PARALLEL_LOOP(collapse=3)
15191519
do l = 0, p
15201520
do k = 0, n
15211521
do j = 0, m
1522-
$:GPU_LOOP(parallelism='[seq]')
1523-
do i = momxb, E_idx
1524-
rhs_vf(i)%sf(j, k, l) = &
1525-
rhs_vf(i)%sf(j, k, l) + 1._wp/dx(j)* &
1526-
(flux_src_n_in(i)%sf(j - 1, k, l) &
1527-
- flux_src_n_in(i)%sf(j, k, l))
1528-
end do
1529-
end do
1530-
end do
1531-
end do
1532-
end if
1522+
if (surface_tension .or. viscous) then
1523+
$:GPU_LOOP(parallelism='[seq]')
1524+
do i = momxb, E_idx
1525+
rhs_vf(i)%sf(j, k, l) = &
1526+
rhs_vf(i)%sf(j, k, l) + 1._wp/dx(j)* &
1527+
(flux_src_n_in(i)%sf(j - 1, k, l) &
1528+
- flux_src_n_in(i)%sf(j, k, l))
1529+
end do
1530+
end if
15331531

1534-
if (chem_params%diffusion) then
1535-
$:GPU_PARALLEL_LOOP(collapse=3)
1536-
do l = 0, p
1537-
do k = 0, n
1538-
do j = 0, m
1539-
$:GPU_LOOP(parallelism='[seq]')
1540-
do i = chemxb, chemxe
1541-
rhs_vf(i)%sf(j, k, l) = &
1542-
rhs_vf(i)%sf(j, k, l) + 1._wp/dx(j)* &
1543-
(flux_src_n_in(i)%sf(j - 1, k, l) &
1544-
- flux_src_n_in(i)%sf(j, k, l))
1545-
end do
1532+
if (chem_params%diffusion) then
1533+
$:GPU_LOOP(parallelism='[seq]')
1534+
do i = chemxb, chemxe
1535+
rhs_vf(i)%sf(j, k, l) = &
1536+
rhs_vf(i)%sf(j, k, l) + 1._wp/dx(j)* &
1537+
(flux_src_n_in(i)%sf(j - 1, k, l) &
1538+
- flux_src_n_in(i)%sf(j, k, l))
1539+
end do
15461540

1547-
if (.not. viscous) then
1548-
rhs_vf(E_idx)%sf(j, k, l) = &
1549-
rhs_vf(E_idx)%sf(j, k, l) + 1._wp/dx(j)* &
1550-
(flux_src_n_in(E_idx)%sf(j - 1, k, l) &
1551-
- flux_src_n_in(E_idx)%sf(j, k, l))
1541+
if (.not. viscous) then
1542+
rhs_vf(E_idx)%sf(j, k, l) = &
1543+
rhs_vf(E_idx)%sf(j, k, l) + 1._wp/dx(j)* &
1544+
(flux_src_n_in(E_idx)%sf(j - 1, k, l) &
1545+
- flux_src_n_in(E_idx)%sf(j, k, l))
1546+
end if
15521547
end if
15531548
end do
15541549
end do
@@ -1622,40 +1617,35 @@ contains
16221617

16231618
else
16241619

1625-
if (viscous .or. surface_tension) then
1620+
if ((surface_tension .or. viscous) .or. chem_params%diffusion) then
16261621
$:GPU_PARALLEL_LOOP(collapse=3)
16271622
do l = 0, p
16281623
do k = 0, n
16291624
do j = 0, m
1630-
$:GPU_LOOP(parallelism='[seq]')
1631-
do i = momxb, E_idx
1632-
rhs_vf(i)%sf(j, k, l) = &
1633-
rhs_vf(i)%sf(j, k, l) + 1._wp/dy(k)* &
1634-
(flux_src_n_in(i)%sf(j, k - 1, l) &
1635-
- flux_src_n_in(i)%sf(j, k, l))
1636-
end do
1637-
end do
1638-
end do
1639-
end do
1640-
end if
1625+
if (surface_tension .or. viscous) then
1626+
$:GPU_LOOP(parallelism='[seq]')
1627+
do i = momxb, E_idx
1628+
rhs_vf(i)%sf(j, k, l) = &
1629+
rhs_vf(i)%sf(j, k, l) + 1._wp/dy(k)* &
1630+
(flux_src_n_in(i)%sf(j, k - 1, l) &
1631+
- flux_src_n_in(i)%sf(j, k, l))
1632+
end do
1633+
end if
16411634

1642-
if (chem_params%diffusion) then
1643-
$:GPU_PARALLEL_LOOP(collapse=3)
1644-
do l = 0, p
1645-
do k = 0, n
1646-
do j = 0, m
1647-
$:GPU_LOOP(parallelism='[seq]')
1648-
do i = chemxb, chemxe
1649-
rhs_vf(i)%sf(j, k, l) = &
1650-
rhs_vf(i)%sf(j, k, l) + 1._wp/dy(k)* &
1651-
(flux_src_n_in(i)%sf(j, k - 1, l) &
1652-
- flux_src_n_in(i)%sf(j, k, l))
1653-
end do
1654-
if (.not. viscous) then
1655-
rhs_vf(E_idx)%sf(j, k, l) = &
1656-
rhs_vf(E_idx)%sf(j, k, l) + 1._wp/dy(k)* &
1657-
(flux_src_n_in(E_idx)%sf(j, k - 1, l) &
1658-
- flux_src_n_in(E_idx)%sf(j, k, l))
1635+
if (chem_params%diffusion) then
1636+
$:GPU_LOOP(parallelism='[seq]')
1637+
do i = chemxb, chemxe
1638+
rhs_vf(i)%sf(j, k, l) = &
1639+
rhs_vf(i)%sf(j, k, l) + 1._wp/dy(k)* &
1640+
(flux_src_n_in(i)%sf(j, k - 1, l) &
1641+
- flux_src_n_in(i)%sf(j, k, l))
1642+
end do
1643+
if (.not. viscous) then
1644+
rhs_vf(E_idx)%sf(j, k, l) = &
1645+
rhs_vf(E_idx)%sf(j, k, l) + 1._wp/dy(k)* &
1646+
(flux_src_n_in(E_idx)%sf(j, k - 1, l) &
1647+
- flux_src_n_in(E_idx)%sf(j, k, l))
1648+
end if
16591649
end if
16601650
end do
16611651
end do
@@ -1733,40 +1723,35 @@ contains
17331723
end do
17341724
end if
17351725

1736-
if (viscous .or. surface_tension) then
1726+
if ((surface_tension .or. viscous) .or. chem_params%diffusion) then
17371727
$:GPU_PARALLEL_LOOP(collapse=3)
17381728
do l = 0, p
17391729
do k = 0, n
17401730
do j = 0, m
1741-
$:GPU_LOOP(parallelism='[seq]')
1742-
do i = momxb, E_idx
1743-
rhs_vf(i)%sf(j, k, l) = &
1744-
rhs_vf(i)%sf(j, k, l) + 1._wp/dz(l)* &
1745-
(flux_src_n_in(i)%sf(j, k, l - 1) &
1746-
- flux_src_n_in(i)%sf(j, k, l))
1747-
end do
1748-
end do
1749-
end do
1750-
end do
1751-
end if
1731+
if (surface_tension .or. viscous) then
1732+
$:GPU_LOOP(parallelism='[seq]')
1733+
do i = momxb, E_idx
1734+
rhs_vf(i)%sf(j, k, l) = &
1735+
rhs_vf(i)%sf(j, k, l) + 1._wp/dz(l)* &
1736+
(flux_src_n_in(i)%sf(j, k, l - 1) &
1737+
- flux_src_n_in(i)%sf(j, k, l))
1738+
end do
1739+
end if
17521740

1753-
if (chem_params%diffusion) then
1754-
$:GPU_PARALLEL_LOOP(collapse=3)
1755-
do l = 0, p
1756-
do k = 0, n
1757-
do j = 0, m
1758-
$:GPU_LOOP(parallelism='[seq]')
1759-
do i = chemxb, chemxe
1760-
rhs_vf(i)%sf(j, k, l) = &
1761-
rhs_vf(i)%sf(j, k, l) + 1._wp/dz(l)* &
1762-
(flux_src_n_in(i)%sf(j, k, l - 1) &
1763-
- flux_src_n_in(i)%sf(j, k, l))
1764-
end do
1765-
if (.not. viscous) then
1766-
rhs_vf(E_idx)%sf(j, k, l) = &
1767-
rhs_vf(E_idx)%sf(j, k, l) + 1._wp/dz(l)* &
1768-
(flux_src_n_in(E_idx)%sf(j, k, l - 1) &
1769-
- flux_src_n_in(E_idx)%sf(j, k, l))
1741+
if (chem_params%diffusion) then
1742+
$:GPU_LOOP(parallelism='[seq]')
1743+
do i = chemxb, chemxe
1744+
rhs_vf(i)%sf(j, k, l) = &
1745+
rhs_vf(i)%sf(j, k, l) + 1._wp/dz(l)* &
1746+
(flux_src_n_in(i)%sf(j, k, l - 1) &
1747+
- flux_src_n_in(i)%sf(j, k, l))
1748+
end do
1749+
if (.not. viscous) then
1750+
rhs_vf(E_idx)%sf(j, k, l) = &
1751+
rhs_vf(E_idx)%sf(j, k, l) + 1._wp/dz(l)* &
1752+
(flux_src_n_in(E_idx)%sf(j, k, l - 1) &
1753+
- flux_src_n_in(E_idx)%sf(j, k, l))
1754+
end if
17701755
end if
17711756
end do
17721757
end do

0 commit comments

Comments
 (0)