%% calculate the first starting field at the bottom

function z0=cal_z0(W,P,kn,Ch,H)
kh=W/Ch;
x=1-(kn/kh)^2;
z0=P*sqrt(sqrt(x))/sqrt(H*kh);
