SUBROUTINE UMIXMODEFATIGUE(DADN,GI_MAX,GII_MAX,GIII_MAX,
     $  GI_MIN,GII_MIN,GIII_MIN,TEMP,DTEMP,PREDEF,DPRED,NFIELD,NPROPS,
     $  PROPS,NSTATV,STATEV,NIARRAY,I_ARRAY,NRARRAY,R_ARRAY,NCARRAY,
     $  C_ARRAY, DA)
c     
      INCLUDE 'ABA_PARAM.INC'
c     
      dimension props(nprops),statev(nstatv), predef(nfield),
     $     dpred(nfield),i_array(niarray), r_array(nrarray)
      character*80 c_array(ncarray)
      parameter (tol=1.0d-6)

       gTT = GI_MAX+GII_MAX+GIII_MAX
       gshear=GII_MAX+GIII_MAX
       test=1.0
       if (abs(gTT).gt.abs(gshear)*tol) test=0.0
       gModeR= abs(gshear)*(1.0-test)/(abs(gTT)+test)
       an1=props(1)
       an2=props(2)
       c1=props(3)
       c2=props(4)
       if(gModeR.gt.0.5)Then
          dadN=c1*gTT**an1
       else
          dadN=c2*gTT**an2
       endif
c       write(*,*) 'an1,c1',an1,an2,c1,c2
c       write(*,*) 'gtt','gmodeR','dadn',gtt,gmodeR,dadn
      return
      end