void ComputeTurbulentViscosityLES(const MultiFab &Tau11, const MultiFab &Tau22, const MultiFab &Tau33, const MultiFab &Tau12, const MultiFab &Tau13, const MultiFab &Tau23, const MultiFab &cons_in, MultiFab &eddyViscosity, MultiFab &Hfx1, MultiFab &Hfx2, MultiFab &Hfx3, MultiFab &Diss, const Geometry &geom, const MultiFab &mapfac_u, const MultiFab &mapfac_v, const std::unique_ptr< MultiFab > &z_phys_nd, const TurbChoice &turbChoice, const Real const_grav, std::unique_ptr< ABLMost > &most, const bool &exp_most)
Definition: ERF_ComputeTurbulentViscosity.cpp:32
-
-
-
-
-
+
+
+
+
+
Definition: ERF_TurbStruct.H:29
347 amrex::Print() << "YSU PBL using ERF coriolis frequency: " << turbChoice[lev].pbl_ysu_coriolis_freq << std::endl;
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- 394 amrex::Warning("\n*** WARNING: Moisture may not yet be compatible with PBL models, \n proceed with caution ***");
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ 349 amrex::Print() << "YSU PBL using ERF coriolis frequency: " << turbChoice[lev].pbl_ysu_coriolis_freq << std::endl;
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ 396 amrex::Warning("\n*** WARNING: Moisture may not yet be compatible with PBL models, \n proceed with caution ***");
+
+
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
- 424 if (windfarm_loc_type == WindFarmLocType::lat_lon and (windfarm_x_shift < 0.0 or windfarm_y_shift < 0.0)) {
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ 426 if (windfarm_loc_type == WindFarmLocType::lat_lon and (windfarm_x_shift < 0.0 or windfarm_y_shift < 0.0)) {
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- 525 amrex::Print() << "Coriolis frequency, f = " << coriolis_factor * sinphi << " 1/s" << std::endl;
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- 540 amrex::Print() << "NOTE: abl_geo_wind_table provided, ignoring input abl_geo_wind" << std::endl;
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ 526 amrex::Print() << "Coriolis frequency, f = " << coriolis_factor * sinphi << " 1/s" << std::endl;
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ 541 amrex::Print() << "NOTE: abl_geo_wind_table provided, ignoring input abl_geo_wind" << std::endl;
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
@@ -833,93 +834,93 @@
Definition: ERF_DiffStruct.H:19
Definition: ERF_DataStruct.H:82
void read_int_string(int max_level, const char *string_to_read, amrex::Vector< int > &vec_to_fill, int default_int)
Definition: ERF_DataStruct.H:545
static void set_flat_terrain_flag()
Definition: ERF_DataStruct.H:572
int force_stage1_single_substep
Definition: ERF_DataStruct.H:585
std::string windfarm_spec_table_extra
Definition: ERF_DataStruct.H:688
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcing
Definition: ERF_DataStruct.H:670
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:641
bool custom_geostrophic_profile
Definition: ERF_DataStruct.H:636
amrex::Real grid_stretching_ratio
Definition: ERF_DataStruct.H:621
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_pressure_grad
Definition: ERF_DataStruct.H:669
amrex::Vector< SubsteppingType > substepping_type
Definition: ERF_DataStruct.H:587
void read_int_string(int max_level, const char *string_to_read, amrex::Vector< int > &vec_to_fill, int default_int)
Definition: ERF_DataStruct.H:546
static void set_flat_terrain_flag()
Definition: ERF_DataStruct.H:573
int force_stage1_single_substep
Definition: ERF_DataStruct.H:586
std::string windfarm_spec_table_extra
Definition: ERF_DataStruct.H:689
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcing
Definition: ERF_DataStruct.H:671
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:642
bool custom_geostrophic_profile
Definition: ERF_DataStruct.H:637
amrex::Real grid_stretching_ratio
Definition: ERF_DataStruct.H:622
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_pressure_grad
Definition: ERF_DataStruct.H:670
amrex::Vector< SubsteppingType > substepping_type
Definition: ERF_DataStruct.H:588
amrex::Real sampling_distance_by_D
Definition: ERF_DataStruct.H:690
std::string windfarm_airfoil_tables
Definition: ERF_DataStruct.H:689
amrex::Real sampling_distance_by_D
Definition: ERF_DataStruct.H:691
std::string windfarm_airfoil_tables
Definition: ERF_DataStruct.H:690
Definition: ERF_SpongeStruct.H:15
67 amrex::Print() << "Not using any molecular diffusivity, i.e. using the modeled turbulent diffusivity"
-
-
-
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+ 75 amrex::Print() << "Not using any molecular diffusivity, i.e. using the modeled turbulent diffusivity"
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
@ None
@ ConstantAlpha
@ Constant
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real pp(amrex::Real y)
Definition: ERF_MicrophysicsUtils.H:219
Definition: ERF_DiffStruct.H:19
@ Theta_hDefinition: ERF_IndexDefines.H:152
Definition: ERF_DiffStruct.H:19
-
-
-
-
-
+
+
+
+
+
-
-
-
-
-
+
+
+
+
+
Definition: ERF_TurbStruct.H:29
diff --git a/ERF__DiffusionSrcForState__T_8cpp.html b/ERF__DiffusionSrcForState__T_8cpp.html
index d20cb14e71..153a272449 100644
--- a/ERF__DiffusionSrcForState__T_8cpp.html
+++ b/ERF__DiffusionSrcForState__T_8cpp.html
@@ -1241,17 +1241,17 @@ @ Theta_hDefinition: ERF_IndexDefines.H:152
Definition: ERF_DiffStruct.H:19
-
-
-
-
-
+
+
+
+
+
-
-
-
-
-
+
+
+
+
+
Definition: ERF_TurbStruct.H:29
diff --git a/ERF__Kessler_8H_source.html b/ERF__Kessler_8H_source.html
index 768207ef2d..04efb42283 100644
--- a/ERF__Kessler_8H_source.html
+++ b/ERF__Kessler_8H_source.html
@@ -327,11 +327,11 @@
Definition: ERF_DataStruct.H:82
-
-
-
-
-
+
+
+
+
+
diff --git a/ERF__MakeBuoyancy_8cpp.html b/ERF__MakeBuoyancy_8cpp.html
index b6289ef31e..9cde7d25be 100644
--- a/ERF__MakeBuoyancy_8cpp.html
+++ b/ERF__MakeBuoyancy_8cpp.html
@@ -386,10 +386,10 @@ @ th0_compDefinition: ERF_IndexDefines.H:66
-
-
-
-
+
+
+
+
Referenced by make_mom_sources().
diff --git a/ERF__MakeMomSources_8cpp.html b/ERF__MakeMomSources_8cpp.html
index c87efd5709..46bdd02c95 100644
--- a/ERF__MakeMomSources_8cpp.html
+++ b/ERF__MakeMomSources_8cpp.html
@@ -399,7 +399,7 @@ 90 // 9. Forest canopy
-
+
@@ -788,112 +788,114 @@ 479 // *****************************************************************************
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
void ApplySpongeZoneBCsForMom(const SpongeChoice &spongeChoice, const Geometry geom, const Box &tbx, const Box &tby, const Box &tbz, const Array4< Real > &rho_u_rhs, const Array4< Real > &rho_v_rhs, const Array4< Real > &rho_w_rhs, const Array4< const Real > &rho_u, const Array4< const Real > &rho_v, const Array4< const Real > &rho_w)Definition: ERF_ApplySpongeZoneBCs.cpp:117
void ApplySpongeZoneBCsForMom_ReadFromFile(const SpongeChoice &spongeChoice, const Geometry geom, const Box &tbx, const Box &tby, const Array4< const Real > &cell_data, const Array4< Real > &rho_u_rhs, const Array4< Real > &rho_v_rhs, const Array4< const Real > &rho_u, const Array4< const Real > &rho_v, const Vector< Real * > d_sponge_ptrs_at_lev)Definition: ERF_ApplySpongeZoneBCs_ReadFromFile.cpp:8
@@ -903,6 +905,8 @@ DirectionSelector< 2 > ZDirDefinition: ERF_DirectionSelector.H:38
void make_buoyancy(Vector< MultiFab > &S_data, const MultiFab &S_prim, MultiFab &buoyancy, const amrex::Geometry geom, const SolverChoice &solverChoice, const MultiFab &base_state, const int n_qstate, const int anelastic)Definition: ERF_MakeBuoyancy.cpp:32
+void NumericalDiffusion_Xmom(const Box &bx, const Real dt, const Real num_diff_coeff, const Array4< const Real > &prim_data, const Array4< const Real > &cell_data, const Array4< Real > &rhs, const Array4< const Real > &mf_arr)Definition: ERF_NumericalDiffusion.cpp:84
+void NumericalDiffusion_Ymom(const Box &bx, const Real dt, const Real num_diff_coeff, const Array4< const Real > &prim_data, const Array4< const Real > &cell_data, const Array4< Real > &rhs, const Array4< const Real > &mf_arr)Definition: ERF_NumericalDiffusion.cpp:149
AMREX_FORCE_INLINE IntVect offset(const int face_dir, const int normal)Definition: ERF_ReadBndryPlanes.cpp:28
Definition: ERF_PlaneAverage.H:14
@@ -915,35 +919,35 @@ amrex::Vector< amrex::Real > input_sounding_timeDefinition: ERF_InputSoundingData.H:315
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > U_inp_sound_dDefinition: ERF_InputSoundingData.H:326
-
-
-
-
-
-amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcingDefinition: ERF_DataStruct.H:670
-
-bool nudging_from_input_soundingDefinition: ERF_DataStruct.H:641
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcingDefinition: ERF_DataStruct.H:671
+
+bool nudging_from_input_soundingDefinition: ERF_DataStruct.H:642
+
+
+
+
+
+
+
+
+
+
+
+
+
+
Here is the call graph for this function:
diff --git a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map
index 7e27f65121..45a1190122 100644
--- a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map
+++ b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map
@@ -1,5 +1,5 @@
diff --git a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5 b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5
index 1d45be5335..555b509dc1 100644
--- a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5
+++ b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5
@@ -1 +1 @@
-c219d0622a1b348eda802510569f41c0
\ No newline at end of file
+ec241a01c30958ea7c35f3d4a480678a
\ No newline at end of file
diff --git a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.png b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.png
index bf33e240600cd06a8fc2cb5a564193e80888c5fe..8b2c87b36faf744333512464ce28e96b1d89ba3b 100644
GIT binary patch
literal 70087
zcmce;by(Hg_C35+#4X(-(%q>@gPU#y=`ImLLZnN&K|n;L4kaiZN|%6uhzN*Cr*wDo
zo1Am|-tYV8%e~KYd|>m5wdR_0%rV9!Oie`=8-ol3fk0r(%Sovt5LYY^h|A>2%kUkN
z$fFSW54xG6tQ6w>;(veC=f)!t)ChU0dk@@FR>s_QG_;QdHk>3SsA-c!j*>T1T@-{LOOLAI>001lrY-zKJ(1P;E>T7J2Q3#l+ROKJ
zX(L{YPEfh}ums>q&?Npo6*-+5kB{H4W9l?FUVF+JJR@-NJP5=ecZtdf|NIw%?gmlB
zU*A(zWB>D~&&BSQBLDo>kradQCGLjNde8mI<^T(;OHIb?92@}w0iECAlQNtFWk~w)
z?uY53p08He5eLO5CnvS8>v5@95s1h~4i@+?sOjheR)`Q1NSYfrZj@;t5YJ{u-&t*xz7lzCNEzdAeRH8nML
zb(1(XJ-e8jq!TV5Kk!{5>A6q*OvKhWo6|=4t;*}&yg7v
z%*@Ql_RYP$JwB^%`T1Wx>gn$O@bP1skn>Wiu&bNMcoLu0n>TOD%FBp=AF6$HZMwL~dWCePUC&~myptUZRRI4|$ubbms6m64IL
z)@8MpP>mn@bXbpdl4NS9Tyi@U+)bUiAjNV&dp^)F$)SNp)%^)zn|jV
z%TX6UUMUxJn!}^wjh=y5$VpG1@VjzN?CikY-Tkof{B&LXWFsIypAEIp5liRc;c=U+
z`^%T5sg@A41~DG@;NakI>Eg!s?$K~?VAE7QTU4ivaXmltrXBh5*$~!=s}iP|71Cl_VsZI^&sY8;OXC
zcgI{O90kbl-n}aNZ6cn181Mf6)eU4Wv4-XGl
zS64;*7Vf;w&$rRg=+Bak!N9~EvfX>HnJ;xpNJzN9Hg>%AThZQrT{B;6N+25Vsgsk&
zg9q$vY+VqsDk>z-5D}a#EY*%PGS1GsJB$6PDJdQ=1x70Ex!Bo7e9uIpWa#wJg`YJ_
zK67^dJwJbgh)73QmnJYx)Kg1WcX4&~bnULXfkB<~^3U1XSy&z_-xHoYcRm!at9Ye3
zIy%B?N}8Kz)z%85LQxq9$H!FSKIAknU%s@kw9F|i6m3$|)zuBbq3}4|$_WZ;gIGut
z_ceO(;DL_L>(tcLq$CzbQCPEp=U1rjv6I+PHwtkk=I7+>?yc&Rx4Vm^)zo;{*syK(
z73h|`tc|iP2|}VH(>r*sJUPU;w9iwA3F`>?erH`hu|uA^4NEwTLKYg@wOGNP&Ksw_#f9MP>?T@TAG
zEiAaPH16Mb-<-aB-;!WwcemEwNq~HFmG;(MI}y)4TS-Yt85#5h1;`lG!*_1qj;0XS
zQdj>pvne7fYG_~(`sUr>Yg_(gg&cTbf+*NiBuyr}B=X(+_XM||K7ESJ`1p~ZGf^?O
zU-Q8O543{3JnCS^HO$z!xW{&OMFj=)M)fYMdL||&PoI{(d-smnuNjBK2RV(B^pjB2
z()#?c{Pg4m4;9a(DuZ;a9(8RX%+AeCeDkK}$;ii{a;xLRLs6d-ogMXk8xlxfjAA{M
z*KoqVCMmsbY&=iCYwz(gZCVn4h&B)c==bm6((T5ucWtAw7l}#j4;?hRXfAq)+D=XwABp+>u%GNT_g)3guNHG(IH#9Vagve@Xk)R%Zj=Tw%9P3GNou8b1k(}%a
zIpfx?TafU7{`_fF8`eYEQsO#M<~Z~D>({T?xA=G5oJ%y&5eV+ti&w&N6fTU7jm>Kt
z9zKHjY-N!PRD&a8WMrhM*ag{3fIK=rKK@{Djo)sx%J=+C*mbSIu=drXsi)OAI5-fq
zo~K978X6kUo;`z8S6%IH%(J_@TUc0lTvK{XDmC!@0o1
z!YU{(=HlX-sCHuNEh;QrU07I{nb9yX7=<`hQ=6Qf{n6XoJhS|}rY9A7I~yBYe_!9j?l(brpSV|A!w5@DOKE6m7#JAfuBz!`*x1<0s;Z4oe}5@5
ztW{7{6c!e~dikEdzP|C@SS)_kz!#80qi3kZEaB`+;)F|cPfN+l@+Ey;%+m
zyz}hj^b`VlU|>K)UHzhd2;CUc_HCVvzb
zQb0zeqH1)QLh714jT*CRc-QC?#dwNm?B`ZdT
zhv$C(KHY4mA|WAxI!{SS86F-E8?rsu4bN=GoAl+&7kqsD#l=N4ElG0Qv}y7B3=9A)
zlcFX9(ebHQNM-x_`^QE{#b&jw%*@K1m$jf6T>7~Gy|*``sEG3BP5IIAV;^(x(1;sc
z9`5d8KYcIF(mvjgjKnGjS+W6-*X>HW-eJ+tNhrO!9pIZ`>S>dae6)6~@6Vk1k8pw+kz;z;Lj-@aMo!u{#UM4GQgBE?U4O5VTs>-N#q
z9a1yw_B05^@L(<9XK|6uQBZy(CG5#B
zp{aQ|9Sux>Qrc
zYVjMg1tQW4fb*8-gCDik(oe=yjt=Mf(EJ;5IMv?~kty@D+4TQ-ry-|=3!fzIfy835
zv_d2$W29eVa=tuhUsm229lcjto-CcCKogjq&CJW&x43AmlDvmMWPJh0XKoafBjq4;
zYfUB1EE%{2+%|);G2w&A4Ds`V!;@_)eALsYyt=x^-6_{q2jMpm9cZ#7mjvVD*4T9k
zRB}~2x-hRJ4^r}o@KOD_`9niDyE^4_bA7KTT%c|Vp$iaSSY7Si9=)?YPjz^xyE-&;
zu=Z5Y_e{5cbd=vJT8w4!Vnr{4!cR)nX2Qv2q{#%OVaf{Wv!Eb7y67@@U=6&n2B^Gg|VsR+{J
z@w!DMmeS$wpw%_pn%e4a8exiy%VJ%IU3AmbymN5$`RBl$+_!cm6}PGX_A{K6VDy%S
zv(u;7LL`YWUM1IFWtE2CF}dw;hL4pvJ6p@i*$#lQygcCcorJ$V=lT^D$pBM)$}~U?
zzb7Z#nwu55vFQG~U>y@PNi!ZHGB6_}iG}6Rb6==fZ4+`$j4allgOI
z=gxel11sxK*~|JbR;DTB;*Lpo?S+08g!%5|y=gn@8zrawd9gB)SK;7h+}9qQUcA{o
z^C3leKOr_&=&y^KJi53j#B4x_ED?&CjZIBkJ80$E#r?P~AX2E_S5$m!V^CLCW?``&
z%a(t(nJP&B4Y~&Ucv8fbjGP=xYioHI68Q8}=-9;MnYAdVxiyT#c(Su=zOK+T=9n1vfFK+^>rDrTbMO9rBi?|stu8M3tccHv8{pja#&gcccX@b*KG-Beu@7kGzH_JE`}pI5
zOKe_V9y~+1i=u);-h;q09-0$
z`Q_z6zILO>uEGyr=y{;5W1*W5S6C-<88Z%78x%pG40qt=74q0&4mK5`#HDGS-wyEi
z2SfzD$ja(!{`>bHv+P6|qN1YssLsw#B8=Zq*IkjU2}eU#>D@_uAqOuU$&hDrUd~
z8#EW|`N+m|d6i;&Y2X&ga#v?OF+IQCvu8VV-N}IWFJHc_k|r{AngdI>-f$eo4dQFwsr(Y1B8Th1j*vULf-rLDzK>|_Kn>=J@(KY#>y1)7#bQTzj;Fm
zX;nw3<6*gWp&r9h`s>$KuIuBlj2@nz#d;sYb$mv2P@TU9cwzkua7Tl<)sI|2ENNbjOwV5LPA0aS6eS}*
zus|Dho2<|?IFK>}xR=)Xn_PluU!Zmk;)9T`D-fssD@Car!yT90Vv`^2UXJ5l*)zG-Blqw`0
zC)xh_Gd%6hn>VpqE?ST}w|ftg(7e4**G9;3WW8!TAPjoaL`S}Ve-RrCbPEMJd8i9E
zc9eO0bi&J*8*U=1s;U64GaoM9VNC9RtNsq!BB`v53>s7@a~B>8*rC?oloWDQUSVNW
z*YCc?XT-z_l$Y38Swl0zj#1g!*>94Qp}7dthIJz$C4HTeV*B_ptAYtM-BkWn-_@w`
zUM3`z4%v!}iwAz6?Tp7q8F22WuA>3S&UyRRZsg;Wy_F#!AD|6?y?smHPK*(htCH$*wDT_Stm^4+g?+wWDQqcbF0R2JKZqtc%`7Y|?CjEByuc$S
z{yj58$1&!VDu*6`kB9dMFik@Gr*=vl9EH$+0@IUOT2Vnnlbw?@F!=2D?b`=lVve!2
z_<_GJKtFW0l*Px!_EM)u^Xa}q!ondzK^*!Y3q13n#eRsy#|X+SE?%FSicLre!h6GG
z9$Zvpx9{|OOi(}ovh{eQ?*s(G@Nf;p#h~5vb}RSFSf%Y~6%QYuQ%f+ms1e1=OGw{X
zg7Av#Dyc_}RaMDTQ)?t7=meB+s#1S;bjbZMd-CKIc*$uZdK82yD>NiIIlr6yu8`xj
z1EJ*Q=?WafXru2rK!cw}hKdyq)dttDh1d)h)%%>*ynDwYBveyY);Bs@SN!_u=m?tM
zOP4M|y0hi??0>7Q#HgXHTv)tFN=Bxqr-zP?{&R3JF)v5xJAUnWmEPvRILYk7g08wEtFyBcH=C@N56L%&Jy4UDJ8
z#}fvjbW2PwIid(@vNJMDO^($=s
zk>T~YHnThY{A8r0?1nWz+S-C%DdXr$mo-f
zfZW{NT7c{w#ot<6EtY>40x5CDc0D62OVZF#l$STIig^F@^AR2}YCBt71cm39JEVWC
zLHe*@0w_gyL!=b@Bqot7*~`o9d=nP{fI-JM%Y`Wc=x$ki1^|pUQ|Wjscm5qQB_vv4jE~8jd$i71Xtnx7#MP@s;XS=pWM`fUBE&j
zbt>&BqGa?NyacONLoaYL`c=T4#)gKUoR_DK>!;^uy}C_0VpCPpq5R#>+gn4gP5B~W
zhSQoYA0L{V0Kb69`bH|^fi~W=Ha3re!Ad|lX#I#eJ{41f%ea9^l-Pm3PsLM
z&;+V$Y!ec$lM7C>Ofe8Q%b{PTh3h~EEP)jD()r+I^-%3Qt+23pnbS>g+b2)La-Pio
z{tdeiwISFZ=?{@;5!}(!bE8lO`e2h$e>2Cr{9=3EO3voKK9jCj45$0(^)!?&?K2mp
za~_!3F<>WO$onVYx<$X^S;*wF?Ua>gpTj2h_k`z^Nk*4okZY
zBEf&4pimc2Pfx5Cps&)>jDUoFG|>6x4HTT3u7q`O8a|eC=YA!kI^AW(4Hqww9I`0Ua3$qQ?zN$|7I`^YRwf*Cz()bdcVu
z!D}`~up*!7^Kd_Z@I7VA&kH5IwJaQG?@P%6uoH+n9;{?hPyXKN)m1}v^`|qRqkxeF
zLPl6n5Wh`VMFpr0IJNVDO}&qvK^f?s0BoYA^|7MjHV{fQG#L~4Gt<*P2U8)RKg+nf
zx@v0=0_p#lU|oFl<_4@~>Gf+fnldteHAizm_bJB(l<5@k
z?B(m1%@+6X11BX#6FMW1+h;-O(hOAd
z_V$jU!9zh+5s2=2b02YNZe!;HP7V$L0d_YwewACPDl12+F-L}mE-Wpfkk*eMH@I)V
zEh>r*g4--2=h@D&2FmDBNsx}u_uz}un0x?poy5(FKb#u*5e=4R{ERR1j5
zdS;@?m1m!cAWHyi>2rR%Z$Hs6@PCvS$k2os7nu^42ILl?>smdrv9QEIMJjGMdJhyY
z@MsW$SE#g&jVFeOg{Z0jc>6ZE`=R}CMOdfB`*-iK+sTjL3_3&tkT;=!Zss@|7FF+l
zxYbW5E5f9b9h=%kN<>Nu0wEa{RScejl2QsNhVh*bG&SKiQ0|hrP5q0C*3M3kMMXq_
zrv{Y|xH#!`%TgULv5BUprt$Idkr7X|nKB)Ap;S9tTU$W985yy!UkkJnF25ag$DVM2
zMd%2cI`#Rs3B4@e;bG=GBWN2MW^E4?#Fgo7bbcqcBy|pF;bTe%Ic}%yO4!s2Y9~7!yA;
z$Cc%Ug`=f>3k{9b;zlE>zuney^7hWJ{nP!SJ9nbw<0iH+
zWupsz{CL194CfWt4lyU~1Dl5bkf
z(?S^9?&XJL7`eN9!~b=44Xv%8C@W`4Jf{n^u*ha6Qt7opI+@xB-?hqqNw+6Oi@V4RQo2Gu_r>Ft!BGpiY268`!bSxz={S=kS#`O8+{0`;|;
zf4T0j8A2LyAu(Y;*y&fDQt{<}Flo0?%1A$35RM@#;C6i|3=49?;R2mBWpPfq#`i+}
z__$1Zrs5C<{7AxUxRju0u@C7!FMZx-R$WsAdOs;0gr^2CFGi7p`_GDQCK1NsbaMHQZ(%nRb;ZbP1oO;r@v)1L?qv;ej$IeB@a3cs&j
zo!;woH}l1doU^lSIeFaY_nTW==|Lu6VV78EXJ-d)KpMzAECyzhcK?zRTa2JGtDYMc
zOnou)J;uXB(iT?YVve7V4>mOw6t>p)l6a*Z``-q|^ZVV?~Pbu1+x!ThsMX=FL`DiSXcKoM=^QIObXG2_ob)DxYnicHBPxzkBXe!iw8pM
z?_}dnRWJ}?el{esl|`IWuzF7OJ|b~l#d~|-L_N#Y)wfPh3iehUU%py9mFKv4Z&JU>
zj_Cyw!gy^xJ;)JVohsR?IXTXXi83kT=Odl*`Hjclt^+3<+zR@%DF}up+)(Hs?DWtxS}hOi*OZa1E{6oBmlJsdT!tCxv1BZ^}IYh*Ku$@mX&4p
zFM^Q)I4ewCT<8{|r{3A$zkv}XRcXZJf3_dF8cQd;v$F$5I3ZyGG$wo$Fd-248@s!0
z++RO^$}A>F7q$^gry?TC&dSQl&bF|(pIBZw3Kx7qU#q$_rdn03i{NH=;C}2_U5t(G
z8W=JkmnxU;bF`S@*X-hXEJ^;2t<33?36TSkB`z*&l_MafLl0eS+&KQ@M+G<-&;#yk
zYlqBKftUtzYN%MN`4!-;qGk@DX8?7(%D@g99{;PD@NWllC<>tyWa;?fduIsoVnj*eKFnG5tQi-6atsyfJu5mQiAWvbZ~!Y6P4i0O@;7}DB`@Qqw^6kpzqoCC@}CAxmTJN
zxHrrk;o&pKZt^_@m2$?LDB>S>g44g-2_^+7@$e>p{|@?-VuiMXmyd6PVi_7VJ3G5f
z#|y8)lj%IlkQSQgJRnMT+-4;FUDJCnn3~o9U~2Hl+uM|+9%v^dCW1lW0%dMqGh)Ya
z_}SXt_4{dT@e$f6jBp8&WL
z$Hx~Htp~mjW{RLk
z4{y8wuD;%Bu;?pS!+`~w7*w0KgW=s-x>TgT5fdDV}*C&AfYBu2pPQA8$QSs^1Es!h{*>!;zI|m=lRZNl{
z;*%?3Nk~&LUnwXk0F9!nD@aRAD=M0*l1mr((SE{bv5!Geu)1L_*e6DpiIGtyT`U!V
z{QLL7M;#p;D81egSJhu-oAaa}8f3zHZeen8-WX36NwWNtLB#t@y`yN00^rsq#4h$8ra`3psY+8{;lL&-7
z8Ozmcq)0!BDl64jA9gEx#P8C=n!hi6VpM-blzhr6dfFxEv~%ox1!2od!(h4>;I$`LTNM``e-SIi$
zNkla_oO$pJg;587FVt^qXmCvTc}z-Cuczg*plT9kWo><3HV$zq9AsBgyOL+$7C<0K
ze(mN7jV8Ke^ooI+)=x>*o_M9XufKG}JGH|vK*I}qc(;)C#-1dhT{_&DacyLL_Q8R^
zpI_(rcshkpMVvy=#D!MQ{khqC;r`@!{Z?0^ud2MRYQ_sMds#QPwdR0C*Ht|-ybZzX
zBm}}^a`OO`xJexDKh2ts*)5JI$LmgF(C#+-FU~&P8grdjP5TOf&zYWl(n}l)8}OnO
zaZE9nxRt)8i3&IQ$Dop^IYg!PJi;q>Wc&KwKY0TAc`HO%6wxGr);2Adn0k899PmXk
zkw5r2*GVN;RYwsU|-);
zKx>p-AbXjU%4DbkIQ#VJlN7EqH;9?_o|REi*a6M2`JW7f9u?;D@4(DhQ7nV^e
zMPTz_a|ZOGMzCw47JikQr+aOMWt^ST|lMz3AH3jI6_8yl3Av+m@($PB1)0NNmXsU-1Q
z?k)`!8`Ts1hA8>Js5veo4CyhK)!_?M-9IpLh0T9pY;3qbvI`<{6L
z7Tt85{*x5Emic0dO%zrQ)M#L!pmu><;r=(XPkWCC-fIG4tplSN>_<}45V$Y#QBd?j
z{End%Wo2i#C?A3cA=vHt@&%N~?*Q~U*x5ms-T}Q65E%pnpy97ymF`#mw<5crAfcxR
zVGF2KBzgu&rFEX5O_Oz_8>KOCTW~1!tj#hxO1;2tLxdxfT~qV$O2BGt?4aCDLuMko
ziL|1^`rk_fx!Ls8=%~;(v#|AYi%x>++t0xL3no4g@r@e1!e&7JbyrgQT6+qwu`)L&
zV+$-Jn}|rtFGid~Fc|YEN8wFBt(K~lpf=$q1vnrO?Nn_58WdQzrv<*vZ2lM+fQs%&
z?&|D3R$)yTJ%h?nla&Rz_h!2zk0ZDXq{;6Q#T$X88k&Z7c)`!t0VTleSH{O8pk#5zGP&{`|FKzE-s)|UbMYx?a_uB
zYA<$oQ`rtP9LaI9yOk7cU6)?2&!^v#7xSZW1eZKk3%uAm#8SVixN;J2j%>)kp;x?+
z7X13`_eT@C_F#{u34|^JG%h?885tQMP*O$)hJn#hr4+5`s3^Xd!9Zx&0ax2*
zMic1e?G4oZc!T#b6f$ria2hv|)6g^GpyXgFy*vy!6IF8a^YZY#!I2w=y=giD;{q%U
z43;J)-|J($A8NdNFulQq3KcL_{^5YkfX0FF7K9~`^|#5CD@K0)WW=CV2>Z2m_AdPH
zFP2Cu7ngEUtZK?iL5+7WffM*>iD|9PcLnj87>ZgZFfiHK!pZy`FXQQ^GIycI$t%=&
z?&J+-O8ds2EJy|kly%H+!`BmATaP**+n1-NO5vyg)ei;|e0-4V85rOElaK9t>fnI)
zcagJrrWjx$%{M6!wRN<#_`wH2L(@DqmYOGg94;A2%ZfZwP?eMn=$R+H!9p7d=({a<
zaPS%6Gmy6(neR&}D(1siP0gEM92^`BWY@N=Bk=B|GmFCU3VZ?QFz=5QpIrQ)t-(g|
zfYeU{j(fHl$s=K^V!ihkipzT+uu2;n-{E~qyqLcO?|c2?eFLs>slJig->Fz>Tvi*Y
z-A4i91cH@t434L`etc;{Sio0(8hq78D0oX=z@ZR|A*H3M2*V!zl^+us19Ib
z#KaoFGyDGiV=!rEegh;Ar&NkUq)GMR$-)wW&$;)vx5_MRTCe}&u)uUc6Ua=Ybno8Y
z*Q8yCV?Kn085A|gUn>CJ_x%HIbalvi5Py{n;iM9FS)IJ+%EY8`9hbut$^)2D_BS>l
z4Kz15w+fnV5)cxCI{M9m7mPPou3Q0leokAkAxg16Q+J~sBk$@}yB=I(BHS66?3
ze=@j*4g{{K#{7#^`BJAIvhc&_rnQPHDrqon0p0}&dMXjO5|if3v;zNPyp>Z@TG=;W
zY@V19zkT}wNQMjyzdr?GwstOqLb~ZD0_FvF;-VrcmuN^1Xf18yHrMF6m=}~xO)ZU#
zkO{24>S8ogE{`6~jn{w6*NR$jKIg$QFV%qzZ^49W_?w-Uq>3grEUXj6wfed`dclQX
z$D(y^y$g#s!IGO6H!bdcWGw4)!C+@)ZAVjVlWqHUT|&r^x%qc>s2drjC+tpcn5tga
z4hIA327o23>!r7*5@=0Ibsn%fZ%(&CPbciM;^6AaBuvNsGEFppb0<
Definition: ERF_DiffStruct.H:19
Definition: ERF_TurbStruct.H:29
@ Theta_h
Definition: ERF_IndexDefines.H:152
Definition: ERF_DiffStruct.H:19
Definition: ERF_TurbStruct.H:29
Definition: ERF_DataStruct.H:82
@ th0_compDefinition: ERF_IndexDefines.H:66
-
-
-
-
+
+
+
+
Referenced by make_mom_sources().
diff --git a/ERF__MakeMomSources_8cpp.html b/ERF__MakeMomSources_8cpp.html
index c87efd5709..46bdd02c95 100644
--- a/ERF__MakeMomSources_8cpp.html
+++ b/ERF__MakeMomSources_8cpp.html
@@ -399,7 +399,7 @@ 90 // 9. Forest canopy
-
+
@@ -788,112 +788,114 @@ 479 // *****************************************************************************
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
void ApplySpongeZoneBCsForMom(const SpongeChoice &spongeChoice, const Geometry geom, const Box &tbx, const Box &tby, const Box &tbz, const Array4< Real > &rho_u_rhs, const Array4< Real > &rho_v_rhs, const Array4< Real > &rho_w_rhs, const Array4< const Real > &rho_u, const Array4< const Real > &rho_v, const Array4< const Real > &rho_w)Definition: ERF_ApplySpongeZoneBCs.cpp:117
void ApplySpongeZoneBCsForMom_ReadFromFile(const SpongeChoice &spongeChoice, const Geometry geom, const Box &tbx, const Box &tby, const Array4< const Real > &cell_data, const Array4< Real > &rho_u_rhs, const Array4< Real > &rho_v_rhs, const Array4< const Real > &rho_u, const Array4< const Real > &rho_v, const Vector< Real * > d_sponge_ptrs_at_lev)Definition: ERF_ApplySpongeZoneBCs_ReadFromFile.cpp:8
@@ -903,6 +905,8 @@ DirectionSelector< 2 > ZDirDefinition: ERF_DirectionSelector.H:38
void make_buoyancy(Vector< MultiFab > &S_data, const MultiFab &S_prim, MultiFab &buoyancy, const amrex::Geometry geom, const SolverChoice &solverChoice, const MultiFab &base_state, const int n_qstate, const int anelastic)Definition: ERF_MakeBuoyancy.cpp:32
+void NumericalDiffusion_Xmom(const Box &bx, const Real dt, const Real num_diff_coeff, const Array4< const Real > &prim_data, const Array4< const Real > &cell_data, const Array4< Real > &rhs, const Array4< const Real > &mf_arr)Definition: ERF_NumericalDiffusion.cpp:84
+void NumericalDiffusion_Ymom(const Box &bx, const Real dt, const Real num_diff_coeff, const Array4< const Real > &prim_data, const Array4< const Real > &cell_data, const Array4< Real > &rhs, const Array4< const Real > &mf_arr)Definition: ERF_NumericalDiffusion.cpp:149
AMREX_FORCE_INLINE IntVect offset(const int face_dir, const int normal)Definition: ERF_ReadBndryPlanes.cpp:28
Definition: ERF_PlaneAverage.H:14
@@ -915,35 +919,35 @@ amrex::Vector< amrex::Real > input_sounding_timeDefinition: ERF_InputSoundingData.H:315
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > U_inp_sound_dDefinition: ERF_InputSoundingData.H:326
-
-
-
-
-
-amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcingDefinition: ERF_DataStruct.H:670
-
-bool nudging_from_input_soundingDefinition: ERF_DataStruct.H:641
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcingDefinition: ERF_DataStruct.H:671
+
+bool nudging_from_input_soundingDefinition: ERF_DataStruct.H:642
+
+
+
+
+
+
+
+
+
+
+
+
+
+
Here is the call graph for this function:
diff --git a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map
index 7e27f65121..45a1190122 100644
--- a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map
+++ b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map
@@ -1,5 +1,5 @@
diff --git a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5 b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5
index 1d45be5335..555b509dc1 100644
--- a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5
+++ b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5
@@ -1 +1 @@
-c219d0622a1b348eda802510569f41c0
\ No newline at end of file
+ec241a01c30958ea7c35f3d4a480678a
\ No newline at end of file
diff --git a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.png b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.png
index bf33e240600cd06a8fc2cb5a564193e80888c5fe..8b2c87b36faf744333512464ce28e96b1d89ba3b 100644
GIT binary patch
literal 70087
zcmce;by(Hg_C35+#4X(-(%q>@gPU#y=`ImLLZnN&K|n;L4kaiZN|%6uhzN*Cr*wDo
zo1Am|-tYV8%e~KYd|>m5wdR_0%rV9!Oie`=8-ol3fk0r(%Sovt5LYY^h|A>2%kUkN
z$fFSW54xG6tQ6w>;(veC=f)!t)ChU0dk@@FR>s_QG_;QdHk>3SsA-c!j*>T1T@-{LOOLAI>001lrY-zKJ(1P;E>T7J2Q3#l+ROKJ
zX(L{YPEfh}ums>q&?Npo6*-+5kB{H4W9l?FUVF+JJR@-NJP5=ecZtdf|NIw%?gmlB
zU*A(zWB>D~&&BSQBLDo>kradQCGLjNde8mI<^T(;OHIb?92@}w0iECAlQNtFWk~w)
z?uY53p08He5eLO5CnvS8>v5@95s1h~4i@+?sOjheR)`Q1NSYfrZj@;t5YJ{u-&t*xz7lzCNEzdAeRH8nML
zb(1(XJ-e8jq!TV5Kk!{5>A6q*OvKhWo6|=4t;*}&yg7v
z%*@Ql_RYP$JwB^%`T1Wx>gn$O@bP1skn>Wiu&bNMcoLu0n>TOD%FBp=AF6$HZMwL~dWCePUC&~myptUZRRI4|$ubbms6m64IL
z)@8MpP>mn@bXbpdl4NS9Tyi@U+)bUiAjNV&dp^)F$)SNp)%^)zn|jV
z%TX6UUMUxJn!}^wjh=y5$VpG1@VjzN?CikY-Tkof{B&LXWFsIypAEIp5liRc;c=U+
z`^%T5sg@A41~DG@;NakI>Eg!s?$K~?VAE7QTU4ivaXmltrXBh5*$~!=s}iP|71Cl_VsZI^&sY8;OXC
zcgI{O90kbl-n}aNZ6cn181Mf6)eU4Wv4-XGl
zS64;*7Vf;w&$rRg=+Bak!N9~EvfX>HnJ;xpNJzN9Hg>%AThZQrT{B;6N+25Vsgsk&
zg9q$vY+VqsDk>z-5D}a#EY*%PGS1GsJB$6PDJdQ=1x70Ex!Bo7e9uIpWa#wJg`YJ_
zK67^dJwJbgh)73QmnJYx)Kg1WcX4&~bnULXfkB<~^3U1XSy&z_-xHoYcRm!at9Ye3
zIy%B?N}8Kz)z%85LQxq9$H!FSKIAknU%s@kw9F|i6m3$|)zuBbq3}4|$_WZ;gIGut
z_ceO(;DL_L>(tcLq$CzbQCPEp=U1rjv6I+PHwtkk=I7+>?yc&Rx4Vm^)zo;{*syK(
z73h|`tc|iP2|}VH(>r*sJUPU;w9iwA3F`>?erH`hu|uA^4NEwTLKYg@wOGNP&Ksw_#f9MP>?T@TAG
zEiAaPH16Mb-<-aB-;!WwcemEwNq~HFmG;(MI}y)4TS-Yt85#5h1;`lG!*_1qj;0XS
zQdj>pvne7fYG_~(`sUr>Yg_(gg&cTbf+*NiBuyr}B=X(+_XM||K7ESJ`1p~ZGf^?O
zU-Q8O543{3JnCS^HO$z!xW{&OMFj=)M)fYMdL||&PoI{(d-smnuNjBK2RV(B^pjB2
z()#?c{Pg4m4;9a(DuZ;a9(8RX%+AeCeDkK}$;ii{a;xLRLs6d-ogMXk8xlxfjAA{M
z*KoqVCMmsbY&=iCYwz(gZCVn4h&B)c==bm6((T5ucWtAw7l}#j4;?hRXfAq)+D=XwABp+>u%GNT_g)3guNHG(IH#9Vagve@Xk)R%Zj=Tw%9P3GNou8b1k(}%a
zIpfx?TafU7{`_fF8`eYEQsO#M<~Z~D>({T?xA=G5oJ%y&5eV+ti&w&N6fTU7jm>Kt
z9zKHjY-N!PRD&a8WMrhM*ag{3fIK=rKK@{Djo)sx%J=+C*mbSIu=drXsi)OAI5-fq
zo~K978X6kUo;`z8S6%IH%(J_@TUc0lTvK{XDmC!@0o1
z!YU{(=HlX-sCHuNEh;QrU07I{nb9yX7=<`hQ=6Qf{n6XoJhS|}rY9A7I~yBYe_!9j?l(brpSV|A!w5@DOKE6m7#JAfuBz!`*x1<0s;Z4oe}5@5
ztW{7{6c!e~dikEdzP|C@SS)_kz!#80qi3kZEaB`+;)F|cPfN+l@+Ey;%+m
zyz}hj^b`VlU|>K)UHzhd2;CUc_HCVvzb
zQb0zeqH1)QLh714jT*CRc-QC?#dwNm?B`ZdT
zhv$C(KHY4mA|WAxI!{SS86F-E8?rsu4bN=GoAl+&7kqsD#l=N4ElG0Qv}y7B3=9A)
zlcFX9(ebHQNM-x_`^QE{#b&jw%*@K1m$jf6T>7~Gy|*``sEG3BP5IIAV;^(x(1;sc
z9`5d8KYcIF(mvjgjKnGjS+W6-*X>HW-eJ+tNhrO!9pIZ`>S>dae6)6~@6Vk1k8pw+kz;z;Lj-@aMo!u{#UM4GQgBE?U4O5VTs>-N#q
z9a1yw_B05^@L(<9XK|6uQBZy(CG5#B
zp{aQ|9Sux>Qrc
zYVjMg1tQW4fb*8-gCDik(oe=yjt=Mf(EJ;5IMv?~kty@D+4TQ-ry-|=3!fzIfy835
zv_d2$W29eVa=tuhUsm229lcjto-CcCKogjq&CJW&x43AmlDvmMWPJh0XKoafBjq4;
zYfUB1EE%{2+%|);G2w&A4Ds`V!;@_)eALsYyt=x^-6_{q2jMpm9cZ#7mjvVD*4T9k
zRB}~2x-hRJ4^r}o@KOD_`9niDyE^4_bA7KTT%c|Vp$iaSSY7Si9=)?YPjz^xyE-&;
zu=Z5Y_e{5cbd=vJT8w4!Vnr{4!cR)nX2Qv2q{#%OVaf{Wv!Eb7y67@@U=6&n2B^Gg|VsR+{J
z@w!DMmeS$wpw%_pn%e4a8exiy%VJ%IU3AmbymN5$`RBl$+_!cm6}PGX_A{K6VDy%S
zv(u;7LL`YWUM1IFWtE2CF}dw;hL4pvJ6p@i*$#lQygcCcorJ$V=lT^D$pBM)$}~U?
zzb7Z#nwu55vFQG~U>y@PNi!ZHGB6_}iG}6Rb6==fZ4+`$j4allgOI
z=gxel11sxK*~|JbR;DTB;*Lpo?S+08g!%5|y=gn@8zrawd9gB)SK;7h+}9qQUcA{o
z^C3leKOr_&=&y^KJi53j#B4x_ED?&CjZIBkJ80$E#r?P~AX2E_S5$m!V^CLCW?``&
z%a(t(nJP&B4Y~&Ucv8fbjGP=xYioHI68Q8}=-9;MnYAdVxiyT#c(Su=zOK+T=9n1vfFK+^>rDrTbMO9rBi?|stu8M3tccHv8{pja#&gcccX@b*KG-Beu@7kGzH_JE`}pI5
zOKe_V9y~+1i=u);-h;q09-0$
z`Q_z6zILO>uEGyr=y{;5W1*W5S6C-<88Z%78x%pG40qt=74q0&4mK5`#HDGS-wyEi
z2SfzD$ja(!{`>bHv+P6|qN1YssLsw#B8=Zq*IkjU2}eU#>D@_uAqOuU$&hDrUd~
z8#EW|`N+m|d6i;&Y2X&ga#v?OF+IQCvu8VV-N}IWFJHc_k|r{AngdI>-f$eo4dQFwsr(Y1B8Th1j*vULf-rLDzK>|_Kn>=J@(KY#>y1)7#bQTzj;Fm
zX;nw3<6*gWp&r9h`s>$KuIuBlj2@nz#d;sYb$mv2P@TU9cwzkua7Tl<)sI|2ENNbjOwV5LPA0aS6eS}*
zus|Dho2<|?IFK>}xR=)Xn_PluU!Zmk;)9T`D-fssD@Car!yT90Vv`^2UXJ5l*)zG-Blqw`0
zC)xh_Gd%6hn>VpqE?ST}w|ftg(7e4**G9;3WW8!TAPjoaL`S}Ve-RrCbPEMJd8i9E
zc9eO0bi&J*8*U=1s;U64GaoM9VNC9RtNsq!BB`v53>s7@a~B>8*rC?oloWDQUSVNW
z*YCc?XT-z_l$Y38Swl0zj#1g!*>94Qp}7dthIJz$C4HTeV*B_ptAYtM-BkWn-_@w`
zUM3`z4%v!}iwAz6?Tp7q8F22WuA>3S&UyRRZsg;Wy_F#!AD|6?y?smHPK*(htCH$*wDT_Stm^4+g?+wWDQqcbF0R2JKZqtc%`7Y|?CjEByuc$S
z{yj58$1&!VDu*6`kB9dMFik@Gr*=vl9EH$+0@IUOT2Vnnlbw?@F!=2D?b`=lVve!2
z_<_GJKtFW0l*Px!_EM)u^Xa}q!ondzK^*!Y3q13n#eRsy#|X+SE?%FSicLre!h6GG
z9$Zvpx9{|OOi(}ovh{eQ?*s(G@Nf;p#h~5vb}RSFSf%Y~6%QYuQ%f+ms1e1=OGw{X
zg7Av#Dyc_}RaMDTQ)?t7=meB+s#1S;bjbZMd-CKIc*$uZdK82yD>NiIIlr6yu8`xj
z1EJ*Q=?WafXru2rK!cw}hKdyq)dttDh1d)h)%%>*ynDwYBveyY);Bs@SN!_u=m?tM
zOP4M|y0hi??0>7Q#HgXHTv)tFN=Bxqr-zP?{&R3JF)v5xJAUnWmEPvRILYk7g08wEtFyBcH=C@N56L%&Jy4UDJ8
z#}fvjbW2PwIid(@vNJMDO^($=s
zk>T~YHnThY{A8r0?1nWz+S-C%DdXr$mo-f
zfZW{NT7c{w#ot<6EtY>40x5CDc0D62OVZF#l$STIig^F@^AR2}YCBt71cm39JEVWC
zLHe*@0w_gyL!=b@Bqot7*~`o9d=nP{fI-JM%Y`Wc=x$ki1^|pUQ|Wjscm5qQB_vv4jE~8jd$i71Xtnx7#MP@s;XS=pWM`fUBE&j
zbt>&BqGa?NyacONLoaYL`c=T4#)gKUoR_DK>!;^uy}C_0VpCPpq5R#>+gn4gP5B~W
zhSQoYA0L{V0Kb69`bH|^fi~W=Ha3re!Ad|lX#I#eJ{41f%ea9^l-Pm3PsLM
z&;+V$Y!ec$lM7C>Ofe8Q%b{PTh3h~EEP)jD()r+I^-%3Qt+23pnbS>g+b2)La-Pio
z{tdeiwISFZ=?{@;5!}(!bE8lO`e2h$e>2Cr{9=3EO3voKK9jCj45$0(^)!?&?K2mp
za~_!3F<>WO$onVYx<$X^S;*wF?Ua>gpTj2h_k`z^Nk*4okZY
zBEf&4pimc2Pfx5Cps&)>jDUoFG|>6x4HTT3u7q`O8a|eC=YA!kI^AW(4Hqww9I`0Ua3$qQ?zN$|7I`^YRwf*Cz()bdcVu
z!D}`~up*!7^Kd_Z@I7VA&kH5IwJaQG?@P%6uoH+n9;{?hPyXKN)m1}v^`|qRqkxeF
zLPl6n5Wh`VMFpr0IJNVDO}&qvK^f?s0BoYA^|7MjHV{fQG#L~4Gt<*P2U8)RKg+nf
zx@v0=0_p#lU|oFl<_4@~>Gf+fnldteHAizm_bJB(l<5@k
z?B(m1%@+6X11BX#6FMW1+h;-O(hOAd
z_V$jU!9zh+5s2=2b02YNZe!;HP7V$L0d_YwewACPDl12+F-L}mE-Wpfkk*eMH@I)V
zEh>r*g4--2=h@D&2FmDBNsx}u_uz}un0x?poy5(FKb#u*5e=4R{ERR1j5
zdS;@?m1m!cAWHyi>2rR%Z$Hs6@PCvS$k2os7nu^42ILl?>smdrv9QEIMJjGMdJhyY
z@MsW$SE#g&jVFeOg{Z0jc>6ZE`=R}CMOdfB`*-iK+sTjL3_3&tkT;=!Zss@|7FF+l
zxYbW5E5f9b9h=%kN<>Nu0wEa{RScejl2QsNhVh*bG&SKiQ0|hrP5q0C*3M3kMMXq_
zrv{Y|xH#!`%TgULv5BUprt$Idkr7X|nKB)Ap;S9tTU$W985yy!UkkJnF25ag$DVM2
zMd%2cI`#Rs3B4@e;bG=GBWN2MW^E4?#Fgo7bbcqcBy|pF;bTe%Ic}%yO4!s2Y9~7!yA;
z$Cc%Ug`=f>3k{9b;zlE>zuney^7hWJ{nP!SJ9nbw<0iH+
zWupsz{CL194CfWt4lyU~1Dl5bkf
z(?S^9?&XJL7`eN9!~b=44Xv%8C@W`4Jf{n^u*ha6Qt7opI+@xB-?hqqNw+6Oi@V4RQo2Gu_r>Ft!BGpiY268`!bSxz={S=kS#`O8+{0`;|;
zf4T0j8A2LyAu(Y;*y&fDQt{<}Flo0?%1A$35RM@#;C6i|3=49?;R2mBWpPfq#`i+}
z__$1Zrs5C<{7AxUxRju0u@C7!FMZx-R$WsAdOs;0gr^2CFGi7p`_GDQCK1NsbaMHQZ(%nRb;ZbP1oO;r@v)1L?qv;ej$IeB@a3cs&j
zo!;woH}l1doU^lSIeFaY_nTW==|Lu6VV78EXJ-d)KpMzAECyzhcK?zRTa2JGtDYMc
zOnou)J;uXB(iT?YVve7V4>mOw6t>p)l6a*Z``-q|^ZVV?~Pbu1+x!ThsMX=FL`DiSXcKoM=^QIObXG2_ob)DxYnicHBPxzkBXe!iw8pM
z?_}dnRWJ}?el{esl|`IWuzF7OJ|b~l#d~|-L_N#Y)wfPh3iehUU%py9mFKv4Z&JU>
zj_Cyw!gy^xJ;)JVohsR?IXTXXi83kT=Odl*`Hjclt^+3<+zR@%DF}up+)(Hs?DWtxS}hOi*OZa1E{6oBmlJsdT!tCxv1BZ^}IYh*Ku$@mX&4p
zFM^Q)I4ewCT<8{|r{3A$zkv}XRcXZJf3_dF8cQd;v$F$5I3ZyGG$wo$Fd-248@s!0
z++RO^$}A>F7q$^gry?TC&dSQl&bF|(pIBZw3Kx7qU#q$_rdn03i{NH=;C}2_U5t(G
z8W=JkmnxU;bF`S@*X-hXEJ^;2t<33?36TSkB`z*&l_MafLl0eS+&KQ@M+G<-&;#yk
zYlqBKftUtzYN%MN`4!-;qGk@DX8?7(%D@g99{;PD@NWllC<>tyWa;?fduIsoVnj*eKFnG5tQi-6atsyfJu5mQiAWvbZ~!Y6P4i0O@;7}DB`@Qqw^6kpzqoCC@}CAxmTJN
zxHrrk;o&pKZt^_@m2$?LDB>S>g44g-2_^+7@$e>p{|@?-VuiMXmyd6PVi_7VJ3G5f
z#|y8)lj%IlkQSQgJRnMT+-4;FUDJCnn3~o9U~2Hl+uM|+9%v^dCW1lW0%dMqGh)Ya
z_}SXt_4{dT@e$f6jBp8&WL
z$Hx~Htp~mjW{RLk
z4{y8wuD;%Bu;?pS!+`~w7*w0KgW=s-x>TgT5fdDV}*C&AfYBu2pPQA8$QSs^1Es!h{*>!;zI|m=lRZNl{
z;*%?3Nk~&LUnwXk0F9!nD@aRAD=M0*l1mr((SE{bv5!Geu)1L_*e6DpiIGtyT`U!V
z{QLL7M;#p;D81egSJhu-oAaa}8f3zHZeen8-WX36NwWNtLB#t@y`yN00^rsq#4h$8ra`3psY+8{;lL&-7
z8Ozmcq)0!BDl64jA9gEx#P8C=n!hi6VpM-blzhr6dfFxEv~%ox1!2od!(h4>;I$`LTNM``e-SIi$
zNkla_oO$pJg;587FVt^qXmCvTc}z-Cuczg*plT9kWo><3HV$zq9AsBgyOL+$7C<0K
ze(mN7jV8Ke^ooI+)=x>*o_M9XufKG}JGH|vK*I}qc(;)C#-1dhT{_&DacyLL_Q8R^
zpI_(rcshkpMVvy=#D!MQ{khqC;r`@!{Z?0^ud2MRYQ_sMds#QPwdR0C*Ht|-ybZzX
zBm}}^a`OO`xJexDKh2ts*)5JI$LmgF(C#+-FU~&P8grdjP5TOf&zYWl(n}l)8}OnO
zaZE9nxRt)8i3&IQ$Dop^IYg!PJi;q>Wc&KwKY0TAc`HO%6wxGr);2Adn0k899PmXk
zkw5r2*GVN;RYwsU|-);
zKx>p-AbXjU%4DbkIQ#VJlN7EqH;9?_o|REi*a6M2`JW7f9u?;D@4(DhQ7nV^e
zMPTz_a|ZOGMzCw47JikQr+aOMWt^ST|lMz3AH3jI6_8yl3Av+m@($PB1)0NNmXsU-1Q
z?k)`!8`Ts1hA8>Js5veo4CyhK)!_?M-9IpLh0T9pY;3qbvI`<{6L
z7Tt85{*x5Emic0dO%zrQ)M#L!pmu><;r=(XPkWCC-fIG4tplSN>_<}45V$Y#QBd?j
z{End%Wo2i#C?A3cA=vHt@&%N~?*Q~U*x5ms-T}Q65E%pnpy97ymF`#mw<5crAfcxR
zVGF2KBzgu&rFEX5O_Oz_8>KOCTW~1!tj#hxO1;2tLxdxfT~qV$O2BGt?4aCDLuMko
ziL|1^`rk_fx!Ls8=%~;(v#|AYi%x>++t0xL3no4g@r@e1!e&7JbyrgQT6+qwu`)L&
zV+$-Jn}|rtFGid~Fc|YEN8wFBt(K~lpf=$q1vnrO?Nn_58WdQzrv<*vZ2lM+fQs%&
z?&|D3R$)yTJ%h?nla&Rz_h!2zk0ZDXq{;6Q#T$X88k&Z7c)`!t0VTleSH{O8pk#5zGP&{`|FKzE-s)|UbMYx?a_uB
zYA<$oQ`rtP9LaI9yOk7cU6)?2&!^v#7xSZW1eZKk3%uAm#8SVixN;J2j%>)kp;x?+
z7X13`_eT@C_F#{u34|^JG%h?885tQMP*O$)hJn#hr4+5`s3^Xd!9Zx&0ax2*
zMic1e?G4oZc!T#b6f$ria2hv|)6g^GpyXgFy*vy!6IF8a^YZY#!I2w=y=giD;{q%U
z43;J)-|J($A8NdNFulQq3KcL_{^5YkfX0FF7K9~`^|#5CD@K0)WW=CV2>Z2m_AdPH
zFP2Cu7ngEUtZK?iL5+7WffM*>iD|9PcLnj87>ZgZFfiHK!pZy`FXQQ^GIycI$t%=&
z?&J+-O8ds2EJy|kly%H+!`BmATaP**+n1-NO5vyg)ei;|e0-4V85rOElaK9t>fnI)
zcagJrrWjx$%{M6!wRN<#_`wH2L(@DqmYOGg94;A2%ZfZwP?eMn=$R+H!9p7d=({a<
zaPS%6Gmy6(neR&}D(1siP0gEM92^`BWY@N=Bk=B|GmFCU3VZ?QFz=5QpIrQ)t-(g|
zfYeU{j(fHl$s=K^V!ihkipzT+uu2;n-{E~qyqLcO?|c2?eFLs>slJig->Fz>Tvi*Y
z-A4i91cH@t434L`etc;{Sio0(8hq78D0oX=z@ZR|A*H3M2*V!zl^+us19Ib
z#KaoFGyDGiV=!rEegh;Ar&NkUq)GMR$-)wW&$;)vx5_MRTCe}&u)uUc6Ua=Ybno8Y
z*Q8yCV?Kn085A|gUn>CJ_x%HIbalvi5Py{n;iM9FS)IJ+%EY8`9hbut$^)2D_BS>l
z4Kz15w+fnV5)cxCI{M9m7mPPou3Q0leokAkAxg16Q+J~sBk$@}yB=I(BHS66?3
ze=@j*4g{{K#{7#^`BJAIvhc&_rnQPHDrqon0p0}&dMXjO5|if3v;zNPyp>Z@TG=;W
zY@V19zkT}wNQMjyzdr?GwstOqLb~ZD0_FvF;-VrcmuN^1Xf18yHrMF6m=}~xO)ZU#
zkO{24>S8ogE{`6~jn{w6*NR$jKIg$QFV%qzZ^49W_?w-Uq>3grEUXj6wfed`dclQX
z$D(y^y$g#s!IGO6H!bdcWGw4)!C+@)ZAVjVlWqHUT|&r^x%qc>s2drjC+tpcn5tga
z4hIA327o23>!r7*5@=0Ibsn%fZ%(&CPbciM;^6AaBuvNsGEFppb0<
90 // 9. Forest canopy
void ApplySpongeZoneBCsForMom(const SpongeChoice &spongeChoice, const Geometry geom, const Box &tbx, const Box &tby, const Box &tbz, const Array4< Real > &rho_u_rhs, const Array4< Real > &rho_v_rhs, const Array4< Real > &rho_w_rhs, const Array4< const Real > &rho_u, const Array4< const Real > &rho_v, const Array4< const Real > &rho_w)
Definition: ERF_ApplySpongeZoneBCs.cpp:117
void ApplySpongeZoneBCsForMom_ReadFromFile(const SpongeChoice &spongeChoice, const Geometry geom, const Box &tbx, const Box &tby, const Array4< const Real > &cell_data, const Array4< Real > &rho_u_rhs, const Array4< Real > &rho_v_rhs, const Array4< const Real > &rho_u, const Array4< const Real > &rho_v, const Vector< Real * > d_sponge_ptrs_at_lev)
Definition: ERF_ApplySpongeZoneBCs_ReadFromFile.cpp:8
DirectionSelector< 2 > ZDirDefinition: ERF_DirectionSelector.H:38
void make_buoyancy(Vector< MultiFab > &S_data, const MultiFab &S_prim, MultiFab &buoyancy, const amrex::Geometry geom, const SolverChoice &solverChoice, const MultiFab &base_state, const int n_qstate, const int anelastic)Definition: ERF_MakeBuoyancy.cpp:32
+void NumericalDiffusion_Xmom(const Box &bx, const Real dt, const Real num_diff_coeff, const Array4< const Real > &prim_data, const Array4< const Real > &cell_data, const Array4< Real > &rhs, const Array4< const Real > &mf_arr)Definition: ERF_NumericalDiffusion.cpp:84
+void NumericalDiffusion_Ymom(const Box &bx, const Real dt, const Real num_diff_coeff, const Array4< const Real > &prim_data, const Array4< const Real > &cell_data, const Array4< Real > &rhs, const Array4< const Real > &mf_arr)Definition: ERF_NumericalDiffusion.cpp:149
AMREX_FORCE_INLINE IntVect offset(const int face_dir, const int normal)Definition: ERF_ReadBndryPlanes.cpp:28
Definition: ERF_PlaneAverage.H:14
@@ -915,35 +919,35 @@ amrex::Vector< amrex::Real > input_sounding_timeDefinition: ERF_InputSoundingData.H:315
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > U_inp_sound_dDefinition: ERF_InputSoundingData.H:326
-
-
-
-
-
-amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcingDefinition: ERF_DataStruct.H:670
-
-bool nudging_from_input_soundingDefinition: ERF_DataStruct.H:641
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcingDefinition: ERF_DataStruct.H:671
+
+bool nudging_from_input_soundingDefinition: ERF_DataStruct.H:642
+
+
+
+
+
+
+
+
+
+
+
+
+
+
Here is the call graph for this function:
diff --git a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map
index 7e27f65121..45a1190122 100644
--- a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map
+++ b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.map
@@ -1,5 +1,5 @@
diff --git a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5 b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5
index 1d45be5335..555b509dc1 100644
--- a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5
+++ b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.md5
@@ -1 +1 @@
-c219d0622a1b348eda802510569f41c0
\ No newline at end of file
+ec241a01c30958ea7c35f3d4a480678a
\ No newline at end of file
diff --git a/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.png b/ERF__MakeMomSources_8cpp_a186d85f6035b61f034bac8995f82862d_cgraph.png
index bf33e240600cd06a8fc2cb5a564193e80888c5fe..8b2c87b36faf744333512464ce28e96b1d89ba3b 100644
GIT binary patch
literal 70087
zcmce;by(Hg_C35+#4X(-(%q>@gPU#y=`ImLLZnN&K|n;L4kaiZN|%6uhzN*Cr*wDo
zo1Am|-tYV8%e~KYd|>m5wdR_0%rV9!Oie`=8-ol3fk0r(%Sovt5LYY^h|A>2%kUkN
z$fFSW54xG6tQ6w>;(veC=f)!t)ChU0dk@@FR>s_QG_;QdHk>3SsA-c!j*>T1T@-{LOOLAI>001lrY-zKJ(1P;E>T7J2Q3#l+ROKJ
zX(L{YPEfh}ums>q&?Npo6*-+5kB{H4W9l?FUVF+JJR@-NJP5=ecZtdf|NIw%?gmlB
zU*A(zWB>D~&&BSQBLDo>kradQCGLjNde8mI<^T(;OHIb?92@}w0iECAlQNtFWk~w)
z?uY53p08He5eLO5CnvS8>v5@95s1h~4i@+?sOjheR)`Q1NSYfrZj@;t5YJ{u-&t*xz7lzCNEzdAeRH8nML
zb(1(XJ-e8jq!TV5Kk!{5>A6q*OvKhWo6|=4t;*}&yg7v
z%*@Ql_RYP$JwB^%`T1Wx>gn$O@bP1skn>Wiu&bNMcoLu0n>TOD%FBp=AF6$HZMwL~dWCePUC&~myptUZRRI4|$ubbms6m64IL
z)@8MpP>mn@bXbpdl4NS9Tyi@U+)bUiAjNV&dp^)F$)SNp)%^)zn|jV
z%TX6UUMUxJn!}^wjh=y5$VpG1@VjzN?CikY-Tkof{B&LXWFsIypAEIp5liRc;c=U+
z`^%T5sg@A41~DG@;NakI>Eg!s?$K~?VAE7QTU4ivaXmltrXBh5*$~!=s}iP|71Cl_VsZI^&sY8;OXC
zcgI{O90kbl-n}aNZ6cn181Mf6)eU4Wv4-XGl
zS64;*7Vf;w&$rRg=+Bak!N9~EvfX>HnJ;xpNJzN9Hg>%AThZQrT{B;6N+25Vsgsk&
zg9q$vY+VqsDk>z-5D}a#EY*%PGS1GsJB$6PDJdQ=1x70Ex!Bo7e9uIpWa#wJg`YJ_
zK67^dJwJbgh)73QmnJYx)Kg1WcX4&~bnULXfkB<~^3U1XSy&z_-xHoYcRm!at9Ye3
zIy%B?N}8Kz)z%85LQxq9$H!FSKIAknU%s@kw9F|i6m3$|)zuBbq3}4|$_WZ;gIGut
z_ceO(;DL_L>(tcLq$CzbQCPEp=U1rjv6I+PHwtkk=I7+>?yc&Rx4Vm^)zo;{*syK(
z73h|`tc|iP2|}VH(>r*sJUPU;w9iwA3F`>?erH`hu|uA^4NEwTLKYg@wOGNP&Ksw_#f9MP>?T@TAG
zEiAaPH16Mb-<-aB-;!WwcemEwNq~HFmG;(MI}y)4TS-Yt85#5h1;`lG!*_1qj;0XS
zQdj>pvne7fYG_~(`sUr>Yg_(gg&cTbf+*NiBuyr}B=X(+_XM||K7ESJ`1p~ZGf^?O
zU-Q8O543{3JnCS^HO$z!xW{&OMFj=)M)fYMdL||&PoI{(d-smnuNjBK2RV(B^pjB2
z()#?c{Pg4m4;9a(DuZ;a9(8RX%+AeCeDkK}$;ii{a;xLRLs6d-ogMXk8xlxfjAA{M
z*KoqVCMmsbY&=iCYwz(gZCVn4h&B)c==bm6((T5ucWtAw7l}#j4;?hRXfAq)+D=XwABp+>u%GNT_g)3guNHG(IH#9Vagve@Xk)R%Zj=Tw%9P3GNou8b1k(}%a
zIpfx?TafU7{`_fF8`eYEQsO#M<~Z~D>({T?xA=G5oJ%y&5eV+ti&w&N6fTU7jm>Kt
z9zKHjY-N!PRD&a8WMrhM*ag{3fIK=rKK@{Djo)sx%J=+C*mbSIu=drXsi)OAI5-fq
zo~K978X6kUo;`z8S6%IH%(J_@TUc0lTvK{XDmC!@0o1
z!YU{(=HlX-sCHuNEh;QrU07I{nb9yX7=<`hQ=6Qf{n6XoJhS|}rY9A7I~yBYe_!9j?l(brpSV|A!w5@DOKE6m7#JAfuBz!`*x1<0s;Z4oe}5@5
ztW{7{6c!e~dikEdzP|C@SS)_kz!#80qi3kZEaB`+;)F|cPfN+l@+Ey;%+m
zyz}hj^b`VlU|>K)UHzhd2;CUc_HCVvzb
zQb0zeqH1)QLh714jT*CRc-QC?#dwNm?B`ZdT
zhv$C(KHY4mA|WAxI!{SS86F-E8?rsu4bN=GoAl+&7kqsD#l=N4ElG0Qv}y7B3=9A)
zlcFX9(ebHQNM-x_`^QE{#b&jw%*@K1m$jf6T>7~Gy|*``sEG3BP5IIAV;^(x(1;sc
z9`5d8KYcIF(mvjgjKnGjS+W6-*X>HW-eJ+tNhrO!9pIZ`>S>dae6)6~@6Vk1k8pw+kz;z;Lj-@aMo!u{#UM4GQgBE?U4O5VTs>-N#q
z9a1yw_B05^@L(<9XK|6uQBZy(CG5#B
zp{aQ|9Sux>Qrc
zYVjMg1tQW4fb*8-gCDik(oe=yjt=Mf(EJ;5IMv?~kty@D+4TQ-ry-|=3!fzIfy835
zv_d2$W29eVa=tuhUsm229lcjto-CcCKogjq&CJW&x43AmlDvmMWPJh0XKoafBjq4;
zYfUB1EE%{2+%|);G2w&A4Ds`V!;@_)eALsYyt=x^-6_{q2jMpm9cZ#7mjvVD*4T9k
zRB}~2x-hRJ4^r}o@KOD_`9niDyE^4_bA7KTT%c|Vp$iaSSY7Si9=)?YPjz^xyE-&;
zu=Z5Y_e{5cbd=vJT8w4!Vnr{4!cR)nX2Qv2q{#%OVaf{Wv!Eb7y67@@U=6&n2B^Gg|VsR+{J
z@w!DMmeS$wpw%_pn%e4a8exiy%VJ%IU3AmbymN5$`RBl$+_!cm6}PGX_A{K6VDy%S
zv(u;7LL`YWUM1IFWtE2CF}dw;hL4pvJ6p@i*$#lQygcCcorJ$V=lT^D$pBM)$}~U?
zzb7Z#nwu55vFQG~U>y@PNi!ZHGB6_}iG}6Rb6==fZ4+`$j4allgOI
z=gxel11sxK*~|JbR;DTB;*Lpo?S+08g!%5|y=gn@8zrawd9gB)SK;7h+}9qQUcA{o
z^C3leKOr_&=&y^KJi53j#B4x_ED?&CjZIBkJ80$E#r?P~AX2E_S5$m!V^CLCW?``&
z%a(t(nJP&B4Y~&Ucv8fbjGP=xYioHI68Q8}=-9;MnYAdVxiyT#c(Su=zOK+T=9n1vfFK+^>rDrTbMO9rBi?|stu8M3tccHv8{pja#&gcccX@b*KG-Beu@7kGzH_JE`}pI5
zOKe_V9y~+1i=u);-h;q09-0$
z`Q_z6zILO>uEGyr=y{;5W1*W5S6C-<88Z%78x%pG40qt=74q0&4mK5`#HDGS-wyEi
z2SfzD$ja(!{`>bHv+P6|qN1YssLsw#B8=Zq*IkjU2}eU#>D@_uAqOuU$&hDrUd~
z8#EW|`N+m|d6i;&Y2X&ga#v?OF+IQCvu8VV-N}IWFJHc_k|r{AngdI>-f$eo4dQFwsr(Y1B8Th1j*vULf-rLDzK>|_Kn>=J@(KY#>y1)7#bQTzj;Fm
zX;nw3<6*gWp&r9h`s>$KuIuBlj2@nz#d;sYb$mv2P@TU9cwzkua7Tl<)sI|2ENNbjOwV5LPA0aS6eS}*
zus|Dho2<|?IFK>}xR=)Xn_PluU!Zmk;)9T`D-fssD@Car!yT90Vv`^2UXJ5l*)zG-Blqw`0
zC)xh_Gd%6hn>VpqE?ST}w|ftg(7e4**G9;3WW8!TAPjoaL`S}Ve-RrCbPEMJd8i9E
zc9eO0bi&J*8*U=1s;U64GaoM9VNC9RtNsq!BB`v53>s7@a~B>8*rC?oloWDQUSVNW
z*YCc?XT-z_l$Y38Swl0zj#1g!*>94Qp}7dthIJz$C4HTeV*B_ptAYtM-BkWn-_@w`
zUM3`z4%v!}iwAz6?Tp7q8F22WuA>3S&UyRRZsg;Wy_F#!AD|6?y?smHPK*(htCH$*wDT_Stm^4+g?+wWDQqcbF0R2JKZqtc%`7Y|?CjEByuc$S
z{yj58$1&!VDu*6`kB9dMFik@Gr*=vl9EH$+0@IUOT2Vnnlbw?@F!=2D?b`=lVve!2
z_<_GJKtFW0l*Px!_EM)u^Xa}q!ondzK^*!Y3q13n#eRsy#|X+SE?%FSicLre!h6GG
z9$Zvpx9{|OOi(}ovh{eQ?*s(G@Nf;p#h~5vb}RSFSf%Y~6%QYuQ%f+ms1e1=OGw{X
zg7Av#Dyc_}RaMDTQ)?t7=meB+s#1S;bjbZMd-CKIc*$uZdK82yD>NiIIlr6yu8`xj
z1EJ*Q=?WafXru2rK!cw}hKdyq)dttDh1d)h)%%>*ynDwYBveyY);Bs@SN!_u=m?tM
zOP4M|y0hi??0>7Q#HgXHTv)tFN=Bxqr-zP?{&R3JF)v5xJAUnWmEPvRILYk7g08wEtFyBcH=C@N56L%&Jy4UDJ8
z#}fvjbW2PwIid(@vNJMDO^($=s
zk>T~YHnThY{A8r0?1nWz+S-C%DdXr$mo-f
zfZW{NT7c{w#ot<6EtY>40x5CDc0D62OVZF#l$STIig^F@^AR2}YCBt71cm39JEVWC
zLHe*@0w_gyL!=b@Bqot7*~`o9d=nP{fI-JM%Y`Wc=x$ki1^|pUQ|Wjscm5qQB_vv4jE~8jd$i71Xtnx7#MP@s;XS=pWM`fUBE&j
zbt>&BqGa?NyacONLoaYL`c=T4#)gKUoR_DK>!;^uy}C_0VpCPpq5R#>+gn4gP5B~W
zhSQoYA0L{V0Kb69`bH|^fi~W=Ha3re!Ad|lX#I#eJ{41f%ea9^l-Pm3PsLM
z&;+V$Y!ec$lM7C>Ofe8Q%b{PTh3h~EEP)jD()r+I^-%3Qt+23pnbS>g+b2)La-Pio
z{tdeiwISFZ=?{@;5!}(!bE8lO`e2h$e>2Cr{9=3EO3voKK9jCj45$0(^)!?&?K2mp
za~_!3F<>WO$onVYx<$X^S;*wF?Ua>gpTj2h_k`z^Nk*4okZY
zBEf&4pimc2Pfx5Cps&)>jDUoFG|>6x4HTT3u7q`O8a|eC=YA!kI^AW(4Hqww9I`0Ua3$qQ?zN$|7I`^YRwf*Cz()bdcVu
z!D}`~up*!7^Kd_Z@I7VA&kH5IwJaQG?@P%6uoH+n9;{?hPyXKN)m1}v^`|qRqkxeF
zLPl6n5Wh`VMFpr0IJNVDO}&qvK^f?s0BoYA^|7MjHV{fQG#L~4Gt<*P2U8)RKg+nf
zx@v0=0_p#lU|oFl<_4@~>Gf+fnldteHAizm_bJB(l<5@k
z?B(m1%@+6X11BX#6FMW1+h;-O(hOAd
z_V$jU!9zh+5s2=2b02YNZe!;HP7V$L0d_YwewACPDl12+F-L}mE-Wpfkk*eMH@I)V
zEh>r*g4--2=h@D&2FmDBNsx}u_uz}un0x?poy5(FKb#u*5e=4R{ERR1j5
zdS;@?m1m!cAWHyi>2rR%Z$Hs6@PCvS$k2os7nu^42ILl?>smdrv9QEIMJjGMdJhyY
z@MsW$SE#g&jVFeOg{Z0jc>6ZE`=R}CMOdfB`*-iK+sTjL3_3&tkT;=!Zss@|7FF+l
zxYbW5E5f9b9h=%kN<>Nu0wEa{RScejl2QsNhVh*bG&SKiQ0|hrP5q0C*3M3kMMXq_
zrv{Y|xH#!`%TgULv5BUprt$Idkr7X|nKB)Ap;S9tTU$W985yy!UkkJnF25ag$DVM2
zMd%2cI`#Rs3B4@e;bG=GBWN2MW^E4?#Fgo7bbcqcBy|pF;bTe%Ic}%yO4!s2Y9~7!yA;
z$Cc%Ug`=f>3k{9b;zlE>zuney^7hWJ{nP!SJ9nbw<0iH+
zWupsz{CL194CfWt4lyU~1Dl5bkf
z(?S^9?&XJL7`eN9!~b=44Xv%8C@W`4Jf{n^u*ha6Qt7opI+@xB-?hqqNw+6Oi@V4RQo2Gu_r>Ft!BGpiY268`!bSxz={S=kS#`O8+{0`;|;
zf4T0j8A2LyAu(Y;*y&fDQt{<}Flo0?%1A$35RM@#;C6i|3=49?;R2mBWpPfq#`i+}
z__$1Zrs5C<{7AxUxRju0u@C7!FMZx-R$WsAdOs;0gr^2CFGi7p`_GDQCK1NsbaMHQZ(%nRb;ZbP1oO;r@v)1L?qv;ej$IeB@a3cs&j
zo!;woH}l1doU^lSIeFaY_nTW==|Lu6VV78EXJ-d)KpMzAECyzhcK?zRTa2JGtDYMc
zOnou)J;uXB(iT?YVve7V4>mOw6t>p)l6a*Z``-q|^ZVV?~Pbu1+x!ThsMX=FL`DiSXcKoM=^QIObXG2_ob)DxYnicHBPxzkBXe!iw8pM
z?_}dnRWJ}?el{esl|`IWuzF7OJ|b~l#d~|-L_N#Y)wfPh3iehUU%py9mFKv4Z&JU>
zj_Cyw!gy^xJ;)JVohsR?IXTXXi83kT=Odl*`Hjclt^+3<+zR@%DF}up+)(Hs?DWtxS}hOi*OZa1E{6oBmlJsdT!tCxv1BZ^}IYh*Ku$@mX&4p
zFM^Q)I4ewCT<8{|r{3A$zkv}XRcXZJf3_dF8cQd;v$F$5I3ZyGG$wo$Fd-248@s!0
z++RO^$}A>F7q$^gry?TC&dSQl&bF|(pIBZw3Kx7qU#q$_rdn03i{NH=;C}2_U5t(G
z8W=JkmnxU;bF`S@*X-hXEJ^;2t<33?36TSkB`z*&l_MafLl0eS+&KQ@M+G<-&;#yk
zYlqBKftUtzYN%MN`4!-;qGk@DX8?7(%D@g99{;PD@NWllC<>tyWa;?fduIsoVnj*eKFnG5tQi-6atsyfJu5mQiAWvbZ~!Y6P4i0O@;7}DB`@Qqw^6kpzqoCC@}CAxmTJN
zxHrrk;o&pKZt^_@m2$?LDB>S>g44g-2_^+7@$e>p{|@?-VuiMXmyd6PVi_7VJ3G5f
z#|y8)lj%IlkQSQgJRnMT+-4;FUDJCnn3~o9U~2Hl+uM|+9%v^dCW1lW0%dMqGh)Ya
z_}SXt_4{dT@e$f6jBp8&WL
z$Hx~Htp~mjW{RLk
z4{y8wuD;%Bu;?pS!+`~w7*w0KgW=s-x>TgT5fdDV}*C&AfYBu2pPQA8$QSs^1Es!h{*>!;zI|m=lRZNl{
z;*%?3Nk~&LUnwXk0F9!nD@aRAD=M0*l1mr((SE{bv5!Geu)1L_*e6DpiIGtyT`U!V
z{QLL7M;#p;D81egSJhu-oAaa}8f3zHZeen8-WX36NwWNtLB#t@y`yN00^rsq#4h$8ra`3psY+8{;lL&-7
z8Ozmcq)0!BDl64jA9gEx#P8C=n!hi6VpM-blzhr6dfFxEv~%ox1!2od!(h4>;I$`LTNM``e-SIi$
zNkla_oO$pJg;587FVt^qXmCvTc}z-Cuczg*plT9kWo><3HV$zq9AsBgyOL+$7C<0K
ze(mN7jV8Ke^ooI+)=x>*o_M9XufKG}JGH|vK*I}qc(;)C#-1dhT{_&DacyLL_Q8R^
zpI_(rcshkpMVvy=#D!MQ{khqC;r`@!{Z?0^ud2MRYQ_sMds#QPwdR0C*Ht|-ybZzX
zBm}}^a`OO`xJexDKh2ts*)5JI$LmgF(C#+-FU~&P8grdjP5TOf&zYWl(n}l)8}OnO
zaZE9nxRt)8i3&IQ$Dop^IYg!PJi;q>Wc&KwKY0TAc`HO%6wxGr);2Adn0k899PmXk
zkw5r2*GVN;RYwsU|-);
zKx>p-AbXjU%4DbkIQ#VJlN7EqH;9?_o|REi*a6M2`JW7f9u?;D@4(DhQ7nV^e
zMPTz_a|ZOGMzCw47JikQr+aOMWt^ST|lMz3AH3jI6_8yl3Av+m@($PB1)0NNmXsU-1Q
z?k)`!8`Ts1hA8>Js5veo4CyhK)!_?M-9IpLh0T9pY;3qbvI`<{6L
z7Tt85{*x5Emic0dO%zrQ)M#L!pmu><;r=(XPkWCC-fIG4tplSN>_<}45V$Y#QBd?j
z{End%Wo2i#C?A3cA=vHt@&%N~?*Q~U*x5ms-T}Q65E%pnpy97ymF`#mw<5crAfcxR
zVGF2KBzgu&rFEX5O_Oz_8>KOCTW~1!tj#hxO1;2tLxdxfT~qV$O2BGt?4aCDLuMko
ziL|1^`rk_fx!Ls8=%~;(v#|AYi%x>++t0xL3no4g@r@e1!e&7JbyrgQT6+qwu`)L&
zV+$-Jn}|rtFGid~Fc|YEN8wFBt(K~lpf=$q1vnrO?Nn_58WdQzrv<*vZ2lM+fQs%&
z?&|D3R$)yTJ%h?nla&Rz_h!2zk0ZDXq{;6Q#T$X88k&Z7c)`!t0VTleSH{O8pk#5zGP&{`|FKzE-s)|UbMYx?a_uB
zYA<$oQ`rtP9LaI9yOk7cU6)?2&!^v#7xSZW1eZKk3%uAm#8SVixN;J2j%>)kp;x?+
z7X13`_eT@C_F#{u34|^JG%h?885tQMP*O$)hJn#hr4+5`s3^Xd!9Zx&0ax2*
zMic1e?G4oZc!T#b6f$ria2hv|)6g^GpyXgFy*vy!6IF8a^YZY#!I2w=y=giD;{q%U
z43;J)-|J($A8NdNFulQq3KcL_{^5YkfX0FF7K9~`^|#5CD@K0)WW=CV2>Z2m_AdPH
zFP2Cu7ngEUtZK?iL5+7WffM*>iD|9PcLnj87>ZgZFfiHK!pZy`FXQQ^GIycI$t%=&
z?&J+-O8ds2EJy|kly%H+!`BmATaP**+n1-NO5vyg)ei;|e0-4V85rOElaK9t>fnI)
zcagJrrWjx$%{M6!wRN<#_`wH2L(@DqmYOGg94;A2%ZfZwP?eMn=$R+H!9p7d=({a<
zaPS%6Gmy6(neR&}D(1siP0gEM92^`BWY@N=Bk=B|GmFCU3VZ?QFz=5QpIrQ)t-(g|
zfYeU{j(fHl$s=K^V!ihkipzT+uu2;n-{E~qyqLcO?|c2?eFLs>slJig->Fz>Tvi*Y
z-A4i91cH@t434L`etc;{Sio0(8hq78D0oX=z@ZR|A*H3M2*V!zl^+us19Ib
z#KaoFGyDGiV=!rEegh;Ar&NkUq)GMR$-)wW&$;)vx5_MRTCe}&u)uUc6Ua=Ybno8Y
z*Q8yCV?Kn085A|gUn>CJ_x%HIbalvi5Py{n;iM9FS)IJ+%EY8`9hbut$^)2D_BS>l
z4Kz15w+fnV5)cxCI{M9m7mPPou3Q0leokAkAxg16Q+J~sBk$@}yB=I(BHS66?3
ze=@j*4g{{K#{7#^`BJAIvhc&_rnQPHDrqon0p0}&dMXjO5|if3v;zNPyp>Z@TG=;W
zY@V19zkT}wNQMjyzdr?GwstOqLb~ZD0_FvF;-VrcmuN^1Xf18yHrMF6m=}~xO)ZU#
zkO{24>S8ogE{`6~jn{w6*NR$jKIg$QFV%qzZ^49W_?w-Uq>3grEUXj6wfed`dclQX
z$D(y^y$g#s!IGO6H!bdcWGw4)!C+@)ZAVjVlWqHUT|&r^x%qc>s2drjC+tpcn5tga
z4hIA327o23>!r7*5@=0Ibsn%fZ%(&CPbciM;^6AaBuvNsGEFppb0<
void make_buoyancy(Vector< MultiFab > &S_data, const MultiFab &S_prim, MultiFab &buoyancy, const amrex::Geometry geom, const SolverChoice &solverChoice, const MultiFab &base_state, const int n_qstate, const int anelastic)
Definition: ERF_MakeBuoyancy.cpp:32
void NumericalDiffusion_Xmom(const Box &bx, const Real dt, const Real num_diff_coeff, const Array4< const Real > &prim_data, const Array4< const Real > &cell_data, const Array4< Real > &rhs, const Array4< const Real > &mf_arr)
Definition: ERF_NumericalDiffusion.cpp:84
void NumericalDiffusion_Ymom(const Box &bx, const Real dt, const Real num_diff_coeff, const Array4< const Real > &prim_data, const Array4< const Real > &cell_data, const Array4< Real > &rhs, const Array4< const Real > &mf_arr)
Definition: ERF_NumericalDiffusion.cpp:149
AMREX_FORCE_INLINE IntVect offset(const int face_dir, const int normal)
Definition: ERF_ReadBndryPlanes.cpp:28
Definition: ERF_PlaneAverage.H:14
amrex::Vector< amrex::Real > input_sounding_time
Definition: ERF_InputSoundingData.H:315
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > U_inp_sound_d
Definition: ERF_InputSoundingData.H:326
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcing
Definition: ERF_DataStruct.H:670
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:641
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcing
Definition: ERF_DataStruct.H:671
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:642
Here is the call graph for this function: