--- MITgcm/pkg/timeave/timeave_statv_write.F	2003/05/09 20:53:27	1.7
+++ MITgcm/pkg/timeave/timeave_statv_write.F	2003/08/03 02:58:29	1.8
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/timeave/timeave_statv_write.F,v 1.7 2003/05/09 20:53:27 jmc Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/timeave/timeave_statv_write.F,v 1.8 2003/08/03 02:58:29 jmc Exp $
 C $Name:  $
 #include "CPP_OPTIONS.h"
 
@@ -80,6 +80,10 @@
 
           CALL TIMEAVE_NORMALIZ(TdiffRtave,TimeAve_full,Nr,
      &                            bi,bj,myThid)
+#ifndef DISABLE_MOM_VECINV
+          CALL TIMEAVE_NORMALIZ(uZetatave,TimeAve_full,Nr,bi,bj,myThid)
+          CALL TIMEAVE_NORMALIZ(vZetatave,TimeAve_full,Nr,bi,bj,myThid)
+#endif
           CALL TIMEAVE_NORMALIZ(phiHydtave,TimeAve_full,Nr,
      &                            bi,bj,myThid)
           CALL TIMEAVE_NORMALIZ(phiHydLowtave,TimeAve_full,1,
@@ -123,6 +127,12 @@
 
        IF (useVariableK)
      & CALL WRITE_FLD_XYZ_RL('Tdiftave.',suff,TdiffRtave,myIter,myThid)
+#ifndef DISABLE_MOM_VECINV
+       IF (vectorInvariantMomentum) THEN
+         CALL WRITE_FLD_XYZ_RL('uZtave.',suff,uZetatave,myIter,myThid)
+         CALL WRITE_FLD_XYZ_RL('vZtave.',suff,vZetatave,myIter,myThid)
+       ENDIF
+#endif
        CALL WRITE_FLD_XYZ_RL('PhHytave.',suff,phiHydtave,myIter,myThid)
        IF ( buoyancyRelation .NE. 'ATMOSPHERIC' ) THEN
        CALL WRITE_FLD_XY_RL('PHLtave.',suff,phiHydLowtave,