*-- Author : S.Burke / J.V. Morris
      SUBROUTINE FKQADD(C,Q)
**********************************************************************                                        
*                                                                    *                                        
* Add a multiple scattering matrix to a covariance matrix            *                                        
*                                                                    *                                        
* Output is C = C + Q, lower half only, where Q(3,x) = Q(x,3) = 0    *                                        
*                                                                    *                                        
**********************************************************************                                        
                                                                        
      DOUBLE PRECISION C(5,5),Q(5,5)                                    
                                                                        
      C(1,1) = C(1,1) + Q(1,1)                                          
      C(2,1) = C(2,1) + Q(2,1)                                          
      C(4,1) = C(4,1) + Q(4,1)                                          
      C(5,1) = C(5,1) + Q(5,1)                                          
      C(2,2) = C(2,2) + Q(2,2)                                          
      C(4,2) = C(4,2) + Q(4,2)                                          
      C(5,2) = C(5,2) + Q(5,2)                                          
      C(4,4) = C(4,4) + Q(4,4)                                          
      C(5,4) = C(5,4) + Q(5,4)                                          
      C(5,5) = C(5,5) + Q(5,5)                                          
                                                                        
      RETURN                                                            
      END                                                               
*