Model Diagnostics - Revision history

Robertson at 13:07, 5 August 2015

2015-08-05T13:07:58Z

Robertson at 13:27, 5 February 2010

2010-02-05T13:27:23Z

Robertson at 20:39, 4 February 2010

2010-02-04T20:39:03Z

@@ Line 99: / Line 99: @@
 | <span class="blue">v_vvisc</span>
 | Vertical viscosity contribution*
 |}
@@ Line 154: / Line 136: @@
 | Net time rate-of-change
 |}
 where <span class="blue">Tname</span> represents the name of the tracer (''i.e.'' temp, salt, etc.).

Robertson at 20:34, 4 February 2010

2010-02-04T20:34:58Z

← Older revision		Revision as of 20:34, 4 February 2010
Line 49:		Line 49:
	\| <span class="blue">vbar_strs</span>		\| <span class="blue">vbar_strs</span>
	\| Stress divergence contribution		\| Stress divergence contribution
	\|-
	~~\| <span class="blue">pvorticity_bar</span>~~
	~~\| <span class="blue">-</span>~~
	~~\| Depth integrated 2D potential vorticity~~
	\|-
	~~\| <span class="blue">rvorticity_bar</span>~~
	~~\| <span class="blue">-</span>~~
	~~\| Depth integrated 2D relative vorticity~~
	\|}		\|}

Line 107:		Line 99:
	\| <span class="blue">v_vvisc</span>		\| <span class="blue">v_vvisc</span>
	\| Vertical viscosity contribution*		\| Vertical viscosity contribution*
			\|}


			:{\| border="1" cellpadding="5" cellspacing="0"
			\|-
			! colspan="2" \| Vorticity Diagnostics
			\|-
			\| <span class="blue">pvorticity_bar</span>
			\| Depth integrated 2D potential vorticity (1/m/s)
			\|-
			\| <span class="blue">rvorticity_bar</span>
			\| Depth integrated 2D relative vorticity (1/s)
	\|-		\|-
	\| <span class="blue">pvorticity</span>		\| <span class="blue">pvorticity</span>
	~~\| <span class="blue">-</span>~~		\| 3D potential vorticity (1/m/s)
	\| 3D potential vorticity
	\|-		\|-
	\| <span class="blue">rvorticity</span>		\| <span class="blue">rvorticity</span>
	~~\| <span class="blue">-</span>~~		\| 3D relative vorticity (1/s)
	\| 3D relative vorticity
	\|}		\|}

Robertson at 20:25, 4 February 2010

2010-02-04T20:25:47Z

@@ Line 4: / Line 4: @@
 ==Output Variables==
