/*forlin radiative acceleration due to line absorption of incident continuum */
double forlin(void);