PRO SUB_READ_AUX_AK19, RAWDATA, SHOT, TIME, LAT, LON, ALT, CO_GAIN, X_GAIN, TILT, temp, inst_temp, V0, V1, ice
	;reads auxilary data at end of shot
	;  time is decimal hours
	; 
	; RAWDATA is data file
	; SHOT is shot number
	; CO_GAIN is gain of co pol channel; uses log-log fit
	; X_GAIN is gain of x pol channel
	; tilt is pitch angle of the lidar
	; temp is not used for this deployment
	
	; GPS time
	hr = double(string(rawdata[shot,4104])+string(rawdata[shot,4105]))
	mn = double(string(rawdata[shot,4106])+string(rawdata[shot,4107]))
	fs = double('0'+string(rawdata[shot,4108])+string(rawdata[shot,4109]) $
	  + string(rawdata[shot,4110])+string(rawdata[shot,4111])$
	  + string(rawdata[shot,4112])+string(rawdata[shot,4113]))
	time = hr + mn/60 + fs/3600
		
	;computer time time
	hr = double(string(rawdata[shot,4122])+string(rawdata[shot,4123]))
  mn = double(string(rawdata[shot,4124])+string(rawdata[shot,4125]))
  fs = double('0'+string(rawdata[shot,4126])+string(rawdata[shot,4127]) $
          + string(rawdata[shot,4128])+string(rawdata[shot,4129])$
          + string(rawdata[shot,4130])+string(rawdata[shot,4131]))
  ;time = hr + mn/60 + fs/3600
          
	; lat, lon, and altitude
	lat = string(rawdata[shot,4001])+string(rawdata[shot,4002])+string(rawdata[shot,4003])+string(rawdata[shot,4004]) $
	    +string(rawdata[shot,4005])+string(rawdata[shot,4006])+string(rawdata[shot,4007])+string(rawdata[shot,4008]) $
	    +string(rawdata[shot,4009])
	lat = double(lat)
  lon = string(rawdata[shot,4012])+string(rawdata[shot,4013])+string(rawdata[shot,4014])+string(rawdata[shot,4015]) $
      +string(rawdata[shot,4016])+string(rawdata[shot,4017])+string(rawdata[shot,4018])+string(rawdata[shot,4019]) $
      +string(rawdata[shot,4020])
  lon = double(lon)
  alt = string(rawdata[shot,4022])+string(rawdata[shot,4023])+string(rawdata[shot,4024])+string(rawdata[shot,4025]) $
      +string(rawdata[shot,4026])+string(rawdata[shot,4027])+string(rawdata[shot,4027])+string(rawdata[shot,4029]) $
      +string(rawdata[shot,4030])+string(rawdata[shot,4031])
  alt = double(alt)
  
  ; PMT gains
  AD0V = string(rawdata[shot,4049])+string(rawdata[shot,4050])+string(rawdata[shot,4051])+string(rawdata[shot,4052]) $
        +string(rawdata[shot,4053])+string(rawdata[shot,4054])+string(rawdata[shot,4055]) 
  V0 = double(AD0V)
  logV0 = alog10(V0)
  CO_GAIN = 10^(-3.2826 + 24.239*logV0 - 14.773*logV0*logV0)
  
  AD1V = string(rawdata[shot,4057])+string(rawdata[shot,4058])+string(rawdata[shot,4059])+string(rawdata[shot,4060]) $
        +string(rawdata[shot,4061])+string(rawdata[shot,4062])+string(rawdata[shot,4063])
  V1 = double(AD1V)
  logV1 = alog10(V1)
  X_GAIN = 10^(-4.2331 + 28.36*logV1 - 18.388*logV1*logV1)
  
    ; tilt
  AD3V = double(string(rawdata[shot,4072])+string(rawdata[shot,4073])+string(rawdata[shot,4074])+string(rawdata[shot,4075]) $
        +string(rawdata[shot,4076])+string(rawdata[shot,4077])+string(rawdata[shot,4078])+string(rawdata[shot,4079]))
  tilt = 100*AD3V
  
  ; temp
  AD4V = double(string(rawdata[shot,4080])+string(rawdata[shot,4081])+string(rawdata[shot,4082])+string(rawdata[shot,4083]) $
        +string(rawdata[shot,4084])+string(rawdata[shot,4085])+string(rawdata[shot,4086])+string(rawdata[shot,4087]))
  temp = 5.096*AD4V - 0.3129

  ; instrument temperature
  AD5V = double(string(rawdata[shot,4088])+string(rawdata[shot,4089])+string(rawdata[shot,4090])+string(rawdata[shot,4091]) $
    +string(rawdata[shot,4092])+string(rawdata[shot,4093])+string(rawdata[shot,4094])+string(rawdata[shot,4095]))
  inst_temp = 100*AD5V
  
  ice = rawdata[shot,4199]
        
RETURN
END