-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathERF__AddTKESources_8H_source.html
More file actions
151 lines (149 loc) · 15.5 KB
/
ERF__AddTKESources_8H_source.html
File metadata and controls
151 lines (149 loc) · 15.5 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "https://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.9.1"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>ERF: Source/Diffusion/ERF_AddTKESources.H Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="navtree.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="resize.js"></script>
<script type="text/javascript" src="navtreedata.js"></script>
<script type="text/javascript" src="navtree.js"></script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/searchdata.js"></script>
<script type="text/javascript" src="search/search.js"></script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 56px;">
<td id="projectalign" style="padding-left: 0.5em;">
<div id="projectname">ERF
</div>
<div id="projectbrief">Energy Research and Forecasting: An Atmospheric Modeling Code</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.9.1 -->
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */
var searchBox = new SearchBox("searchBox", "search",false,'Search','.html');
/* @license-end */
</script>
<script type="text/javascript" src="menudata.js"></script>
<script type="text/javascript" src="menu.js"></script>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */
$(function() {
initMenu('',true,false,'search.php','Search');
$(document).ready(function() { init_search(); });
});
/* @license-end */</script>
<div id="main-nav"></div>
</div><!-- top -->
<div id="side-nav" class="ui-resizable side-nav-resizable">
<div id="nav-tree">
<div id="nav-tree-contents">
<div id="nav-sync" class="sync"></div>
</div>
</div>
<div id="splitbar" style="-moz-user-select:none;"
class="ui-resizable-handle">
</div>
</div>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */
$(document).ready(function(){initNavTree('ERF__AddTKESources_8H_source.html',''); initResizable(); });
/* @license-end */
</script>
<div id="doc-content">
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
</div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<div class="header">
<div class="headertitle">
<div class="title">ERF_AddTKESources.H</div> </div>
</div><!--header-->
<div class="contents">
<a href="ERF__AddTKESources_8H.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span>  <span class="comment">// Using Deardorff (see Sullivan et al 1994)</span></div>
<div class="line"><a name="l00002"></a><span class="lineno"> 2</span>  <span class="comment">// or k-eqn RANS (see Axell & Liungman 2001)</span></div>
<div class="line"><a name="l00003"></a><span class="lineno"> 3</span>  <span class="comment">//</span></div>
<div class="line"><a name="l00004"></a><span class="lineno"> 4</span>  <span class="comment">// Note: At this point, the thermal diffusivity ("Khv" field in ERF), the</span></div>
<div class="line"><a name="l00005"></a><span class="lineno"> 5</span>  <span class="comment">// subgrid heat flux ("hfx_z" here), and the subgrid dissipation</span></div>
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span>  <span class="comment">// ("diss" here) have been updated by ComputeTurbulentViscosityLES --</span></div>
<div class="line"><a name="l00007"></a><span class="lineno"> 7</span>  <span class="comment">// at the beginning of each timestep.</span></div>
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span>  <span class="comment">// The strain rate magnitude is updated at the beginning of the first</span></div>
<div class="line"><a name="l00009"></a><span class="lineno"> 9</span>  <span class="comment">// RK stage only, therefore the shear production term also does not</span></div>
<div class="line"><a name="l00010"></a><span class="lineno"> 10</span>  <span class="comment">// change between RK stages.</span></div>
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span>  <span class="comment">// The surface heat flux hfx_z(i,j,-1) is updated in MOSTStress at</span></div>
<div class="line"><a name="l00012"></a><span class="lineno"> 12</span>  <span class="comment">// each RK stage, but that does not change the buoyancy production term here.</span></div>
<div class="line"><a name="l00013"></a><span class="lineno"><a class="line" href="ERF__AddTKESources_8H.html#a15559095c80f50ca836f275efe183f4f"> 13</a></span>  <span class="keywordflow">if</span> (<a class="code" href="ERF__SetupDiff_8H.html#a4b0a122076eb996c681980a57e58bed4">l_use_keqn</a> && (start_comp <= <a class="code" href="ERF__IndexDefines_8H.html#ae15721db2cb4c2ab4732afe226bbba64">RhoKE_comp</a>) && (<a class="code" href="ERF__SetupDiff_8H.html#a29fc1da4eb54fbf2bb041845ce70d50d">end_comp</a> >= <a class="code" href="ERF__IndexDefines_8H.html#ae15721db2cb4c2ab4732afe226bbba64">RhoKE_comp</a>)) {</div>
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span>  <span class="keywordtype">int</span> qty_index = <a class="code" href="ERF__IndexDefines_8H.html#ae15721db2cb4c2ab4732afe226bbba64">RhoKE_comp</a>;</div>
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span>  <a class="code" href="ERF__InitCustomPert__DataAssimilation__ISV_8H.html#aada27a5057153b6463a47f463ea98742">ParallelFor</a>(bx,[=] AMREX_GPU_DEVICE (<span class="keywordtype">int</span> i, <span class="keywordtype">int</span> j, <span class="keywordtype">int</span> k) noexcept</div>
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span>  {</div>
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span>  <span class="comment">// Add Buoyancy Source</span></div>
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span>  <span class="comment">// where the SGS buoyancy flux tau_{theta,i} = -KH * dtheta/dx_i,</span></div>
<div class="line"><a name="l00019"></a><span class="lineno"> 19</span>  <span class="comment">// such that for dtheta/dz < 0, there is a positive (upward) heat</span></div>
<div class="line"><a name="l00020"></a><span class="lineno"> 20</span>  <span class="comment">// flux; the TKE buoyancy production is then</span></div>
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span>  <span class="comment">// B = g/theta_0 * tau_{theta,w}</span></div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span>  <span class="comment">// for a dry atmosphere.</span></div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span>  <span class="comment">// TODO: To account for moisture, the Brunt-Vaisala frequency,</span></div>
<div class="line"><a name="l00024"></a><span class="lineno"> 24</span>  <span class="comment">// N^2 = g[1/theta * dtheta/dz + ...]</span></div>
<div class="line"><a name="l00025"></a><span class="lineno"> 25</span>  <span class="comment">// **should** be a function of the water vapor and total water</span></div>
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span>  <span class="comment">// mixing ratios, depending on whether conditions are saturated or</span></div>
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span>  <span class="comment">// not (see the WRF model description, Skamarock et al 2019).</span></div>
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span>  cell_rhs(i,j,k,qty_index) += l_abs_g * <a class="code" href="ERF__SetupDiff_8H.html#af95a27e0d3fb9866f74398b99e1b1bd8">l_inv_theta0</a> * hfx_z(i,j,k);</div>
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span>  <span class="keywordflow">if</span> (!<a class="code" href="ERF__SetupDiff_8H.html#a5d0e7ac164885f8b767a91f1e26474a5">use_ref_theta</a>) {</div>
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span>  <span class="comment">// l_inv_theta0 == 1, divide by actual theta</span></div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span>  cell_rhs(i,j,k,qty_index) /= cell_prim(i,j,k,<a class="code" href="ERF__IndexDefines_8H.html#a8b9c5b827b61f426a3879bc5d2883040">PrimTheta_comp</a>);</div>
<div class="line"><a name="l00032"></a><span class="lineno"> 32</span>  }</div>
<div class="line"><a name="l00033"></a><span class="lineno"> 33</span>  </div>
<div class="line"><a name="l00034"></a><span class="lineno"> 34</span>  <span class="comment">// TKE shear production</span></div>
<div class="line"><a name="l00035"></a><span class="lineno"> 35</span>  <span class="comment">// P = -tau_ij * S_ij = 2 * mu_turb * S_ij * S_ij</span></div>
<div class="line"><a name="l00036"></a><span class="lineno"> 36</span>  <span class="comment">// Note: This assumes that the horizontal and vertical diffusivities</span></div>
<div class="line"><a name="l00037"></a><span class="lineno"> 37</span>  <span class="comment">// of momentum are equal</span></div>
<div class="line"><a name="l00038"></a><span class="lineno"> 38</span>  cell_rhs(i,j,k,qty_index) += <a class="code" href="ERF__Constants_8H.html#a2d515b2cda500f9935528386cf6c707a">two</a>*mu_turb(i,j,k,<a class="code" href="namespaceEddyDiff.html#ab21f065008f527d980eaad08e1352205a61ed124b350a1b26da1bd26aebfa5f82">EddyDiff::Mom_v</a>) * SmnSmn_a(i,j,k);</div>
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span>  </div>
<div class="line"><a name="l00040"></a><span class="lineno"> 40</span>  <span class="comment">// TKE dissipation</span></div>
<div class="line"><a name="l00041"></a><span class="lineno"> 41</span>  cell_rhs(i,j,k,qty_index) -= diss(i,j,k);</div>
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span>  });</div>
<div class="line"><a name="l00043"></a><span class="lineno"> 43</span>  }</div>
<div class="ttc" id="aERF__Constants_8H_html_a2d515b2cda500f9935528386cf6c707a"><div class="ttname"><a href="ERF__Constants_8H.html#a2d515b2cda500f9935528386cf6c707a">two</a></div><div class="ttdeci">constexpr amrex::Real two</div><div class="ttdef"><b>Definition:</b> ERF_Constants.H:8</div></div>
<div class="ttc" id="aERF__IndexDefines_8H_html_a8b9c5b827b61f426a3879bc5d2883040"><div class="ttname"><a href="ERF__IndexDefines_8H.html#a8b9c5b827b61f426a3879bc5d2883040">PrimTheta_comp</a></div><div class="ttdeci">#define PrimTheta_comp</div><div class="ttdef"><b>Definition:</b> ERF_IndexDefines.H:50</div></div>
<div class="ttc" id="aERF__IndexDefines_8H_html_ae15721db2cb4c2ab4732afe226bbba64"><div class="ttname"><a href="ERF__IndexDefines_8H.html#ae15721db2cb4c2ab4732afe226bbba64">RhoKE_comp</a></div><div class="ttdeci">#define RhoKE_comp</div><div class="ttdef"><b>Definition:</b> ERF_IndexDefines.H:38</div></div>
<div class="ttc" id="aERF__InitCustomPert__DataAssimilation__ISV_8H_html_aada27a5057153b6463a47f463ea98742"><div class="ttname"><a href="ERF__InitCustomPert__DataAssimilation__ISV_8H.html#aada27a5057153b6463a47f463ea98742">ParallelFor</a></div><div class="ttdeci">ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const Real *dx=geomdata.CellSize();const Real x=(i+0.5) *dx[0];const Real y=(j+0.5) *dx[1];const Real Omg=erf_vortex_Gaussian(x, y, xc, yc, R, beta, sigma);const Real deltaT=-(gamma - 1.0)/(2.0 *sigma *sigma) *Omg *Omg;const Real rho_norm=std::pow(1.0+deltaT, inv_gm1);const Real T=(1.0+deltaT) *T_inf;const Real p=std::pow(rho_norm, Gamma)/Gamma *rho_0 *a_inf *a_inf;const Real rho_theta=rho_0 *rho_norm *(T *std::pow(p_0/p, rdOcp));state_pert(i, j, k, RhoTheta_comp)=rho_theta - getRhoThetagivenP(p_hse(i, j, k));const Real r2d_xy=std::sqrt((x-xc) *(x-xc)+(y-yc) *(y-yc));state_pert(i, j, k, RhoScalar_comp)=0.25 *(1.0+std::cos(PI *std::min(r2d_xy, R)/R));})</div></div>
<div class="ttc" id="aERF__SetupDiff_8H_html_a29fc1da4eb54fbf2bb041845ce70d50d"><div class="ttname"><a href="ERF__SetupDiff_8H.html#a29fc1da4eb54fbf2bb041845ce70d50d">end_comp</a></div><div class="ttdeci">const int end_comp</div><div class="ttdef"><b>Definition:</b> ERF_SetupDiff.H:11</div></div>
<div class="ttc" id="aERF__SetupDiff_8H_html_a4b0a122076eb996c681980a57e58bed4"><div class="ttname"><a href="ERF__SetupDiff_8H.html#a4b0a122076eb996c681980a57e58bed4">l_use_keqn</a></div><div class="ttdeci">bool l_use_keqn</div><div class="ttdef"><b>Definition:</b> ERF_SetupDiff.H:12</div></div>
<div class="ttc" id="aERF__SetupDiff_8H_html_a5d0e7ac164885f8b767a91f1e26474a5"><div class="ttname"><a href="ERF__SetupDiff_8H.html#a5d0e7ac164885f8b767a91f1e26474a5">use_ref_theta</a></div><div class="ttdeci">const bool use_ref_theta</div><div class="ttdef"><b>Definition:</b> ERF_SetupDiff.H:16</div></div>
<div class="ttc" id="aERF__SetupDiff_8H_html_af95a27e0d3fb9866f74398b99e1b1bd8"><div class="ttname"><a href="ERF__SetupDiff_8H.html#af95a27e0d3fb9866f74398b99e1b1bd8">l_inv_theta0</a></div><div class="ttdeci">const Real l_inv_theta0</div><div class="ttdef"><b>Definition:</b> ERF_SetupDiff.H:17</div></div>
<div class="ttc" id="anamespaceEddyDiff_html_ab21f065008f527d980eaad08e1352205a61ed124b350a1b26da1bd26aebfa5f82"><div class="ttname"><a href="namespaceEddyDiff.html#ab21f065008f527d980eaad08e1352205a61ed124b350a1b26da1bd26aebfa5f82">EddyDiff::Mom_v</a></div><div class="ttdeci">@ Mom_v</div><div class="ttdef"><b>Definition:</b> ERF_IndexDefines.H:193</div></div>
</div><!-- fragment --></div><!-- contents -->
</div><!-- doc-content -->
<!-- start footer part -->
<div id="nav-path" class="navpath"><!-- id is needed for treeview function! -->
<ul>
<li class="navelem"><a class="el" href="dir_74389ed8173ad57b461b9d623a1f3867.html">Source</a></li><li class="navelem"><a class="el" href="dir_fc0c15590968603ac439df7b0a572228.html">Diffusion</a></li><li class="navelem"><a class="el" href="ERF__AddTKESources_8H.html">ERF_AddTKESources.H</a></li>
<li class="footer">Generated by <a href="https://www.doxygen.org/index.html"><img class="footer" src="doxygen.svg" width="104" height="31" alt="doxygen"/></a> 1.9.1 </li>
</ul>
</div>
</body>
</html>