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

Adjust length and simulation time

Browse source
This commit is contained in:
Jorge Gonzalez 2024-11-20 13:58:16 +01:00
commit 17021f8194
1 changed files with 3 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

6
plasmaExpansion.f90
View file

@ -99,7 +99,7 @@ program plasmaExpansion
! Set domain boundaries (non-dimensional units) ! Set domain boundaries (non-dimensional units)
r0 = 10.0e-6_dp / L_ref r0 = 10.0e-6_dp / L_ref
rf = 1.0e-3_dp / L_ref rf = 2.0e-3_dp / L_ref
dr = 5.0e2_dp dr = 5.0e2_dp
nr = nint((rf - r0) / dr) + 1 nr = nint((rf - r0) / dr) + 1
dr = (rf - r0) / float(nr-1) dr = (rf - r0) / float(nr-1)
@ -109,7 +109,7 @@ program plasmaExpansion
end do end do
! Set position to calculate cumulative sum of f (non-dimensional units) ! Set position to calculate cumulative sum of f (non-dimensional units)
rCum = 5.0e-4 / L_ref rCum = 1.0e-3 / L_ref
! Index for cumulative sum ! Index for cumulative sum
rCum_index = minloc(abs(r - rCum), 1) rCum_index = minloc(abs(r - rCum), 1)
@ -131,7 +131,7 @@ program plasmaExpansion
end if end if
t0 = 0.0_dp t0 = 0.0_dp
tf = 5.0e-7_dp / t_ref tf = 2.0e-7_dp / t_ref
! tf = 1.0e1_dp * (rf - r0) / c_s ! tf = 1.0e1_dp * (rf - r0) / c_s
dt = 1.0e-2_dp*dr/c_s dt = 1.0e-2_dp*dr/c_s
nt = nint((tf - t0) / dt) nt = nint((tf - t0) / dt)