Added bcdef to read_input_file and define_bc calls

src/DOUAR.f90

@@ -195,7 +195,7 @@ end if
 !------------------------------------------------------------------------------|
 !------------------------------------------------------------------------------|
 call read_input_file (params,threadinfo,material0,mat,surface,boxes,sections,  &
-                     cube_faces,nest)
+                     cube_faces,nest,bcdef)
 
 if (params%irestart==0) then
    current_level=params%initial_refine_level