Rocstar-legacy/SPEC__RFLU__ModChemistry_8F90_source.html

 ! *********************************************************************

 ! * Rocstar Simulation Suite                                          *

 ! * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *

 ! *                                                                   *

 ! * Illinois Rocstar LLC                                              *

 ! * Champaign, IL                                                     *

 ! * www.illinoisrocstar.com                                           *

 ! * sales@illinoisrocstar.com                                         *

 ! *                                                                   *

 ! * License: See LICENSE file in top level of distribution package or *

 ! * http://opensource.org/licenses/NCSA                               *

 ! *********************************************************************

 ! *********************************************************************

 ! * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *

 ! * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *

 ! * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *

 ! * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *

 ! * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *

 ! * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *

 ! * Arising FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *

 ! * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *

 ! *********************************************************************

 ! ******************************************************************************

 !

 ! Purpose: Simple chemistry model for burning crack application.

 !

 ! Description: None.

 !

 ! Notes:

 !   1. Based on collection of routines developed by Luca Massa.

 !

 ! ******************************************************************************

 !

 ! $Id: SPEC_RFLU_ModChemistry.F90,v 1.6 2008/12/06 08:44:40 mtcampbe Exp $

 !

 ! Copyright: (c) 2004-2005 by the University of Illinois

 !

 ! ******************************************************************************


 MODULE spec_rflu_modchemistry


   USE moddatatypes

   USE modparameters

   USE moderror

   USE modglobal, ONLY: t_global

   USE modgrid, ONLY: t_grid

   USE moddatastruct, ONLY: t_region


   IMPLICIT NONE


   PRIVATE

   PUBLIC :: spec_rflu_integratechemsrcterm


 ! ******************************************************************************

 ! Declarations and definitions

 ! ******************************************************************************


   INTEGER, PARAMETER :: MAXEQN = 2,MAXRCT = 1

   INTEGER :: nreaction,ntemp_spec,maxnwtn

   REAL(RFREAL) :: eps,epy,beta,cpGas

   REAL(RFREAL) :: damkholer(MAXRCT),mx(MAXEQN),mn(MAXEQN),nn(MAXRCT), &

                   pre(MAXRCT),qg(MAXRCT),thetay(MAXRCT),VEC(MAXEQN)

   REAL(RFREAL) :: MAT(MAXEQN,MAXEQN),yf(MAXEQN,MAXRCT)

   REAL(RFREAL) :: ratey(MAXRCT)


 ! ******************************************************************************

 ! Routines

 ! ******************************************************************************


   CONTAINS


 ! ******************************************************************************

 !

 ! Purpose: Integrate source term for mixture and species using ODE solver.

 !

 ! Description: None.

 !

 ! Input:

 !   pRegion     Pointer to region

 !

 ! Output: None.

 !

 ! Notes:

 !   1. Numerical method is due to Luca Massa and Mark Short.

 !

 ! ******************************************************************************


   SUBROUTINE spec_rflu_integratechemsrcterm(pRegion,CALLFLAG)


     USE modinterfaces, ONLY: mixtperf_eo_grtuvw, &

                              mixtperf_g_cpr, &

                              mixtperf_r_m, &

                              mixtperf_d_prt


 ! ******************************************************************************

 !   Definitions and declarations

 ! ******************************************************************************


 ! ==============================================================================

 !   Arguments

 ! ==============================================================================


     TYPE(t_region), POINTER :: pregion


 ! ==============================================================================

 !   Locals

 ! ==============================================================================


     CHARACTER(CHRLEN) :: errorstring,rcsidentstring

     INTEGER :: icg,indcp,indmol,ispec,nout,ncumul,ibc,iec,callflag

     REAL(RFREAL) :: fact,ggas,ir,molmass,p,r,rgas,t,u,v,w,rvcumul,x,y,aaa,bbb

     REAL(RFREAL), DIMENSION(MAXEQN) :: f,f_old

     REAL(RFREAL), DIMENSION(:,:), POINTER :: pcvmixt,pcvspec,pdvmixt,pgvmixt,mrhs,srhs

     REAL(RFREAL), DIMENSION(:), POINTER :: vol

     TYPE(t_grid), POINTER :: pgrid

     TYPE(t_global), POINTER :: global


 ! ******************************************************************************

 !   Start

 ! ******************************************************************************


     rcsidentstring = '$RCSfile: SPEC_RFLU_ModChemistry.F90,v $ $Revision: 1.6 $'


     global => pregion%global


     CALL registerfunction(global,'SPEC_RFLU_IntegrateChemSrcTerm',&

   'SPEC_RFLU_ModChemistry.F90')


 ! ******************************************************************************

 !   Set pointers and variables

 ! ******************************************************************************


     pgrid => pregion%grid


     pcvmixt => pregion%mixt%cv

     pdvmixt => pregion%mixt%dv

     pgvmixt => pregion%mixt%gv

     pcvspec => pregion%spec%cv

     mrhs    => pregion%mixt%rhs

     srhs    => pregion%spec%rhs

     vol     => pregion%grid%vol


     indcp  = pregion%mixtInput%indCp

     indmol = pregion%mixtInput%indMol


     fact = 1.0_rfreal

     nout = 100


     CALL set_rate(pregion) ! Set rate constants


 ! ******************************************************************************

 !   Checks

 ! ******************************************************************************


     IF ( pregion%mixt%cvState /= cv_mixt_state_cons ) THEN

       CALL errorstop(global,err_cv_state_invalid,__line__)

     END IF ! pRegion%mixt%cvState


     IF ( pregion%specInput%nSpecies /= maxeqn-1 ) THEN

       WRITE(errorstring,'(I3)') maxeqn

       CALL errorstop(global,err_spec_maxeqn,__line__,errorstring)

     END IF ! pRegion%specInput%nSpecies


 ! ******************************************************************************

 !   Loop limits and sum initialization

 ! ******************************************************************************


     ibc = 1

     iec = merge(pgrid%nCellsTot,pgrid%nCells,callflag == 0)

     ncumul = 0

     rvcumul = 0.0_rfreal


 ! ******************************************************************************

 !   Integrate source term

 ! ******************************************************************************


     DO icg =  ibc,iec

       r =  pcvmixt(cv_mixt_dens,icg)

       ir = 1.0_rfreal/r


       x = pgrid%cofg(xcoord,icg)

       y = pgrid%cofg(ycoord,icg)

       t = pdvmixt(dv_mixt_temp,icg)

       p = pdvmixt(dv_mixt_pres,icg)


       molmass = pgvmixt(gv_mixt_mol,indmol*icg)

       cpgas   = pgvmixt(gv_mixt_cp ,indcp *icg)


       rgas = mixtperf_r_m(molmass)

       ggas = mixtperf_g_cpr(cpgas,rgas)


 ! ==============================================================================

 !     Call integration routine

 ! ==============================================================================


       !CALL rate_split(f,f_old,r,p,fact,global%dtMin,nout)

       pre(1) = damkholer(1)*p**nn(1)


 ! ==============================================================================

 !     Update state vectors. NOTE have conserved state vector for the mixture,

 !     so need to update total internal energy because it is affected by the

 !     temperature. The density is assumed to stay constant during this update.

 !     The dependent variables will be updated outside this routine.

 ! ==============================================================================


       if(callflag == 1) then

          ratey(1) = pre(1) * pcvspec(1,icg) * exp(-thetay(1)/t) / r * global%dtMin

          pcvmixt(cv_mixt_ener,icg) = pcvmixt(cv_mixt_ener,icg) + ratey(1)*yf(1,1)

          pcvspec(           1,icg) = pcvspec(           1,icg) + ratey(1)*yf(2,1)

       elseif(callflag == 0) then

          ratey(1) = pre(1) * pcvspec(1,icg) * exp(-thetay(1)/t) / r * vol(icg)


          mrhs(cv_mixt_ener,icg) = mrhs(cv_mixt_ener,icg) - ratey(1)*yf(1,1)

          srhs(1,icg) = srhs(1,icg) - ratey(1)*yf(2,1)

 !         if (abs(pRegion%grid%cofg(YCOORD,icg)) > 4.5200d-3) Mrhs(CV_MIXT_ENER,icg) = 0_RFREAL

 !         if (abs(pRegion%grid%cofg(YCOORD,icg)) > 4.5200d-3) Srhs(:,icg) = 0_RFREAL

       endif


     END DO ! icg


 ! ******************************************************************************

 !   End

 ! ******************************************************************************


     CALL deregisterfunction(global)


   END SUBROUTINE spec_rflu_integratechemsrcterm


 ! *******************************************************************************

 !

 ! Purpose: Set rate constants.

 !

 ! Description: None.

 !

 ! Input:

 !  pRegion      Pointer to region data

 !

 ! Output: None.

 !

 ! Notes: None.

 !

 ! *******************************************************************************


   SUBROUTINE set_rate(pRegion)


     IMPLICIT NONE


 ! ******************************************************************************

 !   Definitions and declarations

 ! ******************************************************************************


 ! ==============================================================================

 !   Arguments

 ! ==============================================================================


     TYPE(t_region), POINTER :: pregion


 ! ==============================================================================

 !   Locals

 ! ==============================================================================


     INTEGER :: n


 ! ******************************************************************************

 !   Start

 ! ******************************************************************************


     maxnwtn = 1

     eps = 1.0e-8_rfreal

     epy = 2.0e-4_rfreal


     nreaction = 1

     ntemp_spec = 1+pregion%specInput%nSpecies


     beta = 1.0_rfreal


     thetay(1) =  15000_rfreal


     qg(1) = 2511600.0_rfreal


     nn(1) = 0.0_rfreal


     damkholer(1) = 1174775.7803_rfreal


 !!!one reaction, one species

     yf(1,1) = qg(1)


     yf(2,1) = -1.0_rfreal


 !TEMPERATURE AND SPECIES BOUNDS !IN REACTION STEP INTEGRATION


     mx(1) =   6000.0_rfreal

     mn(1) = 0.0_rfreal


     DO n = 2,ntemp_spec

        mn(n) = -10.0_rfreal

        mx(n) = 1.0_rfreal - mn(n)

     END DO ! n


 ! ******************************************************************************

 !   End

 ! ******************************************************************************


   END SUBROUTINE set_rate


 ! ******************************************************************************

 !

 ! Purpose: Integrate source terms using ODE solver.

 !

 ! Description: None.

 !

 ! Input:

 !   old_f

 !   rho

 !   time_step

 !   nout

 !

 ! Output:

 !   f

 !

 ! Notes: None.

 !

 ! ******************************************************************************


   SUBROUTINE rate_split(f,old_f,rho,press,fact,time_step,nout)


     IMPLICIT NONE


 ! ******************************************************************************

 !   Definitions and declarations

 ! ******************************************************************************


 ! ==============================================================================

 !   Arguments

 ! ==============================================================================


     REAL(RFREAL), INTENT(IN) :: fact,press,rho,time_step

     REAL(RFREAL), DIMENSION(MAXEQN) :: old_f,f


 ! ==============================================================================

 !   Locals

 ! ==============================================================================


     LOGICAL :: outcome

     INTEGER :: eqn,i,icnvrgd,nout

     REAL(RFREAL) :: toll,dt_over_rho


 ! ******************************************************************************

 !   Start

 ! ******************************************************************************


     dt_over_rho = time_step/rho*fact


     DO eqn = 1, ntemp_spec

       old_f(eqn) = f(eqn)

     END DO ! eqn


 ! ******************************************************************************

 !   Integrate

 ! ******************************************************************************


     CALL odesolve(old_f(1:ntemp_spec),f(1:ntemp_spec),press, &

                   ntemp_spec,dt_over_rho,nout,outcome)


     IF ( outcome ) THEN

       WRITE(*,*)'DIVERGENCE IN RATE',old_f(1:ntemp_spec),f(1:ntemp_spec), &

                                      dt_over_rho

       f(1:ntemp_spec)=old_f(1:ntemp_spec)

       icnvrgd = 0

     ELSE

       icnvrgd = 1

     END IF ! outcome


 ! ******************************************************************************

 !   End

 ! ******************************************************************************


   END SUBROUTINE rate_split


 ! ******************************************************************************

 !

 ! Purpose: None.

 !

 ! Description: None.

 !

 ! Input:

 !   yin         Old solution

 !   press       Pressure

 !   neqin       Number of equations

 !   dtout       Time step

 !   nout        Number of subiterations

 !

 ! Output:

 !   y           New solution

 !   outcome     Success flag (TRUE if diverged!!!)

 !

 ! Notes: None.

 !

 ! ******************************************************************************


   SUBROUTINE odesolve(yin,y,press,neqin,dtout,nout,outcome)


     IMPLICIT NONE


 ! ******************************************************************************

 !   Definitions and declarations

 ! ******************************************************************************


 ! ==============================================================================

 !   Arguments

 ! ==============================================================================


     LOGICAL :: outcome

     INTEGER, INTENT(IN) :: neqin,nout

     REAL(RFREAL) :: dtout,press

     REAL(RFREAL) :: y(neqin),yin(neqin)


 ! ==============================================================================

 !   Locals

 ! ==============================================================================


     INTEGER :: icol,inwtn,iout,irow

     REAL(RFREAL) :: t,tout

     REAL(RFREAL) :: pp(maxeqn),ytmp(maxeqn)


 ! ******************************************************************************

 !   Start

 ! ******************************************************************************


     dtout = dtout/REAL(nout,kind=rfreal)


     t    = 0.0_rfreal

     tout = t

     y    = yin

     ytmp = yin

     pp   = yin


     DO iout = 1,nout

       DO inwtn = 1, maxnwtn

         call rate_stiff(neqin, tout, y, press, vec)

         call drate_stiff(neqin, tout, y, press, 1,1,mat,neqin)


         DO irow = 1,neqin

           vec(irow) = ytmp(irow) - y(irow) + vec(irow)*dtout


           DO icol = 1,neqin

             IF ( icol == irow ) THEN

               mat(irow,irow) =  mat(irow,irow)*dtout - 1.0_rfreal

             ELSE

               mat(irow,icol) = mat(irow,icol)*dtout

             END IF ! icol

           END DO ! icol

         END DO ! irow


         CALL lusolve(mat(1:neqin,1:neqin),vec(1:neqin),neqin)

         y(1:neqin) = y(1:neqin) - vec(1:neqin)

       END DO ! inwtn


       tout = tout + dtout

       ytmp = y


 !      WRITE(23,*)iout,tout,y(1)

 !      WRITE(6,*)iout,tout,y(1)

     END DO ! iout


     outcome = (.NOT. y(1) < mx(1)) .OR. (.NOT. y(1) > mn(1))


     IF ( outcome ) THEN

        DO irow = 1,neqin

           WRITE (stdout,*)irow,'TXY',yin(irow),y(irow)

        END DO ! irow

     END IF ! outcome


 300 format(////' STATE FLAG =',i3,/,'FAILED INEGRATION')


 ! ******************************************************************************

 !   End

 ! ******************************************************************************


   END SUBROUTINE odesolve


 ! ******************************************************************************

 !

 ! Purpose: Set up rate for two step kinetics and for fine AP/binder mixture

 !

 ! Description: None.

 !

 ! Input: None.

 !

 ! Output: None.

 !

 ! Notes: None.

 !

 ! ******************************************************************************


 ! ***************************************************

 !     ff(1)  = T    (gas phase temperature)

 !     ff(2)  = F    (gas phase fuel Y)

 !     ff(3)  = OX   (gas phase oxidizer X)

 !     ff(4)  = Z    (gas phase intermediate Z)

 ! ********************************************************************


   SUBROUTINE rate_stiff(neqin,time_stiff,ff,press,rrate)


     IMPLICIT NONE


 ! ******************************************************************************

 !   Definitions and declarations

 ! ******************************************************************************


 ! ==============================================================================

 !   Arguments

 ! ==============================================================================


     INTEGER :: neqin

     REAL(RFREAL) :: press,time_stiff

     REAL(RFREAL) :: ff(neqin),rrate(neqin)


 ! ==============================================================================

 !   Locals

 ! ==============================================================================


     INTEGER :: m,n

     REAL(RFREAL) :: ratey(maxrct)

     REAL(RFREAL) :: dratey(maxrct,maxeqn)


 ! ******************************************************************************

 !   Start

 ! ******************************************************************************


     pre(1) = damkholer(1)*press**nn(1)


     ratey(1) = ff(2)*exp(-thetay(1)/ff(1))


     ratey(1) = pre(1)*ratey(1)


 !   yf(species,reaction), m species_eq, n=species_var


     DO m = 1,neqin

       rrate(m) = sum(yf(m,1:nreaction)*ratey(1:nreaction))

     END DO ! m


 ! ******************************************************************************

 !   End

 ! ******************************************************************************


   END SUBROUTINE rate_stiff


 ! ******************************************************************************

 !

 ! Purpose: None.

 !

 ! Description: None.

 !

 ! Input: None.

 !

 ! Output: None.

 !

 ! Notes: None.

 !

 ! ******************************************************************************


   SUBROUTINE drate_stiff(neqin,time_stiff,ff,press,ml,imu,drate,nrowpd)


     IMPLICIT NONE


 ! ******************************************************************************

 !   Definitions and declarations

 ! ******************************************************************************


 ! ==============================================================================

 !   Arguments

 ! ==============================================================================


     INTEGER :: neqin,ml,imu

     REAL(RFREAL) :: press,time_stiff

     REAL(RFREAL) :: ff(neqin),drate(nrowpd,neqin)

     REAL(RFREAL) :: dratey(maxrct,maxeqn)


 ! ==============================================================================

 !   Locals

 ! ==============================================================================


     INTEGER :: m,n,nrowpd

     REAL(RFREAL) :: tsqinv

     REAL(RFREAL) :: ratey(maxrct)


 ! ******************************************************************************

 !   Start

 ! ******************************************************************************


     pre(1) = damkholer(1)*press**nn(1)


     tsqinv = 1.0_rfreal/(ff(1)*ff(1))


     ratey(1) = ff(2)*exp(-thetay(1)/ff(1))


     ratey(1) = pre(1)*ratey(1)


     dratey(1,1) = thetay(1)*ratey(1)*tsqinv

     dratey(1,2) = ratey(1)/merge(eps,ff(2),ff(2)==0.0_rfreal)


 !   yf(species,reaction), m species_eq, n=species_var


     DO m = 1,neqin

        DO n = 1,ntemp_spec

           drate(m,n) = sum(yf(m,1:nreaction)*dratey(1:nreaction,n))

        END DO ! n

     END DO ! m


 ! ******************************************************************************

 !   End

 ! ******************************************************************************


   END SUBROUTINE drate_stiff


 ! ******************************************************************************

 !

 ! Purpose: Solve system of linear equations Ax = f by LU decomposition.

 !

 ! Description: None.

 !

 ! Input:

 !   A           Matrix of size n by n

 !   f           Vector of size n

 !   n           Size of matrix and vector

 !

 ! Output:

 !   f           Solution to Ax=f

 !

 ! Notes: None.

 !

 ! ******************************************************************************


   SUBROUTINE lusolve(A,f,n)


 ! ******************************************************************************

 !   Definitions and declarations

 ! ******************************************************************************


 ! ==============================================================================

 !   Arguments

 ! ==============================================================================


     INTEGER :: n

     REAL(RFREAL) :: f(n)

     REAL(RFREAL) :: a(n,n)


 ! ==============================================================================

 !   Locals

 ! ==============================================================================


     INTEGER :: i,j

     REAL(RFREAL) :: d

     REAL(RFREAL), DIMENSION(MAXEQN) :: b,y

     REAL(RFREAL), DIMENSION(MAXEQN,MAXEQN) :: u,l


 ! ******************************************************************************

 !   Start

 ! ******************************************************************************


     l = 1.0_rfreal

     u = 0.0_rfreal


     u(1,1)   = a(1,1)

     u(1,2:n) = a(1,2:n)

     l(2:n,1) = a(2:n,1)/a(1,1)


     DO i = 2,n-1

        u(i,i) = a(i,i) - sum(l(i,1:i-1)*u(1:i-1,i))

        d = 1.0_rfreal/u(i,i)


        DO j = i+1,n

           u(i,j) =  a(i,j) - sum(l(i,1:i-1)*u(1:i-1,j))

           l(j,i) = (a(j,i) - sum(l(j,1:i-1)*u(1:i-1,i)))*d

        END DO ! j

     END DO ! i


     i=n

     u(i,i) = a(i,i) - sum(l(i,1:i-1)*u(1:i-1,i))


     y(1) = f(1)


     DO i = 2,n

        y(i) = f(i) - sum(l(i,1:i-1)*y(1:i-1))

     END DO ! i


     f(n) = y(n)/u(n,n)

     DO i = n-1,1,-1

        f(i) = (y(i) - sum(u(i,i+1:n)*f(i+1:n)))/u(i,i)

     END DO ! i


 ! ******************************************************************************

 !   End

 ! ******************************************************************************


   END SUBROUTINE lusolve


 END MODULE spec_rflu_modchemistry


 ! ******************************************************************************

 !

 ! RCS Revision history:

 !

 ! $Log: SPEC_RFLU_ModChemistry.F90,v $

 ! Revision 1.6  2008/12/06 08:44:40  mtcampbe

 ! Updated license.

 !

 ! Revision 1.5  2008/11/19 22:17:53  mtcampbe

 ! Added Illinois Open Source License/Copyright

 !

 ! Revision 1.4  2006/04/07 15:19:25  haselbac

 ! Removed tabs

 !

 ! Revision 1.3  2005/06/08 16:06:23  haselbac

 ! Bug fix: Correcting incomplete check-in

 !

 ! Revision 1.2  2005/06/06 14:24:11  haselbac

 ! Adapted to Lucas changes

 !

 ! Revision 1.1  2004/04/01 21:22:41  haselbac

 ! Initial revision

 !

 ! ******************************************************************************


Color::r
unsigned char r() const
Definition: Color.h:68

modglobal
Definition: ModGlobal.F90:39

Aff_transformation_2::m
FT m(int i, int j) const
Definition: Aff_transformation_2.h:228

cimg_library::CImg::sum
Tfloat sum() const
Return the sum of all the pixel values in an image.
Definition: CImg.h:13022

spec_rflu_modchemistry::rate_split
subroutine rate_split(f, old_f, rho, press, fact, time_step, nout)
Definition: SPEC_RFLU_ModChemistry.F90:337

mixtperf_r_m
real(rfreal) function mixtperf_r_m(M)
Definition: MixtPerf_R.F90:54

modparameters
Definition: ModParameters.F90:42

d
const NT & d
Definition: rational_rotation.h:73

y
void int int REAL REAL * y
Definition: read.cpp:74

moderror::registerfunction
subroutine registerfunction(global, funName, fileName)
Definition: ModError.F90:449

Color::b
unsigned char b() const
Definition: Color.h:70

modgrid::t_grid
Definition: ModGrid.F90:56

mixtperf_d_prt
real(rfreal) function mixtperf_d_prt(P, R, T)
Definition: MixtPerf_D.F90:71

v
*********************************************************************Illinois Open Source License ****University of Illinois NCSA **Open Source License University of Illinois All rights reserved ****Developed free of to any person **obtaining a copy of this software and associated documentation to deal with the Software without including without limitation the rights to and or **sell copies of the and to permit persons to whom the **Software is furnished to do subject to the following this list of conditions and the following disclaimers ****Redistributions in binary form must reproduce the above **copyright this list of conditions and the following **disclaimers in the documentation and or other materials **provided with the distribution ****Neither the names of the Center for Simulation of Advanced the University of nor the names of its **contributors may be used to endorse or promote products derived **from this Software without specific prior written permission ****THE SOFTWARE IS PROVIDED AS WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **ARISING OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE v
Definition: roccomf90.h:20

spec_rflu_modchemistry::set_rate
subroutine set_rate(pRegion)
Definition: SPEC_RFLU_ModChemistry.F90:252

i
blockLoc i
Definition: read.cpp:79

spec_rflu_modchemistry::rate_stiff
subroutine rate_stiff(neqin, time_stiff, ff, press, rrate)
Definition: SPEC_RFLU_ModChemistry.F90:522

x
void int int REAL * x
Definition: read.cpp:74

spec_rflu_modchemistry::lusolve
subroutine lusolve(A, f, n)
Definition: SPEC_RFLU_ModChemistry.F90:663

n
const NT & n
Definition: rational_rotation.h:72

spec_rflu_modchemistry::drate_stiff
subroutine drate_stiff(neqin, time_stiff, ff, press, ml, imu, drate, nrowpd)
Definition: SPEC_RFLU_ModChemistry.F90:587

modinterfaces
Definition: ModInterfaces.F90:39

spec_rflu_modchemistry
Definition: SPEC_RFLU_ModChemistry.F90:40

moddatastruct
Definition: ModDataStruct.F90:40

j
j indices j
Definition: Indexing.h:6

spec_rflu_modchemistry::spec_rflu_integratechemsrcterm
subroutine, public spec_rflu_integratechemsrcterm(pRegion, CALLFLAG)
Definition: SPEC_RFLU_ModChemistry.F90:91

modglobal::t_global
Definition: ModGlobal.F90:50

moderror::errorstop
subroutine errorstop(global, errorCode, errorLine, addMessage)
Definition: ModError.F90:483

moderror::deregisterfunction
subroutine deregisterfunction(global)
Definition: ModError.F90:469

merge
*********************************************************************Illinois Open Source License ****University of Illinois NCSA **Open Source License University of Illinois All rights reserved ****Developed free of to any person **obtaining a copy of this software and associated documentation to deal with the Software without including without limitation the rights to merge
Definition: roccomf90.h:20

mixtperf_g_cpr
real(rfreal) function mixtperf_g_cpr(Cp, R)
Definition: MixtPerf_G.F90:39

mixtperf_eo_grtuvw
real(rfreal) function mixtperf_eo_grtuvw(G, R, T, U, V, W)
Definition: MixtPerf_E.F90:70

moddatatypes
Definition: ModDataTypes.F90:40

moderror
Definition: ModError.F90:39

spec_rflu_modchemistry::odesolve
subroutine odesolve(yin, y, press, neqin, dtout, nout, outcome)
Definition: SPEC_RFLU_ModChemistry.F90:416

Line_2::a
RT a() const
Definition: Line_2.h:140

modgrid
Definition: ModGrid.F90:39