subroutine Ampvirt_gggg(j1,j2,j3,j4,A)
      implicit none
      include 'constants.f'
      include 'zprods_com.f'
      double complex A(3,2,2,2,2)
      double complex A1Hggggpppp,A1Hggggpmmm,A1Hggggmmpp,A1Hggggmpmp
      integer j1,j2,j3,j4,i1(3),i2(3),i3(3),i4(4),j

      i1(1)=j1      
      i2(1)=j2      
      i3(1)=j3      
      i4(1)=j4      

      i1(2)=j1      
      i2(2)=j2      
      i3(2)=j4      
      i4(2)=j3      

      i1(3)=j1      
      i2(3)=j4      
      i3(3)=j2      
      i4(3)=j3      

      do j=1,3
      A(j,2,2,2,2)=A1Hggggpppp(i1(j),i2(j),i3(j),i4(j),za,zb)
      A(j,1,1,1,1)=A1Hggggpppp(i1(j),i2(j),i3(j),i4(j),zb,za)

      A(j,2,1,1,1)=A1Hggggpmmm(i1(j),i2(j),i3(j),i4(j),za,zb)
      A(j,1,2,1,1)=A1Hggggpmmm(i2(j),i3(j),i4(j),i1(j),za,zb)
      A(j,1,1,2,1)=A1Hggggpmmm(i3(j),i4(j),i1(j),i2(j),za,zb)
      A(j,1,1,1,2)=A1Hggggpmmm(i4(j),i1(j),i2(j),i3(j),za,zb)

      A(j,1,2,2,2)=A1Hggggpmmm(i1(j),i2(j),i3(j),i4(j),zb,za)
      A(j,2,1,2,2)=A1Hggggpmmm(i2(j),i3(j),i4(j),i1(j),zb,za)
      A(j,2,2,1,2)=A1Hggggpmmm(i3(j),i4(j),i1(j),i2(j),zb,za)
      A(j,2,2,2,1)=A1Hggggpmmm(i4(j),i1(j),i2(j),i3(j),zb,za)

      A(j,1,1,2,2)=A1Hggggmmpp(i1(j),i2(j),i3(j),i4(j),za,zb)
      A(j,1,2,2,1)=A1Hggggmmpp(i4(j),i1(j),i2(j),i3(j),za,zb)
      A(j,2,2,1,1)=A1Hggggmmpp(i3(j),i4(j),i1(j),i2(j),za,zb) 
      A(j,2,1,1,2)=A1Hggggmmpp(i2(j),i3(j),i4(j),i1(j),za,zb)

      A(j,1,2,1,2)=A1Hggggmpmp(i1(j),i2(j),i3(j),i4(j),za,zb)
      A(j,2,1,2,1)=A1Hggggmpmp(i2(j),i3(j),i4(j),i1(j),za,zb)
      enddo
      
      return
      end