fix a minor bug in flux_steger_warming

Author: fangjian19

src/riemann.F90

@@ -51,18 +51,12 @@ subroutine flux_steger_warming(fplus,fmius,rho,vel,prs,tmp,spc,q,dxi,jacob)
       gpd(2)=dxi(i,2)*var0
       gpd(3)=dxi(i,3)*var0
       !
-      if(nondimen) then 
-        gm2=0.5d0/gamma
-        css=sos(tmp(i))
-      else
-        !
 #ifdef COMB
-        gamma = gammarmix(tmp(i),spc(i,:))
-        gm2=0.5d0/gamma
-        call aceval(tmp(i),spc(i,:),css)
+      gamma = gammarmix(tmp(i),spc(i,:))
 #endif
-        !
-      endif
+      gm2=0.5d0/gamma
+      css=sos(tmp(i),spc(i,:))
+      !
       csa=css/var0
       lmach=uu/csa
       !