Hi Omowumi,
The idea is to create the obs object with this flopy class and sample arguments:
obs = flopy.mf6.ModflowUtlobs(gwf, digits=6, continuous={("simulated_values.csv"): data)})
The variable "data" is a list of tuple where each tuple is:
(obsname, obstype, cellid)
for example, assuming disv grid, extracting result on layer 1, cell 1:
("well_1", "head", (1, 1))
Regards,
Mart