Visualize single-grid hyperelastic results

SetDirectory[NotebookDirectory[]];
gFile="./out/global.hdf";
globalDS=Import[gFile];
GlobalInfo=Import[gFile,{"Datasets","GlobalInfo"}];
nFrames=IntegerPart[GlobalInfo[[1]]]+1; nq = IntegerPart[GlobalInfo[[2]]];
qmin=Import[gFile,{"Datasets","qmin"}][[1;;nq]];
qmax=Import[gFile,{"Datasets","qmax"}][[1;;nq]];
qFiles=FileNames["./out/q*.hdf"];
qNames={"(Subscript[v, 1])","(Subscript[v, 2])","(Subscript[F, 11])","(Subscript[F, 12])","(Subscript[F, 21])","(Subscript[F, 22])"};

qPlot[n_]:=Module[{q,qDS},
qDS=Import[qFiles[[n]]];
q=Table[Import[qFiles[[n]],{"Datasets",qDS[[k]]}],{k,1,nq}];
p=GraphicsGrid[Transpose[ ArrayReshape[
Table[ MatrixPlot[q[[k,;;,;;]],DataReversed->True,PlotLegends->Automatic,PlotLabel->qNames[[k]], PlotRange->{qmin[[k]],qmax[[k]]}],{k,nq}],{3,2}]],ImageSize->Full];
Return[p]
]

Manipulate[qPlot[n],{n,1,nFrames,1}]

Export["anim.gif",Manipulate[qPlot[n],{n,1,nFrames,1}]]