From ce4ba4b212f0cca4a369adf799e40b32753309aa Mon Sep 17 00:00:00 2001 From: JHendrikx Date: Wed, 5 Feb 2025 16:09:03 +0100 Subject: [PATCH] Erase Zave array --- vlaplex.f90 | 8 +------- 1 file changed, 1 insertion(+), 7 deletions(-) diff --git a/vlaplex.f90 b/vlaplex.f90 index a73d139..fe6b4e1 100644 --- a/vlaplex.f90 +++ b/vlaplex.f90 @@ -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