&input flfrc = 'Si.fc' fldos = 'Si.phdos' asr = 'crystal' dos = .true. nk1 = 16, nk2 = 16, nk3 = 16, /