source: CIVL/examples/fortran/nek5000/core/cmt/CMTDATA

C> \file CMTDATA include file for CMT variables and common blocks
C> \var U(lxyz,5,lelt) contains nodal values of the five conserved unknowns at
C> each of lx1*ly1*lz1 grid points for all elements
      real                 u(lx1,ly1,lz1,toteq,lelt) 
      common /solnconsvar/ u 
 
      logical            iffltr,    outflsub
     $                , ifcntfilt, ifrestart
      common /input10/   iffltr,    outflsub
     $                , ifcntfilt, ifrestart

      real               phig(lx1,ly1,lz1,lelt)
!     real               ja(lx1*ly1*lz1,ldim*ldim,lelv)
! I really want to set lxd=1 and declare a new rx and call it ja
      common /otherpvar/ phig!,ja

      real res1(lx1,ly1,lz1,lelt,toteq)
     $   , res2(lx1,ly1,lz1,lelt,toteq)
     $   , res3(lx1,ly1,lz1,toteq,lelt) 
      common /cmtres/ res1, res2, res3
      
      real              gradu(lx1*ly1*lz1,toteq,3)
      common /cmtgradu/ gradu     

      real convh (lxd*lyd*lzd,3)
     $   , diffh (lx1*ly1*lz1,3) ! dealias someday?
     $   , totalh(lxd*lyd*lzd,3)
      common /convectdata/ convh, diffh, totalh

      real               usrf(lx1,ly1,lz1,5)
      common /forcedata/ usrf 

      real c_sub_e,c_max,meshh(lelt),maxdiff,gridh(lx1*ly1*lz1,lelt),
     $     nu_s,epsebdg(lelt),maxres(lelt),se0(lelt)
      common /planned_Additions/ c_sub_e,c_max,meshh,maxdiff,
     $        gridh,nu_s,epsebdg,maxres,se0

c We assume we will be working with one gas. Therefore Rgas need not be 
c full filed. To begin with we assume calorically perfect gas.
c In near future we relax the constraint to thermally perfect gas.
c------------------------------------------------- -------------------
c NOTE --- CVG need not be stored on full field... removed in future.
c       THis because, Temp will be stored fullfield and we expect Cv to
c       have explicit temp dependence. At the begining of the time step
c       a scratch space will be use to iteratively compute Cv and T from 
c       internal energy.
c--------------------------------------------------------------------- 
      real                csound(lx1,ly1,lz1,lelt)
      common /cmtgasprop/ csound

      real            cpgref,cvgref,cvjref,gmaref,rgasref,molmass,
     $   MAref,MBref,CAref,CBref,   ! NTN for mixed model
     $    rhoidref,rho0ref,EN_IDref,EN_Oref,
     $    MWref,CWref,MCref,CCref,
     $    AAref,
     $    BBref,R1ref,R2ref,OMref,
     $    AA,BB,OM,CC,CW,se0const

      common /cmtgasref/cpgref,cvgref,cvjref,gmaref,rgasref,molmass,
     $    MAref,MBref,CAref,CBref,   ! NTN for mixed model
     $    rhoidref,rho0ref,EN_IDref,EN_Oref,
     $    MWref,CWref,MCref,CCref,
     $    AAref,
     $    BBref,R1ref,R2ref,OMref,
     $    AA,BB,OM,CC,CW,se0const

      integer              stage,nstage,iostep2
      common /tstepstage/ stage,nstage,iostep2

      real                  tcoef(3,3),time_cmt
      common /timestepcoef/ tcoef,time_cmt
      
      real              rho,phi,pres,cv,e_internal
     $                 ,t0in,p0in,asnd,rgas,cp,mu,lambda,molarmass
      common /nekuscmt/ rho,phi,pres,cv,e_internal
     $                 ,t0in,p0in,asnd,rgas,cp,mu,lambda,molarmass
! let kond=udiff

! index for gas density in vtrans, Q+-
      integer irho
      parameter (irho=1) ! better always be 1
! other VTRANS indices
      integer icp, icv
      parameter (icp =2) ! rho*cp
      parameter (icv =3) ! rho*cv
! indices in VDIFF
      integer imu, iknd, ilam, inus
      parameter (imu =1) ! viscosity
      parameter (iknd=2) ! thermal conductivity
      parameter (ilam=3) ! lambda, 2nd viscosity
      parameter (inus=4) ! nu_s, mass diffusion (EVM)
! indices for fields in U
      integer irg, irpu, irpv, irpw, iret
      parameter (irg =1)     ! phi_g*rho_g
      parameter (irpu=2)     ! phi_g*rho_g*u
      parameter (irpv=3)     ! phi_g*rho_g*v
      parameter (irpw=4)     ! phi_g*rho_g*w
      parameter (iret=toteq) ! phi_g*rho_g*E
! indices for fields in Q-+
      integer iux, iuy, iuz, ipr, ithm, isnd, iph, icvf, icpf, imuf, 
     $        ikndf, ilamf, iu1, iu2, iu3, iu4, iu5, nqq
      parameter (iux =2)     ! u
      parameter (iuy =3)     ! v
      parameter (iuz =4)     ! w
      parameter (ipr =5)     ! P
      parameter (ithm=6)     ! T
      parameter (isnd=7)     ! a
      parameter (iph =8)     ! phi_g
      parameter (icvf=9)     ! rho*cv
      parameter (icpf=10)    ! rho*cp
      parameter (imuf=11)    ! mu
      parameter (ikndf=12)   ! K
      parameter (ilamf=13)   ! lambda
      parameter (iu1=14)   ! U1
      parameter (iu2=15)   ! U2
      parameter (iu3=16)   ! U3
      parameter (iu4=17)   ! U4
      parameter (iu5=18)   ! U5
      parameter (nqq=18)

      logical time4av
      common /cmtflags/ time4av