@@ -22,7 +22,6 @@ void Mixer_TMATS_body(double* y, const double* u, const MixerStruct* prm)
22
22
23
23
/*--------Define Constants-------*/
24
24
double WOut , FARcOut , TtOut , Aphyout , PtOut ;
25
- double Test ;
26
25
double E1 , E2 , Eout1 , Eout2 , A1calc , A2calc , Aoutcalc , Ps1g , Ps2g , Psoutg ;
27
26
double Rt1 , ht1in , Ps1 , S1in , Ts1 , hs1 , rhos1 , V1 , Ps1g_new , Ps1g_old , E1_old , MN1g , gammat1g , Ts1g ;
28
27
double Rt2 , ht2in , Ps2 , S2in , Ts2 , hs2 , rhos2 , V2 , Ps2g_new , Ps2g_old , E2_old , MN2g , gammat2g , Ts2g ;
@@ -32,7 +31,7 @@ void Mixer_TMATS_body(double* y, const double* u, const MixerStruct* prm)
32
31
double Aphy1 , Aphy2 ;
33
32
double Psp , Pss , Ttp , Tts , Ptp , Pts , Tsp , Tss , rhosp , rhoss , Vp , Vs , Ap , As , hsp , hss , htp , hts ;
34
33
double Ws , Wp , MNg , Sin ;
35
- double FARp , FARs , Rtp , Rts , Rsp , Rss ;
34
+ double FARp , FARs , Rtp , Rts , Rsp ;
36
35
double rhosg , Vg , hsg ;
37
36
int iter , iter1 , iter2 , iter3a , iter3b , maxiter ;
38
37
int interpErr = 0 ;
@@ -255,7 +254,7 @@ void Mixer_TMATS_body(double* y, const double* u, const MixerStruct* prm)
255
254
256
255
if (V1 <= 0 ) {
257
256
V1 = 0.00001 ;
258
- if (iter >= maxiter && * (prm -> IWork + Er9 )== 0 ){
257
+ if (iter1 >= maxiter && * (prm -> IWork + Er9 )== 0 ){
259
258
#ifdef MATLAB_MEX_FILE
260
259
printf ("Warning in %s, Input 1 flow velocity is zero\n" , prm -> BlkNm );
261
260
#endif
@@ -316,7 +315,7 @@ void Mixer_TMATS_body(double* y, const double* u, const MixerStruct* prm)
316
315
PcalcStat (Pt2In , Ps2g , Tt2In , ht2in , FARc2In , Rt2 , & S2in , & Ts2 , & hs2 , & rhos2 , & V2 );
317
316
if (V2 <= 0 ) {
318
317
V2 = 0.00001 ;
319
- if (iter >= maxiter && * (prm -> IWork + Er11 )== 0 ){
318
+ if (iter2 >= maxiter && * (prm -> IWork + Er11 )== 0 ){
320
319
#ifdef MATLAB_MEX_FILE
321
320
printf ("Warning in %s, Input 2 flow velocity is zero\n" , prm -> BlkNm );
322
321
#endif
0 commit comments