Divided write(0,*) grain_a into grain-axis-x,y,z because printing the array in the header bleeds into the next line, which is not commented out.