The result of running out is ppmv, how do you convert ppmv to molecule/m^2? I see a formula, t=273/t;
p=P/101325;
Vm=22.4/(t.*p);
for n = 1:(nlayer)
ready_nh3(:,:,n,:)=60.2.*heigth(:,:,n,:).*nh3(:,:,n,:)./Vm;
end, may I ask where the 60.2 in this one comes from