! $Id$
!
! IO of init and run parameters. Subroutines here are `at the end of the
! food chain', i.e. depend on all physics modules plus possibly others.
! Using this module is also a compact way of referring to all physics
! modules at once.
!
module Param_IO
!
use Cdata
use Chemistry
use Chiral
use Cosmicray
use CosmicrayFlux
use Cparam
! use Conductivity
use Density
use Detonate
use Dustdensity
use Dustvelocity
use Energy
use EquationOfState
use Forcing
use General
use Gravity
use Heatflux
use Hydro
use ImplicitDiffusion
use InitialCondition
use Interstellar
use IO, only: lcollective_IO, IO_strategy
use Lorenz_gauge
use Magnetic
use Messages
use NeutralDensity
use NeutralVelocity
use NSCBC
use Opacity, only: read_opacity_run_pars, write_opacity_run_pars
use Poisson
use PointMasses
use Polymer
use Power_spectrum
use Pscalar
use Ascalar
use Radiation
use Selfgravity
use Shear
use Signal_handling
use Shock
use Solid_Cells
use Special
use Streamlines
use Sub
use Testfield
use Testflow
use TestPerturb
use Testscalar
use Timeavg
use Viscosity
!
implicit none
!
private
!
public :: get_datadir, get_snapdir
public :: read_all_init_pars, read_all_run_pars
public :: write_all_init_pars, write_all_run_pars
public :: write_pencil_info
public :: get_downpars
!
logical :: lforce_shear_bc = .true.
!
! local quantities
!
real, dimension(mcom) :: fbcx1=0., fbcx2=0., fbcx1_2=0., fbcx2_2=0., &
fbcy1=0., fbcy2=0., fbcy1_1=0., fbcy1_2=0., fbcy2_1=0., fbcy2_2=0., &
fbcz1=0., fbcz2=0., fbcz1_1=0., fbcz1_2=0., fbcz2_1=0., fbcz2_2=0.
integer :: niter_poisson ! dummy
!
namelist /init_pars/ &
cvsid, ip, xyz0, xyz1, Lxyz, lperi, lshift_origin, lshift_origin_lower,&
xyz_units, wav1, wav1z, coord_system, lpole, ncoarse, lfix_unit_std, &
lequidist, coeff_grid, zeta_grid0, grid_func, xyz_star, lwrite_ic, lwrite_avg1d_binary, &
lnowrite, luniform_z_mesh_aspect_ratio, unit_system, unit_length, &
lmodify,modify_filename, dvid, ldivu_perp, &
unit_velocity, unit_density, unit_temperature, unit_magnetic, c_light, &
G_Newton, hbar, random_gen, seed0, lseed_global, lseed_procdependent, nfilter, lserial_io, der2_type, &
lread_oldsnap, lwrite_var_anyway, lwrite_last_powersnap, &
lread_oldsnap_nomag, lread_oldsnap_nopscalar, lread_oldsnap_notestflow, lread_oldsnap_noisothmhd, &
lread_oldsnap_notestfield, lread_oldsnap_notestscalar, lread_oldsnap_noshear, &
lread_oldsnap_nohydro, lread_oldsnap_nohydro_nomu5, &
lread_oldsnap_nohydro_efield, lread_oldsnap_nohydro_ekfield, &
lread_oldsnap_onlyA, &
ireset_tstart, tstart, lghostfold_usebspline, &
lread_aux, lwrite_aux, lkinflow_as_aux, lenforce_maux_check, &
pretend_lnTT, lprocz_slowest, lzorder, &
lcopysnapshots_exp, bcx, bcy, bcz, r_int, r_ext, r_ref, rsmooth, &
r_int_border, r_ext_border, mu0, force_lower_bound, force_upper_bound, &
lseparate_persist, ldistribute_persist, lpersist, lomit_add_data, &
fbcx1, fbcx2, fbcx1_2, fbcx2_2, &
fbcy1, fbcy2, fbcy1_1, fbcy1_2, fbcy2_1, fbcy2_2, &
fbcz1, fbcz2, fbcz1_1, fbcz1_2, fbcz2_1, fbcz2_2, &
fbcx_bot, fbcx_top, fbcy_bot, fbcy_top, fbcz_bot, fbcz_top, bc_slc_dir, &
vel_spec, mag_spec, &
uxy_spec, bxy_spec, jxbxy_spec, xy_spec, oo_spec, &
uxj_spec, vec_spec, ou_spec, oun_spec, ab_spec, azbz_spec, uzs_spec, ub_spec, &
bb2_spec, jj2_spec, ele_spec, a0_spec, pot_spec, &
Lor_spec, EMF_spec, Tra_spec, GWs_spec, GWh_spec, GWm_spec, Str_spec, Stg_spec, &
SCL_spec, VCT_spec, Tpq_spec, TGW_spec, GWd_spec, GWe_spec, GWf_spec, GWg_spec, &
StT_spec, StX_spec, &
vel_phispec, mag_phispec, &
uxj_phispec, vec_phispec, ou_phispec, ab_phispec, EP_spec, ro_spec, &
nd_spec, ud_spec, ux_spec, uy_spec, uz_spec, phiuu_spec, &
TT_spec, ss_spec, cc_spec, cr_spec, mu_spec, sp_spec, ssp_spec, sssp_spec, &
isaveglobal, lr_spec, r2u_spec, &
np_spec, np_ap_spec, rhop_spec, &
r3u_spec, rhocc_pdf, cc_pdf, lncc_pdf, gcc_pdf, lngcc_pdf, &
xyz_step, xi_step_frac, xi_step_width, dxi_fact, trans_width, &
lcylinder_in_a_box, lsphere_in_a_box, llocal_iso, init_loops, lwrite_2d, &
lcylindrical_gravity, &
border_frac_x, border_frac_y, border_frac_z, lborder_hyper_diff, &
luse_latitude, lshift_datacube_x, lfargo_advection, yequator, lequatory, &
lequatorz, zequator, lav_smallx, xav_max, niter_poisson, &
lforce_shear_bc,lread_from_other_prec, &
pipe_func, glnCrossSec0, CrossSec_x1, CrossSec_x2, CrossSec_w,&
lcorotational_frame, rcorot, lproper_averages, &
ldirect_access, ltolerate_namelist_errors, &
lyinyang, cyinyang_intpol_type, yy_biquad_weights, &
lcutoff_corners, nycut, nzcut, rel_dang, &
sigmaSB_set, c_light_set, cp_set, k_B_set, m_u_set, &
lnoghost_strati, ichannel1, ichannel2, tag_foreign, &
lfatal_num_vector_369, density_scale_factor
!
namelist /run_pars/ &
cvsid, ip, xyz0, xyz1, Lxyz, lperi, lpole, ncoarse, &
lshift_origin, lshift_origin_lower, coord_system, &
nt, it1, it1start, it1d, it_rmv, dt, cdt, ddt, dt_incr, &
lfractional_tstep_advance, lfractional_tstep_negative, &
cdtv, cdtv2, cdtv3, cdtsrc, cdts, cdtr, &
cdtc, isave, itorder, dsnap, dsnap_down, mvar_down, maux_down, &
d1davg, d2davg, dvid, dsound, dtmin, dspec, tmax, toutoff, &
iwig, ldivu_perp, &
dtracers, dfixed_points, unit_system, unit_length, &
unit_velocity, unit_density, unit_temperature, unit_magnetic, &
awig, ialive, max_walltime, dtmax, ldt_paronly, &
lspec_start, lspec_at_tplusdt, vel_spec, mag_spec, &
uxy_spec, bxy_spec, jxbxy_spec, xy_spec, oo_spec, &
uxj_spec, vec_spec, ou_spec, oun_spec, ab_spec, azbz_spec, uzs_spec, ub_spec, &
bb2_spec, jj2_spec, ele_spec, a0_spec, pot_spec, &
Lor_spec, EMF_spec, Tra_spec, GWs_spec, GWh_spec, GWm_spec, Str_spec, Stg_spec, &
SCL_spec, VCT_spec, Tpq_spec, TGW_spec, GWd_spec, GWe_spec, GWf_spec, GWg_spec, &
StT_spec, StX_spec, &
vel_phispec, mag_phispec, &
uxj_phispec, vec_phispec, ou_phispec, ab_phispec, EP_spec, ro_spec, &
nd_spec, ud_spec, ux_spec, uy_spec, uz_spec, phiuu_spec, &
TT_spec, ss_spec, cc_spec, cr_spec, mu_spec, sp_spec, ssp_spec, sssp_spec, &
isaveglobal, lr_spec, r2u_spec, &
np_spec, np_ap_spec, rhop_spec, &
r3u_spec, rhocc_pdf, cc_pdf, lncc_pdf, gcc_pdf, lngcc_pdf, &
lnspecial_pdf, special_pdf, ang_jb_pdf1d, ang_ub_pdf1d, &
kinflow, ladv_der_as_aux, lkinflow_as_aux, &
ampl_kinflow_x, ampl_kinflow_y, ampl_kinflow_z, &
kx_kinflow, ky_kinflow, kz_kinflow, dtphase_kinflow, &
random_gen, der2_type, lrmwig_rho, lrmwig_full, lrmwig_xyaverage, &
lnowrite, noghost_for_isave, nghost_read_fewer, &
test_nonblocking, lwrite_tracers, lwrite_fsum, lwrite_fixed_points, lwrite_ts_hdf5, &
lread_oldsnap_lnrho2rho, lread_oldsnap_nomag, lread_oldsnap_notestflow, lread_oldsnap_nopscalar, &
lread_oldsnap_notestfield, lread_oldsnap_notestscalar, lread_oldsnap_noshear, lrepair_snap, &
lread_oldsnap_nohydro, lread_oldsnap_nohydro_efield, lread_oldsnap_nohydro_ekfield, &
lread_oldsnap_noisothmhd, lread_oldsnap_onlyA, &
lread_oldsnap_rho2lnrho, lread_oldsnap_nosink, lwrite_dim_again, lwrite_last_powersnap, &
lread_aux, comment_char, ix, iy, iy2, iz, iz2, iz3, iz4, slice_position, &
xbot_slice, xtop_slice, ybot_slice, ytop_slice, zbot_slice, ztop_slice, &
r_rslice, nth_rslice, nph_rslice, &
bcx, bcy, bcz, r_int, r_ext, r_int_border, &
r_ext_border, lfreeze_varsquare, lfreeze_varint, lfreeze_varext, &
xfreeze_square, yfreeze_square, rfreeze_int, rfreeze_ext, wfreeze, &
wfreeze_int, wfreeze_ext, wborder, wborder_int, wborder_ext, tborder, &
luse_oldgrid, luse_xyz1, fshift_int, fshift_ext, lpersist, &
ireset_tstart, tstart, lseparate_persist, ldistribute_persist, lomit_add_data, &
fbcx1, fbcx2, fbcx1_2, fbcx2_2, &
fbcy1, fbcy2, fbcy1_1, fbcy1_2, fbcy2_1, fbcy2_2, &
fbcz1, fbcz2, fbcz1_1, fbcz1_2, fbcz2_1, fbcz2_2, &
fbcx_bot, fbcx_top, fbcy_bot, fbcy_top, fbcz_bot, fbcz_top, &
Udrift_bc, ttransient, tavg, idx_tavg, lserial_io, nr_directions, &
lsfu, lsfb, lsfz1, lsfz2, lsfflux, lpdfu, lpdfb, lpdfz1, lpdfz2, &
onedall, pretend_lnTT, old_cdtv, lmaxadvec_sum, save_lastsnap, &
lwrite_aux, lwrite_dvar, force_lower_bound, force_upper_bound, &
lenforce_maux_check, nv1_capitalvar, mskipvar, lread_oldsnap_mskipvar, &
oned, twod, lpoint, mpoint, npoint, lpoint2, mpoint2, npoint2, &
iproc_pt, iproc_p2, &
border_frac_x, border_frac_y, border_frac_z, lborder_hyper_diff, &
lcylinder_in_a_box, lsphere_in_a_box, ipencil_swap, &
lpencil_requested_swap, lpencil_diagnos_swap, lpencil_check, &
lpencil_check_small, lpencil_check_no_zeros, lpencil_check_diagnos_opti, &
lpencil_init, penc0, lwrite_2d, lbidiagonal_derij, lisotropic_advection, &
crash_file_dtmin_factor, ltestperturb, eps_rkf, &
eps_stiff, timestep_scaling, lequatory, lequatorz, zequator, &
lini_t_eq_zero, lini_t_eq_zero_once, &
lav_smallx, xav_max, ldt_paronly, lweno_transport, &
it_timing, har_spec, hav_spec, j_spec, jb_spec, ja_spec, b2_spec, &
lread_less, lread_nogrid, lformat, ltec, &
llsode, lsplit_second, nu_sts, permute_sts, lfargo_advection, &
ldynamical_diffusion, ldyndiff_useumax, re_mesh, lghostfold_usebspline, &
lreset_seed, loutput_varn_at_exact_tsnap, lstop_on_ioerror, mailaddress, &
theta_lower_border, wborder_theta_lower, theta_upper_border, &
wborder_theta_upper, fraction_tborder, lmeridional_border_drive, &
lread_from_other_prec, downsampl, lfullvar_in_slices, ivar_omit, &
lsubstract_reference_state, lzaver_on_input, &
ldirect_access, lproper_averages, lmaximal_cdt, lmaximal_cdtv, &
pipe_func, glnCrossSec0, CrossSec_x1, CrossSec_x2, CrossSec_w, &
cyinyang_intpol_type, yy_biquad_weights, lcutoff_corners, nycut, nzcut, rel_dang, &
lignore_nonequi, tag_foreign, tau_aver1, fmt_avgs, &
ou_omega, ou_kzspec, ab_kzspec, cor_uu, ou_polar, ab_polar, jb_polar, &
uut_spec, uut_polar, ouout_spec, ouout2_spec, ouout_polar, out_spec, uot_spec, &
saffman_ub, saffman_mag, saffman_mag_c, saffman_aa, saffman_aa_c, saffman_bb, &
uu_fft3d, oo_fft3d, bb_fft3d, jj_fft3d, uu_xkyz, oo_xkyz, bb_xkyz, jj_xkyz, &
uu_kx0z, oo_kx0z, bb_kx0z, jj_kx0z, bb_k00z, ee_k00z, gwT_fft3d, &
mh_specflux, mE_specflux, density_scale_factor
!
namelist /IO_pars/ &
lcollective_IO, IO_strategy
!
contains
!***********************************************************************
subroutine get_datadir(dir)
!
! Overwrite datadir from datadir.in, if that exists.
!
! 2-oct-02/wolf: coded
! 25-oct-02/axel: default is taken from cdata.f90 where it's defined
! 14-jan-15/MR : corrected call of mpibcast_char
!
use Mpicomm, only: mpibcast_logical, mpibcast_char,MPI_COMM_WORLD
!
character (len=*) :: dir
logical :: exists
integer, parameter :: unit=1
!
! Let root processor check for existence of datadir.in.
!
if (lroot) then
inquire(FILE='datadir.in',EXIST=exists)
if (exists) then
open(unit,FILE='datadir.in',FORM='formatted')
read(unit,'(A)') dir
close(unit)
endif
endif
!
! Tell other processors whether we need to communicate dir (i.e. datadir).
!
call mpibcast_logical(exists,comm=MPI_COMM_WORLD)
!
! Let root processor communicate dir (i.e. datadir) to all other processors.
!
if (exists) call mpibcast_char(dir,comm=MPI_COMM_WORLD)
!
endsubroutine get_datadir
!***********************************************************************
subroutine get_snapdir(dir)
!
! Read directory_snap from data/directory_snap, if that exists
! wd: I think we should unify these into a subroutine
! `overwrite_string_from_file(dir,file,label[optional])'
!
! 2-nov-02/axel: adapted from get_datadir
!
character (len=*) :: dir
character (len=fnlen) :: tdir
logical :: exists
integer :: unit=1
!
! check for existence of `data/directory_snap'
!
inquire(FILE=trim(datadir)//'/directory_snap',EXIST=exists)
if (exists) then
open(unit,FILE=trim(datadir)//'/directory_snap',FORM='formatted')
read(unit,'(A)') tdir
close(unit)
if (len(trim(tdir)) > 0) call parse_shell(tdir,dir)
endif
if (lroot.and.ip<10) print*,'get_snapdir: dir=',trim(dir)
!
endsubroutine get_snapdir
!***********************************************************************
subroutine read_init_pars(iostat)
!
use File_io, only: parallel_unit
!
integer, intent(out) :: iostat
!
read(parallel_unit, NML=init_pars, IOSTAT=iostat)
!
endsubroutine read_init_pars
!***********************************************************************
subroutine read_run_pars(iostat)
!
use File_io, only: parallel_unit
!
integer, intent(out) :: iostat
!
read(parallel_unit, NML=run_pars, IOSTAT=iostat)
!
endsubroutine read_run_pars
!***********************************************************************
subroutine read_all_init_pars(print)
!
! read input parameters (done by each processor)
!
! Now no warning is produced when somebody sets lcylindrical.
! "lcylindrical=T" is now outdated, use instead: "lcylinder_in_a_box=T"
! This renaming became necessary with the development of
! cylindrical coordinates which led to very similar names
! (coord_system="cylindrical_coords")
!
! 6-jul-02/axel: in case of error, print sample namelist
! 21-oct-03/tony: moved sample namelist stuff to a separate procedure
! 18-aug-15/PABourdin: reworked to simplify code and display all errors at once
! 19-aug-15/PABourdin: renamed from read_startpars to read_all_init_pars
!
use File_io, only: parallel_open, parallel_close, read_namelist
use Mpicomm, only: stop_it_if_any
use Particles_main, only: read_all_particles_init_pars
!
logical, optional, intent(IN) :: print
!
character(len=fnlen) :: file
!
! Set default to shearing sheet if lshear=.true. (even when Sshear==0.).
!
if (lshear .and. lforce_shear_bc) bcx(:)='she'
!
! Open namelist file and read through all items that *may* be present in the various modules.
!
if (lstart) then
file = 'start.in'
call parallel_open(file, remove_comments=.true.)
else
file = trim(datadir)//'/param.nml'
lparam_nml = .true.
call parallel_open(file)
endif
!
call read_namelist(read_init_pars ,'')
call read_namelist(read_initial_condition_pars ,'initial_condition_pars',linitial_condition)
call read_namelist(read_streamlines_init_pars ,'streamlines' ,lstreamlines)
call read_namelist(read_eos_init_pars ,'eos' ,leos)
call read_namelist(read_hydro_init_pars ,'hydro' ,lhydro)
call read_namelist(read_density_init_pars ,'density' ,ldensity)
call read_namelist(read_gravity_init_pars ,'grav' ,lgrav)
call read_namelist(read_selfgravity_init_pars ,'selfgrav' ,lselfgravity)
call read_namelist(read_poisson_init_pars ,'poisson' ,lpoisson)
call read_namelist(read_energy_init_pars ,'entropy' ,lenergy)
call read_namelist(read_magnetic_init_pars ,'magnetic' ,lmagnetic)
call read_namelist(read_lorenz_gauge_init_pars ,'lorenz_gauge' ,llorenz_gauge)
call read_namelist(read_testscalar_init_pars ,'testscalar' ,ltestscalar)
call read_namelist(read_testfield_init_pars ,'testfield' ,ltestfield)
call read_namelist(read_testflow_init_pars ,'testflow' ,ltestflow)
call read_namelist(read_radiation_init_pars ,'radiation' ,lradiation)
call read_namelist(read_pscalar_init_pars ,'pscalar' ,lpscalar)
call read_namelist(read_ascalar_init_pars ,'ascalar' ,lascalar)
call read_namelist(read_chiral_init_pars ,'chiral' ,lchiral)
call read_namelist(read_chemistry_init_pars ,'chemistry' ,lchemistry)
call read_namelist(read_signal_init_pars ,'signal' ,lsignal)
call read_namelist(read_dustvelocity_init_pars ,'dustvelocity' ,ldustvelocity)
call read_namelist(read_dustdensity_init_pars ,'dustdensity' ,ldustdensity)
call read_namelist(read_neutralvelocity_init_pars,'neutralvelocity',lneutralvelocity)
call read_namelist(read_neutraldensity_init_pars ,'neutraldensity' ,lneutraldensity)
call read_namelist(read_cosmicray_init_pars ,'cosmicray' ,lcosmicray)
call read_namelist(read_cosmicrayflux_init_pars ,'cosmicrayflux' ,lcosmicrayflux)
call read_namelist(read_interstellar_init_pars ,'interstellar' ,linterstellar)
call read_namelist(read_shear_init_pars ,'shear' ,lshear)
call read_namelist(read_special_init_pars ,'special' ,lspecial)
call read_namelist(read_solid_cells_init_pars ,'solid_cells' ,lsolid_cells)
call read_namelist(read_NSCBC_init_pars ,'NSCBC' ,lnscbc)
call read_namelist(read_polymer_init_pars ,'polymer' ,lpolymer)
call read_namelist(read_pointmasses_init_pars ,'pointmasses' ,lpointmasses)
!
call read_all_particles_init_pars
!
call parallel_close
!
if (lnamelist_error .and. .not. ltolerate_namelist_errors) then
call sample_pars
call stop_it_if_any (.true., 'read_all_init_pars: Please fix all above WARNINGs for file "'//trim(file)//'"')
endif
call stop_it_if_any (.false., '')
lparam_nml = .false.
if (lyinyang.and.coord_system/='spherical'.and.coord_system/='spherical_coords') then
if (lroot) call warning('read_all_init_pars', 'Yin-Yang grid only implemented for spherical coordinates')
lyinyang=.false.
endif
!
! Parse boundary conditions; compound conditions of the form `a:s' allow
! to have different variables at the lower and upper boundaries.
!
bcx=adjustl(bcx); bcy=adjustl(bcy); bcz=adjustl(bcz)
call parse_bc(bcx,bcx12)
call parse_bc(bcy,bcy12)
call parse_bc(bcz,bcz12)
!
fbcx (:,1) = fbcx1; fbcx (:,2) = fbcx2
fbcx_2(:,1) = fbcx1_2; fbcx_2(:,2) = fbcx2_2
fbcy (:,1) = fbcy1; fbcy (:,2) = fbcy2
fbcy_1(:,1) = fbcy1_1; fbcy_1(:,2) = fbcy2_1
fbcy_2(:,1) = fbcy1_2; fbcy_2(:,2) = fbcy2_2
fbcz (:,1) = fbcz1; fbcz (:,2) = fbcz2
fbcz_1(:,1) = fbcz1_1; fbcz_1(:,2) = fbcz2_1
fbcz_2(:,1) = fbcz1_2; fbcz_2(:,2) = fbcz2_2
!
if (lroot.and.ip<14) then
print*, 'bcx1,bcx2= ', bcx12(:,1)," : ",bcx12(:,2)
print*, 'bcy1,bcy2= ', bcy12(:,1)," : ",bcy12(:,2)
print*, 'bcz1,bcz2= ', bcz12(:,1)," : ",bcz12(:,2)
print*, 'lperi= ', lperi
endif
!
! Option to use maximal rather than total distance for courant time
!
if (old_cdtv) then
lmaximal_cdtv = old_cdtv
if (lroot) call warning('read_all_init_pars', &
'obsolete old_cdtv now replaced by more general lmaximal_cdtv')
endif
!
! Print SVN id from first line.
!
if (lroot) call svn_id(cvsid)
!
! Give online feedback if called with the PRINT optional argument.
!
if (loptest(print)) call write_all_init_pars
!
! Write parameters to log file.
!
if (lstart) call write_all_init_pars(FILE=trim(datadir)//'/params.log')
!
endsubroutine read_all_init_pars
!***********************************************************************
subroutine read_all_run_pars
!
! Read input parameters.
!
! 14-sep-01/axel: inserted from run.f90
! 31-may-02/wolf: renamed from 'cread' to 'read_runpars'
! 6-jul-02/axel: in case of error, print sample namelist
! 21-oct-03/tony: moved sample namelist stuff to a separate procedure
! 18-aug-15/PABourdin: reworked to simplify code and display all errors at once
! 19-aug-15/PABourdin: renamed from 'read_runpars' to 'read_all_run_pars'
! 22-mar-20/MR: removed unneeded call to write_all_run_pars and corresp. param.
! llogging
!
use Dustvelocity, only: copy_bcs_dust
use File_io, only: parallel_open, parallel_close, read_namelist
use General, only: loptest
use Mpicomm, only: stop_it_if_any
use Particles_main, only: read_all_particles_run_pars
use Sub, only: parse_bc
!
character(len=fnlen) :: file = 'run.in'
integer :: idum
!
!AB: putting it to impossible is not correct; it should be inherited from the previous run.
!AB: enabled again...
tstart=impossible
!
! Open namelist file.
!
call parallel_open(file, remove_comments=.true.)
!
! Read through all items that *may* be present in the various modules.
! AB: at some point the sgi_fix stuff should probably be removed (see sgi bug)
!
call read_namelist(read_run_pars ,'')
call read_namelist(read_streamlines_run_pars ,'streamlines' ,lstreamlines)
call read_namelist(read_eos_run_pars ,'eos' ,leos)
call read_namelist(read_hydro_run_pars ,'hydro' ,lhydro.or.lhydro_kinematic)
call read_namelist(read_density_run_pars ,'density' ,ldensity)
call read_namelist(read_forcing_run_pars ,'forcing' ,lforcing)
call read_namelist(read_gravity_run_pars ,'grav' ,lgrav)
call read_namelist(read_selfgravity_run_pars ,'selfgrav' ,lselfgravity)
call read_namelist(read_poisson_run_pars ,'poisson' ,lpoisson)
call read_namelist(read_energy_run_pars ,'entropy' ,lenergy)
! call read_namelist(read_conductivity_run_pars ,'conductivity')
call read_namelist(read_detonate_run_pars ,'detonate' ,ldetonate)
call read_namelist(read_magnetic_run_pars ,'magnetic' ,lmagnetic)
call read_namelist(read_lorenz_gauge_run_pars ,'lorenz_gauge' ,llorenz_gauge)
call read_namelist(read_testscalar_run_pars ,'testscalar' ,ltestscalar)
call read_namelist(read_testfield_run_pars ,'testfield' ,ltestfield)
call read_namelist(read_testflow_run_pars ,'testflow' ,ltestflow)
call read_namelist(read_radiation_run_pars ,'radiation' ,lradiation)
call read_namelist(read_pscalar_run_pars ,'pscalar' ,lpscalar)
call read_namelist(read_ascalar_run_pars ,'ascalar' ,lascalar)
call read_namelist(read_chiral_run_pars ,'chiral' ,lchiral)
call read_namelist(read_chemistry_run_pars ,'chemistry' ,lchemistry)
call read_namelist(read_dustvelocity_run_pars ,'dustvelocity' ,ldustvelocity)
call read_namelist(read_dustdensity_run_pars ,'dustdensity' ,ldustdensity)
call read_namelist(read_neutralvelocity_run_pars,'neutralvelocity' ,lneutralvelocity)
call read_namelist(read_neutraldensity_run_pars ,'neutraldensity' ,lneutraldensity)
call read_namelist(read_cosmicray_run_pars ,'cosmicray' ,lcosmicray)
call read_namelist(read_cosmicrayflux_run_pars ,'cosmicrayflux' ,lcosmicrayflux)
call read_namelist(read_heatflux_run_pars ,'heatflux' ,lheatflux)
call read_namelist(read_interstellar_run_pars ,'interstellar' ,linterstellar)
call read_namelist(read_shear_run_pars ,'shear' ,lshear)
call read_namelist(read_testperturb_run_pars ,'testperturb' ,ltestperturb)
call read_namelist(read_viscosity_run_pars ,'viscosity' ,lviscosity)
call read_namelist(read_special_run_pars ,'special' ,lspecial)
call read_namelist(read_shock_run_pars ,'shock' ,lshock)
call read_namelist(read_solid_cells_run_pars ,'solid_cells' ,lsolid_cells)
call read_namelist(read_NSCBC_run_pars ,'NSCBC' ,lnscbc)
call read_namelist(read_opacity_run_pars ,'opacity' ,lopacity)
call read_namelist(read_polymer_run_pars ,'polymer' ,lpolymer)
call read_namelist(read_pointmasses_run_pars ,'pointmasses' ,lpointmasses)
call read_namelist(read_power_spectrum_run_pars ,'power_spectrum' ,lpower_spectrum)
call read_namelist(read_implicit_diff_run_pars ,'implicit_diffusion',limplicit_diffusion)
!
call read_all_particles_run_pars
!
call parallel_close
!
if (lnamelist_error .and. .not. ltolerate_namelist_errors) then
call sample_pars
call stop_it_if_any (.true., 'read_all_run_pars: Please fix all above WARNINGs for file "'//trim(file)//'"')
endif
call stop_it_if_any (.false., '')
!
! Print SVN id from first line.
!
if (lroot) call svn_id(cvsid)
!
! Set debug logical (easier to use than the combination of ip and lroot).
!
ldebug = lroot .and. (ip < 7)
!
! Parse boundary conditions; compound conditions of the form `a:s' allow
! to have different variables at the lower and upper boundaries.
!
bcx=adjustl(bcx); bcy=adjustl(bcy); bcz=adjustl(bcz)
call parse_bc(bcx,bcx12)
call parse_bc(bcy,bcy12)
call parse_bc(bcz,bcz12)
!
fbcx (:,1) = fbcx1; fbcx (:,2) = fbcx2
fbcx_2(:,1) = fbcx1_2; fbcx_2(:,2) = fbcx2_2
fbcy (:,1) = fbcy1; fbcy (:,2) = fbcy2
fbcy_1(:,1) = fbcy1_1; fbcy_1(:,2) = fbcy2_1
fbcy_2(:,1) = fbcy1_2; fbcy_2(:,2) = fbcy2_2
fbcz (:,1) = fbcz1; fbcz (:,2) = fbcz2
fbcz_1(:,1) = fbcz1_1; fbcz_1(:,2) = fbcz2_1
fbcz_2(:,1) = fbcz1_2; fbcz_2(:,2) = fbcz2_2
!
if (lroot.and.ip<14) then
print*, 'bcx1,bcx2= ', bcx12(:,1)," : ",bcx12(:,2)
print*, 'bcy1,bcy2= ', bcy12(:,1)," : ",bcy12(:,2)
print*, 'bcz1,bcz2= ', bcz12(:,1)," : ",bcz12(:,2)
endif
!
! Ensure that right precision information is written in dim.dat.
!
lwrite_dim_again = lwrite_dim_again .or. lread_from_other_prec
!
! ivar_omit(1:2) defines a range of variables which are omitted when reading the snapshot.
!
if (ivar_omit(1)>0) then
if ((ivar_omit(2)>0).and.(ivar_omit(2)<ivar_omit(1))) then
idum=ivar_omit(1); ivar_omit(1)=ivar_omit(2); ivar_omit(2)=idum
endif
endif
if (lactive_dimension(3)) lzaver_on_input=.false.
!
it_rmv=max(it_rmv,0)
endsubroutine read_all_run_pars
!***********************************************************************
subroutine get_downpars(ind,n,ip)
!
! Calculates start indices & lengths for downsampled output
! Parameters: coordinate direction, number of inner grid points, processor number
!
! 13-feb-14/MR: coded
! 19-aug-15/PABourdin: moved, please do not use 'contains' in subroutines
! MR: Why not?
!
integer, intent(IN) :: ind, n, ip
!
if ( downsampl(ind)>n ) then
print*, 'get_downpars: Warning - stepsize for downsampling in '// &
coornames(ind)//'-direction ', downsampl(ind), 'greater than grid size ', n, &
'! Set to grid size.'
downsampl(ind)=n
endif
!
! first index in direction ind in local farray for output
!
firstind(ind) = downsampl(ind) - modulo(ip*n-1,downsampl(ind))
!
! number of output items in direction ind *without* ghost zones
!
ndown(ind) = get_range_no((/firstind(ind),n,downsampl(ind)/),1)
!
firstind(ind) = firstind(ind) + nghost
!
endsubroutine get_downpars
!***********************************************************************
subroutine write_stub(namelist, needed, omit_suffix)
!
! prints one sample namelist
!
! 1-sep-15/PABourin: coded
!
character(len=*) :: namelist
logical :: needed
logical, optional :: omit_suffix
!
character(len=labellen) :: type
!
if (.not. needed) return
!
if (lstart) then
type = 'init_pars'
else
type = 'run_pars'
endif
if (namelist /= '') type = '_'//trim (type)
if (loptest (omit_suffix)) type = ''
!
write (*,'(A)') '&'//trim (namelist)//trim (type)
write (*,'(A)') '/'
!
endsubroutine write_stub
!***********************************************************************
subroutine sample_pars
!
! standardises handling of errors in namelist, both for init and run
!
! 31-oct-13/MR: coded
! 1-sep-15/PABourin: reworked
!
if (lroot) then
!
write (*,'(A)') ''
if (lstart) then
write (*,'(A)') '=== BEGIN SAMPLE "start.in" ===>'
else
write (*,'(A)') '=== BEGIN SAMPLE "run.in" ===>'
endif
call write_stub ('', .true.)
!
if (lstart) then
call write_stub ('signal', lsignal)
call write_stub ('initial_condition_pars', linitial_condition, .true.)
endif
!
call write_stub ('streamlines', ltracers) !! questionable wg. ltracers
call write_stub ('eos', leos)
call write_stub ('hydro', lhydro .or. lhydro_kinematic)
call write_stub ('density', ldensity)
call write_stub ('grav', lgrav)
call write_stub ('selfgrav', lselfgravity)
call write_stub ('poisson', lpoisson)
call write_stub ('entropy', lentropy .or. ltemperature)
call write_stub ('energy', lthermal_energy)
! call write_stub ('conductivity', lconductivity)
call write_stub ('detonate', ldetonate)
call write_stub ('magnetic', lmagnetic)
call write_stub ('magn_mf', lmagn_mf)
call write_stub ('lorenz_gauge', llorenz_gauge)
call write_stub ('testscalar', ltestscalar)
call write_stub ('testfield', ltestfield)
call write_stub ('testflow', ltestflow)
call write_stub ('radiation', lradiation)
call write_stub ('pscalar', lpscalar)
call write_stub ('ascalar', lascalar)
call write_stub ('chiral', lchiral)
call write_stub ('chemistry', lchemistry)
call write_stub ('dustvelocity', ldustvelocity)
call write_stub ('dustdensity', ldustdensity)
call write_stub ('neutralvelocity', lneutralvelocity)
call write_stub ('neutraldensity', lneutraldensity)
call write_stub ('cosmicray', lcosmicray)
call write_stub ('cosmicrayflux', lcosmicrayflux)
call write_stub ('heatflux', lheatflux)
call write_stub ('interstellar', linterstellar)
call write_stub ('shear', lshear)
call write_stub ('testperturb', ltestperturb)
call write_stub ('special', lspecial)
call write_stub ('solid_cells', lsolid_cells)
call write_stub ('NSCBC', lnscbc)
call write_stub ('polymer', lpolymer)
call write_stub ('pointmasses', lpointmasses)
!
if (.not. lstart) then
call write_stub ('forcing', lforcing)
call write_stub ('shock', lshock)
call write_stub ('viscosity', lviscosity)
call write_stub ('power_spectrum', lpower_spectrum)
call write_stub ('implicit_diffusion', limplicit_diffusion)
endif
!
! Particles section.
!
if (lstart) then
call write_stub ('particles_dens', lparticles_density)
endif
!
call write_stub ('particles', lparticles)
call write_stub ('particles_radius', lparticles_radius)
call write_stub ('particles_cond', lparticles_condensation)
! call write_stub ('particles_potential', lparticles_potential)
call write_stub ('particles_spin', lparticles_spin)
call write_stub ('particles_sink', lparticles_sink)
call write_stub ('particles_number', lparticles_number)
call write_stub ('particles_selfgrav', lparticles_selfgravity)
call write_stub ('particles_stalker', lparticles_stalker)
call write_stub ('particles_mass', lparticles_mass)
call write_stub ('particles_drag', lparticles_drag)
call write_stub ('particles_TT', lparticles_temperature)
call write_stub ('particles_ads', lparticles_adsorbed)
call write_stub ('particles_surf', lparticles_surfspec)
call write_stub ('particles_chem', lparticles_chemistry)
!
if (.not. lstart) then
call write_stub ('particles_adapt', lparticles_adaptation)
call write_stub ('particles_coag', lparticles_coagulation)
call write_stub ('particles_cond', lparticles_condensation)
call write_stub ('particles_coll', lparticles_collisions)
call write_stub ('particles_stirring', lparticles_stirring)
call write_stub ('particles_diagnos_dv', lparticles_diagnos_dv)
call write_stub ('particles_diagnos_state', lparticles_diagnos_state)
endif
!
write (*,'(A)') '<=== END SAMPLE ==='
write (*,'(A)') ''
!
endif
!
endsubroutine sample_pars
!***********************************************************************
subroutine write_all_init_pars(file)
!
! Print input parameters.
!
! 4-oct-02/wolf: adapted
! 19-aug-15/PABourdin: renamed from 'print_startpars' to 'write_all_init_pars'
!
use Particles_main, only: write_all_particles_init_pars
!
character (len=*), optional, intent(in) :: file
!
character (len=datelen) :: date
integer :: unit
logical :: lidl_output
!
lidl_output = .false.
if (lroot) then
if (present(file)) then
unit = 1
if (file == 'IDL') then
open(unit,FILE=trim(datadir)//'/param.nml',DELIM='apostrophe',STATUS='unknown')
lidl_output = .true.
else
open(unit,FILE=file)
call date_time_string(date)
write(unit,*) '! -------------------------------------------------------------'
write(unit,'(A,A)') ' ! ', 'Initializing'
write(unit,'(A,A)') ' ! Date: ', trim(date)
write(unit,*) '! t=', t
endif
else
unit = 6 ! default unit is 6=stdout
endif
!
! If the boundary condition has been changed by subroutines then
! corresponding bc arrays need to be changed too. At present (june 2012)
! it is done only for the bcx array because parker_wind.f90 is the only
! subroutine that changes boundary conditions and it does it only for
! the density in the x direction.
!
if (lidl_output .and. lreset_boundary_values) call inverse_parse_bc(bcx,bcx12)
!
write(unit,NML=init_pars)
!
call write_initial_condition_pars(unit)
call write_streamlines_init_pars(unit)
call write_eos_init_pars(unit)
call write_hydro_init_pars(unit)
call write_density_init_pars(unit)
call write_gravity_init_pars(unit)
call write_selfgravity_init_pars(unit)
call write_poisson_init_pars(unit)
call write_energy_init_pars(unit)
call write_magnetic_init_pars(unit)
call write_lorenz_gauge_init_pars(unit)
call write_testscalar_init_pars(unit)
call write_testfield_init_pars(unit)
call write_testflow_init_pars(unit)
call write_radiation_init_pars(unit)
call write_pscalar_init_pars(unit)
call write_ascalar_init_pars(unit)
call write_chiral_init_pars(unit)
call write_chemistry_init_pars(unit)
call write_signal_init_pars(unit)
call write_dustvelocity_init_pars(unit)
call write_dustdensity_init_pars(unit)
call write_neutralvelocity_init_pars(unit)
call write_neutraldensity_init_pars(unit)
call write_cosmicray_init_pars(unit)
call write_cosmicrayflux_init_pars(unit)
call write_interstellar_init_pars(unit)
call write_shear_init_pars(unit)
call write_special_init_pars(unit)
call write_solid_cells_init_pars(unit)
call write_NSCBC_init_pars(unit)
call write_polymer_init_pars(unit)
call write_pointmasses_init_pars(unit)
!
call write_all_particles_init_pars(unit)
!
if (lidl_output) call write_IDL_logicals(unit)
!
if (present(file)) close(unit)
endif
!
endsubroutine write_all_init_pars
!***********************************************************************
subroutine write_all_run_pars(file)
!
! Print input parameters.
!
! 14-sep-01/axel: inserted from run.f90
! 31-may-02/wolf: renamed from 'cprint' to 'print_runpars'
! 4-oct-02/wolf: added log file stuff
! 19-aug-15/PABourdin: renamed from 'print_runpars' to 'write_all_run_pars'
! 11-jun-19/MR: full list of runpars only output into param2.nml, otherwise
! (in particular when outout in params.log) only the difference
! to the runpars of the preceding run.
! 22-mar-20/MR: added missing full output of run pars at first call of run.x
!
use Particles_main, only: write_all_particles_run_pars
use Syscalls, only: system_cmd
use File_io, only: file_exists,delete_file,flush_file
!
character (len=*), optional :: file
!
integer :: unit=1
character (len=linelen) :: line
character (len=datelen) :: date
logical :: lidl_output
!
if (lroot) then
lidl_output = .false.
if (lreloading) then
! get first line from file RELOAD
line = read_line_from_file('RELOAD')
if ((trim(line) == '') .or. (line == char(0))) line = 'Reloading'
else
line = 'Running'
endif
!
if (present(file)) then
if (file == 'stdout') then
unit = 6
elseif (file == 'IDL') then
!
! Save old param2.nml for later comparison.
!
if (file_exists('data/param2.nml')) &
call system_cmd("sed -e's/[eE]+*000*//g' -e's/ *, *$//' < data/param2.nml > data/param2.nml.sv")
open(unit,FILE=trim(datadir)//'/param2.nml',DELIM='apostrophe')
lidl_output = .true.
else
open(unit,FILE=file,position='append')
endif
else
open(unit,FILE=trim(datadir)//'/params.log',position='append')
endif
!
if (lidl_output) then
!
write(unit,NML=run_pars)
!
call write_streamlines_run_pars(unit)
call write_eos_run_pars(unit)
call write_hydro_run_pars(unit)
call write_density_run_pars(unit)
call write_forcing_run_pars(unit)
call write_gravity_run_pars(unit)
call write_selfgravity_run_pars(unit)
call write_poisson_run_pars(unit)
call write_energy_run_pars(unit)
call write_opacity_run_pars(unit)
! call write_conductivity_run_pars(unit)
call write_detonate_run_pars(unit)
call write_magnetic_run_pars(unit)
call write_lorenz_gauge_run_pars(unit)
call write_testscalar_run_pars(unit)
call write_testfield_run_pars(unit)
call write_testflow_run_pars(unit)
call write_radiation_run_pars(unit)
call write_pscalar_run_pars(unit)
call write_ascalar_run_pars(unit)
call write_chiral_run_pars(unit)
call write_chemistry_run_pars(unit)
call write_dustvelocity_run_pars(unit)
call write_dustdensity_run_pars(unit)
call write_neutralvelocity_run_pars(unit)
call write_neutraldensity_run_pars(unit)
call write_cosmicray_run_pars(unit)
call write_cosmicrayflux_run_pars(unit)
call write_heatflux_run_pars(unit)
call write_interstellar_run_pars(unit)
call write_shear_run_pars(unit)
call write_testperturb_run_pars(unit)
call write_viscosity_run_pars(unit)
call write_special_run_pars(unit)
call write_shock_run_pars(unit)
call write_solid_cells_run_pars(unit)
call write_NSCBC_run_pars(unit)
call write_power_spectrum_run_pars(unit)
call write_polymer_run_pars(unit)
call write_pointmasses_run_pars(unit)
call write_implicit_diff_run_pars(unit)
!
call write_all_particles_run_pars(unit)
!
write(unit,NML=IO_pars)
if (unit /= 6) close(unit)
else ! output in params.log, stdout or other file
! Add separator, comment, and time.
call date_time_string(date)
write(unit,*) '! -------------------------------------------------------------'
write(unit,'(A,A)') ' ! ', trim(line)
write(unit,'(A,A)') ' ! Date: ', trim(date)
write(unit,*) '! t=', t
call flush_file(unit)
if (unit /= 6) close(unit)
!
! Diff to old param2.nml.
!
if (file_exists('data/param2.nml.sv')) then
call system_cmd("sed -e's/[eE]+*000*//g' -e's/ *, *$//' < data/param2.nml |"// &
" diff -w data/param2.nml.sv - | grep '^[><]' >> data/params.log")
call delete_file('data/param2.nml.sv')
else
!
! Append param2.log to params.log for first time run.x is executed.
!
call system_cmd("cat data/param2.nml >> data/params.log")
endif
endif
!
endif
!
endsubroutine write_all_run_pars
!***********************************************************************
subroutine write_IDL_logicals(unit)
!
! Output logicals in namelist format for IDL.
!
integer, intent(in) :: unit
!
write(unit,'(A)') "&lphysics"
write(unit,'(A,L1,A)') " lhydro=", lhydro, ","
write(unit,'(A,L1,A)') " ldensity=", ldensity, ","
write(unit,'(A,L1,A)') " lentropy=", lentropy, ","
write(unit,'(A,L1,A)') " ltemperature=", ltemperature, ","
write(unit,'(A,L1,A)') " lgrav=", lgrav, ","
write(unit,'(A,L1,A)') " lshock=", lshock, ","
write(unit,'(A,L1,A)') " lmagnetic=", lmagnetic, ","
write(unit,'(A,L1,A)') " lforcing=", lforcing, ","
write(unit,'(A,L1,A)') " llorenz_gauge=", llorenz_gauge, ","
write(unit,'(A,L1,A)') " ldustvelocity=", ldustvelocity, ","
write(unit,'(A,L1,A)') " ldustdensity=", ldustdensity, ","
write(unit,'(A,L1,A)') " ltestscalar=", ltestscalar, ","
write(unit,'(A,L1,A)') " ltestfield=", ltestfield, ","
write(unit,'(A,L1,A)') " ltestflow=", ltestflow, ","
write(unit,'(A,L1,A)') " linterstellar=", linterstellar, ","
write(unit,'(A,L1,A)') " lcosmicray=", lcosmicray, ","
write(unit,'(A,L1,A)') " lcosmicrayflux=", lcosmicrayflux, ","
write(unit,'(A,L1,A)') " lheatflux=", lheatflux, ","
write(unit,'(A,L1,A)') " lshear=", lshear, ","
write(unit,'(A,L1,A)') " lpscalar=", lpscalar, ","
write(unit,'(A,L1,A)') " lascalar=", lascalar, ","
write(unit,'(A,L1,A)') " lradiation=", lradiation, ","
write(unit,'(A,L1,A)') " leos=", leos, ","
write(unit,'(A,L1,A)') " lchiral=", lchiral, ","
write(unit,'(A,L1,A)') " lneutralvelocity=", lneutralvelocity, ","
write(unit,'(A,L1,A)') " lneutraldensity=", lneutraldensity, ","
write(unit,'(A,L1,A)') " lpolymer=", lpolymer, ","
write(unit,'(A,L1,A)') " lpointmasses=", lpointmasses, ","
write(unit,'(A,L1,A)') " lsolid_cells=", lsolid_cells, ","
write(unit,'(A,L1,A)') " lpower_spectrum=", lpower_spectrum, ","
write(unit,'(A,L1,A)') " lparticles=", lparticles, ","
write(unit,'(A,L1,A)') " lparticles_drag=", lparticles_drag, ","
write(unit,'(A)') " /"
!
write(unit,NML=IO_pars)
!
endsubroutine write_IDL_logicals
!***********************************************************************
subroutine write_pencil_info
!
! Write information about requested and diagnostic pencils.
! Do this only when on root processor.
!
integer :: i
integer :: unit=1
!
if (lroot) then
open(unit,FILE=trim(datadir)//'/pencils.list')
write(unit,*) 'Pencils requested:'
do i=1,npencils
if (lpenc_requested(i)) write(unit,*) i, pencil_names(i)
enddo
write(unit,*) ''
write(unit,*) 'Pencils requested for diagnostics:'
do i=1,npencils
if (lpenc_diagnos(i)) write(unit,*) i, pencil_names(i)
enddo
write(unit,*) ''
write(unit,*) 'Pencils requested for 2-D diagnostics:'
do i=1,npencils
if (lpenc_diagnos2d(i)) write(unit,*) i, pencil_names(i)
enddo
write(unit,*) ''
write(unit,*) 'Pencils requested video output'
do i=1,npencils
if (lpenc_video(i)) write(unit,*) i, pencil_names(i)
enddo
if (ip<14) call information('write_pencil_info','pencil information written to the file pencils.list')
close(unit)
endif
!
endsubroutine write_pencil_info
!***********************************************************************
endmodule Param_IO