Skip to content

Commit

Permalink
ic changes
Browse files Browse the repository at this point in the history
  • Loading branch information
@wilfonba
wilfonba committed Jun 17, 2024
1 parent c400cb6 commit e8bc25e
Showing 1 changed file with 64 additions and 64 deletions.
128 changes: 64 additions & 64 deletions src/pre_process/include/2dHardcodedIC.fpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -139,13 +139,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (22)
Expand All @@ -163,13 +163,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (23)
Expand All @@ -187,13 +187,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (24)
Expand All @@ -211,13 +211,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (241)
Expand All @@ -235,13 +235,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (242)
Expand All @@ -259,13 +259,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (243)
Expand All @@ -283,13 +283,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (244)
Expand All @@ -307,13 +307,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (11)
Expand All @@ -331,13 +331,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (12)
Expand All @@ -355,13 +355,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (13)
Expand All @@ -379,13 +379,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (14)
Expand All @@ -403,13 +403,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (31)
Expand All @@ -427,13 +427,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (32)
Expand All @@ -451,13 +451,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (33)
Expand All @@ -475,13 +475,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case (34)
Expand All @@ -499,13 +499,13 @@
if (sigma .ne. dflt_real) q_prim_vf(c_idx)%sf(i, j, k) = alph
q_prim_vf(advxb)%sf(i, j, 0) = 1 - alph
q_prim_vf(advxe)%sf(i, j, 0) = alph
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*1d0
q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*1d0
q_prim_vf(contxe)%sf(i, j, 9) = alph*0.1784d0

if (y_cc(j) > ih) then
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 1d0*acc*(y_cc(j) - ih)
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 - 0.1784d0*acc*(y_cc(j) - ih)
else
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*acc*(ih - y_cc(j))
q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 1d0*acc*(ih - y_cc(j))
end if

case default
Expand Down

0 comments on commit e8bc25e

Please sign in to comment.