subroutine fnc_cor(e_raw,t_raw,e_tru,t_tru,ierr)
*
*  apply gain correction, t0 from end of spring '99 run
*
	implicit none
*
	integer ierr
	real e_raw(96),t_raw(96),e_tru(96),t_tru(96)
*
	integer i
*
	real gain(96),t0(96)
	data gain/
     +  1.3490,
     +  1.7760,
     +  1.3431,
     +  1.6535,
     +  1.3822,
     +  0.8921,
     +  1.1266,
     +  1.1075,
     +  0.6885,
     +  1.3780,
     +  1.4963,
     +  1.6854,
     +  1.3282,
     +  1.1951,
     +  1.6311,
     +  1.2355,
     +  0.9265,
     +  1.1017,
     +  1.0387,
     +  1.1235,
     +  3.0631,
     +  1.4574,
     +  2.9403,
     +  2.3105,
     +  0.7018,
     +  0.5689,
     +  1.3023,
     +  0.4308,
     +  7.2744,
     +  4.8130,
     +  6.3464,
     +  4.3965,
     +  6.4672,
     +  2.1223,
     +  5.7190,
     +  3.8270,
     +  5.5933,
     +  3.5573,
     +  7.2953,
     +  4.5155,
     + -1.0000,
     + -1.0000,
     +  0.5685,
     +  4.0755,
     +  0.4863,
     +  1.1021,
     +  1.9790,
     +  1.3107,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000,
     + -1.0000/
	data t0/
     + 11.7093,
     +  5.3815,
     + 11.0336,
     +  6.4755,
     + 11.5943,
     +  7.4655,
     +  8.6438,
     +  7.8608,
     +  5.3950,
     +  8.0930,
     +  2.8922,
     +  4.1906,
     + -1.7038,
     +  1.4062,
     +  0.9303,
     + -0.0107,
     +  3.0629,
     +  1.0491,
     +  2.9087,
     +  3.1353,
     + 11.8983,
     + 10.5449,
     +  8.5466,
     + 10.5997,
     +  9.9274,
     +  9.4582,
     + 14.2622,
     +  8.7502,
     +  8.7058,
     +  7.4473,
     +  6.7529,
     + 10.7423,
     +  5.8481,
     +  4.7975,
     +  7.5267,
     +  3.0696,
     +  5.3129,
     +  8.9596,
     + 16.1758,
     +  7.1330,
     + 10.9356,
     +  9.7191,
     + -8.5086,
     + -0.0884,
     + -3.8019,
     + -1.1336,
     +  4.7208,
     +  2.5485,
     +  0.0000,
     +  0.0000,
     +  8.4536,
     +  7.3361,
     +  8.9089,
     +  9.6676,
     +  9.4657,
     +  6.5749,
     +  7.3859,
     +  7.3731,
     +  6.9575,
     +  6.7292,
     +  9.0566,
     +  6.5508,
     +  6.7255,
     +  7.9918,
     +  8.6127,
     +  7.6008,
     +  5.9447,
     +  9.5933,
     +  7.7858,
     +  7.7879,
     + 15.4427,
     +  7.3080,
     +  9.7525,
     +  8.8299,
     +  8.4940,
     +  8.0453,
     +  9.1049,
     +  6.3088,
     + 12.3727,
     +  5.8498,
     +  9.1602,
     +  8.2551,
     +  4.5354,
     +  2.9856,
     +  1.7248,
     +  1.7374,
     +  2.0582,
     + -0.4280,
     +  2.9363,
     +  4.9999,
     +  8.9026,
     + 11.1938,
     +-10.8447,
     +-10.3728,
     + -9.2904,
     + -9.6806/
*
        real bg_2_ep_scale
        parameter (bg_2_ep_scale = 1.05)
*
*
	ierr = 0
*
	do i=1,96
	   if(gain(i).gt.0.) then
	      e_tru(i) = e_raw(i)/gain(i)
	      t_tru(i) = t_raw(i)-t0(i)
	   else
	      e_tru(i) = 0.
	      t_tru(i) = 0.
	   endif
	end do
*
*  put in the b.g.->e-p scale
*
	do i=1,42
	   e_tru(i) = e_tru(i) / bg_2_ep_scale
	end do
*
	return
	end