SUBROUTINE FKZTOR
*-- Author : S.Burke / J.V. Morris
      SUBROUTINE FKZTOR(SIN,CIN,BZ,COUT)
**********************************************************************                                        
*                                                                    *                                        
* Transform a covariance matrix from fixed z to fixed R              *                                        
*                                                                    *                                        
* Assumed to be in terms of (R,PHI)                                  *                                        
*                                                                    *                                        
*                        *** NOT TESTED ***                          *                                        
*                                                                    *                                        
**********************************************************************                                        
                                                                        
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)                               
                                                                        
      DIMENSION SIN(5),CIN(5,5),COUT(5,5)                               
                                                                        
**********************************************************************                                        
                                                                        
      CPMP = COS(SIN(5) - SIN(2))                                       
      SPMP = (1.D0 - CPMP*CPMP)                                         
                                                                        
      DR = SIN(4)*CPMP                                                  
      DPHI = SIN(4)*SPMP/SIN(1)                                         
      DP = SIN(3)*BZ*SQRT(1.D0 + SIN(4)*SIN(4))                         
                                                                        
      DZ2 = CIN(1,1)/(DR*DR)                                            
      DPHI2 = DPHI*DPHI                                                 
      DP2 = DP*DP                                                       
                                                                        
      CALL UCOPY(CIN(2,2),COUT(2,2),38)                                                                
                                                                        
      COUT(1,1) = DZ2                                                   
      COUT(2,1) = DPHI*DZ2                                              
      COUT(3,1) = 0.D0                                                  
      COUT(4,1) = 0.D0                                                  
      COUT(5,1) = DP*DZ2                                                
                                                                        
      COUT(2,2) = COUT(2,2) + DPHI2*DZ2                                 
      COUT(5,2) = COUT(5,2) + DPHI*DP*DZ2                               
                                                                        
      COUT(5,5) = COUT(5,5) + DP2*DZ2                                   
                                                                        
      RETURN                                                            
      END                                                               
*