55 integer,
public,
protected ::
rho_
58 integer,
allocatable,
public,
protected ::
mom(:)
61 integer,
public,
protected ::
e_
64 integer,
public,
protected ::
p_
67 integer,
public,
protected ::
te_
72 integer,
public,
protected ::
q_
81 double precision,
protected :: small_e
90 double precision,
public,
protected ::
h_ion_fr=1d0
93 double precision,
public,
protected ::
he_ion_fr=1d0
100 double precision,
public,
protected ::
rr=1d0
107 double precision :: gamma_1, inv_gamma_1
112 function fun_kin_en(w, ixI^L, ixO^L, inv_rho)
result(ke)
114 integer,
intent(in) :: ixi^
l, ixo^
l
115 double precision,
intent(in) :: w(ixi^s, nw)
116 double precision :: ke(ixo^s)
117 double precision,
intent(in),
optional :: inv_rho(ixo^s)
118 end function fun_kin_en
146 subroutine ffhd_read_params(files)
149 character(len=*),
intent(in) :: files(:)
157 do n = 1,
size(files)
158 open(
unitpar, file=trim(files(n)), status=
"old")
159 read(
unitpar, ffhd_list,
end=111)
162 end subroutine ffhd_read_params
165 subroutine ffhd_write_info(fh)
167 integer,
intent(in) :: fh
168 integer,
parameter :: n_par = 1
169 double precision :: values(n_par)
170 character(len=name_len) :: names(n_par)
171 integer,
dimension(MPI_STATUS_SIZE) :: st
174 call mpi_file_write(fh, n_par, 1, mpi_integer, st, er)
178 call mpi_file_write(fh, values, n_par, mpi_double_precision, st, er)
179 call mpi_file_write(fh, names, n_par * name_len, mpi_character, st, er)
180 end subroutine ffhd_write_info
199 if(
mype==0)
write(*,*)
'WARNING: set ffhd_thermal_conduction=F when ffhd_energy=F'
203 if(
mype==0)
write(*,*)
'WARNING: set ffhd_hyperbolic_thermal_conduction=F when ffhd_energy=F'
207 if(
mype==0)
write(*,*)
'WARNING: set ffhd_radiative_cooling=F when ffhd_energy=F'
211 if(
mype==0)
write(*,*)
'WARNING: set ffhd_trac=F when ffhd_energy=F'
215 if(
mype==0)
write(*,*)
'WARNING: set ffhd_partial_ionization=F when ffhd_energy=F'
221 if(
mype==0)
write(*,*)
'WARNING: set ffhd_partial_ionization=F when eq_state_units=F'
227 if(
mype==0)
write(*,*)
'WARNING: turn off parabolic TC when using hyperbolic TC'
230 physics_type =
"ffhd"
232 phys_internal_e=.false.
244 if(
mype==0)
write(*,*)
'WARNING: reset ffhd_trac_type=1 for 1D simulation'
249 if(
mype==0)
write(*,*)
'WARNING: set ffhd_trac_mask==bigdouble for global TRAC method'
253 allocate(start_indices(number_species),stop_indices(number_species))
260 mom(:) = var_set_momentum(1)
264 e_ = var_set_energy()
282 stop_indices(1)=nwflux
286 te_ = var_set_auxvar(
'Te',
'Te')
307 if(.not.
allocated(flux_type))
then
308 allocate(flux_type(
ndir, nwflux))
309 flux_type = flux_default
310 else if(any(shape(flux_type) /= [
ndir, nwflux]))
then
311 call mpistop(
"phys_check error: flux_type has wrong shape")
314 phys_get_dt => ffhd_get_dt
315 phys_get_cmax => ffhd_get_cmax_origin
316 phys_get_a2max => ffhd_get_a2max
317 phys_get_cs2max => ffhd_get_cs2max
318 phys_get_tcutoff => ffhd_get_tcutoff
319 phys_get_cbounds => ffhd_get_cbounds
320 phys_to_primitive => ffhd_to_primitive_origin
322 phys_to_conserved => ffhd_to_conserved_origin
324 phys_get_flux => ffhd_get_flux
325 phys_get_v => ffhd_get_v_origin
331 phys_add_source => ffhd_add_source
332 phys_check_params => ffhd_check_params
333 phys_write_info => ffhd_write_info
334 phys_handle_small_values => ffhd_handle_small_values_origin
335 ffhd_handle_small_values => ffhd_handle_small_values_origin
336 phys_check_w => ffhd_check_w_origin
339 phys_get_pthermal => ffhd_get_pthermal_iso
342 phys_get_pthermal => ffhd_get_pthermal_origin
348 ffhd_get_rfactor=>rfactor_from_temperature_ionization
349 phys_update_temperature => ffhd_update_temperature
353 ffhd_get_rfactor=>rfactor_from_constant_ionization
363 call ffhd_physical_units()
369 call mpistop(
"thermal conduction needs ffhd_energy=T")
372 call mpistop(
"hyperbolic thermal conduction needs ffhd_energy=T")
375 call mpistop(
"radiative cooling needs ffhd_energy=T")
385 call add_sts_method(ffhd_get_tc_dt_ffhd,ffhd_sts_set_source_tc_ffhd,
e_,1,
e_,1,.false.)
386 tc_fl%get_temperature_from_conserved => ffhd_get_temperature_from_etot
387 tc_fl%get_temperature_from_eint => ffhd_get_temperature_from_eint
402 rc_fl%get_var_Rfactor => ffhd_get_rfactor
405 rc_fl%has_equi = .false.
410 te_fl_ffhd%get_var_Rfactor => ffhd_get_rfactor
412 phys_te_images => ffhd_te_images
427 subroutine ffhd_te_images
432 case(
'EIvtiCCmpi',
'EIvtuCCmpi')
434 case(
'ESvtiCCmpi',
'ESvtuCCmpi')
436 case(
'SIvtiCCmpi',
'SIvtuCCmpi')
438 case(
'WIvtiCCmpi',
'WIvtuCCmpi')
441 call mpistop(
"Error in synthesize emission: Unknown convert_type")
443 end subroutine ffhd_te_images
446 subroutine ffhd_sts_set_source_tc_ffhd(ixI^L,ixO^L,w,x,wres,fix_conserve_at_step,my_dt,igrid,nflux)
450 integer,
intent(in) :: ixi^
l, ixo^
l, igrid, nflux
451 double precision,
intent(in) :: x(ixi^s,1:
ndim)
452 double precision,
intent(inout) :: wres(ixi^s,1:nw), w(ixi^s,1:nw)
453 double precision,
intent(in) :: my_dt
454 logical,
intent(in) :: fix_conserve_at_step
456 end subroutine ffhd_sts_set_source_tc_ffhd
458 function ffhd_get_tc_dt_ffhd(w,ixI^L,ixO^L,dx^D,x)
result(dtnew)
465 integer,
intent(in) :: ixi^
l, ixo^
l
466 double precision,
intent(in) ::
dx^
d, x(ixi^s,1:
ndim)
467 double precision,
intent(in) :: w(ixi^s,1:nw)
468 double precision :: dtnew
471 end function ffhd_get_tc_dt_ffhd
473 subroutine ffhd_tc_handle_small_e(w, x, ixI^L, ixO^L, step)
476 integer,
intent(in) :: ixi^
l,ixo^
l
477 double precision,
intent(inout) :: w(ixi^s,1:nw)
478 double precision,
intent(in) :: x(ixi^s,1:
ndim)
479 integer,
intent(in) :: step
480 character(len=140) :: error_msg
482 write(error_msg,
"(a,i3)")
"Thermal conduction step ", step
483 call ffhd_handle_small_ei(w,x,ixi^
l,ixo^
l,
e_,error_msg)
484 end subroutine ffhd_tc_handle_small_e
486 subroutine tc_params_read_ffhd(fl)
488 type(tc_fluid),
intent(inout) :: fl
491 logical :: tc_saturate=.false.
492 double precision :: tc_k_para=0d0
493 character(len=std_len) :: tc_slope_limiter=
"MC"
495 namelist /tc_list/ tc_saturate, tc_slope_limiter, tc_k_para
499 read(
unitpar, tc_list,
end=111)
503 fl%tc_saturate = tc_saturate
504 fl%tc_k_para = tc_k_para
505 select case(tc_slope_limiter)
507 fl%tc_slope_limiter = 0
510 fl%tc_slope_limiter = 1
513 fl%tc_slope_limiter = 2
516 fl%tc_slope_limiter = 3
519 fl%tc_slope_limiter = 4
521 call mpistop(
"Unknown tc_slope_limiter, choose MC, minmod")
523 end subroutine tc_params_read_ffhd
525 subroutine rc_params_read(fl)
528 type(rc_fluid),
intent(inout) :: fl
530 integer :: ncool = 4000
531 double precision :: cfrac=0.1d0
534 character(len=std_len) :: coolcurve=
'JCcorona'
537 character(len=std_len) :: coolmethod=
'exact'
540 logical :: tfix=.false.
546 logical :: rc_split=.false.
547 logical :: rad_cut=.false.
548 double precision :: rad_cut_hgt=0.5d0
549 double precision :: rad_cut_dey=0.15d0
551 namelist /rc_list/ coolcurve, coolmethod, ncool, cfrac, tlow, tfix, rc_split, rad_cut, rad_cut_hgt, rad_cut_dey
555 read(
unitpar, rc_list,
end=111)
560 fl%coolcurve=coolcurve
561 fl%coolmethod=coolmethod
567 fl%rad_cut_hgt=rad_cut_hgt
568 fl%rad_cut_dey=rad_cut_dey
569 end subroutine rc_params_read
571 subroutine ffhd_check_params
578 if (
ffhd_gamma <= 0.0d0)
call mpistop (
"Error: ffhd_gamma <= 0")
579 if (
ffhd_adiab < 0.0d0)
call mpistop (
"Error: ffhd_adiab < 0")
583 call mpistop (
"Error: ffhd_gamma <= 0 or ffhd_gamma == 1")
584 inv_gamma_1=1.d0/gamma_1
589 call mpistop(
"usr_set_equi_vars has to be implemented in the user file")
591 end subroutine ffhd_check_params
593 subroutine ffhd_physical_units()
595 double precision :: mp,kb
596 double precision :: a,b
683 end subroutine ffhd_physical_units
685 subroutine ffhd_check_w_origin(primitive,ixI^L,ixO^L,w,flag)
687 logical,
intent(in) :: primitive
688 integer,
intent(in) :: ixi^
l, ixo^
l
689 double precision,
intent(in) :: w(ixi^s,nw)
690 double precision :: tmp(ixi^s)
691 logical,
intent(inout) :: flag(ixi^s,1:nw)
701 where(tmp(ixo^s) < small_e) flag(ixo^s,
e_) = .true.
704 end subroutine ffhd_check_w_origin
706 subroutine ffhd_to_conserved_origin(ixI^L,ixO^L,w,x)
708 integer,
intent(in) :: ixi^
l, ixo^
l
709 double precision,
intent(inout) :: w(ixi^s, nw)
710 double precision,
intent(in) :: x(ixi^s, 1:
ndim)
711 double precision :: inv_gamma2(ixo^s)
715 w(ixo^s,
e_)=w(ixo^s,
p_)*inv_gamma_1+half*w(ixo^s,
mom(1))**2*w(ixo^s,
rho_)
718 end subroutine ffhd_to_conserved_origin
720 subroutine ffhd_to_primitive_origin(ixI^L,ixO^L,w,x)
722 integer,
intent(in) :: ixi^
l, ixo^
l
723 double precision,
intent(inout) :: w(ixi^s, nw)
724 double precision,
intent(in) :: x(ixi^s, 1:
ndim)
725 double precision :: inv_rho(ixo^s), gamma2(ixo^s)
729 call ffhd_handle_small_values(.false., w, x, ixi^
l, ixo^
l,
'ffhd_to_primitive_origin')
732 w(ixo^s,
mom(1)) = w(ixo^s,
mom(1))/w(ixo^s,
rho_)
734 w(ixo^s,
p_)=gamma_1*(w(ixo^s,
e_)-half*w(ixo^s,
rho_)*w(ixo^s,
mom(1))**2)
736 end subroutine ffhd_to_primitive_origin
740 integer,
intent(in) :: ixi^
l, ixo^
l
741 double precision,
intent(inout) :: w(ixi^s, nw)
742 double precision,
intent(in) :: x(ixi^s, 1:
ndim)
749 integer,
intent(in) :: ixi^
l, ixo^
l
750 double precision,
intent(inout) :: w(ixi^s, nw)
751 double precision,
intent(in) :: x(ixi^s, 1:
ndim)
755 call ffhd_handle_small_ei(w,x,ixi^
l,ixi^
l,
e_,
'ffhd_e_to_ei')
759 subroutine ffhd_handle_small_values_origin(primitive, w, x, ixI^L, ixO^L, subname)
762 logical,
intent(in) :: primitive
763 integer,
intent(in) :: ixi^
l,ixo^
l
764 double precision,
intent(inout) :: w(ixi^s,1:nw)
765 double precision,
intent(in) :: x(ixi^s,1:
ndim)
766 character(len=*),
intent(in) :: subname
768 logical :: flag(ixi^s,1:nw)
769 double precision :: tmp2(ixi^s)
771 call phys_check_w(primitive, ixi^
l, ixi^
l, w, flag)
778 where(flag(ixo^s,
rho_)) w(ixo^s,
mom(1)) = 0.0d0
784 where(flag(ixo^s,
e_))
801 if(.not.primitive)
then
810 end subroutine ffhd_handle_small_values_origin
812 subroutine ffhd_get_v_origin(w,x,ixI^L,ixO^L,v)
814 integer,
intent(in) :: ixi^
l, ixo^
l
815 double precision,
intent(in) :: w(ixi^s,nw), x(ixi^s,1:
ndim)
816 double precision,
intent(out) :: v(ixi^s,
ndir)
817 double precision :: rho(ixi^s)
821 rho(ixo^s)=1.d0/rho(ixo^s)
823 v(ixo^s,
ndir) = w(ixo^s,
mom(1))*
block%B0(ixo^s,idir,0)*rho(ixo^s)
825 end subroutine ffhd_get_v_origin
829 integer,
intent(in) :: ixi^
l, ixo^
l, idim
830 double precision,
intent(in) :: w(ixi^s,nw), x(ixi^s,1:
ndim)
831 double precision,
intent(out) :: v(ixi^s)
832 double precision :: rho(ixi^s)
835 v(ixo^s) = (w(ixo^s,
mom(1))*
block%B0(ixo^s,idim,0)) / rho(ixo^s)
838 subroutine ffhd_get_cmax_origin(w,x,ixI^L,ixO^L,idim,cmax)
840 integer,
intent(in) :: ixi^
l, ixo^
l, idim
842 double precision,
intent(in) :: w(ixi^s, nw), x(ixi^s,1:
ndim)
843 double precision,
intent(inout) :: cmax(ixi^s)
850 cmax(ixo^s)=abs(w(ixo^s,
mom(1))*
block%B0(ixo^s,idim,0))+cmax(ixo^s)
852 end subroutine ffhd_get_cmax_origin
854 subroutine ffhd_get_cs2max(w,x,ixI^L,ixO^L,cs2max)
856 integer,
intent(in) :: ixi^
l, ixo^
l
857 double precision,
intent(in) :: w(ixi^s, nw), x(ixi^s,1:
ndim)
858 double precision,
intent(inout) :: cs2max
859 double precision :: cs2(ixi^s)
862 cs2max=maxval(cs2(ixo^s))
863 end subroutine ffhd_get_cs2max
865 subroutine ffhd_get_a2max(w,x,ixI^L,ixO^L,a2max)
867 integer,
intent(in) :: ixi^
l, ixo^
l
868 double precision,
intent(in) :: w(ixi^s, nw), x(ixi^s,1:
ndim)
869 double precision,
intent(inout) :: a2max(
ndim)
870 double precision :: a2(ixi^s,
ndim,nw)
871 integer :: gxo^
l,hxo^
l,jxo^
l,kxo^
l,i,j
876 hxo^
l=ixo^
l-
kr(i,^
d);
877 gxo^
l=hxo^
l-
kr(i,^
d);
878 jxo^
l=ixo^
l+
kr(i,^
d);
879 kxo^
l=jxo^
l+
kr(i,^
d);
880 a2(ixo^s,i,1:nw)=abs(-w(kxo^s,1:nw)+16.d0*w(jxo^s,1:nw)&
881 -30.d0*w(ixo^s,1:nw)+16.d0*w(hxo^s,1:nw)-w(gxo^s,1:nw))
882 a2max(i)=maxval(a2(ixo^s,i,1:nw))/12.d0/
dxlevel(i)**2
884 end subroutine ffhd_get_a2max
886 subroutine ffhd_get_tcutoff(ixI^L,ixO^L,w,x,Tco_local,Tmax_local)
889 integer,
intent(in) :: ixi^
l,ixo^
l
890 double precision,
intent(in) :: x(ixi^s,1:
ndim)
891 double precision,
intent(inout) :: w(ixi^s,1:nw)
892 double precision,
intent(out) :: tco_local,tmax_local
893 double precision,
parameter :: trac_delta=0.25d0
894 double precision :: tmp1(ixi^s),te(ixi^s),lts(ixi^s)
895 double precision,
dimension(ixI^S,1:ndir) :: bunitvec
896 double precision,
dimension(ixI^S,1:ndim) :: gradt
897 double precision :: bdir(
ndim)
898 double precision :: ltrc,ltrp,altr(ixi^s)
899 integer :: idims,jxo^
l,hxo^
l,ixa^
d,ixb^
d
900 integer :: jxp^
l,hxp^
l,ixp^
l,ixq^
l
901 logical :: lrlt(ixi^s)
905 tmax_local=maxval(te(ixo^s))
915 lts(ixo^s)=0.5d0*abs(te(jxo^s)-te(hxo^s))/te(ixo^s)
917 where(lts(ixo^s) > trac_delta)
920 if(any(lrlt(ixo^s)))
then
921 tco_local=maxval(te(ixo^s), mask=lrlt(ixo^s))
932 lts(ixp^s)=0.5d0*abs(te(jxp^s)-te(hxp^s))/te(ixp^s)
933 lts(ixp^s)=max(one, (exp(lts(ixp^s))/ltrc)**ltrp)
934 lts(ixo^s)=0.25d0*(lts(jxo^s)+two*lts(ixo^s)+lts(hxo^s))
935 block%wextra(ixo^s,
tcoff_)=te(ixo^s)*lts(ixo^s)**0.4d0
937 call mpistop(
"ffhd_trac_type not allowed for 1D simulation")
953 gradt(ixo^s,idims)=tmp1(ixo^s)
955 bunitvec(ixo^s,:)=
block%B0(ixo^s,:,0)
960 ixb^
d=(ixomin^
d+ixomax^
d-1)/2+ixa^
d;
963 if(sum(bdir(:)**2) .gt. zero)
then
964 bdir(1:ndim)=bdir(1:ndim)/dsqrt(sum(bdir(:)**2))
966 block%special_values(3:ndim+2)=bdir(1:ndim)
968 tmp1(ixo^s)=dsqrt(sum(bunitvec(ixo^s,:)**2,dim=ndim+1))
969 where(tmp1(ixo^s)/=0.d0)
970 tmp1(ixo^s)=1.d0/tmp1(ixo^s)
972 tmp1(ixo^s)=bigdouble
976 bunitvec(ixo^s,idims)=bunitvec(ixo^s,idims)*tmp1(ixo^s)
979 lts(ixo^s)=abs(sum(gradt(ixo^s,1:ndim)*bunitvec(ixo^s,1:ndim),dim=ndim+1))/te(ixo^s)
981 if(slab_uniform)
then
982 lts(ixo^s)=minval(dxlevel)*lts(ixo^s)
984 lts(ixo^s)=minval(block%ds(ixo^s,:),dim=ndim+1)*lts(ixo^s)
987 where(lts(ixo^s) > trac_delta)
990 if(any(lrlt(ixo^s)))
then
991 block%special_values(1)=maxval(te(ixo^s), mask=lrlt(ixo^s))
993 block%special_values(1)=zero
995 block%special_values(2)=tmax_local
1013 call gradient(te,ixi^l,ixq^l,idims,gradt(ixi^s,idims))
1014 call gradientf(te,x,ixi^l,hxp^l,idims,gradt(ixi^s,idims),nghostcells,.true.)
1015 call gradientf(te,x,ixi^l,jxp^l,idims,gradt(ixi^s,idims),nghostcells,.false.)
1017 bunitvec(ixp^s,:)=block%B0(ixp^s,:,0)
1018 lts(ixp^s)=abs(sum(gradt(ixp^s,1:ndim)*bunitvec(ixp^s,1:ndim),dim=ndim+1))/te(ixp^s)
1019 if(slab_uniform)
then
1020 lts(ixp^s)=minval(dxlevel)*lts(ixp^s)
1022 lts(ixp^s)=minval(block%ds(ixp^s,:),dim=ndim+1)*lts(ixp^s)
1024 lts(ixp^s)=max(one, (exp(lts(ixp^s))/ltrc)**ltrp)
1029 hxo^l=ixp^l-kr(idims,^d);
1030 jxo^l=ixp^l+kr(idims,^d);
1031 altr(ixp^s)=altr(ixp^s)+0.25d0*(lts(hxo^s)+two*lts(ixp^s)+lts(jxo^s))*bunitvec(ixp^s,idims)**2
1033 block%wextra(ixp^s,
tcoff_)=te(ixp^s)*altr(ixp^s)**0.4d0
1037 call mpistop(
"unknown ffhd_trac_type")
1040 end subroutine ffhd_get_tcutoff
1042 subroutine ffhd_get_cbounds(wLC,wRC,wLp,wRp,x,ixI^L,ixO^L,idim,Hspeed,cmax,cmin)
1044 integer,
intent(in) :: ixi^
l, ixo^
l, idim
1045 double precision,
intent(in) :: wlc(ixi^s, nw), wrc(ixi^s, nw)
1046 double precision,
intent(in) :: wlp(ixi^s, nw), wrp(ixi^s, nw)
1047 double precision,
intent(in) :: x(ixi^s,1:
ndim)
1048 double precision,
intent(inout) :: cmax(ixi^s,1:number_species)
1049 double precision,
intent(inout),
optional :: cmin(ixi^s,1:number_species)
1050 double precision,
intent(in) :: hspeed(ixi^s,1:number_species)
1051 double precision :: wmean(ixi^s,nw)
1052 double precision,
dimension(ixI^S) :: umean, dmean, csoundl, csoundr, tmp1,tmp2,tmp3
1059 tmp1(ixo^s)=sqrt(wlp(ixo^s,
rho_))
1060 tmp2(ixo^s)=sqrt(wrp(ixo^s,
rho_))
1061 tmp3(ixo^s)=1.d0/(tmp1(ixo^s)+tmp2(ixo^s))
1062 umean(ixo^s)=(wlp(ixo^s,
mom(1))*
block%B0(ixo^s,idim,idim)*tmp1(ixo^s)&
1063 +wrp(ixo^s,
mom(1))*
block%B0(ixo^s,idim,idim)*tmp2(ixo^s))*tmp3(ixo^s)
1064 call ffhd_get_csound_prim(wlp,x,ixi^
l,ixo^
l,idim,csoundl)
1065 call ffhd_get_csound_prim(wrp,x,ixi^
l,ixo^
l,idim,csoundr)
1066 dmean(ixo^s)=(tmp1(ixo^s)*csoundl(ixo^s)+tmp2(ixo^s)*csoundr(ixo^s)) * &
1067 tmp3(ixo^s) + 0.5d0*tmp1(ixo^s)*tmp2(ixo^s)*tmp3(ixo^s)**2 * &
1068 (wrp(ixo^s,
mom(1))*
block%B0(ixo^s,idim,idim)-wlp(ixo^s,
mom(1))*
block%B0(ixo^s,idim,idim))**2
1069 dmean(ixo^s)=dsqrt(dmean(ixo^s))
1070 if(
present(cmin))
then
1071 cmin(ixo^s,1)=umean(ixo^s)-dmean(ixo^s)
1072 cmax(ixo^s,1)=umean(ixo^s)+dmean(ixo^s)
1074 cmax(ixo^s,1)=abs(umean(ixo^s))+dmean(ixo^s)
1077 wmean(ixo^s,1:nwflux)=0.5d0*(wlc(ixo^s,1:nwflux)+wrc(ixo^s,1:nwflux))
1078 tmp1(ixo^s)=wmean(ixo^s,
mom(1))*
block%B0(ixo^s,idim,idim)/wmean(ixo^s,
rho_)
1079 call ffhd_get_csound(wmean,x,ixi^
l,ixo^
l,idim,csoundr)
1080 if(
present(cmin))
then
1081 cmax(ixo^s,1)=max(tmp1(ixo^s)+csoundr(ixo^s),zero)
1082 cmin(ixo^s,1)=min(tmp1(ixo^s)-csoundr(ixo^s),zero)
1084 cmax(ixo^s,1)=abs(tmp1(ixo^s))+csoundr(ixo^s)
1088 call ffhd_get_csound_prim(wlp,x,ixi^
l,ixo^
l,idim,csoundl)
1089 call ffhd_get_csound_prim(wrp,x,ixi^
l,ixo^
l,idim,csoundr)
1090 csoundl(ixo^s)=max(csoundl(ixo^s),csoundr(ixo^s))
1091 if(
present(cmin))
then
1092 cmin(ixo^s,1)=min(wlp(ixo^s,
mom(1))*
block%B0(ixo^s,idim,idim),&
1093 wrp(ixo^s,
mom(1))*
block%B0(ixo^s,idim,idim))-csoundl(ixo^s)
1094 cmax(ixo^s,1)=max(wlp(ixo^s,
mom(1))*
block%B0(ixo^s,idim,idim),&
1095 wrp(ixo^s,
mom(1))*
block%B0(ixo^s,idim,idim))+csoundl(ixo^s)
1097 cmax(ixo^s,1)=max(wlp(ixo^s,
mom(1))*
block%B0(ixo^s,idim,idim),&
1098 wrp(ixo^s,
mom(1))*
block%B0(ixo^s,idim,idim))+csoundl(ixo^s)
1101 end subroutine ffhd_get_cbounds
1103 subroutine ffhd_get_csound(w,x,ixI^L,ixO^L,idim,csound)
1105 integer,
intent(in) :: ixi^
l, ixo^
l, idim
1106 double precision,
intent(in) :: w(ixi^s, nw), x(ixi^s,1:
ndim)
1107 double precision,
intent(out):: csound(ixi^s)
1110 csound(ixo^s) = dsqrt(csound(ixo^s))*abs(
block%B0(ixo^s,idim,idim))
1111 end subroutine ffhd_get_csound
1114 subroutine ffhd_get_csound_prim(w,x,ixI^L,ixO^L,idim,csound)
1117 integer,
intent(in) :: ixi^
l, ixo^
l, idim
1118 double precision,
intent(in) :: w(ixi^s, nw), x(ixi^s,1:
ndim)
1119 double precision,
intent(out):: csound(ixi^s)
1126 csound(ixo^s) = dsqrt(csound(ixo^s))
1127 end subroutine ffhd_get_csound_prim
1129 subroutine ffhd_get_pthermal_iso(w,x,ixI^L,ixO^L,pth)
1132 integer,
intent(in) :: ixi^
l, ixo^
l
1133 double precision,
intent(in) :: w(ixi^s,nw)
1134 double precision,
intent(in) :: x(ixi^s,1:
ndim)
1135 double precision,
intent(out):: pth(ixi^s)
1139 end subroutine ffhd_get_pthermal_iso
1141 subroutine ffhd_get_pthermal_origin(w,x,ixI^L,ixO^L,pth)
1144 integer,
intent(in) :: ixi^
l, ixo^
l
1145 double precision,
intent(in) :: w(ixi^s,nw)
1146 double precision,
intent(in) :: x(ixi^s,1:
ndim)
1147 double precision,
intent(out):: pth(ixi^s)
1152 {
do ix^db= ixo^lim^db\}
1157 elseif(check_small_values)
then
1158 {
do ix^db= ixo^lim^db\}
1159 if(pth(ix^d)<small_pressure)
then
1160 write(*,*)
"Error: small value of gas pressure",pth(ix^d),&
1161 " encountered when call ffhd_get_pthermal"
1162 write(*,*)
"Iteration: ", it,
" Time: ", global_time
1163 write(*,*)
"Location: ", x(ix^d,:)
1164 write(*,*)
"Cell number: ", ix^d
1166 write(*,*) trim(cons_wnames(iw)),
": ",w(ix^d,iw)
1169 if(trace_small_values)
write(*,*) sqrt(pth(ix^d)-bigdouble)
1170 write(*,*)
"Saving status at the previous time step"
1175 end subroutine ffhd_get_pthermal_origin
1177 subroutine ffhd_get_temperature_from_te(w, x, ixI^L, ixO^L, res)
1179 integer,
intent(in) :: ixi^
l, ixo^
l
1180 double precision,
intent(in) :: w(ixi^s, 1:nw)
1181 double precision,
intent(in) :: x(ixi^s, 1:
ndim)
1182 double precision,
intent(out):: res(ixi^s)
1184 res(ixo^s) = w(ixo^s,
te_)
1185 end subroutine ffhd_get_temperature_from_te
1187 subroutine ffhd_get_temperature_from_eint(w, x, ixI^L, ixO^L, res)
1189 integer,
intent(in) :: ixi^
l, ixo^
l
1190 double precision,
intent(in) :: w(ixi^s, 1:nw)
1191 double precision,
intent(in) :: x(ixi^s, 1:
ndim)
1192 double precision,
intent(out):: res(ixi^s)
1193 double precision :: r(ixi^s)
1195 call ffhd_get_rfactor(w,x,ixi^
l,ixo^
l,r)
1196 res(ixo^s) = gamma_1 * w(ixo^s,
e_)/(w(ixo^s,
rho_)*r(ixo^s))
1197 end subroutine ffhd_get_temperature_from_eint
1199 subroutine ffhd_get_temperature_from_etot(w, x, ixI^L, ixO^L, res)
1201 integer,
intent(in) :: ixi^
l, ixo^
l
1202 double precision,
intent(in) :: w(ixi^s, 1:nw)
1203 double precision,
intent(in) :: x(ixi^s, 1:
ndim)
1204 double precision,
intent(out):: res(ixi^s)
1206 double precision :: r(ixi^s)
1208 call ffhd_get_rfactor(w,x,ixi^
l,ixo^
l,r)
1210 res(ixo^s)=res(ixo^s)/(r(ixo^s)*w(ixo^s,
rho_))
1211 end subroutine ffhd_get_temperature_from_etot
1215 integer,
intent(in) :: ixi^
l, ixo^
l
1216 double precision,
intent(in) :: w(ixi^s,nw)
1217 double precision,
intent(in) :: x(ixi^s,1:
ndim)
1218 double precision,
intent(out) :: csound2(ixi^s)
1219 double precision :: rho(ixi^s)
1224 csound2(ixo^s)=
ffhd_gamma*csound2(ixo^s)/rho(ixo^s)
1230 subroutine ffhd_get_flux(wC,w,x,ixI^L,ixO^L,idim,f)
1233 integer,
intent(in) :: ixi^
l, ixo^
l, idim
1235 double precision,
intent(in) :: wc(ixi^s,nw)
1237 double precision,
intent(in) :: w(ixi^s,nw)
1238 double precision,
intent(in) :: x(ixi^s,1:
ndim)
1239 double precision,
intent(out) :: f(ixi^s,nwflux)
1240 double precision :: ptotal(ixo^s)
1241 double precision :: tmp(ixi^s)
1242 integer :: idirmin, iw, idir, jdir, kdir
1243 double precision,
dimension(ixI^S) :: te,tau,sigt
1248 ptotal(ixo^s)=w(ixo^s,
p_)
1254 f(ixo^s,
mom(1))=(wc(ixo^s,
mom(1))*w(ixo^s,
mom(1))+ptotal(ixo^s))*
block%B0(ixo^s,idim,idim)
1258 f(ixo^s,
e_)=w(ixo^s,
mom(1))*(wc(ixo^s,
e_)+ptotal(ixo^s))*
block%B0(ixo^s,idim,idim)
1260 f(ixo^s,
e_)=f(ixo^s,
e_)+w(ixo^s,
q_)*
block%B0(ixo^s,idim,idim)
1264 end subroutine ffhd_get_flux
1266 subroutine ffhd_add_source(qdt,dtfactor,ixI^L,ixO^L,wCT,wCTprim,w,x,qsourcesplit,active)
1271 integer,
intent(in) :: ixi^
l, ixo^
l
1272 double precision,
intent(in) :: qdt,dtfactor
1273 double precision,
intent(in) :: wct(ixi^s,1:nw),wctprim(ixi^s,1:nw), x(ixi^s,1:
ndim)
1274 double precision,
intent(inout) :: w(ixi^s,1:nw)
1275 logical,
intent(in) :: qsourcesplit
1276 logical,
intent(inout) :: active
1278 if (.not. qsourcesplit)
then
1280 call add_punitb(qdt,ixi^
l,ixo^
l,wct,w,x,wctprim)
1282 call add_hypertc_source(qdt,ixi^
l,ixo^
l,wct,w,x,wctprim)
1288 w,x,qsourcesplit,active,
rc_fl)
1303 if(.not.qsourcesplit)
then
1305 call ffhd_update_temperature(ixi^
l,ixo^
l,wct,w,x)
1308 end subroutine ffhd_add_source
1310 subroutine add_punitb(qdt,ixI^L,ixO^L,wCT,w,x,wCTprim)
1313 integer,
intent(in) :: ixi^
l,ixo^
l
1314 double precision,
intent(in) :: qdt
1315 double precision,
intent(in) :: wct(ixi^s,1:nw),x(ixi^s,1:
ndim)
1316 double precision,
intent(in) :: wctprim(ixi^s,1:nw)
1317 double precision,
intent(inout) :: w(ixi^s,1:nw)
1319 integer :: idims,hxo^
l
1320 double precision :: divb(ixi^s)
1325 hxo^
l=ixo^
l-
kr(idims,^
d);
1326 divb(ixo^s)=divb(ixo^s)+(
block%B0(ixo^s,idims,idims)-
block%B0(hxo^s,idims,idims))/
dxlevel(idims)
1331 w(ixo^s,
mom(1))=w(ixo^s,
mom(1))+qdt*wctprim(ixo^s,
p_)*divb(ixo^s)
1332 end subroutine add_punitb
1336 integer,
intent(in) :: ixi^
l, ixo^
l
1337 double precision,
intent(in) :: w(ixi^s,1:nw),x(ixi^s,1:
ndim)
1338 double precision,
intent(out) :: rho(ixi^s)
1340 rho(ixo^s) = w(ixo^s,
rho_)
1343 subroutine ffhd_handle_small_ei(w, x, ixI^L, ixO^L, ie, subname)
1346 integer,
intent(in) :: ixi^
l,ixo^
l, ie
1347 double precision,
intent(inout) :: w(ixi^s,1:nw)
1348 double precision,
intent(in) :: x(ixi^s,1:
ndim)
1349 character(len=*),
intent(in) :: subname
1351 logical :: flag(ixi^s,1:nw)
1352 double precision :: rho(ixi^s)
1355 where(w(ixo^s,ie)<small_e) flag(ixo^s,ie)=.true.
1356 if(any(flag(ixo^s,ie)))
then
1359 where(flag(ixo^s,ie)) w(ixo^s,ie)=small_e
1363 w(ixo^s,
e_)=w(ixo^s,
e_)*gamma_1
1365 w(ixo^s,
mom(1)) = w(ixo^s,
mom(1))/rho(ixo^s)
1369 end subroutine ffhd_handle_small_ei
1371 subroutine ffhd_update_temperature(ixI^L,ixO^L,wCT,w,x)
1374 integer,
intent(in) :: ixi^
l, ixo^
l
1375 double precision,
intent(in) :: wct(ixi^s,1:nw), x(ixi^s,1:
ndim)
1376 double precision,
intent(inout) :: w(ixi^s,1:nw)
1377 double precision :: iz_h(ixo^s),iz_he(ixo^s), pth(ixi^s)
1383 end subroutine ffhd_update_temperature
1385 subroutine ffhd_get_dt(w,ixI^L,ixO^L,dtnew,dx^D,x)
1391 integer,
intent(in) :: ixi^
l, ixo^
l
1392 double precision,
intent(inout) :: dtnew
1393 double precision,
intent(in) ::
dx^
d
1394 double precision,
intent(in) :: w(ixi^s,1:nw)
1395 double precision,
intent(in) :: x(ixi^s,1:
ndim)
1396 integer :: idirmin,idim
1397 double precision :: dxarr(
ndim)
1398 double precision :: current(ixi^s,7-2*
ndir:3),eta(ixi^s)
1413 end subroutine ffhd_get_dt
1527 function ffhd_kin_en_origin(w, ixI^L, ixO^L, inv_rho)
result(ke)
1529 integer,
intent(in) :: ixi^
l, ixo^
l
1530 double precision,
intent(in) :: w(ixi^s, nw)
1531 double precision :: ke(ixo^s)
1532 double precision,
intent(in),
optional :: inv_rho(ixo^s)
1534 if(
present(inv_rho))
then
1535 ke(ixo^s)=0.5d0*w(ixo^s,
mom(1))**2*inv_rho(ixo^s)
1537 ke(ixo^s)=0.5d0*w(ixo^s,
mom(1))**2/w(ixo^s,
rho_)
1539 end function ffhd_kin_en_origin
1541 subroutine rfactor_from_temperature_ionization(w,x,ixI^L,ixO^L,Rfactor)
1544 integer,
intent(in) :: ixi^
l, ixo^
l
1545 double precision,
intent(in) :: w(ixi^s,1:nw)
1546 double precision,
intent(in) :: x(ixi^s,1:
ndim)
1547 double precision,
intent(out):: rfactor(ixi^s)
1548 double precision :: iz_h(ixo^s),iz_he(ixo^s)
1551 rfactor(ixo^s)=(1.d0+iz_h(ixo^s)+0.1d0*(1.d0+iz_he(ixo^s)*(1.d0+iz_he(ixo^s))))/(2.d0+3.d0*
he_abundance)
1552 end subroutine rfactor_from_temperature_ionization
1554 subroutine rfactor_from_constant_ionization(w,x,ixI^L,ixO^L,Rfactor)
1556 integer,
intent(in) :: ixi^
l, ixo^
l
1557 double precision,
intent(in) :: w(ixi^s,1:nw)
1558 double precision,
intent(in) :: x(ixi^s,1:
ndim)
1559 double precision,
intent(out):: rfactor(ixi^s)
1562 end subroutine rfactor_from_constant_ionization
1564 subroutine get_tau(ixI^L,ixO^L,w,Te,tau,sigT5)
1566 integer,
intent(in) :: ixi^
l, ixo^
l
1567 double precision,
dimension(ixI^S,1:nw),
intent(in) :: w
1568 double precision,
dimension(ixI^S),
intent(in) :: te
1569 double precision,
dimension(ixI^S),
intent(out) :: tau,sigt5
1571 double precision :: dxmin,taumin
1572 double precision,
dimension(ixI^S) :: sigt7,eint
1579 sigt7(ixo^s)=sigt5(ixo^s)*
block%wextra(ixo^s,
tcoff_)
1582 sigt7(ixo^s)=sigt5(ixo^s)*te(ixo^s)
1586 sigt7(ixo^s)=sigt5(ixo^s)*te(ixo^s)
1590 end subroutine get_tau
1592 subroutine add_hypertc_source(qdt,ixI^L,ixO^L,wCT,w,x,wCTprim)
1594 integer,
intent(in) :: ixi^
l,ixo^
l
1595 double precision,
intent(in) :: qdt
1596 double precision,
dimension(ixI^S,1:ndim),
intent(in) :: x
1597 double precision,
dimension(ixI^S,1:nw),
intent(in) :: wct,wctprim
1598 double precision,
dimension(ixI^S,1:nw),
intent(inout) :: w
1600 integer :: hxc^
l,hxo^
l,ixc^
l,jxc^
l,jxo^
l,kxc^
l
1601 double precision :: invdx
1602 double precision,
dimension(ixI^S) :: te,tau,sigt,htc_qsrc,tface
1603 double precision,
dimension(ixI^S) :: htc_esrc
1605 te(ixi^s)=wctprim(ixi^s,
p_)/wct(ixi^s,
rho_)
1606 call get_tau(ixi^
l,ixo^
l,wctprim,te,tau,sigt)
1611 ixcmin^
d=ixomin^
d-
kr(idims,^
d);ixcmax^
d=ixomax^
d;
1612 jxc^
l=ixc^
l+
kr(idims,^
d);
1613 kxc^
l=jxc^
l+
kr(idims,^
d);
1614 hxc^
l=ixc^
l-
kr(idims,^
d);
1615 hxo^
l=ixo^
l-
kr(idims,^
d);
1616 tface(ixc^s)=(7.d0*(te(ixc^s)+te(jxc^s))-(te(hxc^s)+te(kxc^s)))/12.d0
1617 htc_qsrc(ixo^s)=htc_qsrc(ixo^s)+sigt(ixo^s)*
block%B0(ixo^s,idims,0)*(tface(ixo^s)-tface(hxo^s))*invdx
1619 htc_qsrc(ixo^s)=(htc_qsrc(ixo^s)+wct(ixo^s,
q_))/tau(ixo^s)
1620 w(ixo^s,
q_)=w(ixo^s,
q_)-qdt*htc_qsrc(ixo^s)
1621 end subroutine add_hypertc_source
subroutine, public mpistop(message)
Exit MPI-AMRVAC with an error message.
Module for physical and numeric constants.
double precision, parameter bigdouble
A very large real number.
subroutine add_convert_method(phys_convert_vars, nwc, dataset_names, file_suffix)
Frozen-field hydrodynamics module.
integer, public, protected te_
Indices of temperature.
integer, public, protected ffhd_trac_type
Which TRAC method is used.
double precision, public hypertc_kappa
The thermal conductivity kappa in hyperbolic thermal conduction.
integer, public, protected e_
Index of the energy density (-1 if not present)
double precision, public ffhd_gamma
The adiabatic index.
logical, public, protected eq_state_units
double precision, public ffhd_adiab
The adiabatic constant.
procedure(sub_get_pthermal), pointer, public ffhd_get_temperature
logical, public, protected ffhd_hyperbolic_thermal_conduction
Whether hyperbolic type thermal conduction is used.
type(rc_fluid), allocatable, public rc_fl
type of fluid for radiative cooling
integer, public, protected rho_
Index of the density (in the w array)
integer, public, protected tcoff_
Index of the cutoff temperature for the TRAC method.
procedure(sub_get_pthermal), pointer, public ffhd_get_pthermal
subroutine, public ffhd_get_rho(w, x, ixil, ixol, rho)
logical, public, protected ffhd_partial_ionization
Whether plasma is partially ionized.
procedure(sub_convert), pointer, public ffhd_to_conserved
integer, dimension(:), allocatable, public, protected mom
Indices of the momentum density.
double precision, public, protected h_ion_fr
Ionization fraction of H H_ion_fr = H+/(H+ + H)
procedure(fun_kin_en), pointer, public ffhd_kin_en
subroutine, public ffhd_phys_init()
procedure(sub_get_v), pointer, public ffhd_get_v
subroutine, public ffhd_get_csound2(w, x, ixil, ixol, csound2)
logical, public, protected ffhd_energy
Whether an energy equation is used.
double precision, public, protected rr
type(tc_fluid), allocatable, public tc_fl
type of fluid for thermal conduction
type(te_fluid), allocatable, public te_fl_ffhd
type of fluid for thermal emission synthesis
double precision, public, protected he_abundance
Helium abundance over Hydrogen.
logical, public, protected ffhd_viscosity
Whether viscosity is added.
logical, public, protected ffhd_radiative_cooling
Whether radiative cooling is added.
subroutine, public ffhd_get_v_idim(w, x, ixil, ixol, idim, v)
integer, public, protected q_
procedure(sub_convert), pointer, public ffhd_to_primitive
integer, public, protected tweight_
double precision, public, protected he_ion_fr2
Ratio of number He2+ / number He+ + He2+ He_ion_fr2 = He2+/(He2+ + He+)
logical, public, protected ffhd_trac
Whether TRAC method is used.
logical, public, protected ffhd_thermal_conduction
Whether thermal conduction is used.
subroutine, public ffhd_ei_to_e(ixil, ixol, w, x)
logical, public, protected ffhd_gravity
Whether gravity is added.
integer, public, protected p_
Index of the gas pressure (-1 if not present) should equal e_.
double precision, public, protected ffhd_trac_mask
Height of the mask used in the TRAC method.
double precision, public, protected he_ion_fr
Ionization fraction of He He_ion_fr = (He2+ + He+)/(He2+ + He+ + He)
integer, public, protected ffhd_trac_finegrid
Distance between two adjacent traced magnetic field lines (in finest cell size)
subroutine, public ffhd_e_to_ei(ixil, ixol, w, x)
Module for flux conservation near refinement boundaries.
Module with geometry-related routines (e.g., divergence, curl)
subroutine divvector(qvec, ixil, ixol, divq, nth_in)
subroutine gradient(q, ixil, ixol, idir, gradq, nth_in)
This module contains definitions of global parameters and variables and some generic functions/subrou...
type(state), pointer block
Block pointer for using one block and its previous state.
double precision small_pressure
double precision unit_time
Physical scaling factor for time.
double precision unit_density
Physical scaling factor for density.
integer, parameter unitpar
file handle for IO
double precision unit_mass
Physical scaling factor for mass.
integer, dimension(3, 3) kr
Kronecker delta tensor.
double precision phys_trac_mask
double precision unit_numberdensity
Physical scaling factor for number density.
character(len=std_len) convert_type
Which format to use when converting.
double precision unit_pressure
Physical scaling factor for pressure.
integer, parameter ndim
Number of spatial dimensions for grid variables.
double precision unit_length
Physical scaling factor for length.
character(len=std_len), dimension(:), allocatable par_files
Which par files are used as input.
integer mype
The rank of the current MPI task.
double precision, dimension(:), allocatable, parameter d
double precision dt
global time step
integer ndir
Number of spatial dimensions (components) for vector variables.
double precision unit_velocity
Physical scaling factor for velocity.
logical b0field
split magnetic field as background B0 field
logical need_global_cs2max
global value for csound speed
double precision unit_temperature
Physical scaling factor for temperature.
logical si_unit
Use SI units (.true.) or use cgs units (.false.)
double precision, dimension(:,:), allocatable dx
logical phys_trac
Use TRAC for MHD or 1D HD.
logical fix_small_values
fix small values with average or replace methods
double precision, dimension(^nd) dxlevel
store unstretched cell size of current level
double precision cs2max_global
global largest cs2 for hyperbolic thermal conduction
logical slab_uniform
uniform Cartesian geometry or not (stretched Cartesian)
double precision small_density
integer boundspeed
bound (left/min and right.max) speed of Riemann fan
integer phys_trac_finegrid
integer, parameter unitconvert
integer number_equi_vars
number of equilibrium set variables, besides the mag field
Module for including gravity in (magneto)hydrodynamics simulations.
subroutine gravity_get_dt(w, ixil, ixol, dtnew, dxd, x)
subroutine gravity_init()
Initialize the module.
subroutine gravity_add_source(qdt, ixil, ixol, wct, wctprim, w, x, energy, rhov, qsourcesplit, active)
w[iw]=w[iw]+qdt*S[wCT,qtC,x] where S is the source based on wCT within ixO
module ionization degree - get ionization degree for given temperature
subroutine ionization_degree_from_temperature(ixil, ixol, te, iz_h, iz_he)
subroutine ionization_degree_init()
Module containing all the particle routines.
This module defines the procedures of a physics module. It contains function pointers for the various...
module radiative cooling – add optically thin radiative cooling for HD and MHD
subroutine radiative_cooling_init_params(phys_gamma, he_abund)
Radiative cooling initialization.
subroutine cooling_get_dt(w, ixil, ixol, dtnew, dxd, x, fl)
subroutine radiative_cooling_init(fl, read_params)
subroutine radiative_cooling_add_source(qdt, ixil, ixol, wct, wctprim, w, x, qsourcesplit, active, fl)
Module for handling problematic values in simulations, such as negative pressures.
subroutine, public small_values_average(ixil, ixol, w, x, w_flag, windex)
logical, public trace_small_values
trace small values in the source file using traceback flag of compiler
subroutine, public small_values_error(wprim, x, ixil, ixol, w_flag, subname)
logical, dimension(:), allocatable, public small_values_fix_iw
Whether to apply small value fixes to certain variables.
character(len=20), public small_values_method
How to handle small values.
Generic supertimestepping method 1) in amrvac.par in sts_list set the following parameters which have...
subroutine, public add_sts_method(sts_getdt, sts_set_sources, startvar, nflux, startwbc, nwbc, evolve_b)
subroutine which added programatically a term to be calculated using STS Params: sts_getdt function c...
subroutine, public set_conversion_methods_to_head(sts_before_first_cycle, sts_after_last_cycle)
Set the hooks called before the first cycle and after the last cycle in the STS update This method sh...
subroutine, public set_error_handling_to_head(sts_error_handling)
Set the hook of error handling in the STS update. This method is called before updating the BC....
subroutine, public sts_init()
Initialize sts module.
Thermal conduction for HD and MHD or RHD and RMHD or twofl (plasma-neutral) module Adaptation of mod_...
subroutine, public tc_get_hd_params(fl, read_hd_params)
Init TC coefficients: HD case.
double precision function, public get_tc_dt_mhd(w, ixil, ixol, dxd, x, fl)
Get the explicut timestep for the TC (mhd implementation)
subroutine tc_init_params(phys_gamma)
subroutine, public sts_set_source_tc_mhd(ixil, ixol, w, x, wres, fix_conserve_at_step, my_dt, igrid, nflux, fl)
anisotropic thermal conduction with slope limited symmetric scheme Sharma 2007 Journal of Computation...
subroutine get_euv_image(qunit, fl)
subroutine get_sxr_image(qunit, fl)
subroutine get_euv_spectrum(qunit, fl)
subroutine get_whitelight_image(qunit, fl)
Module with all the methods that users can customize in AMRVAC.
procedure(rfactor), pointer usr_rfactor
procedure(set_equi_vars), pointer usr_set_equi_vars
The module add viscous source terms and check time step.
subroutine viscosity_add_source(qdt, ixil, ixol, wct, w, x, energy, qsourcesplit, active)
subroutine viscosity_init(phys_wider_stencil)
Initialize the module.
subroutine viscosity_get_dt(w, ixil, ixol, dtnew, dxd, x)