JGonzalez/vlaplex
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Set correct boundary conditions

Browse source
This commit is contained in:
JHendrikx 2025-02-05 16:08:27 +01:00
commit 868c13e494
1 changed files with 2 additions and 17 deletions
Download patch file Download diff file Expand all files Collapse all files

19
vlaplex.f90
View file

@ -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(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) f0 = f0 * n_bc / (sum(f0)*dv)
!f_i_old(:,1,j0:nv) = 0 f_i_old(:,1,j0:nv) = 0.0_dp
!f_i(:,1,j0:nv) = 0
!T_i(:,1) = 0
f_i_old(z_inj,1,j0:nv) = f0 f_i_old(z_inj,1,j0:nv) = f0
f_i(z_inj,1,j0:nv) = f_i_old(z_inj,1,j0:nv) f_i(:,1,j0:nv) = f_i_old(:,1,j0:nv)
T_i(z_inj,1) = Temp_bc
!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 T_e = Temp_bc
print *, 'TtoZ: ', T_to_Z(Temp_bc) print *, 'TtoZ: ', T_to_Z(Temp_bc)
print *, 'Injection species: ', Zave_bc print *, 'Injection species: ', Zave_bc