-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathERF__AddQKESources_8H_source.html
More file actions
134 lines (132 loc) · 15.3 KB
/
ERF__AddQKESources_8H_source.html
File metadata and controls
134 lines (132 loc) · 15.3 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
<!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_AddQKESources.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__AddQKESources_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_AddQKESources.H</div> </div>
</div><!--header-->
<div class="contents">
<a href="ERF__AddQKESources_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 PBL</span></div>
<div class="line"><a name="l00002"></a><span class="lineno"><a class="line" href="ERF__AddQKESources_8H.html#ad212ca7498a026ec3f73cef6a5f29738"> 2</a></span>  <span class="keywordflow">if</span> (<a class="code" href="ERF__SetupDiff_8H.html#a7b02b17a96d3d80f109a252dbf71ca41">l_use_mynn</a> && start_comp <= RhoKE_comp && end_comp >=<a class="code" href="ERF__IndexDefines_8H.html#ae15721db2cb4c2ab4732afe226bbba64">RhoKE_comp</a>) {</div>
<div class="line"><a name="l00003"></a><span class="lineno"> 3</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="l00004"></a><span class="lineno"> 4</span>  <span class="keyword">auto</span> pbl_mynn_B1_l = <a class="code" href="ERF__SetupVertDiff_8H.html#a4a8c0178282c029acc3d37ed429f9f58">turbChoice</a>.<a class="code" href="structTurbChoice.html#a11fd9429567c4dd4afa871170580aed0">pbl_mynn</a>.<a class="code" href="structMYNNLevel25.html#a8a1dab3e8da3715b39e07436759b3e80">B1</a>;</div>
<div class="line"><a name="l00005"></a><span class="lineno"> 5</span>  </div>
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span>  <a class="code" href="ERF__InitCustomPert__DensityCurrent_8H.html#aa73510fc34fbe6bb9e2610178ba9878f">ParallelFor</a>(bx,[=,moist_indices=solverChoice.moisture_indices] 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="l00007"></a><span class="lineno"> 7</span>  {</div>
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span>  <span class="keywordtype">bool</span> c_ext_dir_on_zlo = ( (bc_ptr[<a class="code" href="namespaceBCVars.html#a00367571247ba920989360b94c156f4aaaf601346b01655b32acf66d07b856aec">BCVars::cons_bc</a>].lo(2) == <a class="code" href="namespaceERFBCType.html#a40fbdd5a379a92debb4bcaaccac6e01ca4b16988b22c277c9f812ddb5fb1539af">ERFBCType::ext_dir</a>) );</div>
<div class="line"><a name="l00009"></a><span class="lineno"> 9</span>  <span class="keywordtype">bool</span> c_ext_dir_on_zhi = ( (bc_ptr[<a class="code" href="namespaceBCVars.html#a00367571247ba920989360b94c156f4aaaf601346b01655b32acf66d07b856aec">BCVars::cons_bc</a>].hi(2) == <a class="code" href="namespaceERFBCType.html#a40fbdd5a379a92debb4bcaaccac6e01ca4b16988b22c277c9f812ddb5fb1539af">ERFBCType::ext_dir</a>) );</div>
<div class="line"><a name="l00010"></a><span class="lineno"> 10</span>  <span class="keywordtype">bool</span> u_ext_dir_on_zlo = ( (bc_ptr[<a class="code" href="namespaceBCVars.html#a00367571247ba920989360b94c156f4aab1782753c68aaad8632b0a7d575e57f4">BCVars::xvel_bc</a>].lo(2) == <a class="code" href="namespaceERFBCType.html#a40fbdd5a379a92debb4bcaaccac6e01ca4b16988b22c277c9f812ddb5fb1539af">ERFBCType::ext_dir</a>) );</div>
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span>  <span class="keywordtype">bool</span> u_ext_dir_on_zhi = ( (bc_ptr[<a class="code" href="namespaceBCVars.html#a00367571247ba920989360b94c156f4aab1782753c68aaad8632b0a7d575e57f4">BCVars::xvel_bc</a>].hi(2) == <a class="code" href="namespaceERFBCType.html#a40fbdd5a379a92debb4bcaaccac6e01ca4b16988b22c277c9f812ddb5fb1539af">ERFBCType::ext_dir</a>) );</div>
<div class="line"><a name="l00012"></a><span class="lineno"> 12</span>  <span class="keywordtype">bool</span> v_ext_dir_on_zlo = ( (bc_ptr[<a class="code" href="namespaceBCVars.html#a00367571247ba920989360b94c156f4aa98ad22810fdfe9e585bc6481fc12454e">BCVars::yvel_bc</a>].lo(2) == <a class="code" href="namespaceERFBCType.html#a40fbdd5a379a92debb4bcaaccac6e01ca4b16988b22c277c9f812ddb5fb1539af">ERFBCType::ext_dir</a>) );</div>
<div class="line"><a name="l00013"></a><span class="lineno"> 13</span>  <span class="keywordtype">bool</span> v_ext_dir_on_zhi = ( (bc_ptr[<a class="code" href="namespaceBCVars.html#a00367571247ba920989360b94c156f4aa98ad22810fdfe9e585bc6481fc12454e">BCVars::yvel_bc</a>].hi(2) == <a class="code" href="namespaceERFBCType.html#a40fbdd5a379a92debb4bcaaccac6e01ca4b16988b22c277c9f812ddb5fb1539af">ERFBCType::ext_dir</a>) );</div>
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span>  </div>
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span>  <span class="comment">// This computes shear production, buoyancy production, and dissipation terms only.</span></div>
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span>  cell_rhs(i, j, k, qty_index) += <a class="code" href="ERF__PBLModels_8H.html#aa037b629b806251024dce3d3e714de7b">ComputeQKESourceTerms</a>(i,j,k,u,v,cell_data,cell_prim,</div>
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span>  mu_turb,cellSizeInv,domain,</div>
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span>  pbl_mynn_B1_l,tm_arr(i,j,0),</div>
<div class="line"><a name="l00019"></a><span class="lineno"> 19</span>  moist_indices,</div>
<div class="line"><a name="l00020"></a><span class="lineno"> 20</span>  c_ext_dir_on_zlo, c_ext_dir_on_zhi,</div>
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span>  u_ext_dir_on_zlo, u_ext_dir_on_zhi,</div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span>  v_ext_dir_on_zlo, v_ext_dir_on_zhi);</div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span>  });</div>
<div class="line"><a name="l00024"></a><span class="lineno"> 24</span>  }</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__DensityCurrent_8H_html_aa73510fc34fbe6bb9e2610178ba9878f"><div class="ttname"><a href="ERF__InitCustomPert__DensityCurrent_8H.html#aa73510fc34fbe6bb9e2610178ba9878f">ParallelFor</a></div><div class="ttdeci">ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const auto prob_lo=geomdata.ProbLo();const auto dx=geomdata.CellSize();const Real x=(prob_lo[0]+(i+myhalf) *dx[0])/mf_m(i, j, 0);const Real z=z_cc(i, j, k);Real L=std::sqrt(std::pow((x - x_c)/x_r, 2)+std::pow((z - z_c)/z_r, 2));if(L<=one) { Real dT=T_pert *(std::cos(PI *L)+one)/two;Real Tbar_hse=p_hse(i, j, k)/(R_d *r_hse(i, j, k));Real theta_perturbed=(Tbar_hse+dT) *std::pow(p_0/p_hse(i, j, k), rdOcp);Real theta_0=(Tbar_hse) *std::pow(p_0/p_hse(i, j, k), rdOcp);if(const_rho) { state_pert(i, j, k, RhoTheta_comp)=r_hse(i, j, k) *(theta_perturbed - theta_0);} else { state_pert(i, j, k, Rho_comp)=getRhoThetagivenP(p_hse(i, j, k))/theta_perturbed - r_hse(i, j, k);} } })</div></div>
<div class="ttc" id="aERF__PBLModels_8H_html_aa037b629b806251024dce3d3e714de7b"><div class="ttname"><a href="ERF__PBLModels_8H.html#aa037b629b806251024dce3d3e714de7b">ComputeQKESourceTerms</a></div><div class="ttdeci">AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real ComputeQKESourceTerms(int i, int j, int k, const amrex::Array4< const amrex::Real > &uvel, const amrex::Array4< const amrex::Real > &vvel, const amrex::Array4< const amrex::Real > &cell_data, const amrex::Array4< const amrex::Real > &cell_prim, const amrex::Array4< const amrex::Real > &K_turb, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Box &domain, amrex::Real pbl_mynn_B1_l, const amrex::Real theta_mean, const MoistureComponentIndices &moisture_indices, bool c_ext_dir_on_zlo, bool c_ext_dir_on_zhi, bool u_ext_dir_on_zlo, bool u_ext_dir_on_zhi, bool v_ext_dir_on_zlo, bool v_ext_dir_on_zhi, const amrex::Real met_h_zeta=one)</div><div class="ttdef"><b>Definition:</b> ERF_PBLModels.H:249</div></div>
<div class="ttc" id="aERF__SetupDiff_8H_html_a7b02b17a96d3d80f109a252dbf71ca41"><div class="ttname"><a href="ERF__SetupDiff_8H.html#a7b02b17a96d3d80f109a252dbf71ca41">l_use_mynn</a></div><div class="ttdeci">bool l_use_mynn</div><div class="ttdef"><b>Definition:</b> ERF_SetupDiff.H:13</div></div>
<div class="ttc" id="aERF__SetupVertDiff_8H_html_a4a8c0178282c029acc3d37ed429f9f58"><div class="ttname"><a href="ERF__SetupVertDiff_8H.html#a4a8c0178282c029acc3d37ed429f9f58">turbChoice</a></div><div class="ttdeci">TurbChoice turbChoice</div><div class="ttdef"><b>Definition:</b> ERF_SetupVertDiff.H:6</div></div>
<div class="ttc" id="anamespaceBCVars_html_a00367571247ba920989360b94c156f4aa98ad22810fdfe9e585bc6481fc12454e"><div class="ttname"><a href="namespaceBCVars.html#a00367571247ba920989360b94c156f4aa98ad22810fdfe9e585bc6481fc12454e">BCVars::yvel_bc</a></div><div class="ttdeci">@ yvel_bc</div><div class="ttdef"><b>Definition:</b> ERF_IndexDefines.H:88</div></div>
<div class="ttc" id="anamespaceBCVars_html_a00367571247ba920989360b94c156f4aaaf601346b01655b32acf66d07b856aec"><div class="ttname"><a href="namespaceBCVars.html#a00367571247ba920989360b94c156f4aaaf601346b01655b32acf66d07b856aec">BCVars::cons_bc</a></div><div class="ttdeci">@ cons_bc</div><div class="ttdef"><b>Definition:</b> ERF_IndexDefines.H:76</div></div>
<div class="ttc" id="anamespaceBCVars_html_a00367571247ba920989360b94c156f4aab1782753c68aaad8632b0a7d575e57f4"><div class="ttname"><a href="namespaceBCVars.html#a00367571247ba920989360b94c156f4aab1782753c68aaad8632b0a7d575e57f4">BCVars::xvel_bc</a></div><div class="ttdeci">@ xvel_bc</div><div class="ttdef"><b>Definition:</b> ERF_IndexDefines.H:87</div></div>
<div class="ttc" id="anamespaceERFBCType_html_a40fbdd5a379a92debb4bcaaccac6e01ca4b16988b22c277c9f812ddb5fb1539af"><div class="ttname"><a href="namespaceERFBCType.html#a40fbdd5a379a92debb4bcaaccac6e01ca4b16988b22c277c9f812ddb5fb1539af">ERFBCType::ext_dir</a></div><div class="ttdeci">@ ext_dir</div><div class="ttdef"><b>Definition:</b> ERF_IndexDefines.H:227</div></div>
<div class="ttc" id="astructMYNNLevel25_html_a8a1dab3e8da3715b39e07436759b3e80"><div class="ttname"><a href="structMYNNLevel25.html#a8a1dab3e8da3715b39e07436759b3e80">MYNNLevel25::B1</a></div><div class="ttdeci">amrex::Real B1</div><div class="ttdef"><b>Definition:</b> ERF_MYNNStruct.H:45</div></div>
<div class="ttc" id="astructTurbChoice_html_a11fd9429567c4dd4afa871170580aed0"><div class="ttname"><a href="structTurbChoice.html#a11fd9429567c4dd4afa871170580aed0">TurbChoice::pbl_mynn</a></div><div class="ttdeci">MYNNLevel25 pbl_mynn</div><div class="ttdef"><b>Definition:</b> ERF_TurbStruct.H:420</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__AddQKESources_8H.html">ERF_AddQKESources.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>