Thanks a lot for your prompt reply. Actually, I also stored other things in data, such as
data = Dict()
data["total_time"] = dtime#Tf
data["time"] = 0.
data["dt"] = dtime
data["friction"] = fric
data["dirichlet_grid"] = [("back", (0,0,1)),
("front",(1,1,1)),
("bottom",(1,1,1)),
("left",(1,1,1)),
("right",(1,1,1)),
] # => fix bottom nodes on X/Y/Z dir
That was why I went for data=Dict(). Any suggestion to have both flexibility and efficiency?