read( funit ) filek
read(cunit) filek
if (filek .NE. 1) then
& filek, 1
if (filek .NE. k) then
& filek, k
& filek,
integer filek