C***********************************************************************
  
      SUBROUTINE TYPE60(TIME,XIN,OUT,T,DTDT,PAR,INFO,ICNTRL,*)
  
C ----------------------------------------------------------------------
C
C       TYPE 60 : IDEAL GAS FLOW NOZZLE
C
C***********************************************************************

      DOUBLE PRECISION  XIN,OUT
      REAL    PAR
C     INTEGER IOSTAT
      DIMENSION  XIN(3),OUT(2),PAR(4)
C     DIMENSION  IOSTAT(2)
  
      PU=  XIN(1)
      PD=  XIN(2)
      T=   XIN(3)
  
      AC=  PAR(1)
      PO=  PAR(2)
      GAMA=PAR(3)
      R=   PAR(4)
  
C     T=T+273.15
      PU=PU+PO
      PD=PD+PO
      PCRIT=(2./(1.+GAMA))**(GAMA/(GAMA-1.))
      RATIO=PD/PU
      C3=(RATIO**(2./GAMA)-RATIO**((1.+GAMA)/GAMA))*2.0/(GAMA-1.0)
      IF (RATIO.LT.PCRIT) C3=(2./(1.0+GAMA))**((1.+GAMA)/(GAMA-1.))
      W=1000.*AC*PU*SQRT(GAMA*C3/(R*T))
C     TEXIT=T*RATIO**((GAMA-1.)/GAMA)-273.15
      TEXIT=T*RATIO**((GAMA-1.)/GAMA)
  
      OUT(1)= W
      OUT(2)= TEXIT
  
C     IOSTAT(1)=1
C     IOSTAT(2)=1
  
      RETURN 1
      END