Initial e2d value can now be define in the input file

src/read_input_file.f90

@@ -68,7 +68,7 @@ subroutine read_input_file (params,threadinfo,material0,mat,surface,boxes,     &
 !    default value is 1.d0
 ! surface are surface structure
 !    for each surface, one needs to define a levelt, itype, material and fnme.
-!    levelt is the inital level for the particles on the surface; to be accurate
+!    levelt is the initial level for the particles on the surface; to be accurate
 !    and avoid wholes in the surface during definition of the lsf, one should use
 !    levelt=levelmax_oct+1 for all surfaces as a minimum value; itype should be 1
 !    for foldable surfaces or 0 for nonfoldable surfaces; material is the material
@@ -175,6 +175,11 @@ endif
 call mpi_bcast(params%init_e2d,1,mpi_logical,0,mpi_comm_world,ierr)
 write(threadinfo%Logunit,*) 'init_e2d ',params%init_e2d
 
+params%e2d0=1.d0
+if (iproc==0) call scanfile (params%infile,'e2d0',params%e2d0,ires)
+call mpi_bcast(params%e2d0,1,mpi_double_precision,0,mpi_comm_world,ierr)
+write(threadinfo%Logunit,*) 'e2d0 ',params%e2d0
+
 do i=0,params%nmat
    write(cm,'(i3)') i
    il=1