Here is one specialized solution.
Let $h=g$. The the original equation becomes:
$$(3n-3)(f^2+2g^2)=(9n+1)(2fg+g^2)\tag{1}$$
We solve (1) for $f$ and obtain (assuming $g>0$):
$$f=A(n)g\tag{2}$$
$$A(n)=\frac{9n+1\pm\sqrt{10(9n^2+3n-2)}}{3(n-1)}\tag{3}$$
(a) For $n=2$, We have $A(2)=-1/3,13$. If $g=3 j$ where $j$ is arbitrary positive polynomial with integer coefficients, then $(f,g,h=g)$ is a solution to the original equation.
(b) For $n=81$, We have $A(81)=-1/6,25/4$. If $g=12 j$ where $j$ is arbitrary positive polynomial with integer coefficients, then $(f,g,h=g)$ is another solution to the original equation.
Other solutions can be found out in similar fashion.