# [[file:~/Dropbox/Research/TetradTurbulence/MLCoarseGrainedVGT/docs/Article/article.org::*Code up correct formulation of stochastic forcing][Code up correct formulation of stochastic forcing:1]]
function createNoiseRateForcingHIT(D_a, D_s)
   c = zeros(3,3,3,3);
   delta(i,j) = (i==j) ? 1 : 0;
   for i in 1:3
      for j in 1:3
         for k in 1:3
            for l in 1:3
   	    term1 = -(1/3)*sqrt(D_s/5)*delta(i,j)*delta(k,l);
   	    term2 = (1/2)*(sqrt(D_s/5) + sqrt(D_a/3))*delta(i,k)*delta(j,l);
   	    term3 = (1/2)*(sqrt(D_s/5) - sqrt(D_a/3))*delta(i,l)*delta(j,k);
   	    c[i,j,k,l] = term1 + term2 + term3;
   	    end
         end
      end
   end
   return c;
end
# Code up correct formulation of stochastic forcing:1 ends here
