remove switches of istiexclg and istiinclg_test

bbb/bbb.v

@@ -220,12 +220,6 @@ rup21          real      /1./        +input #rup21=up2/up1 if isup1up2
 isupgon(ngspmx) integer /ngspmx*0/   +input #user:=1 for par neutral vel. eqn.; index igsp
 isteon      integer  /1/             +input #user:turns on (=1) electron energy eqn.
 istion      integer  /1/             +input #user:turns on (=1) ion enegy eqn.
-istiexclg   integer  /0/             # Ti is only solved for ions (=1, cftiexclg=0.0).
-istiinclg_test integer  /0/          # not often used. for one to turn on Tg
-                                     #       from a converged case without Tg turned on.
-			             #       One can start using =1 while cftiexclg=1,
-				     #       then slowly reducing cftiexclg to 0 and
-				     #       turn this switch back to =0.
 isngon(ngspmx) integer  /6*1/        +input #user:turns on (=1) neutral eqn.; index igsp
 istgon(ngspmx) integer  /6*0/        +input #user:turns on (=1) gas enegy eqn.
 isphion     integer  /0/             +input #user:turns on (=1) potential eqn.

bbb/oderhs.m

@@ -740,9 +740,6 @@ call xerrab('Only ismcnon=0 is validated with openmp.
          ixmp2 = ixpt1(1) + nxomit + 3*(ixpt2(1)-ixpt1(1))/4
       endif
 
-c     set switches for adding Tg from a default Ti case
-      if ((istiexclg .eq. 1) .and. (istiinclg_test.ne.1)) cftiexclg = 0.0
-
 c     Set switches for neutrals-related source terms in plasma equations
 c     (MER 1996/10/28)
 c     (IJ  2015/04/06) add ismcnon>=3 for external call to run_neutrals