define custom_plot( wmin, wmax, xm1, xm2, cm_a, cm_b, fm_a, fm_b, mm_a, mm_b, T_exp)
{

variable phot2ergs = 12.39854*1.602e-9;

xrange(wmin, wmax);
variable lwmin = wmin - (wmax-wmin)/30.;

variable ymax = max(cm_a[where( (xm1 > wmin) and (xm2 <= wmax) ) ] ) *1.3;
yrange(0,ymax);


shplot(xm1, xm2, cm_a, 0.005, 2);
oshplot(xm1, xm2, cm_b, 0.005, 3);

oshplot(xm1, xm2, mm_a, 0.005, 5);
oshplot(xm1, xm2, mm_b, 0.005, 8);


color(1);
xylabel(lwmin, ymax*1.10,
 sprintf("[cm\\u-2\\d]   [phot/cm\\u2\\d/s]  [ergs/cm\\u2\\d/s]    [c/s]"));

color(2);xylabel(lwmin, ymax*1.06,
 sprintf("n\\dH\\u=0.0,      f=%10.3e, %10.3e,   %4.2f",
          sum(fm_a),	  sum(fm_a/xm1*phot2ergs),   sum(cm_a)/T_exp)     );

color(3);xylabel(lwmin, ymax*1.02,
 sprintf("n\\dH\\u=1.0e+22, f=%10.3e, %10.3e,   %4.2f",
          sum(fm_b),	  sum(fm_b/xm1*phot2ergs),   sum(cm_b)/T_exp)     );


% plot_group(brightest(20, where(wl(wmin, wmax) )), 1);

p_groups(wmin, wmax);


}