Gaussian curvature has several applications in engineering and manufacturing. Specifically, a surface whose Gaussian curvature is zero is developable, which means that it can be formed from a flat sheet without any stretching or tearing. Planes, cylinders and cones are the most obvious developable surfaces, but there are others, too. All developable surfaces are ruled surfaces, but the converse is not true.
Why does this matter? Here are a few examples:
Suppose I'm making the hull of a ship from flat metal plates. It's not too difficult to bend these plates, but it's very difficult to stretch them, so developable surfaces are better for ship hulls.
Suppose I'm designing a bottle that will have a paper label glued to its outside. The label can be bent, but not stretched, so the area of the bottle where it is afixed needs to be developable.
Suppose I'm going to manufacture a ruled surface by "swarf cutting". This means that the shape is formed by cutting with the side of a cylindrical cutting tool. At each step, the tool will contact the surface along a straight line "ruling", and the surface normal will be the same at all points along this line. It's not obvious, but it turns out that this implies that the surface is developable. So, if you want to manufacture a ruled surface by swarf cutting (which is desirable, because it's cheap and fast), then you should make it developable. There's a nice picture of swarf cutting on this page, in the section entitled "Advanced 5D Milling".
More generally, if Gaussian curvature is not zero, its magnitude gives you some information about "how far from developable" the surface is. This in turn tells you how much stretching will be needed to form it. If the stretching is too much, the material will tear. This is often an issue when stamping sheet metal to form car bodies, for example. Sometimes the body panel has to be redesigned to avoid the tearing.
Regarding what constitutes an "application" ... when mathematicians say "application", they often mean using some theorem to prove other results, in other areas of mathematics. This always seems a bit incestuous, to me; I much prefer "applications" in the real world.
Best Answer
I'm not sure if there's a simpler way to do it, but here's a way that only involves the basic theory of surfaces.
Your first fundamental form is $E = G = \frac{1}{g^2}, F=0.$ Hence you can compute the Christoffel symbols by solving \begin{equation} \begin{pmatrix} 1/g^2 & 0\\ 0 & 1/g^2 \end{pmatrix} \begin{pmatrix} \Gamma_{11}^1\\ \Gamma_{11}^2 \end{pmatrix} = \begin{pmatrix} \frac{1}{2}E_x\\ F_x-\frac{1}{2}E_y \end{pmatrix} = \begin{pmatrix} -g_x/g^3\\ g_y/g^3 \end{pmatrix} \quad \Rightarrow \quad \begin{pmatrix} \Gamma_{11}^1\\ \Gamma_{11}^2 \end{pmatrix} = \begin{pmatrix} -g_x/g\\ g_y/g \end{pmatrix} \end{equation} \begin{equation} \begin{pmatrix} 1/g^2 & 0\\ 0 & 1/g^2 \end{pmatrix} \begin{pmatrix} \Gamma_{12}^1\\ \Gamma_{12}^2 \end{pmatrix} = \begin{pmatrix} \frac{1}{2}E_y\\ \frac{1}{2}G_x \end{pmatrix} = \begin{pmatrix} -g_y/g^3\\ -g_x/g^3 \end{pmatrix} \quad \Rightarrow \quad \begin{pmatrix} \Gamma_{12}^1\\ \Gamma_{12}^2 \end{pmatrix} = \begin{pmatrix} -g_y/g\\ -g_x/g \end{pmatrix} \end{equation} \begin{equation} \begin{pmatrix} 1/g^2 & 0\\ 0 & 1/g^2 \end{pmatrix} \begin{pmatrix} \Gamma_{22}^1\\ \Gamma_{22}^2 \end{pmatrix} = \begin{pmatrix} F_y-\frac{1}{2}G_x\\ \frac{1}{2}G_y \end{pmatrix} = \begin{pmatrix} g_x/g^3\\ -g_y/g^3 \end{pmatrix} \quad \Rightarrow \quad \begin{pmatrix} \Gamma_{22}^1\\ \Gamma_{22}^2 \end{pmatrix} = \begin{pmatrix} g_x/g\\ -g_y/g \end{pmatrix} \end{equation}
Furthermore, $(\Gamma_{12}^2)_x = -g_{xx}/g+g_x^2/g^2,$ and $(\Gamma_{11}^2)_y = g_{yy}/g-g_y^2/g^2$, so that the Gauss formula (see do Carmo's Differential Geometry of Curves and Surfaces, chapter 4-3, eq. 5) gives \begin{align} K &= -\frac{1}{E}\big( (\Gamma_{12}^2)_x - (\Gamma_{11}^2)_y + \Gamma_{12}^1\Gamma_{11}^2-\Gamma_{11}^2\Gamma_{22}^2-\Gamma_{11}^1\Gamma_{12}^2\big)\\ &= -g^2 \big( -g_{xx}/g+g_x^2/g^2 - g_{yy}/g + g_y^2/g^2 +g_y^2/g^2 +g_x^2/g^2 -g_y^2/g^2 -g_x^2/g^2\big)\\ &= - \big( -g_{xx}g-g_{yy}g + g_x^2 + g_y^2 \big)\\ &= g(g_{xx}+g_{yy}) - (g_x^2+g_y^2) . \end{align}