/*hmiopc derive total H- H minus opacity */
double hmiopc(double freq);