#ifndef ICE_H
#define ICE_H
double psub_Ice(double T);
double g_Ice(double T,double p);
double dg_dp_Ice(double T, double p);
double dg2_dp2_Ice(double T, double p);
double IsothermCompress_Ice(double T, double p);
double dg_dT_Ice(double T, double p);
double h_Ice(double T, double p);
double s_Ice(double T, double p);
double rho_Ice(double T, double p);
#endif