Agustin,
The solution to this is to change the read statement in subroutine get_xy : read (line,fmt='(2f30.0)',iostat=iostat) x1,y1 and replace this by a routine to finds each text string that is defined by a separator character ( comma, space, tab or other ) then read each text string to get each number. (Note: to read tab characters, see OPEN help) The coding for this is straight forward, although the code gets a bit longer as more special characteristics of the data file are included.
John