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

Return to conditions of master branch

Browse source
This commit is contained in:
JHendrikx 2025-02-04 17:23:54 +01:00
commit e17b2050c7
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

@ -98,7 +98,7 @@ program VlaPlEx
! Set input parameters (remember these have to be in non-dimensional units) ! Set input parameters (remember these have to be in non-dimensional units)
c_s = sqrt(11.0_dp * gamma_i * 1.0_dp) c_s = sqrt(11.0_dp * gamma_i * 1.0_dp)
bc_file = 'bc.csv' bc_file = 'bc_fa_T30Z11.csv'
call boundaryConditions%init(bc_file) call boundaryConditions%init(bc_file)
! Set domain boundaries (non-dimensional units) ! Set domain boundaries (non-dimensional units)
@ -139,7 +139,7 @@ program VlaPlEx
t0 = 0.0_dp t0 = 0.0_dp
tf = 2.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 = 7.0e2_dp*dr/c_s dt = 1.0e-2_dp*dr/c_s
nt = nint((tf - t0) / dt) nt = nint((tf - t0) / dt)
dt = (tf - t0) / float(nt) dt = (tf - t0) / float(nt)
print *, '#timesteps: ', nt print *, '#timesteps: ', nt