FUNCTION FRDSYS
*-- Author :    Stephen J. Maxfield   19/08/92
      FUNCTION FRDSYS(RAD, LINT)
                                                                        
*     Calculate systematic correction to charge division coordinate                                           
                                                                        
      DIMENSION ALPHA(7), BETA(7), GAMMA(7)                             
      DATA ALPHA/-17.0, -14.0, -12.3, -10.5,  -9.4,  -8.6,  -7.8/       
      DATA BETA/ 0.361, 0.307, 0.298, 0.275, 0.264, 0.270, 0.281/       
      DATA GAMMA/-0.00123, -0.00092, -0.00105, -0.00096,                
     +           -0.00093, -0.00109, -0.00133/                          
                                                                        
      IF(LINT .GT.12) LINT = 12                                         
      IF(LINT .LT. 6) LINT = 6                                          
                                                                        
      J = LINT - 5                                                      
                                                                        
      ERRSYS = ALPHA(J) + BETA(J) * RAD + GAMMA(J) * RAD * RAD          
                                                                        
      FRDSYS = RAD + ERRSYS                                             
                                                                        
      RETURN                                                            
      END                                                               
*