In their comments to my first answer, the OP has clarified that they did not mean to regard the metric $h$ as a given, but, rather, an output of the problem of prescribing coframings by specifying their structure equations relative to certain given data. I don't want to erase my earlier answer, so I will now describe how to set up this different problem, which I would describe as follows:
One starts with a compact $3$-manifold $M$, a collection of smooth functions $f^i_{jk}$ on $M$ for $1\le i, j, k\le 3$, and a pair of integers $i_0$ and $j_0$ in the range $1\le i_0,j_0\le 3$. One seeks to find (or to determine whether there exists) a coframing $\omega = (\omega^i)$ on $M$ that satisfies
$\mathrm{d}\omega^i = \tfrac12 f^i_{jk}\,\omega^j\wedge\omega^k$, and
$f^{i_0}_{j_0k}\omega^k = \mathrm{d} a$ for some $a\in C^\infty(M)$.
Notes: First, the existence of a global coframing on $M$ implies that $M$ is orientable, so we might as well assume this from the outset. Also, it seems natural to assume that $f^i_{jk}=-f^i_{kj}$, but the OP does not specify this. For the purposes of Condition 1, this anti-symmetry assumption would be harmless, but Condition 2, as formulated, could have a different character without the anti-symmetry assumption.
Second, note that, if the three functions $f^{i_0}_{j_0k},\ 1\le k\le 3$ do not have any common zeros on $M$, then Condition 2 cannot ever be satisfied because the compactness of $M$ would ensure that the unknown function $a$ would have to have critical points, i.e., zeros of $\mathrm{d}a$. In practice, it seems better to weaken Condition 2 slightly by, instead, only requiring
2'. $\mathrm{d}\bigl(f^{i_0}_{j_0k}\,\omega^k \bigr)=0$,
since this is more directly expressed as a differential equation.
At any rate, Conditions 1 and 2/2', which are first order differential equations on $\omega$, clearly imply some zeroth order equations on $\omega$, that have to be taken into account:
$$
\begin{aligned}
0 = \mathrm{d}(\mathrm{d}\omega^i) = \mathrm{d}\bigl(\tfrac12 f^i_{jk}\,\omega^j\wedge\omega^k\bigr)\bigr)
&= \tfrac12\bigl(\mathrm{d} f^i_{jk}\wedge\omega^j\wedge\omega^k
+ Q^i(f)\,\omega^1\wedge\omega^2\wedge\omega^3\bigr)\\
0 =\mathrm{d}(\mathrm{d}a) =\mathrm{d}\bigl(f^{i_0}_{j_0k}\,\omega^k\bigr)
&= \mathrm{d}f^{i_0}_{j_0k}\wedge \omega^k + \tfrac12\,P_{jk}(f)\,\omega^j\wedge\omega^k
\end{aligned}
\tag 1
$$
where $Q^i(f)$ and $P_{jk}(f)$ are quadratic polynomials in the given functions $f=(f^i_{jk})$ that it will not be necessary to write out here. This is the vanishing of three $3$-forms and one $2$-form. When we regard $\omega=(\omega^i)$ as the usual tautological forms on the (trivial) principal $\mathrm{GL}(3,\mathbb{R})$-bundle $F\to M$ consisting of coframes of the tangent bundle of $M$, these equations will cut out a locus $B\subset F$, and solutions to the given problem will be coframings that are sections of $B\to M$. (Of course, not every section of $B$ will satisfy Conditions 1 and 2/2'.)
For generic choice of $f$, $B$ will be submanifold of $F$ of codimension 6, and the mapping $B\to M$ will be a submersion with fibers of dimension $3$, i.e., $B$ will be a $6$-manifold. Moreover, again, for generic choice of $f$, it can be shown that, at least locally in $B$, the three $1$-forms $\omega^i$ can be extended to a local coframing of $B$ by choosing three more $1$-forms $\pi_a$ ($1\le a\le 3$) for which there exist functions $A^{ia}_j$ on $B$ so that, pulled back to $B$ one has
$$
\Theta^i = \mathrm{d}\omega^i-\tfrac12\,f^i_{jk}\,\omega^j\wedge\omega^k
= A^{ia}_j\,\pi_a\wedge\omega^j,
$$
and, in addition, for each $u\in B$, the only $3$-plane $E\subset T_uB$
that satisfies $E^*(\Theta^i)=0$ and $E^*(\omega^1\wedge\omega^2\wedge\omega^3)\not=0$ is given by $\pi_1=\pi_2=\pi_3=0$. Consequently, the only sections of $B$ that satisfy Condition 1 are the ones for which the forms $\pi_a$ vanish when pulled back to the graph in $B$. In other words, every coframing $\omega$ that satisfies Conditions 1 and 2/2' must be a section of $B$ that that is an integral manifold of the ideal generated by the $\pi_a$.
However, we must also take into account Condition 2', and one finds that, on $B$ one must have an expansion of the form
$$
\Theta^4 = \mathrm{d}\bigl(f^{i_0}_{j_0k}\,\omega^k\bigr)
= C^a_k\,\pi_a\wedge\omega^k + T_{jk}\,\omega^j\wedge\omega^k
$$
for some functions $C^a_k$ and $T_{jk}=-T_{kj}$ on $B$. Since Condition 2' requires that $0 = E^*(\Theta^4)= E^*\bigl(T_{jk}\,\omega^j\wedge\omega^k\bigr)$ on any $3$-plane $E\subset T_uB$ that is tangent to the graph of a section $\omega$ satisfying Conditions 1 and 2/2', it follows that such sections must actually lie in the locus $B'\subset B$ defined by $T_{jk}=0$ ($1\le j < k\le 3$), which, generically is 3 independent equations on $B$, so that $B'\to M$ has discrete fibers over $M$.
Thus, in the generic situation, any coframing that satisfies Conditions 1 and 2/2' is found simply by differentiating the equations and imposing the zeroth order conditions that they imply.
Of course, it could happen that $B' = B = F$ in some cases, i.e., all of these zeroth order conditions are identities. This happens, for example, when the $f^i_{jk}$ are constants and $Q^i(f) = P_{jk}(f)= 0$ holds identically. In this case, of course, we know that all of the solutions are locally equivalent to the left-invariant coframing on a fixed Lie group $G$ of dimension $3$. (We also know that the compactness assumption implies further that there cannot be solutions on a compact manifold unless $f^{i_0}_{j_0k}=0$ for $1\le k\le 3$.)
There are many intermediate cases between these two extremes, depending on what identities and relations the functions $f = (f^i_{jk})$ satisfy on $M$. It would be very tedious to go through all the possibilities, so, unless one knows more about the functions $f$ on $M$, it is difficult to say much more in general about the problem and whether there is any need for any significant application of techniques from exterior differential systems.