update_z_q0 function

Gibbs sampling for zz for q=0q=0