mat_data_struct f_matpar_H_water(double H) {
	mat_data_struct matnow;
	double hmelt = 333550;
	if (H <= 0) {
		matnow.H = H;
		matnow.rho = 920;
		matnow.cp = 2000;
		matnow.T = H / matnow.cp;
		matnow.k = 2.3;
		matnow.mu = 0;
	}
	else if (H >= hmelt) {
		matnow.H = H;
		matnow.cp = 4180;
		matnow.T = (H - hmelt) / matnow.cp;
		matnow.rho = 1000*(1+0.0007*matnow.T);
		matnow.k = 0.62;
		matnow.mu = 0.0101;
	}
	else {
		matnow.H = H;
		double x = H / hmelt;
		matnow.rho = 920 + x * 80;
		matnow.cp = 2000 + x * 218;
		matnow.T = 0;
		matnow.k = 2.3 - x * 1.68;
		matnow.mu = 0.;
	}
	return matnow;
}