@@ -185,7 +185,7 @@ function ConstantTVGasDomain(;phase::E2,initialconds::Dict{X,X2},constantspecies
185
185
jacobian= zeros (typeof (T),length (phase. species),length (phase. species))
186
186
end
187
187
rxnarray = getreactionindices (phase)
188
- return ConstantTPDomain (phase,[1 ,length (phase. species),length (phase. species)+ 1 ],[1 ,length (phase. species)+ length (phase. reactions)],constspcinds,
188
+ return ConstantTVGasDomain (phase,[1 ,length (phase. species),length (phase. species)+ 1 ],[1 ,length (phase. species)+ length (phase. reactions)],constspcinds,
189
189
T,V,kfs,krevs,efficiencyinds,Gs,rxnarray,mu,diffs,jacobian,sensitivity,false ,MVector (false ),MVector (0.0 ),p, Dict (" P" => length (phase. species)+ 1 )), y0, p
190
190
end
191
191
export ConstantTVGasDomain
@@ -1180,7 +1180,7 @@ end
1180
1180
for ind in d. efficiencyinds # efficiency related rates may have changed
1181
1181
d. kfs[ind],d. krevs[ind] = getkfkrev (d. phase. reactions[ind],d. phase,d. T,P,C,N,ns,d. Gs,d. diffusivity,d. V,0.0 ;f= kfps[ind])
1182
1182
end
1183
- return ns,cs,d. T,d . P, V,C,N,d. mu,d. kfs,d. krevs,Array {Float64,1} (),Array {Float64,1} (),Array {Float64,1} (),Array {Float64,1} (),0.0 ,Array {Float64,1} (),0.0
1183
+ return ns,cs,d. T,P,d . V,C,N,d. mu,d. kfs,d. krevs,Array {Float64,1} (),Array {Float64,1} (),Array {Float64,1} (),Array {Float64,1} (),0.0 ,Array {Float64,1} (),0.0
1184
1184
elseif nothermochg
1185
1185
d. kfs = kfps
1186
1186
for ind in d. efficiencyinds # efficiency related rates may have changed
@@ -3832,4 +3832,4 @@ end
3832
3832
else
3833
3833
return eval (d)(phase= IdealDiluteSolution (sensspcs, sensrxns, domain. phase. solvent, name= " phase" ), initialconds= initialconds)[1 ], sensspcnames, senstooriginspcind, senstooriginrxnind
3834
3834
end
3835
- end
3835
+ end
0 commit comments