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

Erase Zave array

Browse source
This commit is contained in:
JHendrikx 2025-02-05 16:09:03 +01:00
commit ce4ba4b212
1 changed files with 1 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

8
vlaplex.f90
View file

@ -261,9 +261,6 @@ program VlaPlEx
print *, 'Injection species: ', Zave_bc
print *, '-------------------------'
Zave(1) = Zave_bc
! r = rf, v<0
f_i_old(:,nr,1:j0-1) = 0.0_dp
f_i(:,nr,1:j0-1) = f_i_old(:,nr,1:j0-1)
@ -291,12 +288,9 @@ program VlaPlEx
u_i(iz,i) = sum(v(:) *f_i(iz,i,:))*dv / n_i(iz,i)
E_i(i) = sum(v(:)**2*f_i(iz,i,:))*dv / n_i(iz,i)
T_i(iz,i) = 2.0_dp*E_i(i) - 2.0_dp*u_i(iz,i)**2
Zave(i) = Zave_bc
else
u_i(iz,i) = 0.0_dp
T_i(iz,i) = 0.0_dp
Zave(i) = 0.0_dp
end if
end do
@ -420,7 +414,7 @@ program VlaPlEx
if (mod(t,everyOutput) == 0 .or. t == nt) then
! call writeOutputF(t, dt, nr, r, nv, v, f_i_old)
call writeOutputPhi(t, dt, nr, r, phi, E, n_e)
call writeOutputMom(t, dt, nr, r, n_i(1,:), u_i(1,:), T_i(1,:), Zave)
call writeOutputMom(t, dt, nr, r, n_i(1,:), u_i(1,:), T_i(1,:), Z_list(1))
call writeOutputFCum(t, dt, r(rCum_index), nv, v, fCum_i(1,:))
end if