ICs

wilfonba · Feb 14, 2024 · 3c42479 · 3c42479
1 parent 158fc26
commit 3c42479
Show file tree

Hide file tree

Showing 3 changed files with 8 additions and 8 deletions.
diff --git a/src/pre_process/include/2dHardcodedIC.fpp b/src/pre_process/include/2dHardcodedIC.fpp
@@ -144,7 +144,7 @@
 
         case(23004) ! 2D Interface
 
-            ih = 0.0436 - 0.00186/20*(sin((2*pi/1.86d-3)*x_cc(i)+pi/2))
+            ih = 0.01116 - 0.00186/20*(sin((2*pi/1.86d-3)*x_cc(i)+pi/2))
             alph = 5d-1*(1 + tanh((y_cc(j) - ih)/1d-16))
 
             if (alph < 1e-6) alph = 1e-6
@@ -156,9 +156,9 @@
             q_prim_vf(contxb)%sf(i, j, 0) = (1 - alph)*950d0
             q_prim_vf(contxe)%sf(i, j, 0) = alph*1d0
 
-            if (y_cc(j)  > 0.0436) then
-                pInterface = 1d5 + 950*9.81*0.0439
-                q_prim_vf(E_idx)%sf(i, j, 0) = pInterface + 1d0*9.81*(y_cc(j) - 0.0436)
+            if (y_cc(j)  > 0.01116) then
+                pInterface = 1d5 + 950*9.81*0.01116
+                q_prim_vf(E_idx)%sf(i, j, 0) = pInterface + 1d0*9.81*(y_cc(j) - 0.01116)
             else
                 q_prim_vf(E_idx)%sf(i, j, 0) = 1d5 + 950d0*9.81*(y_cc(j))
             end if

diff --git a/src/simulation/m_boundary_conditions.fpp b/src/simulation/m_boundary_conditions.fpp
@@ -1114,7 +1114,7 @@ contains
                                         -q_prim_vf(i)%sf(l, j - 1, k) + 2d0*bc_y%vb1
                                 elseif (i == momxb + 1 .and. num_dims > 1) then
                                     q_prim_vf(i)%sf(l, -j, k) = &
-                                        -q_prim_vf(i)%sf(l, j - 1, k) + 2d0*(0.5*sin(600*pi*t))
+                                        -q_prim_vf(i)%sf(l, j - 1, k) + 2d0*bc_y%vb1
                                 elseif (i == momxb + 2 .and. num_dims > 2) then
                                     q_prim_vf(i)%sf(l, -j, k) = &
                                         -q_prim_vf(i)%sf(l, j - 1, k) + 2d0*bc_y%vb3
@@ -1138,7 +1138,7 @@ contains
                                         -q_prim_vf(i)%sf(l, n - (j - 1), k) + 2d0*bc_y%ve1
                                 elseif (i == momxb + 1 .and. num_dims > 1) then
                                     q_prim_vf(i)%sf(l, n + j, k) = &
-                                        -q_prim_vf(i)%sf(l, n - (j - 1), k) + 2d0*(0.5*sin(600*pi*t))
+                                        -q_prim_vf(i)%sf(l, n - (j - 1), k) + 2d0*bc_y%ve2
                                 elseif (i == momxb + 2 .and. num_dims > 2) then
                                     q_prim_vf(i)%sf(l, n + j, k) = &
                                         -q_prim_vf(i)%sf(l, n - (j - 1), k) + 2d0*bc_y%ve3

diff --git a/toolchain/mfc/bench.py b/toolchain/mfc/bench.py
@@ -1,6 +1,6 @@
 import os, sys, uuid, subprocess, dataclasses
 
-import rich.table
+import rich.table, typing
 
 from .printer import cons
 from .state   import ARG, CFG
@@ -13,7 +13,7 @@
 class BenchCase:
     slug: str
     path: str
-    args: list[str]
+    args: typing.List[str]
 
 
 def bench(targets = None):