MHD1D_shock Program

Uses

  • program~~mhd1d_shock~~UsesGraph program~mhd1d_shock MHD1D_shock module~advec advec program~mhd1d_shock->module~advec module~phys_consts phys_consts module~advec->module~phys_consts iso_fortran_env iso_fortran_env module~phys_consts->iso_fortran_env

Calls

program~~mhd1d_shock~~CallsGraph program~mhd1d_shock MHD1D_shock proc~writearray~2 writearray program~mhd1d_shock->proc~writearray~2 proc~advec1d_mc advec1D_MC program~mhd1d_shock->proc~advec1d_mc proc~derivative~2 derivative program~mhd1d_shock->proc~derivative~2 proc~minmod~2 minmod proc~advec1d_mc->proc~minmod~2

Contents


Variables

Type AttributesNameInitial
integer, parameter:: method =1
integer, parameter:: npts =1001
real(kind=wp), parameter:: tcfl =0.75
real(kind=wp), parameter:: pi =3.14
real(kind=wp), parameter:: amu =1.67e-27
real(kind=wp), parameter:: kb =1.38e-23
real(kind=wp), parameter:: gammap =5.0/3.0
real(kind=wp), parameter:: mu0 =4.0*pi*1e-7
real(kind=wp), parameter:: stride =2e3
real(kind=wp), dimension(npts):: dx1i
real(kind=wp), dimension(-1:npts+2):: x1
real(kind=wp), dimension(-1:npts+2):: rhom
real(kind=wp), dimension(-1:npts+2):: rhou1
real(kind=wp), dimension(-1:npts+2):: rhou2
real(kind=wp), dimension(-1:npts+2):: rhou3
real(kind=wp), dimension(-1:npts+2):: B1
real(kind=wp), dimension(-1:npts+2):: B2
real(kind=wp), dimension(-1:npts+2):: B3
real(kind=wp), dimension(-1:npts+2):: rhoeps
real(kind=wp), dimension(-1:npts+2):: u1
real(kind=wp), dimension(-1:npts+2):: u2
real(kind=wp), dimension(-1:npts+2):: u3
real(kind=wp), dimension(npts+1):: x1i
real(kind=wp), dimension(npts+1):: v1i
real(kind=wp), dimension(0:npts+2):: dx1
integer :: lx1
integer :: it
integer :: ix1
integer :: u
real(kind=wp) :: t =0
real(kind=wp) :: dt =1e-6
real(kind=wp) :: dtout =5
real(kind=wp) :: tdur =55
real(kind=wp) :: xi =2
real(kind=wp) :: toutnext
real(kind=wp), dimension(npts):: Q
real(kind=wp), dimension(npts):: p
real(kind=wp), dimension(npts):: grad1B1u3
real(kind=wp), dimension(npts):: rhoepshalf
real(kind=wp), dimension(npts):: sourceterm
real(kind=wp), dimension(npts):: du1full
real(kind=wp), dimension(npts):: grad1u1
real(kind=wp), dimension(npts):: vA
real(kind=wp), dimension(npts):: vsnd
real(kind=wp) :: Ti =8000

Functions

function derivative(f, dx)

Arguments

Type IntentOptional AttributesName
real(kind=wp), intent(in), dimension(:):: f
real(kind=wp), intent(in), dimension(:):: dx

Return Value real(kind=wp), dimension(1:size(f))


Subroutines

subroutine writearray(fileunit, array)

Arguments

Type IntentOptional AttributesName
integer, intent(in) :: fileunit
real(kind=wp), intent(in), dimension(:):: array