function dy = glaze_97(z,y)

global phi_av phi_ad r_ab c_ab rho_ad rho_ab rho_av wa m_s c_d c_l c_v c_s al r_v ep  omega g  c_b  l rho_l rho_s r_d 

dy=zeros(5,1);

m = y(1)+y(2)+y(3)+m_s; % gesamtmasse in column
c_b = (y(1)*c_d+y(2)*c_v+y(3)*c_l+m_s*c_s)/m; % bulk specific heat

[thetaa,p] = usstand(z);


a1=r_v*y(5)/p;
b1=y(2)+ep*y(1);
c1 = b1/(m^2*c_b);
d1 = m/y(4);
e1=g/y(4);
f1=a1*b1/c_b;
g1=b1/m^2*c_b;
h1=a1*b1/m*c_b;
i1=y(3)/rho_l;
j1=m_s/rho_s;

dy(1)=2*al*rho_ad*phi_ad*(a1*y(4)*c1)^0.5;
dy(2)=2*al*rho_av*phi_av*(a1*y(4)*c1)^0.5-omega*y(2)*d1;
dy(3)=omega*y(2)*d1;
dy(4)=e1*(rho_ab*(f1)-m^2);
dy(5)=2*al*rho_ab*c_ab*thetaa*(a1*y(4)*g1)^0.5+l*omega*y(2)*d1...
    -rho_ab*g*((h1)-(i1+j1));


end