From 868c13e4945de1b5a1bcf4821d5d335d51723d2c Mon Sep 17 00:00:00 2001 From: JHendrikx Date: Wed, 5 Feb 2025 16:08:27 +0100 Subject: [PATCH] Set correct boundary conditions --- vlaplex.f90 | 19 ++----------------- 1 file changed, 2 insertions(+), 17 deletions(-) diff --git a/vlaplex.f90 b/vlaplex.f90 index 734910a..a73d139 100644 --- a/vlaplex.f90 +++ b/vlaplex.f90 @@ -252,25 +252,10 @@ program VlaPlEx f0(j0:nv) = 1.0_dp / sqrt(PI*Temp_bc) * exp(-(v(j0:nv) - u_bc)**2 / Temp_bc) f0 = f0 * n_bc / (sum(f0)*dv) - !f_i_old(:,1,j0:nv) = 0 - !f_i(:,1,j0:nv) = 0 - !T_i(:,1) = 0 - + f_i_old(:,1,j0:nv) = 0.0_dp f_i_old(z_inj,1,j0:nv) = f0 - f_i(z_inj,1,j0:nv) = f_i_old(z_inj,1,j0:nv) - T_i(z_inj,1) = Temp_bc + f_i(:,1,j0:nv) = f_i_old(:,1,j0:nv) - !do iz = 1, nz - ! ! f0(j0:nv) = v(j0:nv)**2 / sqrt(PI*Temp_bc**3) * exp(-(v(j0:nv) - u_bc)**2 / Temp_bc) - ! f0(j0:nv) = 1.0_dp / sqrt(PI*Temp_bc) * exp(-(v(j0:nv) - u_bc)**2 / Temp_bc) - ! f0 = f0 * n_bc / (sum(f0)*dv) - - ! ! Boundary conditions - ! ! r = r0, v>0 - ! f_i_old(iz,1,j0:nv) = f0 - ! f_i(iz,1,j0:nv) = f_i_old(iz,1,j0:nv) - ! T_i(iz,1) = Temp_bc - !end do T_e = Temp_bc print *, 'TtoZ: ', T_to_Z(Temp_bc) print *, 'Injection species: ', Zave_bc