text:
c --------------------------------------------------------------------
      real*8 function BB_GGV(s,t1,ms1,mg,mt,ms,ms2,mx1,mx2,
     &                       s2t,c2t,c2b,susy)

      implicit none

      real*8 Pi,Nc,Co,Ck,Cqed,Nf,Ns,alphas,fac,susy
      real*8 s,t1,u1,ms1,mg,mt,ms,ms2,s2t,c2t,t,u,ms12,mx1,mx2
      real*8 m2(2,8),softin(4),c2b
      real*8 SK3B0A(7),SK3B0B(8),SK3B0C(3),SK3B0D(3,2),SK3BP(4)
      real*8 SK3C0A(10),SK3C0B(6),SK3C0C(8,2),SK3D0(11,2)
      real*8 SOF1(8)
      real*8 mur,dec,TB_GGB,STGGV

      Pi=4.D0*atan(1.D0)
      Nc = 3.D0
      Co = Nc*(Nc**2-1.D0)
      Ck = (Nc**2-1.D0)/Nc
      Cqed = Nc**2 - 1/Nc**2
      Nf = 6.D0
      Ns = 6.D0
      alphas = 1.D0

      u1 = - s - t1
      t  = t1 + ms1**2
      u  = u1 + ms1**2
      ms12 = ms1**2
      fac = 1.D0/(64*Pi*(Nc**2-1.D0)**2)

      m2(1,1) = s
      m2(1,2) = t1 + ms1**2
      m2(1,3) = u1 + ms1**2
      m2(1,4) = t1
      m2(1,5) = u1
      m2(2,1) = ms1**2
      m2(2,2) = ms2**2
      m2(2,3) = ms**2
      m2(2,4) = mg**2
      m2(2,5) = mt**2
      m2(2,6) = ms1**2
      m2(2,7) = mx1**2
      m2(2,8) = mx2**2

      call SCALAR_ARRAY_TB_B_GS(m2,SK3B0A,SK3B0B,SK3B0C,SK3B0D,SK3BP)
      call SCALAR_ARRAY_TB_C_GS(m2,SK3C0A,SK3C0B,SK3C0C,SK3D0)

      softin(1) = s
      softin(2) = t1
      softin(3) = u1
      softin(4) = ms1**2

      call SOFT_ARRAY_LQ(softin,SOF1)

      mur = ms1

      dec = 2.D0*alphas/Pi * TB_GGB(s,t1,ms1)*
     &      (  (Ns-2.D0)/12.D0 * susy * log(mur**2/ms**2)
     &       + 1.D0/24.D0             * log(mur**2/ms1**2)
     &       + 1.D0/24.D0      * susy * log(mur**2/ms2**2)
     &       + 1.D0/24.D0      * susy * log(mur**2/mx1**2)
     &       + 1.D0/24.D0      * susy * log(mur**2/mx2**2)
     &       + 1.D0/6.D0              * log(mur**2/mt**2)
     &       + Nc/6.D0         * susy * log(mur**2/mg**2)   )

      STGGV =
     &  + Cqed*Pi*alphas**3 * ( 4*ms1**2*s**(-2)*t**(-1) + 4*ms1**2*
     &    s**(-2)*u**(-1) - 12*ms1**2*s**(-2)*t1**(-1) - 12*ms1**2*
     &    s**(-2)*u1**(-1) + 4*ms1**2*s**(-1)*t**(-1)*t1**(-1) + 4*
     &    ms1**2*s**(-1)*t**(-1)*u1**(-1) + 4*ms1**2*s**(-1)*u**(-1)*
     &    t1**(-1) + 4*ms1**2*s**(-1)*u**(-1)*u1**(-1) - 4*ms1**2*
     &    s**(-1)*t1**(-2) - 4*ms1**2*s**(-1)*u1**(-2) + 40*ms1**4*
     &    s**(-2)*t**(-1)*t1**(-1) + 20*ms1**4*s**(-2)*t**(-1)*u1**(-1)
     &     + 20*ms1**4*s**(-2)*u**(-1)*t1**(-1) + 40*ms1**4*s**(-2)*
     &    u**(-1)*u1**(-1) - 40*ms1**4*s**(-2)*t1**(-2) - 48*ms1**4*
     &    s**(-2)*t1**(-1)*u1**(-1) - 40*ms1**4*s**(-2)*u1**(-2) + 4*
     &    ms1**4*s**(-1)*t**(-1)*t1**(-2) + 4*ms1**4*s**(-1)*u**(-1)*
     &    u1**(-2) + 72*ms1**6*s**(-2)*t**(-1)*t1**(-2) + 56*ms1**6*
     &    s**(-2)*t**(-1)*t1**(-1)*u1**(-1) + 56*ms1**6*s**(-2)*u**(-1)
     &    *t1**(-1)*u1**(-1) + 72*ms1**6*s**(-2)*u**(-1)*u1**(-2) - 32*
     &    ms1**6*s**(-2)*t1**(-3) - 32*ms1**6*s**(-2)*t1**(-2)*u1**(-1)
     &     - 32*ms1**6*s**(-2)*t1**(-1)*u1**(-2) )
      STGGV = STGGV + Cqed*Pi*alphas**3 * (  - 32*ms1**6*s**(-2)*
     &    u1**(-3) + 32*ms1**8*s**(-2)*t**(-1)*t1**(-3) + 32*ms1**8*
     &    s**(-2)*t**(-1)*t1**(-2)*u1**(-1) + 32*ms1**8*s**(-2)*u**(-1)
     &    *t1**(-1)*u1**(-2) + 32*ms1**8*s**(-2)*u**(-1)*u1**(-3) )
      STGGV = STGGV + Cqed*s2t**2*Pi*alphas**3 * ( 8*ms1**2*s**(-2)*
     &    t1**(-1) + 8*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + Ck*Pi*alphas**3 * ( 16*ms1**2*s**(-2)*t1**(-1) +
     &    16*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + Ck*s2t**2*Pi*alphas**3 * (  - 16*ms1**2*s**(-2)*
     &    t1**(-1) - 16*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + Co*Pi*alphas**3 * ( 16./3.*ms1**2*s**(-3) - 8./3.
     &    *ms1**2*s**(-2)*t1**(-1) - 8./3.*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + Co*s2t**2*Pi*alphas**3 * ( 4*ms1**2*s**(-2)*
     &    t1**(-1) + 4*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + Co*(Nf-1)*Pi*alphas**3 * ( 32./3.*ms1**2*s**(-3)
     &     + 8./3.*ms1**2*s**(-2)*t1**(-1) + 8./3.*ms1**2*s**(-2)*
     &    u1**(-1) )
      STGGV = STGGV + Nc*Ck*Pi*alphas**3 * (  - 4*ms1**2*s**(-3)*
     &    t**(-1)*t1 + 4*ms1**2*s**(-3)*u**(-1)*t1 + 8*ms1**2*s**(-3)
     &     + 6*ms1**2*s**(-2)*t**(-1) + 10*ms1**2*s**(-2)*u**(-1) - 30*
     &    ms1**2*s**(-2)*t1**(-1) - 30*ms1**2*s**(-2)*u1**(-1) - 12*
     &    ms1**2*s**(-1)*t**(-1)*t1**(-1) - 16*ms1**2*s**(-1)*t**(-1)*
     &    u1**(-1) - 16*ms1**2*s**(-1)*u**(-1)*t1**(-1) - 12*ms1**2*
     &    s**(-1)*u**(-1)*u1**(-1) + 12*ms1**2*s**(-1)*t1**(-2) + 12*
     &    ms1**2*s**(-1)*u1**(-2) - 20*ms1**4*s**(-3)*t**(-1) - 20*
     &    ms1**4*s**(-3)*u**(-1) + 16*ms1**4*s**(-3)*t1**(-1) + 16*
     &    ms1**4*s**(-3)*u1**(-1) - 18*ms1**4*s**(-2)*t**(-1)*t1**(-1)
     &     + 24*ms1**4*s**(-2)*t**(-1)*u1**(-1) + 24*ms1**4*s**(-2)*
     &    u**(-1)*t1**(-1) - 18*ms1**4*s**(-2)*u**(-1)*u1**(-1) + 40*
     &    ms1**4*s**(-2)*t1**(-2) - 16*ms1**4*s**(-2)*t1**(-1)*u1**(-1)
     &     + 40*ms1**4*s**(-2)*u1**(-2) - 20*ms1**4*s**(-1)*t**(-1)*
     &    t1**(-2) - 20*ms1**4*s**(-1)*u**(-1)*u1**(-2) + 8*ms1**4*
     &    s**(-1)*t1**(-3) )
      STGGV = STGGV + Nc*Ck*Pi*alphas**3 * ( 8*ms1**4*s**(-1)*u1**(-3)
     &     - 16*ms1**6*s**(-3)*t**(-1)*t1**(-1) - 16*ms1**6*s**(-3)*
     &    u**(-1)*u1**(-1) - 104*ms1**6*s**(-2)*t**(-1)*t1**(-2) - 24*
     &    ms1**6*s**(-2)*t**(-1)*t1**(-1)*u1**(-1) - 24*ms1**6*s**(-2)*
     &    u**(-1)*t1**(-1)*u1**(-1) - 104*ms1**6*s**(-2)*u**(-1)*
     &    u1**(-2) + 64*ms1**6*s**(-2)*t1**(-3) + 32*ms1**6*s**(-2)*
     &    t1**(-2)*u1**(-1) + 32*ms1**6*s**(-2)*t1**(-1)*u1**(-2) + 64*
     &    ms1**6*s**(-2)*u1**(-3) - 8*ms1**6*s**(-1)*t**(-1)*t1**(-3)
     &     - 8*ms1**6*s**(-1)*u**(-1)*u1**(-3) - 64*ms1**8*s**(-2)*
     &    t**(-1)*t1**(-3) - 32*ms1**8*s**(-2)*t**(-1)*t1**(-2)*
     &    u1**(-1) - 32*ms1**8*s**(-2)*u**(-1)*t1**(-1)*u1**(-2) - 64*
     &    ms1**8*s**(-2)*u**(-1)*u1**(-3) )
      STGGV = STGGV + Nc*Ck*s2t**2*Pi*alphas**3 * (  - 4*ms1**2*s**(-2)
     &    *t1**(-1) - 4*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + Nc*Co*Pi*alphas**3 * ( ms1**2*s**(-3)*t**(-1)*t1
     &     - ms1**2*s**(-3)*u**(-1)*t1 - 38./3.*ms1**2*s**(-3) + 3*
     &    ms1**2*s**(-2)*t**(-1) + 2*ms1**2*s**(-2)*u**(-1) - 41./3.*
     &    ms1**2*s**(-2)*t1**(-1) - 41./3.*ms1**2*s**(-2)*u1**(-1) + 5*
     &    ms1**2*s**(-1)*t**(-1)*t1**(-1) + 5*ms1**2*s**(-1)*u**(-1)*
     &    u1**(-1) - 5*ms1**2*s**(-1)*t1**(-2) - 5*ms1**2*s**(-1)*
     &    u1**(-2) + 5*ms1**4*s**(-3)*t**(-1) + 5*ms1**4*s**(-3)*
     &    u**(-1) - 4*ms1**4*s**(-3)*t1**(-1) - 4*ms1**4*s**(-3)*
     &    u1**(-1) + 3*ms1**4*s**(-2)*t**(-1)*t1**(-1) + 3*ms1**4*
     &    s**(-2)*u**(-1)*u1**(-1) + 11*ms1**4*s**(-1)*t**(-1)*t1**(-2)
     &     + 11*ms1**4*s**(-1)*u**(-1)*u1**(-2) - 6*ms1**4*s**(-1)*
     &    t1**(-3) - 6*ms1**4*s**(-1)*u1**(-3) + 4*ms1**6*s**(-3)*
     &    t**(-1)*t1**(-1) + 4*ms1**6*s**(-3)*u**(-1)*u1**(-1) + 6*
     &    ms1**6*s**(-1)*t**(-1)*t1**(-3) + 6*ms1**6*s**(-1)*u**(-1)*
     &    u1**(-3) )
      STGGV = STGGV + susy*Cqed*Pi*alphas**3 * ( 16*ms1**2*s**(-2)*
     &    t1**(-1) + 16*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + susy*Cqed*s2t**2*Pi*alphas**3 * (  - 8*ms1**2*
     &    s**(-2)*t1**(-1) - 8*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + susy*Ck*Pi*alphas**3 * (  - 16*ms1**2*s**(-2)*
     &    t1**(-1) - 16*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + susy*Ck*s2t**2*Pi*alphas**3 * ( 16*ms1**2*s**(-2)
     &    *t1**(-1) + 16*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + susy*Co*Pi*alphas**3 * (  - 16*ms1**2*s**(-3) )
      STGGV = STGGV + susy*Co*s2t**2*Pi*alphas**3 * (  - 4*ms1**2*
     &    s**(-2)*t1**(-1) - 4*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + susy*Co*(Ns-2)*Pi*alphas**3 * (  - 32./3.*ms1**2*
     &    s**(-3) - 8./3.*ms1**2*s**(-2)*t1**(-1) - 8./3.*ms1**2*
     &    s**(-2)*u1**(-1) )
      STGGV = STGGV + susy*Nc*Ck*Pi*alphas**3 * (  - 8*ms1**2*s**(-2)*
     &    t1**(-1) - 8*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + susy*Nc*Ck*s2t**2*Pi*alphas**3 * ( 4*ms1**2*
     &    s**(-2)*t1**(-1) + 4*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + susy*Nc*Co*Pi*alphas**3 * (  - 8*(s-4*ms12)**(-1)
     &    *mg**2*ms1**(-2)*s**(-1) + 4*(s-4*ms12)**(-1)*mg**2*ms1**(-2)
     &    *t1**(-1) + 4*(s-4*ms12)**(-1)*mg**2*ms1**(-2)*u1**(-1) + 32*
     &    (s-4*ms12)**(-1)*mg**2*s**(-2) - 16*(s-4*ms12)**(-1)*mg**2*
     &    s**(-1)*t1**(-1) - 16*(s-4*ms12)**(-1)*mg**2*s**(-1)*u1**(-1)
     &     + 32./3.*ms1**2*s**(-3) + 32./3.*ms1**2*s**(-2)*t1**(-1) +
     &    32./3.*ms1**2*s**(-2)*u1**(-1) + 8*mg**2*ms1**(-2)*s**(-2) -
     &    4*mg**2*ms1**(-2)*s**(-1)*t1**(-1) - 4*mg**2*ms1**(-2)*
     &    s**(-1)*u1**(-1) - 16*mg**2*s**(-3)*t**(-1)*t1 + 16*mg**2*
     &    s**(-3)*u**(-1)*t1 + 32*mg**2*s**(-3) - 40*mg**2*s**(-2)*
     &    t**(-1) - 24*mg**2*s**(-2)*u**(-1) + 40*mg**2*s**(-2)*
     &    t1**(-1) + 40*mg**2*s**(-2)*u1**(-1) - 16*mg**2*s**(-1)*
     &    t**(-1)*t1**(-1) - 16*mg**2*s**(-1)*u**(-1)*u1**(-1) + 16*
     &    mg**2*s**(-1)*t1**(-2) + 16*mg**2*s**(-1)*u1**(-2) - 32*mg**2
     &    *ms1**2*s**(-3)*t**(-1) - 32*mg**2*ms1**2*s**(-3)*u**(-1) +
     &    16*mg**2*ms1**2*s**(-3)*t1**(-1) )
      STGGV = STGGV + susy*Nc*Co*Pi*alphas**3 * ( 16*mg**2*ms1**2*
     &    s**(-3)*u1**(-1) - 72*mg**2*ms1**2*s**(-2)*t**(-1)*t1**(-1)
     &     - 72*mg**2*ms1**2*s**(-2)*u**(-1)*u1**(-1) + 32*mg**2*ms1**2
     &    *s**(-2)*t1**(-2) + 32*mg**2*ms1**2*s**(-2)*u1**(-2) - 32*
     &    mg**2*ms1**2*s**(-1)*t**(-1)*t1**(-2) - 32*mg**2*ms1**2*
     &    s**(-1)*u**(-1)*u1**(-2) + 16*mg**2*ms1**2*s**(-1)*t1**(-3)
     &     + 16*mg**2*ms1**2*s**(-1)*u1**(-3) - 16*mg**2*ms1**4*s**(-3)
     &    *t**(-1)*t1**(-1) - 16*mg**2*ms1**4*s**(-3)*u**(-1)*u1**(-1)
     &     - 32*mg**2*ms1**4*s**(-2)*t**(-1)*t1**(-2) - 32*mg**2*ms1**4
     &    *s**(-2)*u**(-1)*u1**(-2) - 16*mg**2*ms1**4*s**(-1)*t**(-1)*
     &    t1**(-3) - 16*mg**2*ms1**4*s**(-1)*u**(-1)*u1**(-3) )
      STGGV = STGGV + SK3B0A(1)*susy*Nc*Ck*Pi*alphas**3 * ( 32./3.*
     &    s**(-2) + 64./3.*ms1**2*s**(-2)*t1**(-1) + 64./3.*ms1**2*
     &    s**(-2)*u1**(-1) + 64./3.*ms1**4*s**(-2)*t1**(-2) + 128./3.*
     &    ms1**4*s**(-2)*t1**(-1)*u1**(-1) + 64./3.*ms1**4*s**(-2)*
     &    u1**(-2) )
      STGGV = STGGV + SK3B0A(1)*susy*Nc*Co*Pi*alphas**3 * (  - 8*
     &    (s-4*ms12)**(-1)*mg**2*ms1**(-2)*s**(-1) + 4*(s-4*ms12)**(-1)
     &    *mg**2*ms1**(-2)*t1**(-1) + 4*(s-4*ms12)**(-1)*mg**2*
     &    ms1**(-2)*u1**(-1) + 32*(s-4*ms12)**(-1)*mg**2*s**(-2) - 16*
     &    (s-4*ms12)**(-1)*mg**2*s**(-1)*t1**(-1) - 16*(s-4*ms12)**(-1)
     &    *mg**2*s**(-1)*u1**(-1) - 64./3.*s**(-4)*t1**2 - 64./3.*
     &    s**(-3)*t1 - 32./3.*s**(-2) - 128./3.*ms1**2*s**(-3) - 64./3.
     &    *ms1**2*s**(-2)*t1**(-1) - 64./3.*ms1**2*s**(-2)*u1**(-1) -
     &    64./3.*ms1**4*s**(-2)*t1**(-2) - 64./3.*ms1**4*s**(-2)*
     &    u1**(-2) + 8*mg**2*ms1**(-2)*s**(-2) - 4*mg**2*ms1**(-2)*
     &    s**(-1)*t1**(-1) - 4*mg**2*ms1**(-2)*s**(-1)*u1**(-1) - 16*
     &    mg**2*s**(-3)*t**(-1)*t1 + 16*mg**2*s**(-3)*u**(-1)*t1 + 32*
     &    mg**2*s**(-3) - 40*mg**2*s**(-2)*t**(-1) - 24*mg**2*s**(-2)*
     &    u**(-1) + 40*mg**2*s**(-2)*t1**(-1) + 40*mg**2*s**(-2)*
     &    u1**(-1) - 16*mg**2*s**(-1)*t**(-1)*t1**(-1) - 16*mg**2*
     &    s**(-1)*u**(-1)*u1**(-1) )
      STGGV = STGGV + SK3B0A(1)*susy*Nc*Co*Pi*alphas**3 * ( 16*mg**2*
     &    s**(-1)*t1**(-2) + 16*mg**2*s**(-1)*u1**(-2) - 128*mg**2*
     &    ms1**2*s**(-4) - 32*mg**2*ms1**2*s**(-3)*t**(-1) - 32*mg**2*
     &    ms1**2*s**(-3)*u**(-1) - 80./3.*mg**2*ms1**2*s**(-3)*t1**(-1)
     &     - 80./3.*mg**2*ms1**2*s**(-3)*u1**(-1) - 72*mg**2*ms1**2*
     &    s**(-2)*t**(-1)*t1**(-1) - 72*mg**2*ms1**2*s**(-2)*u**(-1)*
     &    u1**(-1) + 32*mg**2*ms1**2*s**(-2)*t1**(-2) - 32./3.*mg**2*
     &    ms1**2*s**(-2)*t1**(-1)*u1**(-1) + 32*mg**2*ms1**2*s**(-2)*
     &    u1**(-2) - 32*mg**2*ms1**2*s**(-1)*t**(-1)*t1**(-2) - 32*
     &    mg**2*ms1**2*s**(-1)*u**(-1)*u1**(-2) + 16*mg**2*ms1**2*
     &    s**(-1)*t1**(-3) + 16*mg**2*ms1**2*s**(-1)*u1**(-3) - 16*
     &    mg**2*ms1**4*s**(-3)*t**(-1)*t1**(-1) - 16*mg**2*ms1**4*
     &    s**(-3)*u**(-1)*u1**(-1) - 32*mg**2*ms1**4*s**(-2)*t**(-1)*
     &    t1**(-2) - 32*mg**2*ms1**4*s**(-2)*u**(-1)*u1**(-2) - 16*
     &    mg**2*ms1**4*s**(-1)*t**(-1)*t1**(-3) - 16*mg**2*ms1**4*
     &    s**(-1)*u**(-1)*u1**(-3) )
      STGGV = STGGV + SK3B0A(2)*Ck*Pi*alphas**3 * ( 32./3.*s**(-2) + 64.
     &    /3.*ms1**2*s**(-2)*t1**(-1) + 64./3.*ms1**2*s**(-2)*u1**(-1)
     &     + 64./3.*ms1**4*s**(-2)*t1**(-2) + 128./3.*ms1**4*s**(-2)*
     &    t1**(-1)*u1**(-1) + 64./3.*ms1**4*s**(-2)*u1**(-2) )
      STGGV = STGGV + SK3B0A(2)*Co*Pi*alphas**3 * (  - 64./3.*s**(-4)*
     &    t1**2 - 64./3.*s**(-3)*t1 - 32./3.*s**(-2) - 128./3.*ms1**2*
     &    s**(-3) - 64./3.*ms1**2*s**(-2)*t1**(-1) - 64./3.*ms1**2*
     &    s**(-2)*u1**(-1) - 64./3.*ms1**4*s**(-2)*t1**(-2) - 64./3.*
     &    ms1**4*s**(-2)*u1**(-2) - 128*mt**2*ms1**2*s**(-4) - 128./3.*
     &    mt**2*ms1**2*s**(-3)*t1**(-1) - 128./3.*mt**2*ms1**2*s**(-3)*
     &    u1**(-1) - 32./3.*mt**2*ms1**2*s**(-2)*t1**(-1)*u1**(-1) )
      STGGV = STGGV + SK3B0A(4)*susy*Ck*(Ns-2)*Pi*alphas**3 * ( 16./3.*
     &    s**(-2) + 32./3.*ms1**2*s**(-2)*t1**(-1) + 32./3.*ms1**2*
     &    s**(-2)*u1**(-1) + 32./3.*ms1**4*s**(-2)*t1**(-2) + 64./3.*
     &    ms1**4*s**(-2)*t1**(-1)*u1**(-1) + 32./3.*ms1**4*s**(-2)*
     &    u1**(-2) )
      STGGV = STGGV + SK3B0A(4)*susy*Co*(Ns-2)*Pi*alphas**3 * (  - 32./
     &    3.*s**(-4)*t1**2 - 32./3.*s**(-3)*t1 - 16./3.*s**(-2) - 64./3.
     &    *ms1**2*s**(-3) - 32./3.*ms1**2*s**(-2)*t1**(-1) - 32./3.*
     &    ms1**2*s**(-2)*u1**(-1) - 32./3.*ms1**4*s**(-2)*t1**(-2) - 32.
     &    /3.*ms1**4*s**(-2)*u1**(-2) + 128*ms**2*ms1**2*s**(-4) + 128./
     &    3.*ms**2*ms1**2*s**(-3)*t1**(-1) + 128./3.*ms**2*ms1**2*
     &    s**(-3)*u1**(-1) + 32./3.*ms**2*ms1**2*s**(-2)*t1**(-1)*
     &    u1**(-1) )
      STGGV = STGGV + SK3B0A(5)*susy*Ck*Pi*alphas**3 * ( 8./3.*s**(-2)
     &     + 16./3.*ms1**2*s**(-2)*t1**(-1) + 16./3.*ms1**2*s**(-2)*
     &    u1**(-1) + 16./3.*ms1**4*s**(-2)*t1**(-2) + 32./3.*ms1**4*
     &    s**(-2)*t1**(-1)*u1**(-1) + 16./3.*ms1**4*s**(-2)*u1**(-2) )
      STGGV = STGGV + SK3B0A(5)*susy*Co*Pi*alphas**3 * (  - 16./3.*
     &    s**(-4)*t1**2 - 16./3.*s**(-3)*t1 - 8./3.*s**(-2) - 32./3.*
     &    ms1**2*s**(-3) - 16./3.*ms1**2*s**(-2)*t1**(-1) - 16./3.*
     &    ms1**2*s**(-2)*u1**(-1) - 16./3.*ms1**4*s**(-2)*t1**(-2) - 16.
     &    /3.*ms1**4*s**(-2)*u1**(-2) + 64*ms2**2*ms1**2*s**(-4) + 64./
     &    3.*ms2**2*ms1**2*s**(-3)*t1**(-1) + 64./3.*ms2**2*ms1**2*
     &    s**(-3)*u1**(-1) + 16./3.*ms2**2*ms1**2*s**(-2)*t1**(-1)*
     &    u1**(-1) )
      STGGV = STGGV + SK3B0A(6)*susy*Ck*Pi*alphas**3 * ( 8./3.*s**(-2)
     &     + 16./3.*ms1**2*s**(-2)*t1**(-1) + 16./3.*ms1**2*s**(-2)*
     &    u1**(-1) + 16./3.*ms1**4*s**(-2)*t1**(-2) + 32./3.*ms1**4*
     &    s**(-2)*t1**(-1)*u1**(-1) + 16./3.*ms1**4*s**(-2)*u1**(-2) )
      STGGV = STGGV + SK3B0A(6)*susy*Co*Pi*alphas**3 * (  - 16./3.*
     &    s**(-4)*t1**2 - 16./3.*s**(-3)*t1 - 8./3.*s**(-2) - 32./3.*
     &    ms1**2*s**(-3) - 16./3.*ms1**2*s**(-2)*t1**(-1) - 16./3.*
     &    ms1**2*s**(-2)*u1**(-1) + 64*ms1**2*mx1**2*s**(-4) + 64./3.*
     &    ms1**2*mx1**2*s**(-3)*t1**(-1) + 64./3.*ms1**2*mx1**2*s**(-3)
     &    *u1**(-1) + 16./3.*ms1**2*mx1**2*s**(-2)*t1**(-1)*u1**(-1) -
     &    16./3.*ms1**4*s**(-2)*t1**(-2) - 16./3.*ms1**4*s**(-2)*
     &    u1**(-2) )
      STGGV = STGGV + SK3B0A(7)*susy*Ck*Pi*alphas**3 * ( 8./3.*s**(-2)
     &     + 16./3.*ms1**2*s**(-2)*t1**(-1) + 16./3.*ms1**2*s**(-2)*
     &    u1**(-1) + 16./3.*ms1**4*s**(-2)*t1**(-2) + 32./3.*ms1**4*
     &    s**(-2)*t1**(-1)*u1**(-1) + 16./3.*ms1**4*s**(-2)*u1**(-2) )
      STGGV = STGGV + SK3B0A(7)*susy*Co*Pi*alphas**3 * (  - 16./3.*
     &    s**(-4)*t1**2 - 16./3.*s**(-3)*t1 - 8./3.*s**(-2) - 32./3.*
     &    ms1**2*s**(-3) - 16./3.*ms1**2*s**(-2)*t1**(-1) - 16./3.*
     &    ms1**2*s**(-2)*u1**(-1) - 16./3.*ms1**4*s**(-2)*t1**(-2) - 16.
     &    /3.*ms1**4*s**(-2)*u1**(-2) + 64*mx2**2*ms1**2*s**(-4) + 64./
     &    3.*mx2**2*ms1**2*s**(-3)*t1**(-1) + 64./3.*mx2**2*ms1**2*
     &    s**(-3)*u1**(-1) + 16./3.*mx2**2*ms1**2*s**(-2)*t1**(-1)*
     &    u1**(-1) )
      STGGV = STGGV + SK3B0B(1)*Nc*Co*Pi*alphas**3 * (  - 32*
     &    (s-4*ms12)**(-1)*s**(-3)*t1**2 - 32*(s-4*ms12)**(-1)*s**(-2)*
     &    t1 - 16*(s-4*ms12)**(-1)*s**(-1) - 2*(s-4*ms12)**(-1)*
     &    t1**(-1) - 2*(s-4*ms12)**(-1)*u1**(-1) - 32*s**(-4)*t1**2 -
     &    32*s**(-3)*t1 + 2*s**(-1)*t1**(-1) + 2*s**(-1)*u1**(-1) - 32*
     &    ms1**2*s**(-3) - 8*ms1**2*s**(-2)*t1**(-1) - 8*ms1**2*s**(-2)
     &    *u1**(-1) )
      STGGV = STGGV + SK3B0B(1)*susy*Nc*Ck*Pi*alphas**3 * (  - 16*
     &    (s-4*ms12)**(-1)*s**(-3)*t1**2 - 16*(s-4*ms12)**(-1)*s**(-2)*
     &    t1 - 8*(s-4*ms12)**(-1)*s**(-1) - (s-4*ms12)**(-1)*t1**(-1)
     &     - (s-4*ms12)**(-1)*u1**(-1) - 64*(s-4*ms12)**(-1)*mg**2*
     &    s**(-4)*t1**2 - 64*(s-4*ms12)**(-1)*mg**2*s**(-3)*t1 - 32*
     &    (s-4*ms12)**(-1)*mg**2*s**(-2) - 4*(s-4*ms12)**(-1)*mg**2*
     &    s**(-1)*t1**(-1) - 4*(s-4*ms12)**(-1)*mg**2*s**(-1)*u1**(-1)
     &     - 16*s**(-4)*t1**2 - 16*s**(-3)*t1 + s**(-1)*t1**(-1) +
     &    s**(-1)*u1**(-1) - 16*ms1**2*s**(-3) - 4*ms1**2*s**(-2)*
     &    t1**(-1) - 4*ms1**2*s**(-2)*u1**(-1) + 16*mg**2*s**(-3) + 4*
     &    mg**2*s**(-2)*t1**(-1) + 4*mg**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3B0B(2)*Co*Pi*alphas**3 * (  - 64*ms1**4*
     &    s**(-4) - 64./3.*ms1**4*s**(-3)*t1**(-1) - 64./3.*ms1**4*
     &    s**(-3)*u1**(-1) - 16./3.*ms1**4*s**(-2)*t1**(-1)*u1**(-1) )
      STGGV = STGGV + SK3B0B(2)*Nc*Ck*Pi*alphas**3 * ( 32*
     &    (s-4*ms12)**(-1)*s**(-3)*t1**2 + 32*(s-4*ms12)**(-1)*s**(-2)*
     &    t1 + 16*(s-4*ms12)**(-1)*s**(-1) + 2*(s-4*ms12)**(-1)*
     &    t1**(-1) + 2*(s-4*ms12)**(-1)*u1**(-1) + 32*s**(-4)*t1**2 +
     &    32*s**(-3)*t1 - 2*s**(-1)*t1**(-1) - 2*s**(-1)*u1**(-1) + 32*
     &    ms1**2*s**(-3) + 8*ms1**2*s**(-2)*t1**(-1) + 8*ms1**2*s**(-2)
     &    *u1**(-1) )
      STGGV = STGGV + SK3B0B(3)*susy*Nc*Co*Pi*alphas**3 * ( 16*
     &    (s-4*ms12)**(-1)*s**(-3)*t1**2 + 16*(s-4*ms12)**(-1)*s**(-2)*
     &    t1 + 8*(s-4*ms12)**(-1)*s**(-1) + (s-4*ms12)**(-1)*t1**(-1)
     &     + (s-4*ms12)**(-1)*u1**(-1) + 64*(s-4*ms12)**(-1)*mg**2*
     &    s**(-4)*t1**2 + 64*(s-4*ms12)**(-1)*mg**2*s**(-3)*t1 + 32*
     &    (s-4*ms12)**(-1)*mg**2*s**(-2) + 4*(s-4*ms12)**(-1)*mg**2*
     &    s**(-1)*t1**(-1) + 4*(s-4*ms12)**(-1)*mg**2*s**(-1)*u1**(-1)
     &     + 16*s**(-4)*t1**2 + 16*s**(-3)*t1 - s**(-1)*t1**(-1) -
     &    s**(-1)*u1**(-1) + 16*ms1**2*s**(-3) + 4*ms1**2*s**(-2)*
     &    t1**(-1) + 4*ms1**2*s**(-2)*u1**(-1) - 16*mg**2*s**(-3) - 4*
     &    mg**2*s**(-2)*t1**(-1) - 4*mg**2*s**(-2)*u1**(-1) + 128*mg**2
     &    *ms1**2*s**(-4) + 128./3.*mg**2*ms1**2*s**(-3)*t1**(-1) + 128.
     &    /3.*mg**2*ms1**2*s**(-3)*u1**(-1) + 32./3.*mg**2*ms1**2*
     &    s**(-2)*t1**(-1)*u1**(-1) )
      STGGV = STGGV + SK3B0B(4)*Co*Pi*alphas**3 * ( 128*mt**2*ms1**2*
     &    s**(-4) + 128./3.*mt**2*ms1**2*s**(-3)*t1**(-1) + 128./3.*
     &    mt**2*ms1**2*s**(-3)*u1**(-1) + 32./3.*mt**2*ms1**2*s**(-2)*
     &    t1**(-1)*u1**(-1) )
      STGGV = STGGV + SK3B0B(5)*susy*Co*(Ns-2)*Pi*alphas**3 * (  - 128*
     &    ms**2*ms1**2*s**(-4) - 128./3.*ms**2*ms1**2*s**(-3)*t1**(-1)
     &     - 128./3.*ms**2*ms1**2*s**(-3)*u1**(-1) - 32./3.*ms**2*
     &    ms1**2*s**(-2)*t1**(-1)*u1**(-1) )
      STGGV = STGGV + SK3B0B(6)*susy*Co*Pi*alphas**3 * (  - 64*ms2**2*
     &    ms1**2*s**(-4) - 64./3.*ms2**2*ms1**2*s**(-3)*t1**(-1) - 64./
     &    3.*ms2**2*ms1**2*s**(-3)*u1**(-1) - 16./3.*ms2**2*ms1**2*
     &    s**(-2)*t1**(-1)*u1**(-1) )
      STGGV = STGGV + SK3B0B(7)*susy*Co*Pi*alphas**3 * (  - 64*ms1**2*
     &    mx1**2*s**(-4) - 64./3.*ms1**2*mx1**2*s**(-3)*t1**(-1) - 64./
     &    3.*ms1**2*mx1**2*s**(-3)*u1**(-1) - 16./3.*ms1**2*mx1**2*
     &    s**(-2)*t1**(-1)*u1**(-1) )
      STGGV = STGGV + SK3B0B(8)*susy*Co*Pi*alphas**3 * (  - 64*mx2**2*
     &    ms1**2*s**(-4) - 64./3.*mx2**2*ms1**2*s**(-3)*t1**(-1) - 64./
     &    3.*mx2**2*ms1**2*s**(-3)*u1**(-1) - 16./3.*mx2**2*ms1**2*
     &    s**(-2)*t1**(-1)*u1**(-1) )
      STGGV = STGGV + SK3B0C(1)*Cqed*Pi*alphas**3 * ( 32*s**(-2) + 64*
     &    ms1**2*s**(-2)*t1**(-1) + 64*ms1**2*s**(-2)*u1**(-1) + 96*
     &    ms1**4*s**(-2)*t1**(-2) + 128*ms1**4*s**(-2)*t1**(-1)*
     &    u1**(-1) + 96*ms1**4*s**(-2)*u1**(-2) + 64*ms1**6*s**(-2)*
     &    t1**(-3) + 64*ms1**6*s**(-2)*t1**(-2)*u1**(-1) + 64*ms1**6*
     &    s**(-2)*t1**(-1)*u1**(-2) + 64*ms1**6*s**(-2)*u1**(-3) )
      STGGV = STGGV + SK3B0C(1)*Nc*Ck*Pi*alphas**3 * (  - 32*
     &    (s-4*ms12)**(-1)*s**(-3)*t1**2 - 32*(s-4*ms12)**(-1)*s**(-2)*
     &    t1 - 16*(s-4*ms12)**(-1)*s**(-1) - 2*(s-4*ms12)**(-1)*
     &    t1**(-1) - 2*(s-4*ms12)**(-1)*u1**(-1) - 32*s**(-4)*t1**2 -
     &    32*s**(-3)*t1 - 80*s**(-2) + 2*s**(-1)*t1**(-1) + 2*s**(-1)*
     &    u1**(-1) - 96*ms1**2*s**(-3) - 184*ms1**2*s**(-2)*t1**(-1) -
     &    184*ms1**2*s**(-2)*u1**(-1) - 32*ms1**4*s**(-3)*t1**(-1) - 32
     &    *ms1**4*s**(-3)*u1**(-1) - 288*ms1**4*s**(-2)*t1**(-2) - 256*
     &    ms1**4*s**(-2)*t1**(-1)*u1**(-1) - 288*ms1**4*s**(-2)*
     &    u1**(-2) - 192*ms1**6*s**(-2)*t1**(-3) - 128*ms1**6*s**(-2)*
     &    t1**(-2)*u1**(-1) - 128*ms1**6*s**(-2)*t1**(-1)*u1**(-2) -
     &    192*ms1**6*s**(-2)*u1**(-3) )
      STGGV = STGGV + SK3B0C(1)*Nc*Co*Pi*alphas**3 * ( 32*
     &    (s-4*ms12)**(-1)*s**(-3)*t1**2 + 32*(s-4*ms12)**(-1)*s**(-2)*
     &    t1 + 16*(s-4*ms12)**(-1)*s**(-1) + 2*(s-4*ms12)**(-1)*
     &    t1**(-1) + 2*(s-4*ms12)**(-1)*u1**(-1) + 32*s**(-4)*t1**2 +
     &    32*s**(-3)*t1 + 16*s**(-2) - 2*s**(-1)*t1**(-1) - 2*s**(-1)*
     &    u1**(-1) + 96*ms1**2*s**(-3) + 56*ms1**2*s**(-2)*t1**(-1) +
     &    56*ms1**2*s**(-2)*u1**(-1) + 32*ms1**4*s**(-3)*t1**(-1) + 32*
     &    ms1**4*s**(-3)*u1**(-1) + 96*ms1**4*s**(-2)*t1**(-2) + 96*
     &    ms1**4*s**(-2)*u1**(-2) + 64*ms1**6*s**(-2)*t1**(-3) + 64*
     &    ms1**6*s**(-2)*u1**(-3) )
      STGGV = STGGV + SK3B0C(3)*susy*Cqed*Pi*alphas**3 * (  - 16*
     &    s**(-2) - 32*ms1**2*s**(-2)*t1**(-1) - 32*ms1**2*s**(-2)*
     &    u1**(-1) - 48*ms1**4*s**(-2)*t1**(-2) - 64*ms1**4*s**(-2)*
     &    t1**(-1)*u1**(-1) - 48*ms1**4*s**(-2)*u1**(-2) - 32*ms1**6*
     &    s**(-2)*t1**(-3) - 32*ms1**6*s**(-2)*t1**(-2)*u1**(-1) - 32*
     &    ms1**6*s**(-2)*t1**(-1)*u1**(-2) - 32*ms1**6*s**(-2)*u1**(-3)
     &     + 16*mg**2*s**(-2)*t1**(-1) + 16*mg**2*s**(-2)*u1**(-1) + 48
     &    *mg**2*ms1**2*s**(-2)*t1**(-2) + 64*mg**2*ms1**2*s**(-2)*
     &    t1**(-1)*u1**(-1) + 48*mg**2*ms1**2*s**(-2)*u1**(-2) + 32*
     &    mg**2*ms1**4*s**(-2)*t1**(-3) + 32*mg**2*ms1**4*s**(-2)*
     &    t1**(-2)*u1**(-1) + 32*mg**2*ms1**4*s**(-2)*t1**(-1)*u1**(-2)
     &     + 32*mg**2*ms1**4*s**(-2)*u1**(-3) )
      STGGV = STGGV + SK3B0C(3)*susy*Nc*Ck*Pi*alphas**3 * ( 16*
     &    (s-4*ms12)**(-1)*s**(-3)*t1**2 + 16*(s-4*ms12)**(-1)*s**(-2)*
     &    t1 + 8*(s-4*ms12)**(-1)*s**(-1) + (s-4*ms12)**(-1)*t1**(-1)
     &     + (s-4*ms12)**(-1)*u1**(-1) + 64*(s-4*ms12)**(-1)*mg**2*
     &    s**(-4)*t1**2 + 64*(s-4*ms12)**(-1)*mg**2*s**(-3)*t1 + 32*
     &    (s-4*ms12)**(-1)*mg**2*s**(-2) + 4*(s-4*ms12)**(-1)*mg**2*
     &    s**(-1)*t1**(-1) + 4*(s-4*ms12)**(-1)*mg**2*s**(-1)*u1**(-1)
     &     + 16*s**(-4)*t1**2 + 16*s**(-3)*t1 + 40*s**(-2) - s**(-1)*
     &    t1**(-1) - s**(-1)*u1**(-1) + 48*ms1**2*s**(-3) + 92*ms1**2*
     &    s**(-2)*t1**(-1) + 92*ms1**2*s**(-2)*u1**(-1) + 16*ms1**4*
     &    s**(-3)*t1**(-1) + 16*ms1**4*s**(-3)*u1**(-1) + 144*ms1**4*
     &    s**(-2)*t1**(-2) + 128*ms1**4*s**(-2)*t1**(-1)*u1**(-1) + 144
     &    *ms1**4*s**(-2)*u1**(-2) + 96*ms1**6*s**(-2)*t1**(-3) + 64*
     &    ms1**6*s**(-2)*t1**(-2)*u1**(-1) + 64*ms1**6*s**(-2)*t1**(-1)
     &    *u1**(-2) + 96*ms1**6*s**(-2)*u1**(-3) - 48*mg**2*s**(-3) -
     &    52*mg**2*s**(-2)*t1**(-1) )
      STGGV = STGGV + SK3B0C(3)*susy*Nc*Ck*Pi*alphas**3 * (  - 52*mg**2
     &    *s**(-2)*u1**(-1) - 16*mg**2*ms1**2*s**(-3)*t1**(-1) - 16*
     &    mg**2*ms1**2*s**(-3)*u1**(-1) - 144*mg**2*ms1**2*s**(-2)*
     &    t1**(-2) - 128*mg**2*ms1**2*s**(-2)*t1**(-1)*u1**(-1) - 144*
     &    mg**2*ms1**2*s**(-2)*u1**(-2) - 96*mg**2*ms1**4*s**(-2)*
     &    t1**(-3) - 64*mg**2*ms1**4*s**(-2)*t1**(-2)*u1**(-1) - 64*
     &    mg**2*ms1**4*s**(-2)*t1**(-1)*u1**(-2) - 96*mg**2*ms1**4*
     &    s**(-2)*u1**(-3) )
      STGGV = STGGV + SK3B0C(3)*susy*Nc*Co*Pi*alphas**3 * (  - 16*
     &    (s-4*ms12)**(-1)*s**(-3)*t1**2 - 16*(s-4*ms12)**(-1)*s**(-2)*
     &    t1 - 8*(s-4*ms12)**(-1)*s**(-1) - (s-4*ms12)**(-1)*t1**(-1)
     &     - (s-4*ms12)**(-1)*u1**(-1) + 8*(s-4*ms12)**(-1)*mg**2*
     &    ms1**(-2)*s**(-1) - 4*(s-4*ms12)**(-1)*mg**2*ms1**(-2)*
     &    t1**(-1) - 4*(s-4*ms12)**(-1)*mg**2*ms1**(-2)*u1**(-1) - 64*
     &    (s-4*ms12)**(-1)*mg**2*s**(-4)*t1**2 - 64*(s-4*ms12)**(-1)*
     &    mg**2*s**(-3)*t1 - 64*(s-4*ms12)**(-1)*mg**2*s**(-2) + 12*
     &    (s-4*ms12)**(-1)*mg**2*s**(-1)*t1**(-1) + 12*(s-4*ms12)**(-1)
     &    *mg**2*s**(-1)*u1**(-1) - 16*s**(-4)*t1**2 - 16*s**(-3)*t1 -
     &    8*s**(-2) + s**(-1)*t1**(-1) + s**(-1)*u1**(-1) - 48*ms1**2*
     &    s**(-3) - 28*ms1**2*s**(-2)*t1**(-1) - 28*ms1**2*s**(-2)*
     &    u1**(-1) - 16*ms1**4*s**(-3)*t1**(-1) - 16*ms1**4*s**(-3)*
     &    u1**(-1) - 48*ms1**4*s**(-2)*t1**(-2) - 48*ms1**4*s**(-2)*
     &    u1**(-2) - 32*ms1**6*s**(-2)*t1**(-3) - 32*ms1**6*s**(-2)*
     &    u1**(-3) )
      STGGV = STGGV + SK3B0C(3)*susy*Nc*Co*Pi*alphas**3 * (  - 8*mg**2*
     &    ms1**(-2)*s**(-2) + 4*mg**2*ms1**(-2)*s**(-1)*t1**(-1) + 4*
     &    mg**2*ms1**(-2)*s**(-1)*u1**(-1) + 48*mg**2*s**(-3) + 20*
     &    mg**2*s**(-2)*t1**(-1) + 20*mg**2*s**(-2)*u1**(-1) + 16*mg**2
     &    *ms1**2*s**(-3)*t1**(-1) + 16*mg**2*ms1**2*s**(-3)*u1**(-1)
     &     + 48*mg**2*ms1**2*s**(-2)*t1**(-2) + 48*mg**2*ms1**2*s**(-2)
     &    *u1**(-2) + 32*mg**2*ms1**4*s**(-2)*t1**(-3) + 32*mg**2*
     &    ms1**4*s**(-2)*u1**(-3) )
      STGGV = STGGV + SK3B0D(1,1)*Cqed*Pi*alphas**3 * ( 4*s**(-2)*
     &    t**(-1)*t1 - 20*s**(-2) - 4*s**(-1)*t1**(-1) - 4*s**(-1)*
     &    u1**(-1) - 4*t**(-1)*u1**(-1) - 4*ms1**2*s**(-2)*t**(-1) - 24
     &    *ms1**2*s**(-2)*t1**(-1) + 4*ms1**2*s**(-1)*t**(-1)*t1**(-1)
     &     + 20*ms1**2*s**(-1)*t**(-1)*u1**(-1) - 120*ms1**4*s**(-2)*
     &    t**(-1)*t1**(-1) - 80*ms1**4*s**(-2)*t**(-1)*u1**(-1) - 160*
     &    ms1**6*s**(-2)*t**(-1)*t1**(-2) - 128*ms1**6*s**(-2)*t**(-1)*
     &    t1**(-1)*u1**(-1) - 64*ms1**8*s**(-2)*t**(-1)*t1**(-3) - 64*
     &    ms1**8*s**(-2)*t**(-1)*t1**(-2)*u1**(-1) )
      STGGV = STGGV + SK3B0D(1,1)*Nc*Ck*Pi*alphas**3 * (  - 4*s**(-3)*
     &    t**(-1)*t1**2 + 20*s**(-3)*t1 + 2*s**(-2)*t**(-1)*t1 + 46*
     &    s**(-2) - 4*s**(-1)*t**(-1) + 12*s**(-1)*t1**(-1) + 8*s**(-1)
     &    *u1**(-1) + 8*t**(-1)*u1**(-1) + 4*ms1**2*s**(-3)*t**(-1)*t1
     &     + 24*ms1**2*s**(-3) + 38*ms1**2*s**(-2)*t**(-1) + 84*ms1**2*
     &    s**(-2)*t1**(-1) + 8*ms1**2*s**(-2)*u1**(-1) - 20*ms1**2*
     &    s**(-1)*t**(-1)*t1**(-1) - 32*ms1**2*s**(-1)*t**(-1)*u1**(-1)
     &     + 8*ms1**2*s**(-1)*t1**(-2) + 40*ms1**4*s**(-3)*t**(-1) +
     &    236*ms1**4*s**(-2)*t**(-1)*t1**(-1) + 72*ms1**4*s**(-2)*
     &    t**(-1)*u1**(-1) + 112*ms1**4*s**(-2)*t1**(-2) + 80*ms1**4*
     &    s**(-2)*t1**(-1)*u1**(-1) - 8*ms1**4*s**(-1)*t**(-1)*t1**(-2)
     &     + 32*ms1**6*s**(-3)*t**(-1)*t1**(-1) + 304*ms1**6*s**(-2)*
     &    t**(-1)*t1**(-2) + 112*ms1**6*s**(-2)*t**(-1)*t1**(-1)*
     &    u1**(-1) + 64*ms1**6*s**(-2)*t1**(-3) + 64*ms1**6*s**(-2)*
     &    t1**(-2)*u1**(-1) + 128*ms1**8*s**(-2)*t**(-1)*t1**(-3) + 64*
     &    ms1**8*s**(-2)*t**(-1)*t1**(-2)*u1**(-1) )
      STGGV = STGGV + SK3B0D(1,1)*Nc*Co*Pi*alphas**3 * ( 7*s**(-3)*
     &    t**(-1)*t1**2 - 23*s**(-3)*t1 + 17*s**(-2)*t**(-1)*t1 - 33*
     &    s**(-2) + 5*s**(-1)*t**(-1) - 5*s**(-1)*t1**(-1) + 19*ms1**2*
     &    s**(-3)*t**(-1)*t1 - 44*ms1**2*s**(-3) + 47*ms1**2*s**(-2)*
     &    t**(-1) - 78*ms1**2*s**(-2)*t1**(-1) + 11*ms1**2*s**(-1)*
     &    t**(-1)*t1**(-1) - 6*ms1**2*s**(-1)*t1**(-2) + 12*ms1**4*
     &    s**(-3)*t**(-1) - 32*ms1**4*s**(-3)*t1**(-1) + 46*ms1**4*
     &    s**(-2)*t**(-1)*t1**(-1) - 112*ms1**4*s**(-2)*t1**(-2) + 6*
     &    ms1**4*s**(-1)*t**(-1)*t1**(-2) + 16*ms1**6*s**(-2)*t**(-1)*
     &    t1**(-2) - 64*ms1**6*s**(-2)*t1**(-3) )
      STGGV = STGGV + SK3B0D(1,2)*Cqed*Pi*alphas**3 * (  - 4*s**(-2)*
     &    u**(-1)*t1 - 20*s**(-2) - 4*s**(-1)*u**(-1) - 4*s**(-1)*
     &    t1**(-1) - 4*s**(-1)*u1**(-1) - 4*u**(-1)*t1**(-1) - 4*ms1**2
     &    *s**(-2)*u**(-1) - 24*ms1**2*s**(-2)*u1**(-1) + 20*ms1**2*
     &    s**(-1)*u**(-1)*t1**(-1) + 4*ms1**2*s**(-1)*u**(-1)*u1**(-1)
     &     - 80*ms1**4*s**(-2)*u**(-1)*t1**(-1) - 120*ms1**4*s**(-2)*
     &    u**(-1)*u1**(-1) - 128*ms1**6*s**(-2)*u**(-1)*t1**(-1)*
     &    u1**(-1) - 160*ms1**6*s**(-2)*u**(-1)*u1**(-2) - 64*ms1**8*
     &    s**(-2)*u**(-1)*t1**(-1)*u1**(-2) - 64*ms1**8*s**(-2)*u**(-1)
     &    *u1**(-3) )
      STGGV = STGGV + SK3B0D(1,2)*Nc*Ck*Pi*alphas**3 * (  - 4*s**(-3)*
     &    u**(-1)*t1**2 - 20*s**(-3)*t1 - 10*s**(-2)*u**(-1)*t1 + 26*
     &    s**(-2) - 10*s**(-1)*u**(-1) + 8*s**(-1)*t1**(-1) + 12*
     &    s**(-1)*u1**(-1) + 8*u**(-1)*t1**(-1) - 4*ms1**2*s**(-3)*
     &    u**(-1)*t1 + 24*ms1**2*s**(-3) + 34*ms1**2*s**(-2)*u**(-1) +
     &    8*ms1**2*s**(-2)*t1**(-1) + 84*ms1**2*s**(-2)*u1**(-1) - 32*
     &    ms1**2*s**(-1)*u**(-1)*t1**(-1) - 20*ms1**2*s**(-1)*u**(-1)*
     &    u1**(-1) + 8*ms1**2*s**(-1)*u1**(-2) + 40*ms1**4*s**(-3)*
     &    u**(-1) + 72*ms1**4*s**(-2)*u**(-1)*t1**(-1) + 236*ms1**4*
     &    s**(-2)*u**(-1)*u1**(-1) + 80*ms1**4*s**(-2)*t1**(-1)*
     &    u1**(-1) + 112*ms1**4*s**(-2)*u1**(-2) - 8*ms1**4*s**(-1)*
     &    u**(-1)*u1**(-2) + 32*ms1**6*s**(-3)*u**(-1)*u1**(-1) + 112*
     &    ms1**6*s**(-2)*u**(-1)*t1**(-1)*u1**(-1) + 304*ms1**6*s**(-2)
     &    *u**(-1)*u1**(-2) + 64*ms1**6*s**(-2)*t1**(-1)*u1**(-2) + 64*
     &    ms1**6*s**(-2)*u1**(-3) + 64*ms1**8*s**(-2)*u**(-1)*t1**(-1)*
     &    u1**(-2) )
      STGGV = STGGV + SK3B0D(1,2)*Nc*Ck*Pi*alphas**3 * ( 128*ms1**8*
     &    s**(-2)*u**(-1)*u1**(-3) )
      STGGV = STGGV + SK3B0D(1,2)*Nc*Co*Pi*alphas**3 * ( 7*s**(-3)*
     &    u**(-1)*t1**2 + 23*s**(-3)*t1 - 3*s**(-2)*u**(-1)*t1 - 10*
     &    s**(-2) - 5*s**(-1)*u**(-1) - 5*s**(-1)*u1**(-1) - 19*ms1**2*
     &    s**(-3)*u**(-1)*t1 - 44*ms1**2*s**(-3) + 28*ms1**2*s**(-2)*
     &    u**(-1) - 78*ms1**2*s**(-2)*u1**(-1) + 11*ms1**2*s**(-1)*
     &    u**(-1)*u1**(-1) - 6*ms1**2*s**(-1)*u1**(-2) + 12*ms1**4*
     &    s**(-3)*u**(-1) - 32*ms1**4*s**(-3)*u1**(-1) + 46*ms1**4*
     &    s**(-2)*u**(-1)*u1**(-1) - 112*ms1**4*s**(-2)*u1**(-2) + 6*
     &    ms1**4*s**(-1)*u**(-1)*u1**(-2) + 16*ms1**6*s**(-2)*u**(-1)*
     &    u1**(-2) - 64*ms1**6*s**(-2)*u1**(-3) )
      STGGV = STGGV + SK3B0D(3,1)*susy*Cqed*Pi*alphas**3 * ( 8*s**(-2)
     &     + 24*ms1**2*s**(-2)*t1**(-1) + 8*ms1**2*s**(-2)*u1**(-1) +
     &    48*ms1**4*s**(-2)*t1**(-2) + 32*ms1**4*s**(-2)*t1**(-1)*
     &    u1**(-1) + 32*ms1**6*s**(-2)*t1**(-3) + 32*ms1**6*s**(-2)*
     &    t1**(-2)*u1**(-1) - 16*mg**2*s**(-2)*t1**(-1) - 48*mg**2*
     &    ms1**2*s**(-2)*t1**(-2) - 32*mg**2*ms1**2*s**(-2)*t1**(-1)*
     &    u1**(-1) - 32*mg**2*ms1**4*s**(-2)*t1**(-3) - 32*mg**2*ms1**4
     &    *s**(-2)*t1**(-2)*u1**(-1) )
      STGGV = STGGV + SK3B0D(3,1)*susy*Nc*Ck*Pi*alphas**3 * (  - 8*
     &    s**(-3)*t1 - 24*s**(-2) - 16*ms1**2*s**(-3) - 72*ms1**2*
     &    s**(-2)*t1**(-1) - 16*ms1**2*s**(-2)*u1**(-1) - 16*ms1**4*
     &    s**(-3)*t1**(-1) - 144*ms1**4*s**(-2)*t1**(-2) - 64*ms1**4*
     &    s**(-2)*t1**(-1)*u1**(-1) - 96*ms1**6*s**(-2)*t1**(-3) - 64*
     &    ms1**6*s**(-2)*t1**(-2)*u1**(-1) + 16*mg**2*s**(-3) + 48*
     &    mg**2*s**(-2)*t1**(-1) + 16*mg**2*ms1**2*s**(-3)*t1**(-1) +
     &    144*mg**2*ms1**2*s**(-2)*t1**(-2) + 64*mg**2*ms1**2*s**(-2)*
     &    t1**(-1)*u1**(-1) + 96*mg**2*ms1**4*s**(-2)*t1**(-3) + 64*
     &    mg**2*ms1**4*s**(-2)*t1**(-2)*u1**(-1) )
      STGGV = STGGV + SK3B0D(3,1)*susy*Nc*Co*Pi*alphas**3 * ( 16*
     &    s**(-3)*t**(-1)*t1**2 - 8*s**(-3)*t1 + 40*s**(-2)*t**(-1)*t1
     &     - 32*s**(-2) + 16*s**(-1)*t**(-1) - 16*s**(-1)*t1**(-1) + 48
     &    *ms1**2*s**(-3)*t**(-1)*t1 - 16*ms1**2*s**(-3) + 112*ms1**2*
     &    s**(-2)*t**(-1) - 48*ms1**2*s**(-2)*t1**(-1) + 48*ms1**2*
     &    s**(-1)*t**(-1)*t1**(-1) - 32*ms1**2*s**(-1)*t1**(-2) + 48*
     &    ms1**4*s**(-3)*t**(-1) + 104*ms1**4*s**(-2)*t**(-1)*t1**(-1)
     &     + 16*ms1**4*s**(-2)*t1**(-2) + 48*ms1**4*s**(-1)*t**(-1)*
     &    t1**(-2) - 16*ms1**4*s**(-1)*t1**(-3) + 16*ms1**6*s**(-3)*
     &    t**(-1)*t1**(-1) + 32*ms1**6*s**(-2)*t**(-1)*t1**(-2) + 32*
     &    ms1**6*s**(-2)*t1**(-3) + 16*ms1**6*s**(-1)*t**(-1)*t1**(-3)
     &     + 16*mg**2*s**(-3)*t**(-1)*t1 - 32*mg**2*s**(-3) + 40*mg**2*
     &    s**(-2)*t**(-1) - 56*mg**2*s**(-2)*t1**(-1) + 16*mg**2*
     &    s**(-1)*t**(-1)*t1**(-1) - 16*mg**2*s**(-1)*t1**(-2) + 32*
     &    mg**2*ms1**2*s**(-3)*t**(-1) - 32*mg**2*ms1**2*s**(-3)*
     &    t1**(-1) )
      STGGV = STGGV + SK3B0D(3,1)*susy*Nc*Co*Pi*alphas**3 * ( 72*mg**2*
     &    ms1**2*s**(-2)*t**(-1)*t1**(-1) - 80*mg**2*ms1**2*s**(-2)*
     &    t1**(-2) + 32*mg**2*ms1**2*s**(-1)*t**(-1)*t1**(-2) - 16*
     &    mg**2*ms1**2*s**(-1)*t1**(-3) + 16*mg**2*ms1**4*s**(-3)*
     &    t**(-1)*t1**(-1) + 32*mg**2*ms1**4*s**(-2)*t**(-1)*t1**(-2)
     &     - 32*mg**2*ms1**4*s**(-2)*t1**(-3) + 16*mg**2*ms1**4*s**(-1)
     &    *t**(-1)*t1**(-3) )
      STGGV = STGGV + SK3B0D(3,2)*susy*Cqed*Pi*alphas**3 * ( 8*s**(-2)
     &     + 8*ms1**2*s**(-2)*t1**(-1) + 24*ms1**2*s**(-2)*u1**(-1) +
     &    32*ms1**4*s**(-2)*t1**(-1)*u1**(-1) + 48*ms1**4*s**(-2)*
     &    u1**(-2) + 32*ms1**6*s**(-2)*t1**(-1)*u1**(-2) + 32*ms1**6*
     &    s**(-2)*u1**(-3) - 16*mg**2*s**(-2)*u1**(-1) - 32*mg**2*
     &    ms1**2*s**(-2)*t1**(-1)*u1**(-1) - 48*mg**2*ms1**2*s**(-2)*
     &    u1**(-2) - 32*mg**2*ms1**4*s**(-2)*t1**(-1)*u1**(-2) - 32*
     &    mg**2*ms1**4*s**(-2)*u1**(-3) )
      STGGV = STGGV + SK3B0D(3,2)*susy*Nc*Ck*Pi*alphas**3 * ( 8*s**(-3)
     &    *t1 - 16*s**(-2) - 16*ms1**2*s**(-3) - 16*ms1**2*s**(-2)*
     &    t1**(-1) - 72*ms1**2*s**(-2)*u1**(-1) - 16*ms1**4*s**(-3)*
     &    u1**(-1) - 64*ms1**4*s**(-2)*t1**(-1)*u1**(-1) - 144*ms1**4*
     &    s**(-2)*u1**(-2) - 64*ms1**6*s**(-2)*t1**(-1)*u1**(-2) - 96*
     &    ms1**6*s**(-2)*u1**(-3) + 16*mg**2*s**(-3) + 48*mg**2*s**(-2)
     &    *u1**(-1) + 16*mg**2*ms1**2*s**(-3)*u1**(-1) + 64*mg**2*
     &    ms1**2*s**(-2)*t1**(-1)*u1**(-1) + 144*mg**2*ms1**2*s**(-2)*
     &    u1**(-2) + 64*mg**2*ms1**4*s**(-2)*t1**(-1)*u1**(-2) + 96*
     &    mg**2*ms1**4*s**(-2)*u1**(-3) )
      STGGV = STGGV + SK3B0D(3,2)*susy*Nc*Co*Pi*alphas**3 * ( 16*
     &    s**(-3)*u**(-1)*t1**2 + 8*s**(-3)*t1 - 8*s**(-2)*u**(-1)*t1
     &     - 24*s**(-2) - 8*s**(-1)*u**(-1) - 16*s**(-1)*u1**(-1) - 48*
     &    ms1**2*s**(-3)*u**(-1)*t1 - 16*ms1**2*s**(-3) + 64*ms1**2*
     &    s**(-2)*u**(-1) - 48*ms1**2*s**(-2)*u1**(-1) + 48*ms1**2*
     &    s**(-1)*u**(-1)*u1**(-1) - 32*ms1**2*s**(-1)*u1**(-2) + 48*
     &    ms1**4*s**(-3)*u**(-1) + 104*ms1**4*s**(-2)*u**(-1)*u1**(-1)
     &     + 16*ms1**4*s**(-2)*u1**(-2) + 48*ms1**4*s**(-1)*u**(-1)*
     &    u1**(-2) - 16*ms1**4*s**(-1)*u1**(-3) + 16*ms1**6*s**(-3)*
     &    u**(-1)*u1**(-1) + 32*ms1**6*s**(-2)*u**(-1)*u1**(-2) + 32*
     &    ms1**6*s**(-2)*u1**(-3) + 16*ms1**6*s**(-1)*u**(-1)*u1**(-3)
     &     - 16*mg**2*s**(-3)*u**(-1)*t1 - 32*mg**2*s**(-3) + 24*mg**2*
     &    s**(-2)*u**(-1) - 56*mg**2*s**(-2)*u1**(-1) + 16*mg**2*
     &    s**(-1)*u**(-1)*u1**(-1) - 16*mg**2*s**(-1)*u1**(-2) + 32*
     &    mg**2*ms1**2*s**(-3)*u**(-1) - 32*mg**2*ms1**2*s**(-3)*
     &    u1**(-1) )
      STGGV = STGGV + SK3B0D(3,2)*susy*Nc*Co*Pi*alphas**3 * ( 72*mg**2*
     &    ms1**2*s**(-2)*u**(-1)*u1**(-1) - 80*mg**2*ms1**2*s**(-2)*
     &    u1**(-2) + 32*mg**2*ms1**2*s**(-1)*u**(-1)*u1**(-2) - 16*
     &    mg**2*ms1**2*s**(-1)*u1**(-3) + 16*mg**2*ms1**4*s**(-3)*
     &    u**(-1)*u1**(-1) + 32*mg**2*ms1**4*s**(-2)*u**(-1)*u1**(-2)
     &     - 32*mg**2*ms1**4*s**(-2)*u1**(-3) + 16*mg**2*ms1**4*s**(-1)
     &    *u**(-1)*u1**(-3) )
      STGGV = STGGV + SK3BP(1)*Cqed*Pi*alphas**3 * ( 32*ms1**2*s**(-2)
     &     + 64*ms1**4*s**(-2)*t1**(-1) + 64*ms1**4*s**(-2)*u1**(-1) +
     &    64*ms1**6*s**(-2)*t1**(-2) + 128*ms1**6*s**(-2)*t1**(-1)*
     &    u1**(-1) + 64*ms1**6*s**(-2)*u1**(-2) )
      STGGV = STGGV + SK3BP(1)*Nc*Ck*Pi*alphas**3 * (  - 64*ms1**2*
     &    s**(-4)*t1**2 - 64*ms1**2*s**(-3)*t1 - 96*ms1**2*s**(-2) -
     &    128*ms1**4*s**(-3) - 192*ms1**4*s**(-2)*t1**(-1) - 192*ms1**4
     &    *s**(-2)*u1**(-1) - 192*ms1**6*s**(-2)*t1**(-2) - 256*ms1**6*
     &    s**(-2)*t1**(-1)*u1**(-1) - 192*ms1**6*s**(-2)*u1**(-2) )
      STGGV = STGGV + SK3BP(1)*Nc*Co*Pi*alphas**3 * ( 64*ms1**2*s**(-4)
     &    *t1**2 + 64*ms1**2*s**(-3)*t1 + 32*ms1**2*s**(-2) + 128*
     &    ms1**4*s**(-3) + 64*ms1**4*s**(-2)*t1**(-1) + 64*ms1**4*
     &    s**(-2)*u1**(-1) + 64*ms1**6*s**(-2)*t1**(-2) + 64*ms1**6*
     &    s**(-2)*u1**(-2) )
      STGGV = STGGV + SK3BP(4)*susy*Cqed*Pi*alphas**3 * (  - 16*ms1**2*
     &    s**(-2) - 32*ms1**4*s**(-2)*t1**(-1) - 32*ms1**4*s**(-2)*
     &    u1**(-1) - 32*ms1**6*s**(-2)*t1**(-2) - 64*ms1**6*s**(-2)*
     &    t1**(-1)*u1**(-1) - 32*ms1**6*s**(-2)*u1**(-2) + 16*mg**2*
     &    s**(-2) + 32*mg**2*ms1**2*s**(-2)*t1**(-1) + 32*mg**2*ms1**2*
     &    s**(-2)*u1**(-1) + 32*mg**2*ms1**4*s**(-2)*t1**(-2) + 64*
     &    mg**2*ms1**4*s**(-2)*t1**(-1)*u1**(-1) + 32*mg**2*ms1**4*
     &    s**(-2)*u1**(-2) )
      STGGV = STGGV + SK3BP(4)*susy*Nc*Ck*Pi*alphas**3 * ( 32*ms1**2*
     &    s**(-4)*t1**2 + 32*ms1**2*s**(-3)*t1 + 48*ms1**2*s**(-2) + 64
     &    *ms1**4*s**(-3) + 96*ms1**4*s**(-2)*t1**(-1) + 96*ms1**4*
     &    s**(-2)*u1**(-1) + 96*ms1**6*s**(-2)*t1**(-2) + 128*ms1**6*
     &    s**(-2)*t1**(-1)*u1**(-1) + 96*ms1**6*s**(-2)*u1**(-2) - 32*
     &    mg**2*s**(-4)*t1**2 - 32*mg**2*s**(-3)*t1 - 48*mg**2*s**(-2)
     &     - 64*mg**2*ms1**2*s**(-3) - 96*mg**2*ms1**2*s**(-2)*t1**(-1)
     &     - 96*mg**2*ms1**2*s**(-2)*u1**(-1) - 96*mg**2*ms1**4*s**(-2)
     &    *t1**(-2) - 128*mg**2*ms1**4*s**(-2)*t1**(-1)*u1**(-1) - 96*
     &    mg**2*ms1**4*s**(-2)*u1**(-2) )
      STGGV = STGGV + SK3BP(4)*susy*Nc*Co*Pi*alphas**3 * (  - 32*ms1**2
     &    *s**(-4)*t1**2 - 32*ms1**2*s**(-3)*t1 - 16*ms1**2*s**(-2) -
     &    64*ms1**4*s**(-3) - 32*ms1**4*s**(-2)*t1**(-1) - 32*ms1**4*
     &    s**(-2)*u1**(-1) - 32*ms1**6*s**(-2)*t1**(-2) - 32*ms1**6*
     &    s**(-2)*u1**(-2) + 32*mg**2*s**(-4)*t1**2 + 32*mg**2*s**(-3)*
     &    t1 + 16*mg**2*s**(-2) + 64*mg**2*ms1**2*s**(-3) + 32*mg**2*
     &    ms1**2*s**(-2)*t1**(-1) + 32*mg**2*ms1**2*s**(-2)*u1**(-1) +
     &    32*mg**2*ms1**4*s**(-2)*t1**(-2) + 32*mg**2*ms1**4*s**(-2)*
     &    u1**(-2) )
      STGGV = STGGV + SK3C0A(1)*Nc*Co*Pi*alphas**3 * (  - 32*ms1**2*
     &    s**(-2) - 16*ms1**2*s**(-1)*t1**(-1) - 16*ms1**2*s**(-1)*
     &    u1**(-1) )
      STGGV = STGGV + SK3C0A(1)*susy*Cqed*Pi*alphas**3 * ( 8*s**(-1) +
     &    4*t1**(-1) + 4*u1**(-1) - 8*ms1**2*s**(-1)*t1**(-1) - 8*
     &    ms1**2*s**(-1)*u1**(-1) + 16*ms1**4*s**(-2)*t1**(-1) + 16*
     &    ms1**4*s**(-2)*u1**(-1) + 16*mg**2*s**(-1)*t1**(-1) + 16*
     &    mg**2*s**(-1)*u1**(-1) - 32*mg**2*ms1**2*s**(-2)*t1**(-1) -
     &    32*mg**2*ms1**2*s**(-2)*u1**(-1) + 16*mg**4*s**(-2)*t1**(-1)
     &     + 16*mg**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(1)*susy*Nc*Ck*Pi*alphas**3 * ( 4*s**(-1)
     &     - 16*ms1**2*s**(-2) - 8*ms1**4*s**(-2)*t1**(-1) - 8*ms1**4*
     &    s**(-2)*u1**(-1) + 16*mg**2*s**(-2) - 4*mg**2*s**(-1)*
     &    t1**(-1) - 4*mg**2*s**(-1)*u1**(-1) + 16*mg**2*ms1**2*s**(-2)
     &    *t1**(-1) + 16*mg**2*ms1**2*s**(-2)*u1**(-1) - 8*mg**4*
     &    s**(-2)*t1**(-1) - 8*mg**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(2)*Cqed*Pi*alphas**3 * (  - 8*t1**(-1) - 8
     &    *u1**(-1) + 48*ms1**2*s**(-1)*t1**(-1) + 48*ms1**2*s**(-1)*
     &    u1**(-1) - 96*ms1**4*s**(-2)*t1**(-1) - 96*ms1**4*s**(-2)*
     &    u1**(-1) )
      STGGV = STGGV + SK3C0A(2)*Cqed*s2t**2*Pi*alphas**3 * ( 16*ms1**4*
     &    s**(-2)*t1**(-1) + 16*ms1**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(2)*Ck*Pi*alphas**3 * ( 32*ms1**4*s**(-2)*
     &    t1**(-1) + 32*ms1**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(2)*Ck*s2t**2*Pi*alphas**3 * (  - 32*ms1**4
     &    *s**(-2)*t1**(-1) - 32*ms1**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(2)*Co*Pi*alphas**3 * (  - 32*ms1**4*
     &    s**(-3) - 16*ms1**4*s**(-2)*t1**(-1) - 16*ms1**4*s**(-2)*
     &    u1**(-1) )
      STGGV = STGGV + SK3C0A(2)*Co*s2t**2*Pi*alphas**3 * ( 8*ms1**4*
     &    s**(-2)*t1**(-1) + 8*ms1**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(2)*Nc*Ck*Pi*alphas**3 * (  - 16*s**(-1) +
     &    32*ms1**2*s**(-2) - 16*ms1**2*s**(-1)*t1**(-1) - 16*ms1**2*
     &    s**(-1)*u1**(-1) + 48*ms1**4*s**(-2)*t1**(-1) + 48*ms1**4*
     &    s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(2)*Nc*Ck*s2t**2*Pi*alphas**3 * (  - 8*
     &    ms1**4*s**(-2)*t1**(-1) - 8*ms1**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(3)*susy*Nc*Co*Pi*alphas**3 * (  - 4*
     &    s**(-1) + 16*ms1**2*s**(-2) + 8*ms1**4*s**(-2)*t1**(-1) + 8*
     &    ms1**4*s**(-2)*u1**(-1) - 16*mg**2*s**(-2) - 4*mg**2*s**(-1)*
     &    t1**(-1) - 4*mg**2*s**(-1)*u1**(-1) + 64*mg**2*ms1**2*s**(-3)
     &     + 32*mg**2*ms1**2*s**(-2)*t1**(-1) + 32*mg**2*ms1**2*s**(-2)
     &    *u1**(-1) - 8*mg**4*s**(-2)*t1**(-1) - 8*mg**4*s**(-2)*
     &    u1**(-1) )
      STGGV = STGGV + SK3C0A(4)*Co*Pi*alphas**3 * ( 64*mt**2*ms1**2*
     &    s**(-3) + 16*mt**2*ms1**2*s**(-2)*t1**(-1) + 16*mt**2*ms1**2*
     &    s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(7)*susy*Co*(Ns-2)*Pi*alphas**3 * (  - 64*
     &    ms**2*ms1**2*s**(-3) - 16*ms**2*ms1**2*s**(-2)*t1**(-1) - 16*
     &    ms**2*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(8)*susy*Cqed*s2t**2*Pi*alphas**3 * (  - 16
     &    *ms2**2*ms1**2*s**(-2)*t1**(-1) - 16*ms2**2*ms1**2*s**(-2)*
     &    u1**(-1) )
      STGGV = STGGV + SK3C0A(8)*susy*Ck*Pi*alphas**3 * (  - 32*ms2**2*
     &    ms1**2*s**(-2)*t1**(-1) - 32*ms2**2*ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(8)*susy*Ck*s2t**2*Pi*alphas**3 * ( 32*
     &    ms2**2*ms1**2*s**(-2)*t1**(-1) + 32*ms2**2*ms1**2*s**(-2)*
     &    u1**(-1) )
      STGGV = STGGV + SK3C0A(8)*susy*Co*Pi*alphas**3 * (  - 32*ms2**2*
     &    ms1**2*s**(-3) )
      STGGV = STGGV + SK3C0A(8)*susy*Co*s2t**2*Pi*alphas**3 * (  - 8*
     &    ms2**2*ms1**2*s**(-2)*t1**(-1) - 8*ms2**2*ms1**2*s**(-2)*
     &    u1**(-1) )
      STGGV = STGGV + SK3C0A(8)*susy*Nc*Ck*s2t**2*Pi*alphas**3 * ( 8*
     &    ms2**2*ms1**2*s**(-2)*t1**(-1) + 8*ms2**2*ms1**2*s**(-2)*
     &    u1**(-1) )
      STGGV = STGGV + SK3C0A(9)*susy*Ck*c2t*Pi*alphas**3 * ( 32*ms1**2*
     &    mx1**2*s**(-2)*t1**(-1)*c2b + 32*ms1**2*mx1**2*s**(-2)*
     &    u1**(-1)*c2b )
      STGGV = STGGV + SK3C0A(9)*susy*Co*Pi*alphas**3 * (  - 32*ms1**2*
     &    mx1**2*s**(-3) - 8*ms1**2*mx1**2*s**(-2)*t1**(-1) - 8*ms1**2*
     &    mx1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(9)*susy*Co*c2t*Pi*alphas**3 * (  - 8*
     &    ms1**2*mx1**2*s**(-2)*t1**(-1)*c2b - 8*ms1**2*mx1**2*s**(-2)*
     &    u1**(-1)*c2b )
      STGGV = STGGV + SK3C0A(10)*susy*Ck*c2t*Pi*alphas**3 * (  - 32*
     &    mx2**2*ms1**2*s**(-2)*t1**(-1)*c2b - 32*mx2**2*ms1**2*s**(-2)
     &    *u1**(-1)*c2b )
      STGGV = STGGV + SK3C0A(10)*susy*Co*Pi*alphas**3 * (  - 32*mx2**2*
     &    ms1**2*s**(-3) - 8*mx2**2*ms1**2*s**(-2)*t1**(-1) - 8*mx2**2*
     &    ms1**2*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0A(10)*susy*Co*c2t*Pi*alphas**3 * ( 8*mx2**2*
     &    ms1**2*s**(-2)*t1**(-1)*c2b + 8*mx2**2*ms1**2*s**(-2)*
     &    u1**(-1)*c2b )
      STGGV = STGGV + SK3C0B(1)*Cqed*Pi*alphas**3 * (  - 16*s**(-1) - 8
     &    *t1**(-1) - 8*u1**(-1) + 32*ms1**2*s**(-2) + 32*ms1**2*
     &    s**(-1)*t1**(-1) + 32*ms1**2*s**(-1)*u1**(-1) - 32*ms1**4*
     &    s**(-2)*t1**(-1) - 32*ms1**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0B(1)*Nc*Ck*Pi*alphas**3 * ( 32*s**(-3)*t1**2
     &     + 32*s**(-2)*t1 - 64*ms1**2*s**(-4)*t1**2 - 64*ms1**2*
     &    s**(-3)*t1 + 96*ms1**2*s**(-2) + 16*ms1**2*s**(-1)*t1**(-1)
     &     + 16*ms1**2*s**(-1)*u1**(-1) - 192*ms1**4*s**(-3) - 32*
     &    ms1**4*s**(-2)*t1**(-1) - 32*ms1**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0B(2)*Nc*Co*Pi*alphas**3 * (  - 16*
     &    (s-4*ms12)**(-1)*s**(-2)*t1**2 - 16*(s-4*ms12)**(-1)*s**(-1)*
     &    t1 + (s-4*ms12)**(-1)*t1**(-1)*u1 - 6*(s-4*ms12)**(-1) +
     &    (s-4*ms12)**(-1)*t1*u1**(-1) + 48*s**(-3)*t1**2 + 48*s**(-2)*
     &    t1 + 24*s**(-1) + t1**(-1) + u1**(-1) + 32*ms1**2*s**(-2) +
     &    20*ms1**2*s**(-1)*t1**(-1) + 20*ms1**2*s**(-1)*u1**(-1) - 16*
     &    ms1**4*s**(-2)*t1**(-1) - 16*ms1**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0B(5)*susy*Nc*Co*Pi*alphas**3 * ( 4*
     &    (s-4*ms12)**(-1)*s**(-2)*t1**2 + 4*(s-4*ms12)**(-1)*s**(-1)*
     &    t1 - 1./4.*(s-4*ms12)**(-1)*t1**(-1)*u1 + 3./2.*
     &    (s-4*ms12)**(-1) - 1./4.*(s-4*ms12)**(-1)*t1*u1**(-1) - 64*
     &    (s-4*ms12)**(-1)*mg**4*s**(-4)*t1**2 - 64*(s-4*ms12)**(-1)*
     &    mg**4*s**(-3)*t1 - 32*(s-4*ms12)**(-1)*mg**4*s**(-2) - 4*
     &    (s-4*ms12)**(-1)*mg**4*s**(-1)*t1**(-1) - 4*(s-4*ms12)**(-1)*
     &    mg**4*s**(-1)*u1**(-1) - 4*s**(-3)*t1**2 - 4*s**(-2)*t1 - 6*
     &    s**(-1) - 1./4.*t1**(-1) - 1./4.*u1**(-1) - 16*ms1**2*s**(-4)
     &    *t1**2 - 16*ms1**2*s**(-3)*t1 + 16*ms1**2*s**(-2) - ms1**2*
     &    s**(-1)*t1**(-1) - ms1**2*s**(-1)*u1**(-1) - 48*ms1**4*
     &    s**(-3) - 4*ms1**4*s**(-2)*t1**(-1) - 4*ms1**4*s**(-2)*
     &    u1**(-1) + 8*mg**2*s**(-2) + 4*mg**2*s**(-1)*t1**(-1) + 4*
     &    mg**2*s**(-1)*u1**(-1) - 8*mg**2*ms1**2*s**(-2)*t1**(-1) - 8*
     &    mg**2*ms1**2*s**(-2)*u1**(-1) + 48*mg**4*s**(-3) + 12*mg**4*
     &    s**(-2)*t1**(-1) )
      STGGV = STGGV + SK3C0B(5)*susy*Nc*Co*Pi*alphas**3 * ( 12*mg**4*
     &    s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0B(6)*susy*Cqed*Pi*alphas**3 * ( 8*s**(-1) +
     &    4*t1**(-1) + 4*u1**(-1) - 16*ms1**2*s**(-2) - 16*ms1**2*
     &    s**(-1)*t1**(-1) - 16*ms1**2*s**(-1)*u1**(-1) + 16*ms1**4*
     &    s**(-2)*t1**(-1) + 16*ms1**4*s**(-2)*u1**(-1) + 16*mg**2*
     &    s**(-2) + 8*mg**2*s**(-1)*t1**(-1) + 8*mg**2*s**(-1)*u1**(-1)
     &     - 16*mg**2*ms1**2*s**(-2)*t1**(-1) - 16*mg**2*ms1**2*s**(-2)
     &    *u1**(-1) )
      STGGV = STGGV + SK3C0B(6)*susy*Nc*Ck*Pi*alphas**3 * (  - 4*
     &    (s-4*ms12)**(-1)*s**(-2)*t1**2 - 4*(s-4*ms12)**(-1)*s**(-1)*
     &    t1 + 1./4.*(s-4*ms12)**(-1)*t1**(-1)*u1 - 3./2.*
     &    (s-4*ms12)**(-1) + 1./4.*(s-4*ms12)**(-1)*t1*u1**(-1) - 32*
     &    (s-4*ms12)**(-1)*mg**2*s**(-3)*t1**2 - 32*(s-4*ms12)**(-1)*
     &    mg**2*s**(-2)*t1 - 16*(s-4*ms12)**(-1)*mg**2*s**(-1) - 2*
     &    (s-4*ms12)**(-1)*mg**2*t1**(-1) - 2*(s-4*ms12)**(-1)*mg**2*
     &    u1**(-1) - 64*(s-4*ms12)**(-1)*mg**4*s**(-4)*t1**2 - 64*
     &    (s-4*ms12)**(-1)*mg**4*s**(-3)*t1 - 32*(s-4*ms12)**(-1)*mg**4
     &    *s**(-2) - 4*(s-4*ms12)**(-1)*mg**4*s**(-1)*t1**(-1) - 4*
     &    (s-4*ms12)**(-1)*mg**4*s**(-1)*u1**(-1) + 4*s**(-3)*t1**2 + 4
     &    *s**(-2)*t1 + 6*s**(-1) + 1./4.*t1**(-1) + 1./4.*u1**(-1) +
     &    16*ms1**2*s**(-4)*t1**2 + 16*ms1**2*s**(-3)*t1 - 16*ms1**2*
     &    s**(-2) + ms1**2*s**(-1)*t1**(-1) + ms1**2*s**(-1)*u1**(-1)
     &     + 48*ms1**4*s**(-3) + 4*ms1**4*s**(-2)*t1**(-1) + 4*ms1**4*
     &    s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0B(6)*susy*Nc*Ck*Pi*alphas**3 * (  - 32*mg**2
     &    *s**(-4)*t1**2 - 32*mg**2*s**(-3)*t1 + 24*mg**2*s**(-2) + 6*
     &    mg**2*s**(-1)*t1**(-1) + 6*mg**2*s**(-1)*u1**(-1) - 96*mg**2*
     &    ms1**2*s**(-3) - 16*mg**2*ms1**2*s**(-2)*t1**(-1) - 16*mg**2*
     &    ms1**2*s**(-2)*u1**(-1) + 48*mg**4*s**(-3) + 12*mg**4*s**(-2)
     &    *t1**(-1) + 12*mg**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0C(1,1)*Nc*Ck*Pi*alphas**3 * ( 16*s**(-2)*t1
     &     - 32*ms1**2*s**(-1)*u1**(-1) + 64*ms1**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0C(1,1)*Nc*Co*Pi*alphas**3 * (  - 32*s**(-3)*
     &    t1**2 - 16*s**(-2)*t1 - 32*ms1**2*s**(-3)*t1 - 32*ms1**2*
     &    s**(-2) )
      STGGV = STGGV + SK3C0C(1,2)*Nc*Ck*Pi*alphas**3 * (  - 16*s**(-2)*
     &    t1 - 16*s**(-1) - 32*ms1**2*s**(-1)*t1**(-1) + 64*ms1**4*
     &    s**(-2)*t1**(-1) )
      STGGV = STGGV + SK3C0C(1,2)*Nc*Co*Pi*alphas**3 * (  - 32*s**(-3)*
     &    t1**2 - 48*s**(-2)*t1 - 16*s**(-1) + 32*ms1**2*s**(-3)*t1 )
      STGGV = STGGV + SK3C0C(3,1)*Cqed*Pi*alphas**3 * ( 16*s**(-1) + 16
     &    *u1**(-1) - 32*ms1**2*s**(-2) - 64*ms1**2*s**(-1)*u1**(-1) +
     &    64*ms1**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0C(3,1)*Nc*Ck*Pi*alphas**3 * (  - 32*s**(-2)*
     &    t1 + 32*ms1**2*s**(-3)*t1 - 32*ms1**2*s**(-2) + 32*ms1**2*
     &    s**(-1)*u1**(-1) + 64*ms1**4*s**(-3) )
      STGGV = STGGV + SK3C0C(3,2)*Cqed*Pi*alphas**3 * ( 16*s**(-1) + 16
     &    *t1**(-1) - 32*ms1**2*s**(-2) - 64*ms1**2*s**(-1)*t1**(-1) +
     &    64*ms1**4*s**(-2)*t1**(-1) )
      STGGV = STGGV + SK3C0C(3,2)*Nc*Ck*Pi*alphas**3 * ( 32*s**(-2)*t1
     &     + 32*s**(-1) - 32*ms1**2*s**(-3)*t1 - 64*ms1**2*s**(-2) + 32
     &    *ms1**2*s**(-1)*t1**(-1) + 64*ms1**4*s**(-3) )
      STGGV = STGGV + SK3C0C(7,1)*susy*Nc*Ck*Pi*alphas**3 * (  - 16*
     &    ms1**4*s**(-2)*u1**(-1) - 8*mg**2*s**(-1)*u1**(-1) + 32*mg**2
     &    *ms1**2*s**(-2)*u1**(-1) - 16*mg**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0C(7,1)*susy*Nc*Co*Pi*alphas**3 * ( 8*s**(-3)
     &    *t1**2 + 16*ms1**2*s**(-3)*t1 + 16*ms1**4*s**(-3) - 16*mg**2*
     &    s**(-3)*t1 - 8*mg**2*s**(-2) - 16*mg**4*s**(-3) )
      STGGV = STGGV + SK3C0C(7,2)*susy*Nc*Ck*Pi*alphas**3 * (  - 16*
     &    ms1**4*s**(-2)*t1**(-1) - 8*mg**2*s**(-1)*t1**(-1) + 32*mg**2
     &    *ms1**2*s**(-2)*t1**(-1) - 16*mg**4*s**(-2)*t1**(-1) )
      STGGV = STGGV + SK3C0C(7,2)*susy*Nc*Co*Pi*alphas**3 * ( 8*s**(-3)
     &    *t1**2 + 16*s**(-2)*t1 + 8*s**(-1) - 16*ms1**2*s**(-3)*t1 -
     &    16*ms1**2*s**(-2) + 16*ms1**4*s**(-3) + 16*mg**2*s**(-3)*t1
     &     + 8*mg**2*s**(-2) - 16*mg**4*s**(-3) )
      STGGV = STGGV + SK3C0C(8,1)*susy*Cqed*Pi*alphas**3 * ( 8*s**(-2)*
     &    t1 - 8*s**(-1) - 8*u1**(-1) + 16*ms1**2*s**(-2) + 16*ms1**2*
     &    s**(-1)*u1**(-1) - 16*ms1**4*s**(-2)*u1**(-1) - 16*mg**2*
     &    s**(-2) - 24*mg**2*s**(-1)*u1**(-1) + 32*mg**2*ms1**2*s**(-2)
     &    *u1**(-1) - 16*mg**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0C(8,1)*susy*Nc*Ck*Pi*alphas**3 * (  - 8*
     &    s**(-3)*t1**2 - 16*ms1**2*s**(-3)*t1 - 16*ms1**4*s**(-3) - 16
     &    *ms1**4*s**(-2)*u1**(-1) + 16*mg**2*s**(-3)*t1 - 8*mg**2*
     &    s**(-2) + 8*mg**2*s**(-1)*u1**(-1) + 32*mg**2*ms1**2*s**(-3)
     &     - 16*mg**4*s**(-3) + 16*mg**4*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3C0C(8,2)*susy*Cqed*Pi*alphas**3 * (  - 8*
     &    s**(-2)*t1 - 16*s**(-1) - 8*t1**(-1) + 16*ms1**2*s**(-2) + 16
     &    *ms1**2*s**(-1)*t1**(-1) - 16*ms1**4*s**(-2)*t1**(-1) - 16*
     &    mg**2*s**(-2) - 24*mg**2*s**(-1)*t1**(-1) + 32*mg**2*ms1**2*
     &    s**(-2)*t1**(-1) - 16*mg**4*s**(-2)*t1**(-1) )
      STGGV = STGGV + SK3C0C(8,2)*susy*Nc*Ck*Pi*alphas**3 * (  - 8*
     &    s**(-3)*t1**2 - 16*s**(-2)*t1 - 8*s**(-1) + 16*ms1**2*s**(-3)
     &    *t1 + 16*ms1**2*s**(-2) - 16*ms1**4*s**(-3) - 16*ms1**4*
     &    s**(-2)*t1**(-1) - 16*mg**2*s**(-3)*t1 - 24*mg**2*s**(-2) + 8
     &    *mg**2*s**(-1)*t1**(-1) + 32*mg**2*ms1**2*s**(-3) - 16*mg**4*
     &    s**(-3) + 16*mg**4*s**(-2)*t1**(-1) )
      STGGV = STGGV + SK3D0(1,1)*Nc*Co*Pi*alphas**3 * (  - 16*s**(-3)*
     &    t1**3 - 16*s**(-2)*t1**2 - 8*s**(-1)*t1 - 16*ms1**2*s**(-2)*
     &    t1 - 16*ms1**2*s**(-1) - 32*ms1**4*s**(-1)*t1**(-1) )
      STGGV = STGGV + SK3D0(1,2)*Nc*Co*Pi*alphas**3 * ( 8 + 16*s**(-3)*
     &    t1**3 + 32*s**(-2)*t1**2 + 24*s**(-1)*t1 + 16*ms1**2*s**(-2)*
     &    t1 - 32*ms1**4*s**(-1)*u1**(-1) )
      STGGV = STGGV + SK3D0(2,1)*Cqed*Pi*alphas**3 * ( 8*t1*u1**(-1) +
     &    32*ms1**2*s**(-1) + 48*ms1**2*u1**(-1) - 32*ms1**4*s**(-2) -
     &    32*ms1**4*s**(-1)*t1**(-1) - 96*ms1**4*s**(-1)*u1**(-1) + 64*
     &    ms1**6*s**(-2)*t1**(-1) + 64*ms1**6*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3D0(2,1)*Nc*Ck*Pi*alphas**3 * ( 8*s**(-1)*t1 -
     &    32*ms1**2*s**(-2)*t1 + 16*ms1**2*s**(-1) + 32*ms1**4*s**(-3)*
     &    t1 - 32*ms1**4*s**(-2) + 32*ms1**4*s**(-1)*t1**(-1) - 64*
     &    ms1**6*s**(-2)*t1**(-1) )
      STGGV = STGGV + SK3D0(2,2)*Cqed*Pi*alphas**3 * ( 8*t1**(-1)*u1 +
     &    32*ms1**2*s**(-1) + 48*ms1**2*t1**(-1) - 32*ms1**4*s**(-2) -
     &    96*ms1**4*s**(-1)*t1**(-1) - 32*ms1**4*s**(-1)*u1**(-1) + 64*
     &    ms1**6*s**(-2)*t1**(-1) + 64*ms1**6*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3D0(2,2)*Nc*Ck*Pi*alphas**3 * (  - 8 - 8*
     &    s**(-1)*t1 + 32*ms1**2*s**(-2)*t1 + 48*ms1**2*s**(-1) - 32*
     &    ms1**4*s**(-3)*t1 - 64*ms1**4*s**(-2) + 32*ms1**4*s**(-1)*
     &    u1**(-1) - 64*ms1**6*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3D0(3,1)*Nc*Ck*Pi*alphas**3 * (  - 8*s**(-2)*
     &    t1**2 - 8*s**(-1)*t1 - 16*ms1**2*s**(-1) - 32*ms1**4*s**(-2)
     &     - 32*ms1**4*s**(-1)*t1**(-1) - 32*ms1**4*s**(-1)*u1**(-1) )
      STGGV = STGGV + SK3D0(3,2)*Nc*Ck*Pi*alphas**3 * (  - 8*s**(-2)*
     &    t1**2 - 8*s**(-1)*t1 - 16*ms1**2*s**(-1) - 32*ms1**4*s**(-2)
     &     - 32*ms1**4*s**(-1)*t1**(-1) - 32*ms1**4*s**(-1)*u1**(-1) )
      STGGV = STGGV + SK3D0(9,1)*susy*Nc*Co*Pi*alphas**3 * (  - 4*
     &    s**(-2)*t1**2 - 12*ms1**2*s**(-2)*t1 - 16*ms1**4*s**(-2) - 8*
     &    ms1**6*s**(-2)*t1**(-1) - 4*mg**2*s**(-2)*t1 - 4*mg**2*
     &    s**(-1) + 16*mg**2*ms1**2*s**(-3)*t1 + 16*mg**2*ms1**2*
     &    s**(-2) - 4*mg**2*ms1**2*s**(-1)*t1**(-1) - 8*mg**2*ms1**4*
     &    s**(-2)*t1**(-1) - 16*mg**4*s**(-3)*t1 + 4*mg**4*s**(-1)*
     &    t1**(-1) + 8*mg**4*ms1**2*s**(-2)*t1**(-1) + 8*mg**6*s**(-2)*
     &    t1**(-1) )
      STGGV = STGGV + SK3D0(9,2)*susy*Nc*Co*Pi*alphas**3 * (  - 4 - 4*
     &    s**(-2)*t1**2 - 8*s**(-1)*t1 + 12*ms1**2*s**(-2)*t1 + 12*
     &    ms1**2*s**(-1) - 16*ms1**4*s**(-2) - 8*ms1**6*s**(-2)*
     &    u1**(-1) + 4*mg**2*s**(-2)*t1 - 16*mg**2*ms1**2*s**(-3)*t1 -
     &    4*mg**2*ms1**2*s**(-1)*u1**(-1) - 8*mg**2*ms1**4*s**(-2)*
     &    u1**(-1) + 16*mg**4*s**(-3)*t1 + 16*mg**4*s**(-2) + 4*mg**4*
     &    s**(-1)*u1**(-1) + 8*mg**4*ms1**2*s**(-2)*u1**(-1) + 8*mg**6*
     &    s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3D0(10,1)*susy*Cqed*Pi*alphas**3 * (  - 4*
     &    s**(-1)*t1 - 4*t1*u1**(-1) - 12*ms1**2*s**(-1) - 12*ms1**2*
     &    u1**(-1) + 16*ms1**4*s**(-1)*u1**(-1) - 8*ms1**6*s**(-2)*
     &    t1**(-1) - 8*ms1**6*s**(-2)*u1**(-1) + 12*mg**2*s**(-1) + 16*
     &    mg**2*u1**(-1) - 4*mg**2*ms1**2*s**(-1)*t1**(-1) - 36*mg**2*
     &    ms1**2*s**(-1)*u1**(-1) + 24*mg**2*ms1**4*s**(-2)*t1**(-1) +
     &    24*mg**2*ms1**4*s**(-2)*u1**(-1) + 4*mg**4*s**(-1)*t1**(-1)
     &     + 20*mg**4*s**(-1)*u1**(-1) - 24*mg**4*ms1**2*s**(-2)*
     &    t1**(-1) - 24*mg**4*ms1**2*s**(-2)*u1**(-1) + 8*mg**6*s**(-2)
     &    *t1**(-1) + 8*mg**6*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3D0(10,1)*susy*Nc*Ck*Pi*alphas**3 * ( 4*s**(-2)
     &    *t1**2 + 12*ms1**2*s**(-2)*t1 + 16*ms1**4*s**(-2) + 8*ms1**6*
     &    s**(-2)*t1**(-1) - 12*mg**2*s**(-2)*t1 + 4*mg**2*s**(-1) - 32
     &    *mg**2*ms1**2*s**(-2) + 4*mg**2*ms1**2*s**(-1)*t1**(-1) - 24*
     &    mg**2*ms1**4*s**(-2)*t1**(-1) + 16*mg**4*s**(-2) - 4*mg**4*
     &    s**(-1)*t1**(-1) + 24*mg**4*ms1**2*s**(-2)*t1**(-1) - 8*mg**6
     &    *s**(-2)*t1**(-1) )
      STGGV = STGGV + SK3D0(10,2)*susy*Cqed*Pi*alphas**3 * ( 4 + 4*
     &    s**(-1)*t1 - 4*t1**(-1)*u1 - 12*ms1**2*s**(-1) - 12*ms1**2*
     &    t1**(-1) + 16*ms1**4*s**(-1)*t1**(-1) - 8*ms1**6*s**(-2)*
     &    t1**(-1) - 8*ms1**6*s**(-2)*u1**(-1) + 12*mg**2*s**(-1) + 16*
     &    mg**2*t1**(-1) - 36*mg**2*ms1**2*s**(-1)*t1**(-1) - 4*mg**2*
     &    ms1**2*s**(-1)*u1**(-1) + 24*mg**2*ms1**4*s**(-2)*t1**(-1) +
     &    24*mg**2*ms1**4*s**(-2)*u1**(-1) + 20*mg**4*s**(-1)*t1**(-1)
     &     + 4*mg**4*s**(-1)*u1**(-1) - 24*mg**4*ms1**2*s**(-2)*
     &    t1**(-1) - 24*mg**4*ms1**2*s**(-2)*u1**(-1) + 8*mg**6*s**(-2)
     &    *t1**(-1) + 8*mg**6*s**(-2)*u1**(-1) )
      STGGV = STGGV + SK3D0(10,2)*susy*Nc*Ck*Pi*alphas**3 * ( 4 + 4*
     &    s**(-2)*t1**2 + 8*s**(-1)*t1 - 12*ms1**2*s**(-2)*t1 - 12*
     &    ms1**2*s**(-1) + 16*ms1**4*s**(-2) + 8*ms1**6*s**(-2)*
     &    u1**(-1) + 12*mg**2*s**(-2)*t1 + 16*mg**2*s**(-1) - 32*mg**2*
     &    ms1**2*s**(-2) + 4*mg**2*ms1**2*s**(-1)*u1**(-1) - 24*mg**2*
     &    ms1**4*s**(-2)*u1**(-1) + 16*mg**4*s**(-2) - 4*mg**4*s**(-1)*
     &    u1**(-1) + 24*mg**4*ms1**2*s**(-2)*u1**(-1) - 8*mg**6*s**(-2)
     &    *u1**(-1) )
      STGGV = STGGV + SK3D0(11,2)*susy*Nc*Ck*Pi*alphas**3 * ( 16*ms1**4
     &    *s**(-2) + 16*ms1**6*s**(-2)*t1**(-1) + 16*ms1**6*s**(-2)*
     &    u1**(-1) - 8*mg**2*s**(-1) - 8*mg**2*ms1**2*s**(-1)*t1**(-1)
     &     - 8*mg**2*ms1**2*s**(-1)*u1**(-1) - 16*mg**2*ms1**4*s**(-2)*
     &    t1**(-1) - 16*mg**2*ms1**4*s**(-2)*u1**(-1) - 16*mg**4*
     &    s**(-2) + 8*mg**4*s**(-1)*t1**(-1) + 8*mg**4*s**(-1)*u1**(-1)
     &     - 16*mg**4*ms1**2*s**(-2)*t1**(-1) - 16*mg**4*ms1**2*s**(-2)
     &    *u1**(-1) + 16*mg**6*s**(-2)*t1**(-1) + 16*mg**6*s**(-2)*
     &    u1**(-1) )
      STGGV = STGGV + SOF1(1)*Cqed*Pi*alphas**3 * ( 64*s**(-1) - 128*
     &    ms1**2*s**(-2) + 128*ms1**2*s**(-1)*t1**(-1) + 128*ms1**2*
     &    s**(-1)*u1**(-1) - 512*ms1**4*s**(-2)*t1**(-1) - 512*ms1**4*
     &    s**(-2)*u1**(-1) + 128*ms1**4*s**(-1)*t1**(-2) + 128*ms1**4*
     &    s**(-1)*u1**(-2) - 256*ms1**6*s**(-2)*t1**(-2) - 512*ms1**6*
     &    s**(-2)*t1**(-1)*u1**(-1) - 256*ms1**6*s**(-2)*u1**(-2) )
      STGGV = STGGV + SOF1(1)*Nc*Ck*Pi*alphas**3 * (  - 128*s**(-3)*
     &    t1**2 - 128*s**(-2)*t1 - 64*s**(-1) + 256*ms1**2*s**(-4)*
     &    t1**2 + 256*ms1**2*s**(-3)*t1 - 128*ms1**2*s**(-2) - 128*
     &    ms1**2*s**(-1)*t1**(-1) - 128*ms1**2*s**(-1)*u1**(-1) + 512*
     &    ms1**4*s**(-3) + 256*ms1**4*s**(-2)*t1**(-1) + 256*ms1**4*
     &    s**(-2)*u1**(-1) - 128*ms1**4*s**(-1)*t1**(-2) - 128*ms1**4*
     &    s**(-1)*u1**(-2) + 256*ms1**6*s**(-2)*t1**(-2) + 256*ms1**6*
     &    s**(-2)*u1**(-2) )
      STGGV = STGGV + SOF1(2)*Cqed*Pi*alphas**3 * (  - 64*ms1**2*
     &    s**(-2) - 128*ms1**4*s**(-2)*t1**(-1) - 128*ms1**4*s**(-2)*
     &    u1**(-1) - 128*ms1**6*s**(-2)*t1**(-2) - 256*ms1**6*s**(-2)*
     &    t1**(-1)*u1**(-1) - 128*ms1**6*s**(-2)*u1**(-2) )
      STGGV = STGGV + SOF1(2)*Nc*Ck*Pi*alphas**3 * ( 128*ms1**2*s**(-4)
     &    *t1**2 + 128*ms1**2*s**(-3)*t1 + 192*ms1**2*s**(-2) + 256*
     &    ms1**4*s**(-3) + 384*ms1**4*s**(-2)*t1**(-1) + 384*ms1**4*
     &    s**(-2)*u1**(-1) + 384*ms1**6*s**(-2)*t1**(-2) + 512*ms1**6*
     &    s**(-2)*t1**(-1)*u1**(-1) + 384*ms1**6*s**(-2)*u1**(-2) )
      STGGV = STGGV + SOF1(2)*Nc*Co*Pi*alphas**3 * (  - 128*ms1**2*
     &    s**(-4)*t1**2 - 128*ms1**2*s**(-3)*t1 - 64*ms1**2*s**(-2) -
     &    256*ms1**4*s**(-3) - 128*ms1**4*s**(-2)*t1**(-1) - 128*ms1**4
     &    *s**(-2)*u1**(-1) - 128*ms1**6*s**(-2)*t1**(-2) - 128*ms1**6*
     &    s**(-2)*u1**(-2) )
      STGGV = STGGV + SOF1(3)*Cqed*Pi*alphas**3 * (  - 64*ms1**2*
     &    s**(-2) - 128*ms1**4*s**(-2)*t1**(-1) - 128*ms1**4*s**(-2)*
     &    u1**(-1) - 128*ms1**6*s**(-2)*t1**(-2) - 256*ms1**6*s**(-2)*
     &    t1**(-1)*u1**(-1) - 128*ms1**6*s**(-2)*u1**(-2) )
      STGGV = STGGV + SOF1(3)*Nc*Ck*Pi*alphas**3 * ( 128*ms1**2*s**(-4)
     &    *t1**2 + 128*ms1**2*s**(-3)*t1 + 192*ms1**2*s**(-2) + 256*
     &    ms1**4*s**(-3) + 384*ms1**4*s**(-2)*t1**(-1) + 384*ms1**4*
     &    s**(-2)*u1**(-1) + 384*ms1**6*s**(-2)*t1**(-2) + 512*ms1**6*
     &    s**(-2)*t1**(-1)*u1**(-1) + 384*ms1**6*s**(-2)*u1**(-2) )
      STGGV = STGGV + SOF1(3)*Nc*Co*Pi*alphas**3 * (  - 128*ms1**2*
     &    s**(-4)*t1**2 - 128*ms1**2*s**(-3)*t1 - 64*ms1**2*s**(-2) -
     &    256*ms1**4*s**(-3) - 128*ms1**4*s**(-2)*t1**(-1) - 128*ms1**4
     &    *s**(-2)*u1**(-1) - 128*ms1**6*s**(-2)*t1**(-2) - 128*ms1**6*
     &    s**(-2)*u1**(-2) )
      STGGV = STGGV + SOF1(4)*Nc*Ck*Pi*alphas**3 * ( 64*s**(-2)*t1 + 64
     &    *s**(-1) + 128*ms1**2*s**(-1)*t1**(-1) - 128*ms1**4*s**(-2)*
     &    t1**(-1) - 128*ms1**4*s**(-2)*u1**(-1) + 128*ms1**4*s**(-1)*
     &    t1**(-2) )
      STGGV = STGGV + SOF1(4)*Nc*Co*Pi*alphas**3 * (  - 64*s**(-4)*
     &    t1**3 - 64*s**(-3)*t1**2 - 128*ms1**2*s**(-3)*t1 + 128*ms1**4
     &    *s**(-2)*u1**(-1) )
      STGGV = STGGV + SOF1(5)*Nc*Ck*Pi*alphas**3 * (  - 64*s**(-2)*t1
     &     + 128*ms1**2*s**(-1)*u1**(-1) - 128*ms1**4*s**(-2)*t1**(-1)
     &     - 128*ms1**4*s**(-2)*u1**(-1) + 128*ms1**4*s**(-1)*u1**(-2)
     &     )
      STGGV = STGGV + SOF1(5)*Nc*Co*Pi*alphas**3 * ( 64*s**(-4)*t1**3
     &     + 128*s**(-3)*t1**2 + 64*s**(-2)*t1 + 128*ms1**2*s**(-3)*t1
     &     + 128*ms1**2*s**(-2) + 128*ms1**4*s**(-2)*t1**(-1) )
      STGGV = STGGV + SOF1(6)*Nc*Ck*Pi*alphas**3 * (  - 64*s**(-2)*t1
     &     + 128*ms1**2*s**(-1)*u1**(-1) - 128*ms1**4*s**(-2)*t1**(-1)
     &     - 128*ms1**4*s**(-2)*u1**(-1) + 128*ms1**4*s**(-1)*u1**(-2)
     &     )
      STGGV = STGGV + SOF1(6)*Nc*Co*Pi*alphas**3 * ( 64*s**(-4)*t1**3
     &     + 128*s**(-3)*t1**2 + 64*s**(-2)*t1 + 128*ms1**2*s**(-3)*t1
     &     + 128*ms1**2*s**(-2) + 128*ms1**4*s**(-2)*t1**(-1) )
      STGGV = STGGV + SOF1(7)*Nc*Ck*Pi*alphas**3 * ( 64*s**(-2)*t1 + 64
     &    *s**(-1) + 128*ms1**2*s**(-1)*t1**(-1) - 128*ms1**4*s**(-2)*
     &    t1**(-1) - 128*ms1**4*s**(-2)*u1**(-1) + 128*ms1**4*s**(-1)*
     &    t1**(-2) )
      STGGV = STGGV + SOF1(7)*Nc*Co*Pi*alphas**3 * (  - 64*s**(-4)*
     &    t1**3 - 64*s**(-3)*t1**2 - 128*ms1**2*s**(-3)*t1 + 128*ms1**4
     &    *s**(-2)*u1**(-1) )
      STGGV = STGGV + SOF1(8)*Nc*Co*Pi*alphas**3 * ( 128*s**(-3)*t1**2
     &     + 128*s**(-2)*t1 + 64*s**(-1) + 256*ms1**2*s**(-2) + 128*
     &    ms1**2*s**(-1)*t1**(-1) + 128*ms1**2*s**(-1)*u1**(-1) + 128*
     &    ms1**4*s**(-1)*t1**(-2) + 128*ms1**4*s**(-1)*u1**(-2) )

      BB_GGV = fac * STGGV  + dec

      end

max_stars_repo_name: sliem/docker-hep

max_issues_repo_name: sliem/docker-hep

max_forks_repo_name: sliem/docker-hep


text:
subroutine nsc_endite(a,itask)
   !-----------------------------------------------------------------------
   ! DESCRIPTION
   !    This routine checks convergence and performs updates at:
   !    - itask=1 The end of an internal iteration
   !    - itask=2 The end of the external loop iteration
   !-----------------------------------------------------------------------
   use typre

   use def_parame
   use Mod_NSCompressible

   implicit none
   class(NSCompressibleProblem) :: a
   integer(ip) :: itask

   select case(itask)

   case(0)

      call a%SpecificNSCompEndite(0)

   case(1)

      call a%SpecificNSCompEndite(1)

   case(2)

      call a%SpecificNSCompEndite(2)

   end select

end subroutine nsc_endite

max_stars_repo_name: ciaid-colombia/InsFEM

max_issues_repo_name: ciaid-colombia/InsFEM

max_forks_repo_name: ciaid-colombia/InsFEM


text:
!     Floating point initialization of integer parameter

!     Make sure partial eval is OK

      program main

      integer*4 i
      parameter ( i = 100 )
      integer*4 j
      parameter ( j = i*0.1 )

      integer*4 array

      common / C1 /
     *     array(j)

      array(1) = 2*j
      print *,array(1)

      end

ext: f


text:
      subroutine coorot(y,z,yz,numnp)
c     implicit double precision (a-h,o-z)                                    dp
c
c     write nodal coordinates on plot tape
c
      common/bk05/ifil,iadd,maxsiz,head(12)
      common/bk14/lfna(15),lfnt(6)
      common/bk20/ntotal
      common/bk34/bb(1)
      dimension y(*),z(*),yz(*)
c
      do 10 j=1,numnp
      yz(2*j-1)=y(j)
      yz(2*j)=z(j)
   10 continue
      nn2=2*numnp
ck      call wrabsg(lfna(7),yz(1),nn2,iadd,bb(ntotal+1),0)
ck      call riosta(lfna(7))
ck      iadd=iadd+nn2
c
      return
c
      end

max_stars_repo_name: imohame/LabCode

max_issues_repo_name: imohame/LabCode

max_forks_repo_name: imohame/LabCode