gfortran fixes in initialize_p.F for src_MD-M