&sml_param
sml_marker_temp_factor=1.5D0
sml_node_file='ma_1.00_rhoi_mesh.eqd.node'
sml_ele_file='ma_1.00_rhoi_mesh.eqd.ele'
sml_surf_file='ma_1.00_rhoi_mesh.eqd.flx.aif'
sml_monte_num=1000000
sml_f0_grid=.true.
sml_pol_decomp=.true.
sml_f0_grid_alpha=0.0000D0
sml_bt_sign=-1D0
sml_special=0             ! 0: normal simulation, 1: single particle simulation
sml_dt=4.0D-4            ! delta-t for one time step - unit of toroidal transit time. 
sml_mstep=4           ! total time step
sml_deltaf=.true.          ! delta-f simulation switch - incomplete
sml_turb_efield=.true.
sml_electron_on=.true.
sml_nphi_total=1
!sml_wedge_n=6
sml_bounce=0              ! Particle motion boundary condition 
sml_sheath_mode=0
sml_sheath_adjust=.true.
sml_sheath_init_pot_factor=2.5
sml_inpsi=0.00            ! inner boundary of simulation - unit of eq_x_psi
sml_outpsi=1.0           ! outter boundary of simulation - unit of eq_x_psi
sml_dwdt_exb_only=.true.
sml_restart_write_period=100
sml_restart=.false.
sml_deltaf_f0_mode=-2
sml_guess_table_size=1000
sml_00_efield=.false.
sml_00_npsi=30
sml_input_file_dir='.'
sml_bd_ext_delta2=0.01
sml_bd_ext_delta1=0.015
sml_bd_ext_delta1H=-0.01
sml_bd_ext_delta2H=0.02
sml_bd_ext_delta3=0.115
sml_bd_ext_delta4=0.05
sml_zero_inner_bd=.true.
sml_bounce_zero_weight=1
sml_rhomax=6.0D-2
sml_grid_nrho=6
sml_ncycle_half=30
sml_f_source_period=2
sml_source=.true.
sml_no_fp_in_f=.false.
sml_poisson_solver_type=0
sml_initial_deltaf_noise=1D-18
sml_f0_1_ln=2.2D0
sml_f0_1_lt=10D0
sml_f0_1_lt_e=10D0
sml_mode_select_mode=1
sml_mode_select_n=1
!sml_bp_mult=0.001
/

&ptl_param
ptl_mass_au=1D0             ! 1 for hydrogen, 2 for deutron
ptl_charge_eu=1D0            ! ion charge
ptl_e_mass_au=5.45D-4      ! previous value of 5.45D-4
ptl_num=2000              ! number of particle for simulation
ptl_maxnum=4000
/


&eq_param ! Initial equilibrium profile - Tanh profile
eq_filename='ma_151x151.eqd'
eq_den_shape=-1
eq_den_file='ma_den_bet02p500.dat'
eq_den_v1=1.0D19      ! inside value of density m^-3
eq_den_v2=1.0D19      ! outside value of density  m^-3
eq_den_x1=0.30D0      ! pedestal center
eq_den_x2=0.1D0       ! pedestal width
eq_den_v3=4.5D19
eq_den_x3=0D0
eq_tempi_shape=-1
eq_tempe_shape=-1
eq_tempi_file='ma_ti.dat'
eq_tempe_file='ma_te.dat'
eq_tempi_x1=0.50D0
eq_tempe_x1=0.50D0
eq_tempi_x2=0.10D0
eq_tempe_x2=0.10D0
eq_tempi_v1=2D3      ! ion temperature (inside) - eV
eq_tempi_v2=4d2     ! ion temperature (outside) - eV
eq_tempe_v1=2D3
eq_tempe_v2=4D2
eq_tempi_v3=5.7D3
eq_tempi_x3=0D0
eq_tempe_v3=5.7D3
eq_tempe_x3=0D0
/


&col_param  ! Collision
col_mode=4 ! 0 : off , 1 monte-carlo (non-conserving) 2: monte-carlo (conserving)
col_f_nthreads=1
/

&diag_param ! diagnosis
diag_1d_period=1
diag_f0_period=500
diag_tracer_period=10000
!diag_3d_on=.true.
/

&neu_param            ! neutral collision
!neu_col_mode=0
/

&src_param
src_narea=2
src_period=10
src_heat_power1=2.5D6
src_pout1=0.05
src_decay_width1=0.15
src_pout2=0.90D0
src_pin2=0.75D0
src_decay_width2=0.06
src_heat_power2=-2.5D6
/

&f0_param
f0_nmu=31
f0_nvp=15
f0_smu_max=3.0
f0_vp_max=3.0
/

&lim_param            ! limiter 
/

/
&smooth_param
smooth_mode_in=0
smooth_H_mode_in=2
smooth_H_n_in=1
smooth_r1_n_in=2
smooth_r1_d0_in=0.0042
smooth_r1_type_in=1
smooth_diag_mode_in=-1
/

&tbl_param
/

&prof_inparam
profile_papi_enable=.true.
/

&mon_param
/

&adios_param
/

&ptb_3db_param
/

&coupling_param
/

&cce_phicoupling
/

&cce_fcoupling
/

&diff_param
/
