#pragma once
double Zeta_Solar(double T, double f, int ABcode);
double Zeta_arbitrary(double T, double f, int ABcode, int N_fZ, int N_TZ, double *lnfZ_arr, double *lnTZ_arr, double *Z_arr);