fixed bug

fixed bug: sub3d_elem_module, after changing plyangle into optional,
should've put if(present(plyangle)) statement bf assigning plyangle to
elem%plyangle

unfixed: sth wrong in coh3d8 and interface type routines

notes:
- print*, .... output at log file
- can combine preproc_libnee and libmat