diff --git a/manual/eqs/output.tex b/manual/eqs/output.tex
index 84ace12526..1f512b16ab 100644
--- a/manual/eqs/output.tex
+++ b/manual/eqs/output.tex
@@ -232,7 +232,7 @@ \subsection{~Output parameters} \label{sub:outpars}
       where
       \begin{equation} n = \frac{1}{2} + \frac{kd}{\sinh 2kd}
       \: . \label{eq:n} \end{equation}
-\item \textbf{TWO} Wave to ocean momentum flux
+\item \textbf{TWO} Wave to ocean momentum flux (m$^2$/s$^2$)
 \item \textbf{BHD} Bernoulli head (m$^2$/s$^2$)
       \begin{equation} J =  g \int \!\!\!\! \int  \frac{k}
         {\sinh 2kd}  \: F(k,\theta) \: dk d\theta
@@ -265,6 +265,9 @@ \subsection{~Output parameters} \label{sub:outpars}
         {\sinh^2 kd}  \: F(k,\theta) \: dk d\theta
       \: . \label{eq:Ussp} \end{equation}
 The number $n$ and the central wavenumbers for the partitions $k_n$ are `run-time' in the sense that they are defined in the model definition file generated from ww3\_grid.inp (see the sample input file for details on setting these quantities).  The wavenumber integral bounds for each $n$ partition divide the model wave spectra to accomodate the prescribed values of $k_n$.
+\item \textbf{TOC} Total ocean stress (Pa)
+      \begin{equation} \tau_{oc} = \tau_a - \tau_{aw} + \tau_{woc}
+      \: . \label{eq:Toc} \end{equation}
 \end{list}
 
 \item{Wave-bottom layer}
diff --git a/model/inp/ww3_shel.inp b/model/inp/ww3_shel.inp
index 775ccf1420..d10e3df652 100644
--- a/model/inp/ww3_shel.inp
+++ b/model/inp/ww3_shel.inp
@@ -197,6 +197,7 @@ $  F  F  6     9   P2SMS      P2L   Micro seism  source term
 $  F  F  6    10   TAUICE     TWI   Wave to sea ice stress
 $  F  F  6    11   PHICE      FIC   Wave to sea ice energy flux
 $  F  F  6    12   USSP       USP   Partitioned surface Stokes drift
+$  F  F  6    13   TAUOC[X,Y] TOC   Total momentum to the ocean
 $   -------------------------------------------------
 $        7                          Wave-bottom layer 
 $   -------------------------------------------------
diff --git a/model/nml/ww3_multi.nml b/model/nml/ww3_multi.nml
index 96b15c5d8d..1e851efe13 100644
--- a/model/nml/ww3_multi.nml
+++ b/model/nml/ww3_multi.nml
@@ -199,7 +199,7 @@
 !  EF TH1M STH1M TH2M STH2M WN
 !  PHS PTP PLP PDIR PSPR PWS PDP PQP PPE PGW PSW PTM10 PT01 PT02 PEP TWS PNR
 !  UST CHA CGE FAW TAW TWA WCC WCF WCH WCM FWS
-!  SXY TWO BHD FOC TUS USS P2S USF P2L TWI FIC
+!  SXY TWO BHD FOC TUS USS P2S USF P2L TWI FIC USP TOC
 !  ABR UBR BED FBB TBB
 !  MSS MSC WL02 AXT AYT AXY
 !  DTD FC CFX CFD CFK
diff --git a/model/nml/ww3_shel.nml b/model/nml/ww3_shel.nml
index c5cd53af45..7989eda08a 100644
--- a/model/nml/ww3_shel.nml
+++ b/model/nml/ww3_shel.nml
@@ -187,6 +187,7 @@
 !  F  F  6    10   TAUICE     TWI   Wave to sea ice stress
 !  F  F  6    11   PHICE      FIC   Wave to sea ice energy flux
 !  F  F  6    12   USSP       USP   Partitioned surface Stokes drift
+!  F  F  6    13   TAUOC[X,Y] TOC   Total momentum to the ocean
 !   -------------------------------------------------
 !        7                          Wave-bottom layer 
 !   -------------------------------------------------
diff --git a/model/tools/bash/ww3_multi_inp2nml.sh b/model/tools/bash/ww3_multi_inp2nml.sh
index 9119b908d9..6059400689 100755
--- a/model/tools/bash/ww3_multi_inp2nml.sh
+++ b/model/tools/bash/ww3_multi_inp2nml.sh
@@ -672,7 +672,7 @@ cat >> $nmlfile << EOF
 !  EF TH1M STH1M TH2M STH2M WN
 !  PHS PTP PLP PDIR PSPR PWS PDP PQP PPE PGW PSW PTM10 PT01 PT02 PEP TWS PNR
 !  UST CHA CGE FAW TAW TWA WCC WCF WCH WCM FWS
-!  SXY TWO BHD FOC TUS USS P2S USF P2L TWI FIC
+!  SXY TWO BHD FOC TUS USS P2S USF P2L TWI FIC USP TOC
 !  ABR UBR BED FBB TBB
 !  MSS MSC WL02 AXT AYT AXY
 !  DTD FC CFX CFD CFK
diff --git a/model/tools/bash/ww3_ounf_inp2nml.sh b/model/tools/bash/ww3_ounf_inp2nml.sh
index 560a4f9649..c0c8ec4d14 100755
--- a/model/tools/bash/ww3_ounf_inp2nml.sh
+++ b/model/tools/bash/ww3_ounf_inp2nml.sh
@@ -183,7 +183,7 @@ cat >> $nmlfile << EOF
 !  EF TH1M STH1M TH2M STH2M WN
 !  PHS PTP PLP PDIR PSPR PWS PDP PQP PPE PGW PSW PTM10 PT01 PT02 PEP TWS PNR
 !  UST CHA CGE FAW TAW TWA WCC WCF WCH WCM FWS
-!  SXY TWO BHD FOC TUS USS P2S USF P2L TWI FIC
+!  SXY TWO BHD FOC TUS USS P2S USF P2L TWI FIC USP TOC
 !  ABR UBR BED FBB TBB
 !  MSS MSC WL02 AXT AYT AXY
 !  DTD FC CFX CFD CFK
diff --git a/model/tools/bash/ww3_shel_inp2nml.sh b/model/tools/bash/ww3_shel_inp2nml.sh
index de95b4a60b..6489da5015 100755
--- a/model/tools/bash/ww3_shel_inp2nml.sh
+++ b/model/tools/bash/ww3_shel_inp2nml.sh
@@ -829,6 +829,8 @@ cat >> $nmlfile << EOF
 !  F  F  6     9   P2SMS      P2L   Micro seism  source term
 !  F  F  6    10   TAUICE     TWI   Wave to sea ice stress
 !  F  F  6    11   PHICE      FIC   Wave to sea ice energy flux
+!  F  F  6    12   USSP       USP   Partitioned surface Stokes drift
+!  F  F  6    13   TAUOC[X,Y] TOC   Total momentum to the ocean
 !   -------------------------------------------------
 !        7                          Wave-bottom layer
 !   -------------------------------------------------
@@ -869,7 +871,7 @@ cat >> $nmlfile << EOF
 !
 ! * coupling fields exchanged list is :
 !   - Sent fields by ww3:
-!       - Ocean model : T0M1 OCHA OHS DIR BHD TWO UBR FOC TAW TUS USS LM DRY
+!       - Ocean model : T0M1 OCHA OHS DIR BHD TWO UBR FOC TAW TUS USS LM DRY TOC
 !       - Atmospheric model : ACHA AHS TP (or FP) FWS
 !       - Ice model : IC5 TWI
 !   - Received fields by ww3: