Tweaked trace function to not use or pass the off-diagonal terms

src/module_invariants.f90

@@ -7,9 +7,16 @@ contains
 !this function computes the trace of a tensor
 !Tr[\sigma]=\sum_{i} \sigma_{ii}
 
-function trace (exx,eyy,ezz,exy,eyz,ezx)
+! Original version pre 05/11
+!function trace (exx,eyy,ezz,exy,eyz,ezx)
+!double precision trace
+!double precision exx,eyy,ezz,exy,eyz,ezx
+!trace=exx+eyy+ezz
+!end function trace
+
+function trace (exx,eyy,ezz)
 double precision trace
-double precision exx,eyy,ezz,exy,eyz,ezx
+double precision exx,eyy,ezz
 trace=exx+eyy+ezz
 end function trace