-The tables below show the variables that are written to the [[Variables#DIAname|diagnostics output file]].
+<wikitex>The tables below show the variables that are written to the [[Variables#DIAname|diagnostics output file]].
 :{| border="1" cellpadding="5" cellspacing="0"
@@ Line 21: / Line 21: @@
 | <span class="blue">vbar_hadv</span>
 | Horizontal advection contribution*
 |-
 | <span class="blue">ubar_prsgrd</span>
@@ Line 29: / Line 49: @@
 | <span class="blue">vbar_strs</span>
 | Stress divergence contribution
 |}
@@ Line 47: / Line 75: @@
 | <span class="blue">v_hadv</span>
 | Horizontal advection contribution*
 |-
 | <span class="blue">u_hvisc</span>
 | <span class="blue">v_hvisc</span>
 | Horizontal viscosity contribution*
 |-
 | <span class="blue">u_prsgrd</span>
 | <span class="blue">v_prsgrd</span>
 | Pressure gradient contribution
 |-
 | <span class="blue">u_vvisc</span>
 | <span class="blue">v_vvisc</span>
 | Vertical viscosity contribution*
 |}
@@ Line 72: / Line 124: @@
 | <span class="blue">Tname_hadv</span>
 | Horizontal advection contribution
 |-
 | <span class="blue">Tname_hdiff</span>
 | Horizontal diffusivity contribution*
 |-
 | <span class="blue">Tname_vadv</span>
@@ Line 93: / Line 160: @@
 Because the diagnostic terms are accelerations (changes in current velocity over time), the terms from the depth-average momentum calculations must be integrated over all short timesteps when the model is configured in 3D.
+</wikitex>
 ==Derivation of the 2D momentum diagnostic terms integrated over short timesteps.==
 <wikitex>

Robertson: /* Derivation of the 2D momentum diagnostic terms integrated over short timesteps. */

2010-01-07T19:22:48Z

Derivation of the 2D momentum diagnostic terms integrated over short timesteps.

← Older revision		Revision as of 19:22, 7 January 2010
Line 95:		Line 95:

	==Derivation of the 2D momentum diagnostic terms integrated over short timesteps.==		==Derivation of the 2D momentum diagnostic terms integrated over short timesteps.==
	:<~~math~~>\bar u_{n+1} = \frac{(\bar u_{n} D_{n} + R_{n+1})}{D_{n+1}}\;\;\;\;and\;\;\;\;\bar u_{n} = \frac{(\bar u_{n-1} D_{n-1} + R_{n})}{D_{n}}~~</math>~~		<wikitex>
			$$\bar u_{n+1} = \frac{(\bar u_{n} D_{n} + R_{n+1})}{D_{n+1}}\;\;\;\;and\;\;\;\;\bar u_{n} = \frac{(\bar u_{n-1} D_{n-1} + R_{n})}{D_{n}}$$


	where ~~<span class="blue">D<sub>n</sub></span>~~ is water column depth at barotropic (short) timestep		where $D_n$ is water column depth at barotropic (short) timestep
	~~<span class="blue">~~n~~</span>~~ and ~~<span class="blue">R<sub>n</sub></span>~~ is the right-hand-side term.~~<br>~~		$n$ and $R_n$ is the right-hand-side term.
	Combining:		Combining:


	~~:<math>~~\bar u_{n+1} = \frac{\left(\frac{(\bar u_{n-1}D_{n-1} + R_{n})}{D_{n}}D_{n} + R_{n+1}\right)}{D_{n+1}}=\frac{(\bar u_{n-1}D_{n-1} + R_{n} + R_{n+1})}{D_{n+1}}~~</math>~~		$$\bar u_{n+1} = \frac{\left(\frac{(\bar u_{n-1}D_{n-1} + R_{n})}{D_{n}}D_{n} + R_{n+1}\right)}{D_{n+1}}=\frac{(\bar u_{n-1}D_{n-1} + R_{n} + R_{n+1})}{D_{n+1}}$$


Line 109:		Line 110:


	~~:<math>~~\~~begin~~{~~align}~~\bar u_{n+1} &= \frac{(\bar u_o D_o + R_1 + R_2 + R_3 + \ldots + R_{n-1} + R_n + R_{n+1})}{D_{n+1}} \\ \bar u_{N+1} &= \frac{\left(\bar u_oD_o + \sum_{iif = 1}^{n+1} R_{iif}\right)}{D_{n+1}} \\ DU\_avg1 &= \sum_{iif=1}^{ndtfast} \frac{1}{n_u} \left\{ weight(1,iif)\; \bar u_{iif}^* \zeta_{iif}^\right\} \\ DiaU2int &= \sum_{iif=1}^{ndtfast} \frac{1}{n_u} \left\{ weight(1,iif)\; \zeta_{iif}^ \; \sum_{n=1}^{iif} R_n^* \right\}~~\end{align~~}~~</math>~~		$$ \eqalign{\bar u_{n+1} &= \frac{(\bar u_o D_o + R_1 + R_2 + R_3 + \ldots + R_{n-1} + R_n + R_{n+1})}
			{D_{n+1}} \cr
			\noalign{\medskip}
			\bar u_{N+1} &= \frac{\left(\bar u_oD_o + \sum_{iif = 1}^{n+1} R_{iif}\right)}{D_{n+1}} \cr
			\noalign{\medskip}
			DU\_avg1 &= \sum_{iif=1}^{ndtfast} \frac{1}{n_u} \left\{ weight(1,iif)\; \bar u_{iif}^* \zeta_{iif}^*\right\} \cr
			\noalign{\medskip}
			DiaU2int &= \sum_{iif=1}^{ndtfast} \frac{1}{n_u} \left\{ weight(1,iif)\; \zeta_{iif}^* \; \sum_{n=1}^{iif} R_n^* \right\} }$$


	where the asterisks ( * ) superscript denotes the value calculated during corrector timestep.		where the asterisks ( * ) superscript denotes the value calculated during corrector timestep.
			</wikitex>

Robertson: Diagnostics moved to Model Diagnostics: Wrong name

2008-04-14T19:28:52Z

Diagnostics moved to Model Diagnostics: Wrong name

← Older revision	Revision as of 19:28, 14 April 2008
(No difference)

Robertson at 19:28, 14 April 2008

2008-04-14T19:28:26Z

← Older revision		Revision as of 19:28, 14 April 2008
Line 1:		Line 1:
	<div class="title">Diagnostics</div>		<div class="title">Model Diagnostics</div>
	Time-averaged values for each of the terms in the momentum and tracer equations can be written out by defining the variables [[Options#DIAGNOSTICS_UV\|DIAGNOSTICS_UV]] and [[Options#DIAGNOSTICS_TS\|DIAGNOSTICS_TS]], respectively, in [[cppdefs.h]].		Time-averaged values for each of the terms in the momentum and tracer equations can be written out by defining the variables [[Options#DIAGNOSTICS_UV\|DIAGNOSTICS_UV]] and [[Options#DIAGNOSTICS_TS\|DIAGNOSTICS_TS]], respectively, in [[cppdefs.h]].

Robertson: New page:
Diagnostics
Time-averaged values for each of the terms in the momentum and tracer equations can be written out by defining the variables [[Options#DIAGNOSTICS_UV|D...

2008-04-14T19:27:00Z

New page: <div class="title">Diagnostics</div> Time-averaged values for each of the terms in the momentum and tracer equations can be written out by defining the variables [[Options#DIAGNOSTICS_UV|D...

New page

<div class="title">Diagnostics</div>
Time-averaged values for each of the terms in the momentum and tracer equations can be written out by defining the variables [[Options#DIAGNOSTICS_UV|DIAGNOSTICS_UV]] and [[Options#DIAGNOSTICS_TS|DIAGNOSTICS_TS]], respectively, in [[cppdefs.h]].

==Output Variables==
The tables below show the variables that are written to the [[Variables#DIAname|diagnostics output file]].

:{| border="1" cellpadding="5" cellspacing="0"
|-
! colspan="3" | 2D Diagnostic momentum terms (m/s2)
|-
| ubar_accel
| vbar_accel
| Net acceleration
|-
| ubar_cor
| vbar_cor
| Coriolis contribution*
|-
| ubar_hadv
| vbar_hadv
| Horizontal advection contribution*
|-
| ubar_prsgrd
| vbar_prsgrd
| Pressure gradient contribution
|-
| ubar_strs
| vbar_strs
| Stress divergence contribution
|}

:{| border="1" cellpadding="5" cellspacing="0"
|-
! colspan="3" | 3D Diagnostic momentum terms (m/s2)
|-
| u_accel
| v_accel
| Net acceleration
|-
| u_cor
| v_cor
| Coriolis contribution*
|-
| u_hadv
| v_hadv
| Horizontal advection contribution*
|-
| u_hvisc
| v_hvisc
| Horizontal viscosity contribution*
|-
| u_prsgrd
| v_prsgrd
| Pressure gradient contribution
|-
| u_vadv
| v_vadv
| Vertical advection contribution*
|-
| u_vvisc
| v_vvisc
| Vertical viscosity contribution*
|}

:{| border="1" cellpadding="5" cellspacing="0"
|-
! colspan="2" | Tracer Diagnostic terms (Tunits/s)
|-
| Tname_hadv
| Horizontal advection contribution
|-
| Tname_hdiff
| Horizontal diffusivity contribution*
|-
| Tname_vadv
| Vertical advection contribution
|-
| Tname_vdiff
| Vertical diffusivity contribution
|-
| Tname_rate
| Net time rate-of-change
|}

where Tname represents the name of the tracer (''i.e.'' temp, salt, etc.).

{{note}}The diagnostics code adapts to the model configuration, the terms marked with * are written to the diagnostics output file when appropriate. Diagnostic terms are reported for all tracers.

Because the diagnostic terms are accelerations (changes in current velocity over time), the terms from the depth-average momentum calculations must be integrated over all short timesteps when the model is configured in 3D.

==Derivation of the 2D momentum diagnostic terms integrated over short timesteps.==
:<math>\bar u_{n+1} = \frac{(\bar u_{n} D_{n} + R_{n+1})}{D_{n+1}}\;\;\;\;and\;\;\;\;\bar u_{n} = \frac{(\bar u_{n-1} D_{n-1} + R_{n})}{D_{n}}</math>

where Dn is water column depth at barotropic (short) timestep
n and Rn is the right-hand-side term. 
Combining:

:<math>\bar u_{n+1} = \frac{\left(\frac{(\bar u_{n-1}D_{n-1} + R_{n})}{D_{n}}D_{n} + R_{n+1}\right)}{D_{n+1}}=\frac{(\bar u_{n-1}D_{n-1} + R_{n} + R_{n+1})}{D_{n+1}}</math>

Expanding:

:<math>\begin{align}\bar u_{n+1} &= \frac{(\bar u_o D_o + R_1 + R_2 + R_3 + \ldots + R_{n-1} + R_n + R_{n+1})}{D_{n+1}} \\ \bar u_{N+1} &= \frac{\left(\bar u_oD_o + \sum_{iif = 1}^{n+1} R_{iif}\right)}{D_{n+1}} \\ DU\_avg1 &= \sum_{iif=1}^{ndtfast} \frac{1}{n_u} \left\{ weight(1,iif)\; \bar u_{iif}^* \zeta_{iif}^*\right\} \\ DiaU2int &= \sum_{iif=1}^{ndtfast} \frac{1}{n_u} \left\{ weight(1,iif)\; \zeta_{iif}^* \; \sum_{n=1}^{iif} R_n^* \right\}\end{align}</math>

where the asterisks ( * ) superscript denotes the value calculated during corrector timestep.

← Older revision		Revision as of 13:27, 5 February 2010
Line 6:		Line 6:
	<wikitex>The tables below show the variables that are written to the [[Variables#DIAname\|diagnostics output file]].		<wikitex>The tables below show the variables that are written to the [[Variables#DIAname\|diagnostics output file]].

	:{\| border="1" cellpadding="5" cellspacing="0"		:{\| border="1" cellpadding="5" cellspacing="0" width="600"
	\|-		\|-
	! colspan="3" \| 2D Diagnostic momentum terms (m/s<sup>2</sup>)		! colspan="3" \| 2D Diagnostic momentum terms (m/s<sup>2</sup>)
	\|-		\|-
	\| <span class="blue">ubar_accel</span>		\| width="110" \| <span class="blue">ubar_accel</span>
	\| <span class="blue">vbar_accel</span>		\| width="110" \| <span class="blue">vbar_accel</span>
	\| Net acceleration		\| Net acceleration
	\|-		\|-
Line 52:		Line 52:


	:{\| border="1" cellpadding="5" cellspacing="0"		:{\| border="1" cellpadding="5" cellspacing="0" width="600"
	\|-		\|-
	! colspan="3" \| 3D Diagnostic momentum terms (m/s<sup>2</sup>)		! colspan="3" \| 3D Diagnostic momentum terms (m/s<sup>2</sup>)
	\|-		\|-
	\| <span class="blue">u_accel</span>		\| width="110" \| <span class="blue">u_accel</span>
	\| <span class="blue">v_accel</span>		\| width="110" \| <span class="blue">v_accel</span>
	\| Net acceleration		\| Net acceleration
	\|-		\|-
Line 102:		Line 102:


	:{\| border="1" cellpadding="5" cellspacing="0"		:{\| border="1" cellpadding="5" cellspacing="0" width="600"
	\|-		\|-
	! colspan="2" \| Tracer Diagnostic terms (Tunits/s)		! colspan="2" \| Tracer Diagnostic terms (Tunits/s)
	\|-		\|-
	\| <span class="blue">Tname_hadv</span>		\| width="110" \| <span class="blue">Tname_hadv</span>
	\| Horizontal advection contribution		\| Horizontal advection contribution
	\|-		\|-
Line 138:		Line 138:


	:{\| border="1" cellpadding="5" cellspacing="0"		:{\| border="1" cellpadding="5" cellspacing="0" width="600"
	\|-		\|-
	! colspan="2" \| Vorticity Diagnostics		! colspan="2" \| Vorticity Diagnostics
	\|-		\|-
	\| <span class="blue">pvorticity_bar</span>		\| width="110" \| <span class="blue">pvorticity_bar</span>
	\| Depth integrated 2D potential vorticity (m<sup>-1</sup> s<sup>-1</sup>)		\| Depth integrated 2D potential vorticity (m<sup>-1</sup> s<sup>-1</sup>)
	\|-		\|-