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

Use easiest case of 2 Zave values

Browse source
This commit is contained in:
JHendrikx 2025-02-04 17:24:32 +01:00
commit 23c5d5fb45
1 changed files with 2 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

4
vlaplex.f90
View file

@ -158,7 +158,7 @@ program VlaPlEx
write(*, '(A,ES0.4e3)') 'CFL: ', dt*vf/dr write(*, '(A,ES0.4e3)') 'CFL: ', dt*vf/dr
nz = 16 nz = 2
! Allocate vectors ! Allocate vectors
allocate(f_i(1:nz,1:nr,1:nv), f_i_old(1:nz,1:nr,1:nv)) allocate(f_i(1:nz,1:nr,1:nv), f_i_old(1:nz,1:nr,1:nv))
allocate(n_i(1:nz,1:nr)) allocate(n_i(1:nz,1:nr))
@ -179,7 +179,7 @@ program VlaPlEx
n_e = 0.0_dp n_e = 0.0_dp
T_e = 0.0_dp T_e = 0.0_dp
Zave = 0.0_dp Zave = 0.0_dp
Z_list = (/(I, I = 1, nz) /) Z_list = (/ 6.0, 12.0 /)
phi = 0.0_dp phi = 0.0_dp
phi_old = 0.0_dp phi_old = 0.0_dp
E = 0.0_dp E = 0.0_dp