c Be sure to include ilut_pormed.h before this file. double precision coef, myxcur(MAXN) common /scratch/ coef, myxcur