FUNCTION CP(T,KC)                                                       
C                                                                              
         COMMON /b_1/TEMPD(5,5),A(5,5),AX(5,5),B(5,5),BX(5,5),                      
     1          C(5,5),CX(5,5),D(5,5),DX(5,5),ND,NTHE(6),FRAC(10,5)            
C                                                                              
C                                                                              
C        SUBROUTINE CALCULATES THE AVERAGE THERMAL CAPACITY PER UNIT VOL       
C        COMPOSITE CONDUCTOR                                                   
C                                                                              
       CA=0.0                                                                  
       DO 2 K=1,ND                                                             
       NTH=NTHE(K)                                                             
       DO 1 L=1,NTH                                                            
       IF(T.GE.TEMPD(K,L))CM=(C(K,L)*(T**CX(K,L)))+(D(K,L)*(T**DX(K,L)))       
    1  CONTINUE                                                                
       CA=CA+FRAC(KC,K)*CM                                                     
    2  CONTINUE                                                                
       CP=CA                                                                   
       RETURN                                                                  
C                                                                              
       END                                                                     
C                                                                              
C     FILE MN.CQUENCH.RHO(NO)                                                  
C                                                                              
C     PART OF THE CULHAM VERSION OF QUENCH. SEE CQMAIN.                        
C                                                                              
       FUNCTION RHO(T,KC)                                                      
C                                                                              
         COMMON /b_1/TEMPD(5,5),A(5,5),AX(5,5),B(5,5),BX(5,5),                      
     1          C(5,5),CX(5,5),D(5,5),DX(5,5),ND,NTHE(6),FRAC(10,5)            
      common/bcu/r_rr(10),b_t(10)
C                                                                              
C                                                                              
C        SUBROUTINE CALCULATES THE AVERAGE ELECTRICAL RESISTIVITY OF A         
C        COMPOSITE CONDUCTOR                                                   
C                                                                              
       RA=0.0                                                                  
       DO 2 K=1,ND                                                             
       NTH=NTHE(K)                                                             
       DO 1 L=1,NTH                                                            
       IF(T.GE.TEMPD(K,L))RM=(A(K,L)*(T**AX(K,L)))+(B(K,L)*(T**BX(K,L)))       
       if(k.eq.2)rm=rcu(r_rr(kc),t,b_t(kc))
    1  CONTINUE                                                                
       RA=RA+FRAC(KC,K)/RM                                                     
    2  CONTINUE                                                                
       RHO=1.0/RA                                                              
       RETURN                                                                  
C                                                                              
       END