!    -*- f90 -*-
! Note: the context of this file is case sensitive.

python module py_cea ! in 
    interface  ! in :py_cea
        block data  ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8, optional,dimension(100) :: atmwt=1.00794d0
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8, optional,dimension(100) :: valnce=1.
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8, optional :: avgdr=6.0221367d0
            real*8, optional :: boltz=1.380658d0
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8, optional :: pi=3.14159265d0
            real*8 :: pp
            real*8 :: r
            real*8, optional :: rr=8314.51d0
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character, optional,dimension(30,4),intent(c) :: fmt='(1x'
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character, optional,dimension(100,2),intent(c) :: symbol='h '
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            integer, parameter,optional :: maxpv=26
            integer, parameter,optional :: iosch=13
            integer, parameter,optional :: maxmat=50
            integer, parameter,optional :: ncol=8
            integer, parameter,optional :: maxng=500
            integer, parameter,optional :: maxr=24
            integer, parameter,optional :: nlines=50
            integer, parameter,optional :: iotrn=18
            integer, parameter,optional :: maxtr=30
            integer, parameter,optional :: maxt=51
            integer, parameter,optional :: maxngc=600
            integer, parameter,optional :: maxnc=300
            integer, parameter,optional :: iothm=14
            integer, parameter,optional :: maxmix=52
            integer, parameter,optional :: maxel=20
            integer, parameter,optional :: ioplt=15
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end block data 
        subroutine setinpline(ninp,line) ! in :py_cea:py_cea.f
            integer :: ninp
            character*132 :: line
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine setinpline
        subroutine py_cea(prefix,myfile,make_out_files,read_data) ! in :py_cea:py_cea.f
            character*196 :: prefix
            character*200 :: myfile
            integer :: make_out_files
            integer :: read_data
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine py_cea
        subroutine cphs ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
            entry allcon()
        end subroutine cphs
        subroutine deton ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine deton
        subroutine efmt(fone,aa,vx) ! in :py_cea:py_cea.f
            character*4 :: fone
            character*15 :: aa
            real*8 dimension(50) :: vx
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine efmt
        subroutine eqlbrm ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine eqlbrm
        subroutine frozen ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine frozen
        subroutine gauss ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine gauss
        subroutine hcalc ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine hcalc
        subroutine infree(readok,cin,ncin,lcin,dpin) ! in :py_cea:py_cea.f
            logical :: readok
            character dimension(600,15),intent(c) :: cin
            integer :: ncin
            integer dimension(600) :: lcin
            real*8 dimension(600) :: dpin
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine infree
        subroutine input(readok,caseok,ensert) ! in :py_cea:py_cea.f
            logical :: readok
            logical :: caseok
            character dimension(20,15),intent(c) :: ensert
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine input
        subroutine matrix ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine matrix
        subroutine newof ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine newof
        subroutine out1 ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
            entry out3()
            entry out2()
            entry out4()
        end subroutine out1
        subroutine react ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine react
        subroutine rktout ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine rktout
        subroutine rocket ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine rocket
        subroutine search ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
            entry readtr()
        end subroutine search
        subroutine seten ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine seten
        subroutine shck ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine shck
        subroutine thermp ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine thermp
        subroutine tranin ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine tranin
        subroutine tranp ! in :py_cea:py_cea.f
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine tranp
        subroutine utherm(readok) ! in :py_cea:py_cea.f
            logical :: readok
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine utherm
        subroutine utran(readok) ! in :py_cea:py_cea.f
            logical :: readok
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine utran
        subroutine varfmt(vx) ! in :py_cea:py_cea.f
            real*8 dimension(8) :: vx
            integer :: ip
            integer :: iplt
            integer :: it
            integer dimension(45) :: jcond
            integer dimension(20) :: jx
            integer :: nc
            integer :: ng
            integer :: ngp1
            integer :: nlm
            integer :: nplt
            integer :: nof
            integer :: nomit
            integer :: nonly
            integer :: np
            integer :: npr
            integer :: npt
            integer :: ngc
            integer :: nsert
            integer :: nspr
            integer :: nspx
            integer :: nt
            integer dimension(24) :: nfla
            integer dimension(300) :: ifz
            real*8 dimension(2) :: am
            real*8 dimension(20,2) :: b0p
            real*8 :: cpmix
            real*8 dimension(2) :: hpp
            real*8 dimension(2) :: vmin
            real*8 dimension(2) :: vpls
            real*8 :: wmix
            real*8 dimension(2) :: wp
            real*8 dimension(100) :: atmwt
            real*8 :: bcheck
            real*8 dimension(52) :: oxf
            real*8 dimension(26) :: p
            real*8 dimension(2) :: rh
            real*8 dimension(51) :: t
            real*8 dimension(26) :: v
            real*8 dimension(100) :: valnce
            real*8 dimension(20,600) :: a
            real*8 dimension(20) :: atwt
            real*8 :: avgdr
            real*8 :: boltz
            real*8 dimension(20) :: b0
            real*8 :: eqrat
            real*8 dimension(50,51) :: g
            real*8 :: hsub0
            real*8 :: oxfl
            real*8 :: pi
            real*8 :: pp
            real*8 :: r
            real*8 :: rr
            real*8 :: size_bn
            real*8 :: s0
            real*8 :: tln
            real*8 :: tm
            real*8 :: trace
            real*8 :: tt
            real*8 :: viscns
            real*8 :: vv
            real*8 dimension(50) :: x
            integer :: imat
            integer :: iq1
            integer :: isv
            integer :: jliq
            integer :: jsol
            integer :: lsave
            integer :: msing
            real*8 dimension(300,9) :: cft
            real*8 dimension(500,9,3) :: coef
            real*8 dimension(600) :: cp
            real*8 :: cpsum
            real*8 dimension(600) :: h0
            real*8 dimension(600) :: mu
            real*8 dimension(600) :: mw
            real*8 dimension(600) :: s
            real*8 dimension(2,300) :: temp
            real*8 dimension(4) :: tg
            integer :: makeio
            integer :: read_d
            integer :: nlinesset
            integer :: nline
            real*8 dimension(600) :: deln
            real*8 dimension(600,8) :: en
            real*8 dimension(600) :: enln
            real*8 :: enn
            real*8 :: ennl
            real*8 :: enlsav
            real*8 :: ensave
            real*8 dimension(600) :: sln
            real*8 :: sumn
            character*15 :: case_bn
            character dimension(20,2),intent(c) :: elmt
            character dimension(24,15),intent(c) :: energy
            character dimension(30,4),intent(c) :: fmt
            character dimension(24,8),intent(c) :: fox
            character dimension(601,15),intent(c) :: omit
            character dimension(20,15),intent(c) :: pltvar
            character dimension(601,15),intent(c) :: prod
            character dimension(24,12,2),intent(c) :: ratom
            character dimension(24,15),intent(c) :: rname
            character dimension(100,2),intent(c) :: symbol
            character*10 :: thdate
            character*200 :: pfile
            character dimension(50,132),intent(c) :: inplines
            real*8 dimension(24) :: dens
            real*8 dimension(24) :: enth
            real*8 dimension(24) :: pecwt
            real*8 dimension(24) :: rmw
            real*8 dimension(24,12) :: rnum
            real*8 dimension(24) :: rtemp
            integer dimension(24) :: jray
            integer :: nreac
            logical :: convg
            logical dimension(8) :: debug
            logical :: detdbg
            logical :: detn
            logical :: eql
            logical :: gonly
            logical :: hp
            logical :: ions
            logical :: massf
            logical :: moles
            logical :: newr
            logical :: pderiv
            logical :: shock
            logical :: short_bn
            logical :: siunit
            logical :: sp
            logical :: tp
            logical :: trnspt
            logical :: vol
            real*8 dimension(8) :: coneql
            real*8 dimension(8) :: confro
            real*8 dimension(8) :: cpeql
            real*8 dimension(8) :: cpfro
            real*8 dimension(8) :: preql
            real*8 dimension(8) :: prfro
            real*8 dimension(8) :: vis
            real*8 dimension(8) :: u1
            real*8 dimension(8) :: mach1
            real*8 :: a1
            real*8 :: gamma1
            logical :: incdeq
            logical :: incdfz
            logical :: refleq
            logical :: reflfz
            logical :: shkdbg
            integer :: nsk
            real*8 dimension(8) :: cpr
            real*8 dimension(8) :: dlvpt
            real*8 dimension(8) :: dlvtp
            real*8 dimension(8) :: gammas
            real*8 dimension(8) :: hsum
            real*8 dimension(8) :: ppp
            real*8 dimension(8) :: ssum
            real*8 dimension(8) :: totn
            real*8 dimension(8) :: ttt
            real*8 dimension(8) :: vlm
            real*8 dimension(8) :: wm
            real*8 dimension(500,20) :: pltout
            real*8 dimension(30) :: cprr
            real*8 dimension(30) :: con
            real*8 dimension(30,30) :: eta
            real*8 dimension(30) :: wmol
            real*8 dimension(30) :: xs
            real*8 dimension(30,30) :: stc
            integer dimension(30) :: ind
            integer dimension(20) :: jcm
            integer :: nm
            integer :: nr
            integer :: ntape
            real*8 :: acat
            real*8 dimension(8) :: aeat
            real*8 dimension(8) :: app
            real*8 :: awt
            real*8 :: cstr
            real*8 :: ma
            real*8 dimension(16) :: pcp
            real*8 dimension(8) :: sonvel
            real*8 dimension(8) :: spim
            real*8 dimension(13) :: subar
            real*8 dimension(8) :: supar
            real*8 dimension(8) :: vmoc
            real*8 dimension(8) :: vaci
            real*8 :: tcest
            logical :: area
            integer :: iopt
            logical :: debugf
            logical :: fac
            logical :: froz
            integer :: isup
            integer :: nfz
            integer :: npp
            integer :: nsub
            integer :: nsup
            logical :: page1
            logical :: rkt
            common /indx/ ip,iplt,it,jcond,jx,nc,ng,ngp1,nlm,nplt,nof,nomit,nonly,np,npr,npt,ngc,nsert,nspr,nspx,nt,nfla,ifz
            common /inpt/ am,b0p,cpmix,hpp,vmin,vpls,wmix,wp,atmwt,bcheck,oxf,p,rh,t,v,valnce
            common /miscr/ a,atwt,avgdr,boltz,b0,eqrat,g,hsub0,oxfl,pi,pp,r,rr,size_bn,s0,tln,tm,trace,tt,viscns,vv,x
            common /misci/ imat,iq1,isv,jliq,jsol,lsave,msing
            common /therm/ cft,coef,cp,cpsum,h0,mu,mw,s,temp,tg
            common /cet_comm/ makeio,read_d,nlinesset,nline
            common /comp/ deln,en,enln,enn,ennl,enlsav,ensave,sln,sumn
            common /cdata/ case_bn,elmt,energy,fmt,fox,omit,pltvar,prod,ratom,rname,symbol,thdate,pfile
            common /cet_inp/ inplines
            common /reactn/ dens,enth,pecwt,rmw,rnum,rtemp,jray,nreac
            common /miscl/ convg,debug,detdbg,detn,eql,gonly,hp,ions,massf,moles,newr,pderiv,shock,short_bn,siunit,sp,tp,trnspt,vol
            common /trpts/ coneql,confro,cpeql,cpfro,preql,prfro,vis
            common /shocks/ u1,mach1,a1,gamma1,incdeq,incdfz,refleq,reflfz,shkdbg,nsk
            common /prtout/ cpr,dlvpt,dlvtp,gammas,hsum,ppp,ssum,totn,ttt,vlm,wm,pltout
            common /trnp/ cprr,con,eta,wmol,xs,stc,ind,jcm,nm,nr,ntape
            common /rockt/ acat,aeat,app,awt,cstr,ma,pcp,sonvel,spim,subar,supar,vmoc,vaci,tcest,area,iopt,debugf,fac,froz,isup,nfz,npp,nsub,nsup,page1,rkt
        end subroutine varfmt
    end interface 
end python module py_cea

! This file was auto-generated with f2py (version:2).
! See http://cens.ioc.ee/projects/f2py2e/
