C      MAIN PROGRAM--QUENCH CALCULATION, MULTIPLE COILS.                       
C   CONCATENATION FILE FOR R.L. VERSION OF QUENCH WITH UP TO A MAXIMUM         
C   OF 40 COILS .  QUENCH INITIATED FROM HOT BORE TUBES IS AVAILABLE           
C                                                                              
C   CHANGED TO DO 10 COILS FOR UP TO 400 STEPS FOR VAX11/780
C   6/11/82
C
c   changed to 800 steps 11/21/03
C     FILE MN.CQUENCH.MAIN(NO)                                                 
C                                                                              
C     MAIN ROUTINE FROM THE CULHAM VERSION OF QUENCH RECEIVED FROM             
C     CULHAM IN JANUARY 1973, AS MODIFIED BY M.NEWMAN IN THE SAME              
C     MONTH TO COMPILE ON THE 360/195.                                         
C                                                                              
C     EACH ROUTINE HAS BEEN MADE INTO AN ELECTRIC FILE. THERE IS A             
C     CONCATENATION FILE FOR ALL THE ROUTINES CALLED CQUENCH.LIST1             
C                                                                              
C     PROGRAM IN FOLLOWING SEGMENTS:
C     BURN - MAIN ROUTINE
C     BD - DATA1,DATA2,DATA3
C     MS - MFLIP,SPLIC,SPLIT
C     CV - CALC,NORVOL,OVLAP,NEWVOL
C     DDD - DA01A,DYBDX,DEKM
C     BF - CP,RHO
C     rcut- replaces RHO for Cu, calculates rho(RRR,B)
C        NC        = NUMBER OF COILS IN MAGNET SYSTEM                          
c      common /b_1/rd(225),nd,nthe(6),frac(50)
c         COMMON TEMPD(5,5),A(5,5),AX(5,5),B(5,5),BX(5,5),                      
c     1          C(5,5),CX(5,5),D(5,5),DX(5,5),ND,NTHE(6),FRAC(10,5)            
c      common /b_2/rd2(240)
c         COMMON AMPSSP(10),AMPSCR(10),OHM(10),RT(10),                          
c     1          FK(10,10),RK(10,10)                                            
c      common /b_3/rd3(4070),ind(10),dt
c         COMMON AMPS(10),VELZ(10),P1(10),P2(10),P3(10),P4(10),                 
c     1          X(10),VOL(10,800),IND(10),DT                                   
      common /b_4/nc,nt,rd4(4051)
c         COMMON NC,NT,AREA(10),AREARS(10),TSP(10),TQ(10),IQ(10),               
c     1          TEMPZP,TEMP(10,800)                                            
c      common /b_5/rd5(60)
c      COMMON TOTVOL(10),ARCON(10),PROVEL(6,5),ALPHA(10)                        
c      common /b_6/rd6(48)
c      COMMON A1(6),B1(6),C1(6),D1(6),A2(6),B2(6),C2(6),D2(6)                   
c      common /b_7/rd7(48)
c      COMMON A3(6),B3(6),C3(6),D3(6),A4(6),B4(6),C4(6),D4(6)                   
c      common /b_8/rd8(48)
c      COMMON A5(6),B5(6),C5(6),D5(6),A6(6),B6(6),C6(6),D6(6)                   
c      common /b_9/rd9(36)
c      COMMON PROV1(6),PROV2(6),PROV3(6),PROV4(6),PROV5(6),PROV6(6)             
c      common /b_10/rd10(26)
c      COMMON FIELD1(6),FIELD2(10),BETA(10)                                     
      DIMENSION CARD(20)   
	CHARACTER*80 FILENAME                                                    
C                                                                              
C   
      CALL GETENV('DAT',FILENAME)
      OPEN(UNIT=5,NAME=FILENAME,STATUS='OLD')
      CALL GETENV('OUT',FILENAME)
      OPEN(UNIT=6,NAME=FILENAME,STATUS='NEW')
       CALL GETENV('VUE',FILENAME)
      OPEN(UNIT=7,NAME=FILENAME,STATUS='NEW')
      CALL GETENV('MAT',FILENAME)
      OPEN(UNIT=8,NAME=FILENAME,STATUS='OLD')
      OPEN(UNIT=99,FILE='TEMP.OUT',STATUS='SCRATCH')
      ITERM=0                                                                  
    1 CONTINUE                                                                 
      READ(5,100)CARD                                                          
  100 FORMAT(20A4)                                                             
      WRITE(6,102)CARD                                                         
      WRITE(7,102)CARD
  102 FORMAT('0'20A4)                                                          
      READ(5,*)NC                                                            
      IF(NC.GT.1)GO TO 3
      CALL GETENV('PLT',FILENAME)
      OPEN(UNIT=1,NAME=FILENAME,STATUS='NEW')
      CALL GETENV('PLA',FILENAME)
      OPEN(UNIT=2,NAME=FILENAME,STATUS='NEW')
      CALL GETENV('PLV',FILENAME)
      OPEN(UNIT=3,NAME=FILENAME,STATUS='NEW')
      CALL GETENV('PLN',FILENAME)
      OPEN(UNIT=4,NAME=FILENAME,STATUS='NEW')
3     IF(ITERM.GE.1)GO TO 2                                                    
       CALL DATA1                                                              
    2  CALL DATA2                                                              
       CALL DATA3                                                              
       CALL CALC                                                               
      READ(5,*)ITERM                                                         
      IF(ITERM.GT.0)GO TO 1                                                    
      CALL EXIT
  999 FORMAT (20F)                                                             
       END