From 4ff2e04c1b00527ecfe394e38629d25b1faa3490 Mon Sep 17 00:00:00 2001 From: Anand Date: Sat, 14 Sep 2024 07:51:20 +0530 Subject: [PATCH] quasi 2d --- src/pre_process/m_assign_variables.f90 | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) diff --git a/src/pre_process/m_assign_variables.f90 b/src/pre_process/m_assign_variables.f90 index f13b082b9..bf59897d5 100644 --- a/src/pre_process/m_assign_variables.f90 +++ b/src/pre_process/m_assign_variables.f90 @@ -191,11 +191,13 @@ subroutine s_perturb_primitive(j, k, l, q_prim_vf) V_in = 0.1*sqrt(1.4*P_in) beta = 0.01d0 - if(p == 0) then - q_prim_vf(momxb)%sf(j, k, l) = V_in*sin(x_cc(j))* cos(y_cc(k)) - q_prim_vf(momxb+1)%sf(j, k, l) = -V_in*cos(x_cc(j))* sin(y_cc(k)) - q_prim_vf(E_idx)%sf(j, k, l) = P_in + (V_in**2/4d0)*(cos(2d0*x_cc(j)) + cos(2d0*y_cc(k))) - + if(p /= -1) then + q_prim_vf(momxb+1)%sf(j, k, l) = V_in*sin(y_cc(k))* cos(z_cc(l)) + q_prim_vf(momxb+2)%sf(j, k, l) = -V_in*cos(y_cc(k))* sin(z_cc(l)) + q_prim_vf(E_idx)%sf(j, k, l) = P_in + (V_in**2/4d0)*(cos(2d0*y_cc(k)) + cos(2d0*z_cc(l))) + if(p > 0) then + q_prim_vf(momxb)%sf(j, k, l) = 0d0 + end if R_in = sqrt((x_cc(j) - pi)**2d0 + (y_cc(k) - pi)**2d0) !q_prim_vf(momxb)%sf(j, k, l) = q_prim_vf(momxb)%sf(j, k, l) - V_in*beta*(y_cc(k) - pi)*exp(-0.5d0*R_in**2d0) / R_in