Solved by tom in post #2 
              
              
                So you have prior = p = N(0,1) and q = N(mu, diag(var)) 
It seems that you compute D_KL(p || q) in the first and the more common D_KL(q || p) in the second. KL Divergence is not symmetric, so these will differ. 
Best regards 
Thomas