12 #include <AMReX_AmrCore.H>
13 #include <AMReX_BCRec.H>
14 #include <AMReX_Print.H>
16 #include <AMReX_ParallelDescriptor.H>
17 #include <AMReX_ParmParse.H>
18 #include <AMReX_MultiFabUtil.H>
19 #include <AMReX_FillPatchUtil.H>
20 #include <AMReX_VisMF.H>
21 #include <AMReX_PhysBCFunct.H>
22 #include <AMReX_YAFluxRegister.H>
23 #include <AMReX_ErrorList.H>
24 #include <AMReX_EBFabFactory.H>
25 #include <AMReX_EBMultiFabUtil.H>
30 #include <AMReX_FFT_Poisson.H>
33 #ifdef AMREX_MEM_PROFILING
34 #include <AMReX_MemProfiler.H>
54 #ifdef ERF_USE_PARTICLES
62 #ifdef ERF_USE_WINDFARM
73 #include "ERF_P3Interface.H"
79 #include <AMReX_Lazy.H>
83 using amrex::MPI_COMM_WORLD;
84 using amrex::MPI_Comm;
87 #ifdef ERF_USE_MULTIBLOCK
88 class MultiBlockContainer;
95 FullState, Perturbational
110 :
public amrex::AmrCore
129 ERF& operator=(
ERF&& other) noexcept = delete;
135 ERF& operator=(const
ERF& other) = delete;
145 const
amrex::MultiFab& U_new,
146 const
amrex::MultiFab& V_new,
147 const
amrex::MultiFab& W_new,
149 const
bool is_track_io,
150 amrex::TagBoxArray* tags =
nullptr);
172 amrex::Geometry const& geom_at_lev);
187 void solve_with_fft (
int lev,
const amrex::Box& subdomain,
188 amrex::MultiFab& rhs, amrex::MultiFab& p,
189 amrex::Array<amrex::MultiFab,AMREX_SPACEDIM>& fluxes);
191 void solve_with_EB_mlmg (
int lev, amrex::Vector<amrex::MultiFab>& rhs, amrex::Vector<amrex::MultiFab>& p,
192 amrex::Vector<amrex::Array<amrex::MultiFab,AMREX_SPACEDIM>>& fluxes);
195 amrex::MultiFab& rhs, amrex::MultiFab& p,
196 amrex::Array<amrex::MultiFab,AMREX_SPACEDIM>& fluxes,
197 amrex::MultiFab& ax_sub, amrex::MultiFab& ay_sub,
198 amrex::MultiFab& znd_sub);
200 void solve_with_mlmg (
int lev, amrex::Vector<amrex::MultiFab>& rhs, amrex::Vector<amrex::MultiFab>& p,
201 amrex::Vector<amrex::Array<amrex::MultiFab,AMREX_SPACEDIM>>& fluxes);
203 void ImposeBCsOnPhi (
int lev, amrex::MultiFab& phi,
const amrex::Box& subdomain);
206 amrex::Array<amrex::LinOpBCType,AMREX_SPACEDIM>
234 void FillBdyCCVels (amrex::Vector<amrex::MultiFab>& mf_cc_vel,
int levc=0);
240 amrex::Gpu::HostVector<amrex::Real>& h_avg_u , amrex::Gpu::HostVector<amrex::Real>& h_avg_v,
241 amrex::Gpu::HostVector<amrex::Real>& h_avg_w , amrex::Gpu::HostVector<amrex::Real>& h_avg_rho,
242 amrex::Gpu::HostVector<amrex::Real>& h_avg_th , amrex::Gpu::HostVector<amrex::Real>& h_avg_ksgs,
243 amrex::Gpu::HostVector<amrex::Real>& h_avg_Kmv , amrex::Gpu::HostVector<amrex::Real>& h_avg_Khv,
244 amrex::Gpu::HostVector<amrex::Real>& h_avg_qv , amrex::Gpu::HostVector<amrex::Real>& h_avg_qc,
245 amrex::Gpu::HostVector<amrex::Real>& h_avg_qr ,
246 amrex::Gpu::HostVector<amrex::Real>& h_avg_wqv , amrex::Gpu::HostVector<amrex::Real>& h_avg_wqc,
247 amrex::Gpu::HostVector<amrex::Real>& h_avg_wqr ,
248 amrex::Gpu::HostVector<amrex::Real>& h_avg_qi , amrex::Gpu::HostVector<amrex::Real>& h_avg_qs,
249 amrex::Gpu::HostVector<amrex::Real>& h_avg_qg ,
250 amrex::Gpu::HostVector<amrex::Real>& h_avg_uu , amrex::Gpu::HostVector<amrex::Real>& h_avg_uv,
251 amrex::Gpu::HostVector<amrex::Real>& h_avg_uw,
252 amrex::Gpu::HostVector<amrex::Real>& h_avg_vv , amrex::Gpu::HostVector<amrex::Real>& h_avg_vw,
253 amrex::Gpu::HostVector<amrex::Real>& h_avg_ww,
254 amrex::Gpu::HostVector<amrex::Real>& h_avg_uth , amrex::Gpu::HostVector<amrex::Real>& h_avg_vth,
255 amrex::Gpu::HostVector<amrex::Real>& h_avg_wth, amrex::Gpu::HostVector<amrex::Real>& h_avg_thth,
256 amrex::Gpu::HostVector<amrex::Real>& h_avg_ku, amrex::Gpu::HostVector<amrex::Real>& h_avg_kv,
257 amrex::Gpu::HostVector<amrex::Real>& h_avg_kw,
258 amrex::Gpu::HostVector<amrex::Real>& h_avg_p,
259 amrex::Gpu::HostVector<amrex::Real>& h_avg_pu, amrex::Gpu::HostVector<amrex::Real>& h_avg_pv,
260 amrex::Gpu::HostVector<amrex::Real>& h_avg_pw, amrex::Gpu::HostVector<amrex::Real>& h_avg_wthv);
262 amrex::Gpu::HostVector<amrex::Real>& h_avg_u , amrex::Gpu::HostVector<amrex::Real>& h_avg_v,
263 amrex::Gpu::HostVector<amrex::Real>& h_avg_w , amrex::Gpu::HostVector<amrex::Real>& h_avg_rho,
264 amrex::Gpu::HostVector<amrex::Real>& h_avg_th , amrex::Gpu::HostVector<amrex::Real>& h_avg_ksgs,
265 amrex::Gpu::HostVector<amrex::Real>& h_avg_Kmv , amrex::Gpu::HostVector<amrex::Real>& h_avg_Khv,
266 amrex::Gpu::HostVector<amrex::Real>& h_avg_qv , amrex::Gpu::HostVector<amrex::Real>& h_avg_qc,
267 amrex::Gpu::HostVector<amrex::Real>& h_avg_qr ,
268 amrex::Gpu::HostVector<amrex::Real>& h_avg_wqv , amrex::Gpu::HostVector<amrex::Real>& h_avg_wqc,
269 amrex::Gpu::HostVector<amrex::Real>& h_avg_wqr ,
270 amrex::Gpu::HostVector<amrex::Real>& h_avg_qi , amrex::Gpu::HostVector<amrex::Real>& h_avg_qs,
271 amrex::Gpu::HostVector<amrex::Real>& h_avg_qg ,
272 amrex::Gpu::HostVector<amrex::Real>& h_avg_uu , amrex::Gpu::HostVector<amrex::Real>& h_avg_uv,
273 amrex::Gpu::HostVector<amrex::Real>& h_avg_uw,
274 amrex::Gpu::HostVector<amrex::Real>& h_avg_vv , amrex::Gpu::HostVector<amrex::Real>& h_avg_vw,
275 amrex::Gpu::HostVector<amrex::Real>& h_avg_ww,
276 amrex::Gpu::HostVector<amrex::Real>& h_avg_uth , amrex::Gpu::HostVector<amrex::Real>& h_avg_vth,
277 amrex::Gpu::HostVector<amrex::Real>& h_avg_wth, amrex::Gpu::HostVector<amrex::Real>& h_avg_thth,
278 amrex::Gpu::HostVector<amrex::Real>& h_avg_ku, amrex::Gpu::HostVector<amrex::Real>& h_avg_kv,
279 amrex::Gpu::HostVector<amrex::Real>& h_avg_kw,
280 amrex::Gpu::HostVector<amrex::Real>& h_avg_p,
281 amrex::Gpu::HostVector<amrex::Real>& h_avg_pu, amrex::Gpu::HostVector<amrex::Real>& h_avg_pv,
282 amrex::Gpu::HostVector<amrex::Real>& h_avg_pw, amrex::Gpu::HostVector<amrex::Real>& h_avg_wthv);
284 void derive_stress_profiles (amrex::Gpu::HostVector<amrex::Real>& h_avg_tau11, amrex::Gpu::HostVector<amrex::Real>& h_avg_tau12,
285 amrex::Gpu::HostVector<amrex::Real>& h_avg_tau13, amrex::Gpu::HostVector<amrex::Real>& h_avg_tau22,
286 amrex::Gpu::HostVector<amrex::Real>& h_avg_tau23, amrex::Gpu::HostVector<amrex::Real>& h_avg_tau33,
287 amrex::Gpu::HostVector<amrex::Real>& h_avg_hfx3, amrex::Gpu::HostVector<amrex::Real>& h_avg_q1fx3,
288 amrex::Gpu::HostVector<amrex::Real>& h_avg_q2fx3, amrex::Gpu::HostVector<amrex::Real>& h_avg_diss);
289 void derive_stress_profiles_stag (amrex::Gpu::HostVector<amrex::Real>& h_avg_tau11, amrex::Gpu::HostVector<amrex::Real>& h_avg_tau12,
290 amrex::Gpu::HostVector<amrex::Real>& h_avg_tau13, amrex::Gpu::HostVector<amrex::Real>& h_avg_tau22,
291 amrex::Gpu::HostVector<amrex::Real>& h_avg_tau23, amrex::Gpu::HostVector<amrex::Real>& h_avg_tau33,
292 amrex::Gpu::HostVector<amrex::Real>& h_avg_hfx3, amrex::Gpu::HostVector<amrex::Real>& h_avg_q1fx3,
293 amrex::Gpu::HostVector<amrex::Real>& h_avg_q2fx3, amrex::Gpu::HostVector<amrex::Real>& h_avg_diss);
297 volWgtSumMF (
int lev,
const amrex::MultiFab& mf,
int comp,
bool finemask);
307 const amrex::DistributionMapping& dm)
override;
313 const amrex::DistributionMapping& dm)
override;
323 const amrex::DistributionMapping& dm)
override;
328 #ifdef ERF_USE_WW3_COUPLING
330 void send_to_ww3(
int lev);
332 void read_waves(
int lev);
341 amrex::Vector<amrex::MultiFab>& state_old,
342 amrex::Vector<amrex::MultiFab>& state_new,
343 amrex::MultiFab& xvel_old, amrex::MultiFab& yvel_old, amrex::MultiFab& zvel_old,
344 amrex::MultiFab& xvel_new, amrex::MultiFab& yvel_new, amrex::MultiFab& zvel_new,
345 amrex::MultiFab& source, amrex::MultiFab& xmom_src,
346 amrex::MultiFab& ymom_src, amrex::MultiFab& zmom_src,
347 amrex::MultiFab& buoyancy, amrex::Geometry fine_geom,
351 amrex::MultiFab& cons_in,
353 const int& iteration,
357 amrex::MultiFab& cons_in,
358 amrex::MultiFab& xvel_in,
359 amrex::MultiFab& yvel_in,
363 amrex::MultiFab& cons_in,
367 void compute_shoc_tendencies (
int lev,
368 amrex::MultiFab& cons_in,
369 amrex::MultiFab& xvel_in,
370 amrex::MultiFab& yvel_in,
371 amrex::MultiFab& zvel_in,
372 amrex::MultiFab& source, amrex::MultiFab& xmom_src,
373 amrex::MultiFab& ymom_src, amrex::MultiFab& zmom_src,
378 void compute_p3_tendencies (
int lev,
379 amrex::MultiFab& cons_in,
380 amrex::MultiFab& source,
391 void Write3DPlotFile (
int which, PlotFileType plotfile_type, amrex::Vector<std::string> plot_var_names);
392 void Write2DPlotFile (
int which, PlotFileType plotfile_type, amrex::Vector<std::string> plot_var_names);
399 const amrex::Vector<const amrex::MultiFab*> &mf,
400 const amrex::Vector<const amrex::MultiFab*> &mf_nd,
401 const amrex::Vector<std::string> &varnames,
402 const amrex::Vector<amrex::Geometry>& my_geom,
404 const amrex::Vector<int> &level_steps,
405 const amrex::Vector<amrex::IntVect>& my_ref_ratio,
406 const std::string &versionName =
"HyperCLaw-V1.1",
407 const std::string &levelPrefix =
"Level_",
408 const std::string &mfPrefix =
"Cell",
409 const amrex::Vector<std::string>& extra_dirs = amrex::Vector<std::string>())
const;
414 const amrex::Vector<amrex::BoxArray> &bArray,
415 const amrex::Vector<std::string> &varnames,
416 const amrex::Vector<amrex::Geometry>& my_geom,
418 const amrex::Vector<int> &level_steps,
419 const amrex::Vector<amrex::IntVect>& my_ref_ratio,
420 const std::string &versionName,
421 const std::string &levelPrefix,
422 const std::string &mfPrefix)
const;
425 amrex::MultiFab& dens, amrex::MultiFab&
pres, amrex::MultiFab&
pi,
426 amrex::MultiFab& th, amrex::MultiFab&
qv,
427 std::unique_ptr<amrex::MultiFab>& z_cc);
429 #ifdef ERF_USE_NETCDF
431 void writeToNCColumnFile (
int lev,
442 void init_thin_body (
int lev,
const amrex::BoxArray& ba,
const amrex::DistributionMapping& dm);
445 amrex::Geometry& geom_weather,
446 amrex::BoxArray& nba,
447 amrex::DistributionMapping& dm);
450 const amrex::Geometry& geom_weather,
451 const amrex::BoxArray& nba,
452 const amrex::DistributionMapping& dm,
453 amrex::Vector<amrex::MultiFab>& weather_forecast_data);
456 amrex::MultiFab& weather_forecast_interp);
460 #ifdef ERF_USE_MULTIBLOCK
463 ERF (
const amrex::RealBox& rb,
int max_level_in,
464 const amrex::Vector<int>& n_cell_in,
int coord,
465 const amrex::Vector<amrex::IntVect>& ref_ratio,
466 const amrex::Array<int,AMREX_SPACEDIM>& is_per,
470 void Evolve_MB (
int MBstep,
int max_block_step);
478 void SetParmParsePrefix (std::string name) {
pp_prefix = name; }
481 void SetMultiBlockPointer (MultiBlockContainer *mbc) { m_mbc = mbc; }
484 std::vector<amrex::Box> domain_p;
487 amrex::Vector<amrex::Vector<amrex::MultiFab> >
vars_new;
488 amrex::Vector<amrex::Vector<amrex::MultiFab> >
vars_old;
491 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
vel_t_avg;
500 #ifdef ERF_USE_NETCDF
502 void fill_from_realbdy (
const amrex::Vector<amrex::MultiFab*>& mfs,
507 amrex::IntVect ngvect_cons,
508 amrex::IntVect ngvect_vels);
511 void fill_from_realbdy_upwind (
const amrex::Vector<amrex::MultiFab*>& mfs,
516 amrex::IntVect ngvect_cons,
517 amrex::IntVect ngvect_vels);
520 #ifdef ERF_USE_NETCDF
521 void init_from_wrfinput (
int lev,
526 void init_from_metgrid (
int lev);
527 void init_from_ncfile (
int lev);
530 #ifdef ERF_USE_WINDFARM
531 void init_windfarm(
int lev);
532 void advance_windfarm (
const amrex::Geometry& a_geom,
534 amrex::MultiFab& cons_in,
535 amrex::MultiFab& U_old, amrex::MultiFab& V_old, amrex::MultiFab& W_old,
536 amrex::MultiFab& mf_vars_windfarm,
537 const amrex::MultiFab& mf_Nturb,
538 const amrex::MultiFab& mf_SMark,
589 void init_stuff (
int lev,
const amrex::BoxArray& ba,
const amrex::DistributionMapping& dm,
590 amrex::Vector<amrex::MultiFab>& lev_new, amrex::Vector<amrex::MultiFab>& lev_old,
591 amrex::MultiFab& tmp_base_state,
592 std::unique_ptr<amrex::MultiFab>& tmp_zphys_nd);
607 #ifdef ERF_USE_WINDFARM
609 void initializeWindFarm (
const int&);
618 const amrex::Vector<amrex::MultiFab*>& mfs_vel,
619 bool cons_only=
false);
623 const amrex::Vector<amrex::MultiFab*>& mfs_vel,
624 const amrex::Vector<amrex::MultiFab*>& mfs_mom,
625 const amrex::MultiFab& old_base_state,
626 const amrex::MultiFab& new_base_state,
627 bool fillset=
true,
bool cons_only=
false);
638 const amrex::Vector<amrex::MultiFab*>& mfs_vel,
639 const amrex::Vector<amrex::MultiFab*>& mfs_mom,
640 int ng_cons,
int ng_vel,
bool cons_only,
int icomp_cons,
int ncomp_cons);
676 static amrex::Vector<std::string>
PlotFileVarNames (amrex::Vector<std::string> plot_var_names) ;
679 void setPlotVariables (
const std::string& pp_plot_var_names, amrex::Vector<std::string>& plot_var_names);
681 void appendPlotVariables (
const std::string& pp_plot_var_names, amrex::Vector<std::string>& plot_var_names);
683 #ifdef ERF_USE_NETCDF
685 void createNCColumnFile (
int lev,
689 void init_from_wrfbdy (amrex::Vector<amrex::FArrayBox*> x_vel_lateral,
690 amrex::Vector<amrex::FArrayBox*> y_vel_lateral,
691 amrex::Vector<amrex::FArrayBox*> z_vel_lateral,
692 amrex::Vector<amrex::FArrayBox*> T_lateral);
710 amrex::Vector<amrex::Vector<amrex::FArrayBox>> bdy_data_xlo;
711 amrex::Vector<amrex::Vector<amrex::FArrayBox>> bdy_data_xhi;
712 amrex::Vector<amrex::Vector<amrex::FArrayBox>> bdy_data_ylo;
713 amrex::Vector<amrex::Vector<amrex::FArrayBox>> bdy_data_yhi;
718 amrex::Vector<amrex::FArrayBox> low_data_zlo;
722 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
lat_m,
lon_m;
736 amrex::Vector<amrex::Gpu::DeviceVector<amrex::Real>>
th_bc_data;
754 void fill_rhs (amrex::MultiFab& rhs_mf,
const amrex::MultiFab& state_mf,
amrex::Real time,
const amrex::Geometry& geom);
759 std::unique_ptr<ProblemBase>
prob =
nullptr;
771 amrex::Vector<amrex::Real>
dt;
774 #ifndef ERF_USE_MULTIBLOCK
777 amrex::Vector<amrex::Vector<amrex::MultiFab> >
vars_new;
778 amrex::Vector<amrex::Vector<amrex::MultiFab> >
vars_old;
781 amrex::Vector<amrex::Vector<amrex::MultiFab> >
gradp;
784 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
vel_t_avg;
787 amrex::Vector<std::unique_ptr<MRISplitIntegrator<amrex::Vector<amrex::MultiFab> > > >
mri_integrator_mem;
793 amrex::Vector<std::unique_ptr<ERFPhysBCFunct_u>>
physbcs_u;
794 amrex::Vector<std::unique_ptr<ERFPhysBCFunct_v>>
physbcs_v;
795 amrex::Vector<std::unique_ptr<ERFPhysBCFunct_w>>
physbcs_w;
800 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
Qv_prim;
801 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
Qr_prim;
815 std::unique_ptr<Microphysics>
micro;
816 amrex::Vector<amrex::Vector<amrex::MultiFab*>>
qmoist;
820 #ifdef ERF_USE_WINDFARM
821 std::unique_ptr<WindFarm> windfarm;
822 amrex::Vector<amrex::MultiFab> Nturb;
823 amrex::Vector<amrex::MultiFab> vars_windfarm;
826 amrex::Vector<amrex::MultiFab> SMark;
832 amrex::Vector<amrex::Vector<amrex::MultiFab*>>
lsm_data;
833 amrex::Vector<amrex::Vector<amrex::MultiFab*>>
lsm_flux;
835 amrex::Vector<std::unique_ptr<IRadiation>>
rad;
839 amrex::Vector<std::unique_ptr<SHOCInterface>> shoc_interface;
843 amrex::Vector<std::unique_ptr<P3Interface>> p3_interface;
856 amrex::Vector<ERFFillPatcher>
FPr_c;
857 amrex::Vector<ERFFillPatcher>
FPr_u;
858 amrex::Vector<ERFFillPatcher>
FPr_v;
859 amrex::Vector<ERFFillPatcher>
FPr_w;
862 amrex::Vector<amrex::Vector<std::unique_ptr<amrex::MultiFab>>>
Tau;
867 amrex::Vector<amrex::Vector<std::unique_ptr<amrex::MultiFab>>>
sst_lev;
868 amrex::Vector<amrex::Vector<std::unique_ptr<amrex::MultiFab>>>
tsk_lev;
869 amrex::Vector<amrex::Vector<std::unique_ptr<amrex::iMultiFab>>>
lmask_lev;
881 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
z_phys_nd;
882 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
z_phys_cc;
884 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
detJ_cc;
885 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
ax;
886 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
ay;
887 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
az;
892 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
ax_src;
893 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
ay_src;
894 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
az_src;
899 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
z_t_rk;
905 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
walldist;
908 amrex::Vector<amrex::Vector<std::unique_ptr<amrex::MultiFab>>>
mapfac;
917 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
Hwave;
918 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
Lwave;
1041 const amrex::Vector<std::string>
cons_names {
"density",
"rhotheta",
"rhoKE",
"rhoadv_0",
1042 "rhoQ1",
"rhoQ2",
"rhoQ3",
1043 "rhoQ4",
"rhoQ5",
"rhoQ6"};
1048 const amrex::Vector<std::string>
derived_names {
"soundspeed",
"temp",
"theta",
"KE",
"scalar",
1049 "vorticity_x",
"vorticity_y",
"vorticity_z",
1051 "pres_hse",
"dens_hse",
"theta_hse",
"pressure",
"pert_pres",
"pert_dens",
1053 #ifdef ERF_USE_WINDFARM
1054 "num_turb",
"SMark0",
"SMark1",
1056 "dpdx",
"dpdy",
"dpdz",
"pres_hse_x",
"pres_hse_y",
1057 "z_phys",
"detJ",
"mapfac",
"lat_m",
"lon_m",
1059 "u_t_avg",
"v_t_avg",
"w_t_avg",
"umag_t_avg",
1073 "moist_density",
"qv",
"qc",
"qi",
"qrain",
"qsnow",
"qgraup",
"qt",
"qn",
"qp",
"qsat",
1074 "rain_accum",
"snow_accum",
"graup_accum",
"reflectivity",
"max_reflectivity",
1079 #ifdef ERF_COMPUTE_ERROR
1081 ,
"xvel_err",
"yvel_err",
"zvel_err",
"pp_err"
1083 ,
"qsrc_sw",
"qsrc_lw"
1092 #ifdef ERF_USE_PARTICLES
1094 ParticleData particleData;
1097 bool m_use_tracer_particles;
1100 void readTracersParams();
1103 void initializeTracers ( amrex::ParGDBBase*,
1104 const amrex::Vector<std::unique_ptr<amrex::MultiFab>>&,
1108 void restartTracers ( amrex::ParGDBBase*,
const std::string& );
1111 void evolveTracers(
int,
1113 amrex::Vector<amrex::Vector<amrex::MultiFab>>&,
1114 const amrex::Vector<std::unique_ptr<amrex::MultiFab>>& );
1118 #ifdef ERF_USE_MULTIBLOCK
1119 MultiBlockContainer *m_mbc =
nullptr;
1169 amrex::Vector<amrex::BoxArray>
ba1d;
1170 amrex::Vector<amrex::BoxArray>
ba2d;
1175 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
mf_PSFC;
1205 amrex::Vector<amrex::Gpu::DeviceVector<amrex::Real> >
d_w_subsid;
1208 amrex::Vector<amrex::Gpu::DeviceVector<amrex::Real> >
d_u_geos;
1211 amrex::Vector<amrex::Gpu::DeviceVector<amrex::Real> >
d_v_geos;
1215 amrex::Vector<amrex::Real>& u_geos,
1216 amrex::Gpu::DeviceVector<amrex::Real>& u_geos_d,
1217 amrex::Vector<amrex::Real>& v_geos,
1218 amrex::Gpu::DeviceVector<amrex::Real>& v_geos_d,
1219 const amrex::Geometry& lgeom,
1220 const amrex::Vector<amrex::Real>& zlev_stag);
1232 amrex::Vector<amrex::Vector<amrex::Gpu::DeviceVector<amrex::Real> > >
d_sponge_ptrs;
1248 std::unique_ptr<WriteBndryPlanes>
m_w2d =
nullptr;
1249 std::unique_ptr<ReadBndryPlanes>
m_r2d =
nullptr;
1268 static AMREX_FORCE_INLINE
1422 int numCores = amrex::ParallelDescriptor::NProcs();
1424 numCores = numCores * omp_get_max_threads();
1428 numCores * (amrex::ParallelDescriptor::second() -
startCPUTime) +
1436 if (amrex::ParallelDescriptor::IOProcessor())
1438 datalog[i] = std::make_unique<std::fstream>();
1439 datalog[i]->open(filename.c_str(),std::ios::out|std::ios::app);
1441 amrex::FileOpenFailed(filename);
1444 amrex::ParallelDescriptor::Barrier(
"ERF::setRecordDataInfo");
1449 if (amrex::ParallelDescriptor::IOProcessor())
1451 der_datalog[i] = std::make_unique<std::fstream>();
1452 der_datalog[i]->open(filename.c_str(),std::ios::out|std::ios::app);
1454 amrex::FileOpenFailed(filename);
1457 amrex::ParallelDescriptor::Barrier(
"ERF::setRecordDerDataInfo");
1462 if (amrex::ParallelDescriptor::IOProcessor())
1465 tot_e_datalog[i]->open(filename.c_str(),std::ios::out|std::ios::app);
1467 amrex::FileOpenFailed(filename);
1470 amrex::ParallelDescriptor::Barrier(
"ERF::setRecordEnergyDataInfo");
1475 amrex::MultiFab dummy(grids[lev],dmap[lev],1,0);
1476 for (amrex::MFIter mfi(dummy); mfi.isValid(); ++mfi)
1478 const amrex::Box& bx = mfi.validbox();
1479 if (bx.contains(cell)) {
1480 sampleptlog[i] = std::make_unique<std::fstream>();
1481 sampleptlog[i]->open(filename.c_str(),std::ios::out|std::ios::app);
1483 amrex::FileOpenFailed(filename);
1487 amrex::ParallelDescriptor::Barrier(
"ERF::setRecordSamplePointInfo");
1492 amrex::MultiFab dummy(grids[lev],dmap[lev],1,0);
1493 for (amrex::MFIter mfi(dummy); mfi.isValid(); ++mfi)
1495 const amrex::Box& bx = mfi.validbox();
1496 if (bx.contains(cell)) {
1498 samplelinelog[i]->open(filename.c_str(),std::ios::out|std::ios::app);
1500 amrex::FileOpenFailed(filename);
1504 amrex::ParallelDescriptor::Barrier(
"ERF::setRecordSampleLineInfo");
1512 amrex::Vector<std::unique_ptr<std::fstream> >
datalog;
1538 amrex::Vector<std::unique_ptr<eb_>>
eb;
1541 AMREX_ASSERT(lev >= 0 && lev <
eb.size() &&
eb[lev] !=
nullptr);
1545 [[nodiscard]] amrex::EBFArrayBoxFactory
const&
1547 return *(
eb[lev]->get_const_factory());
1561 amrex::Vector<std::unique_ptr<amrex::FFT::Poisson<amrex::MultiFab>>> m_3D_poisson;
1562 amrex::Vector<std::unique_ptr<amrex::FFT::PoissonHybrid<amrex::MultiFab>>> m_2D_poisson;
1570 static void print_usage(MPI_Comm , std::ostream& );
1571 static void print_error(MPI_Comm ,
const std::string& msg);
AMREX_ENUM(StateInterpType, FullState, Perturbational)
Contains the Eulerian microphysics class.
#define NBCVAR_max
Definition: ERF_IndexDefines.H:29
Contains the Lagrangian microphysics class.
amrex::Real Real
Definition: ERF_ShocInterface.H:16
void MakeHorizontalAverages()
Definition: ERF.cpp:2219
amrex::Vector< amrex::MultiFab > rU_new
Definition: ERF.H:805
int last_check_file_step
Definition: ERF.H:961
amrex::Vector< std::unique_ptr< amrex::MultiFab > > walldist
Definition: ERF.H:905
void initRayleigh()
Initialize Rayleigh damping profiles.
Definition: ERF_InitRayleigh.cpp:14
static amrex::Real start_time
Definition: ERF.H:977
bool metgrid_basic_linear
Definition: ERF.H:1161
amrex::Vector< amrex::Vector< amrex::Box > > boxes_at_level
Definition: ERF.H:763
amrex::Vector< std::string > samplelinelogname
Definition: ERF.H:1524
int max_step
Definition: ERF.H:974
bool metgrid_debug_msf
Definition: ERF.H:1159
AMREX_FORCE_INLINE std::ostream & DerDataLog(int i)
Definition: ERF.H:1342
amrex::Vector< amrex::Vector< std::unique_ptr< amrex::MultiFab > > > mapfac
Definition: ERF.H:908
void FillBdyCCVels(amrex::Vector< amrex::MultiFab > &mf_cc_vel, int levc=0)
Definition: ERF_FillBdyCCVels.cpp:11
void setRayleighRefFromSounding(bool restarting)
Set Rayleigh mean profiles from input sounding.
Definition: ERF_InitRayleigh.cpp:55
amrex::Vector< std::unique_ptr< MRISplitIntegrator< amrex::Vector< amrex::MultiFab > > > > mri_integrator_mem
Definition: ERF.H:787
void Evolve()
Definition: ERF.cpp:477
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > d_rhotheta_src
Definition: ERF.H:1199
amrex::Vector< amrex::MultiFab > pp_inc
Definition: ERF.H:789
amrex::Vector< ERFFillPatcher > FPr_u
Definition: ERF.H:857
amrex::Real cloud_fraction(amrex::Real time)
Definition: ERF_WriteScalarProfiles.cpp:451
amrex::Vector< amrex::IntVect > sampleline
Definition: ERF.H:1525
amrex::Vector< std::unique_ptr< amrex::MultiFab > > SFS_q1fx3_lev
Definition: ERF.H:874
amrex::Vector< std::unique_ptr< amrex::MultiFab > > Hwave_onegrid
Definition: ERF.H:919
amrex::Vector< amrex::Vector< amrex::MultiFab > > vars_new
Definition: ERF.H:777
ERF(ERF &&) noexcept=delete
AMREX_FORCE_INLINE int NumSampleLineLogs() noexcept
Definition: ERF.H:1385
static void print_tpls(std::ostream &)
Definition: ERF_ConsoleIO.cpp:137
amrex::Vector< amrex::Real > dz_min
Definition: ERF.H:1266
amrex::Real m_plot2d_per_1
Definition: ERF.H:1021
std::string DataLogName(int i) const noexcept
The filename of the ith datalog file.
Definition: ERF.H:1528
std::string plot2d_file_2
Definition: ERF.H:1011
amrex::Vector< std::unique_ptr< amrex::MultiFab > > thin_yforce
Definition: ERF.H:952
void ErrorEst(int lev, amrex::TagBoxArray &tags, amrex::Real time, int ngrow) override
Definition: ERF_Tagging.cpp:20
amrex::Vector< std::unique_ptr< amrex::MultiFab > > SFS_hfx3_lev
Definition: ERF.H:872
amrex::Vector< std::unique_ptr< amrex::MultiFab > > sw_lw_fluxes
Definition: ERF.H:847
amrex::Vector< amrex::Vector< amrex::Real > > h_w_subsid
Definition: ERF.H:1204
static void print_banner(MPI_Comm, std::ostream &)
Definition: ERF_ConsoleIO.cpp:60
amrex::Vector< std::unique_ptr< std::fstream > > datalog
Definition: ERF.H:1512
static amrex::Real sum_per
Definition: ERF.H:1129
amrex::Vector< ERFFillPatcher > FPr_v
Definition: ERF.H:858
int cf_set_width
Definition: ERF.H:855
static amrex::Real previousCPUTimeUsed
Definition: ERF.H:1417
void setPlotVariables(const std::string &pp_plot_var_names, amrex::Vector< std::string > &plot_var_names)
Definition: ERF_Plotfile.cpp:23
amrex::Gpu::DeviceVector< amrex::Real > d_havg_temperature
Definition: ERF.H:1241
amrex::Vector< amrex::Vector< amrex::Vector< amrex::Real > > > h_sponge_ptrs
Definition: ERF.H:1226
int last_plot2d_file_step_2
Definition: ERF.H:958
void fill_from_bndryregs(const amrex::Vector< amrex::MultiFab * > &mfs, amrex::Real time)
Definition: ERF_BoundaryConditionsBndryReg.cpp:13
const int timeprecision
Definition: ERF.H:967
int m_subvol_int
Definition: ERF.H:1018
void setRecordDataInfo(int i, const std::string &filename)
Definition: ERF.H:1434
amrex::Vector< std::unique_ptr< amrex::MultiFab > > SFS_hfx1_lev
Definition: ERF.H:872
amrex::Vector< amrex::BoxArray > ba2d
Definition: ERF.H:1170
void HurricaneTracker(int lev, const amrex::MultiFab &U_new, const amrex::MultiFab &V_new, const amrex::MultiFab &W_new, const amrex::Real velmag_threshold, const bool is_track_io, amrex::TagBoxArray *tags=nullptr)
Definition: ERF_Tagging.cpp:455
amrex::Array< amrex::Array< amrex::Real, AMREX_SPACEDIM *2 >, AMREX_SPACEDIM+NBCVAR_max > m_bc_extdir_vals
Definition: ERF.H:936
void init_from_input_sounding(int lev)
Definition: ERF_InitFromInputSounding.cpp:52
void erf_enforce_hse(int lev, amrex::MultiFab &dens, amrex::MultiFab &pres, amrex::MultiFab &pi, amrex::MultiFab &th, amrex::MultiFab &qv, std::unique_ptr< amrex::MultiFab > &z_cc)
Definition: ERF_Init1D.cpp:160
amrex::Vector< amrex::Vector< amrex::MultiFab > > gradp
Definition: ERF.H:781
std::string plot3d_file_1
Definition: ERF.H:1008
void FillPatchCrseLevel(int lev, amrex::Real time, const amrex::Vector< amrex::MultiFab * > &mfs_vel, bool cons_only=false)
Definition: ERF_FillPatch.cpp:282
eb_ const & get_eb(int lev) const noexcept
Definition: ERF.H:1540
amrex::Gpu::DeviceVector< amrex::Real > d_havg_qv
Definition: ERF.H:1243
static amrex::Real column_loc_y
Definition: ERF.H:1182
static int mg_verbose
Definition: ERF.H:1123
void ReadParameters()
Definition: ERF.cpp:1803
static amrex::Vector< std::string > PlotFileVarNames(amrex::Vector< std::string > plot_var_names)
Definition: ERF_Plotfile.cpp:251
void InitializeFromFile()
amrex::Vector< std::unique_ptr< ERFPhysBCFunct_cons > > physbcs_cons
Definition: ERF.H:792
void project_velocity_tb(int lev, amrex::Real dt, amrex::Vector< amrex::MultiFab > &vars)
Definition: ERF_PoissonSolve_tb.cpp:20
amrex::Vector< std::unique_ptr< amrex::MultiFab > > mf_PSFC
Definition: ERF.H:1175
ERF()
Definition: ERF.cpp:104
void init_Dirichlet_bc_data(const std::string input_file)
Definition: ERF_InitBCs.cpp:652
amrex::Vector< std::unique_ptr< amrex::MultiFab > > z_phys_nd_src
Definition: ERF.H:889
amrex::Vector< std::unique_ptr< amrex::MultiFab > > detJ_cc
Definition: ERF.H:884
void WriteMyEBSurface()
Definition: ERF_EBWriteSurface.cpp:5
amrex::Vector< std::unique_ptr< amrex::MultiFab > > z_phys_cc
Definition: ERF.H:882
amrex::Vector< std::unique_ptr< eb_ > > eb
Definition: ERF.H:1538
amrex::Vector< std::unique_ptr< amrex::MultiFab > > eddyDiffs_lev
Definition: ERF.H:863
static SolverChoice solverChoice
Definition: ERF.H:1087
amrex::Vector< ERFFillPatcher > FPr_c
Definition: ERF.H:856
bool plot_rad
Definition: ERF.H:850
static bool use_fft
Definition: ERF.H:1124
bool m_plot_face_vels
Definition: ERF.H:1024
amrex::Vector< amrex::MultiFab > base_state_new
Definition: ERF.H:914
std::string plot3d_file_2
Definition: ERF.H:1009
void WriteSubvolume()
Definition: ERF_WriteSubvolume.cpp:9
amrex::Vector< std::unique_ptr< amrex::MultiFab > > az
Definition: ERF.H:887
int regrid_int
Definition: ERF.H:1001
amrex::Vector< amrex::Real > fixed_dt
Definition: ERF.H:995
void write_1D_profiles_stag(amrex::Real time)
Definition: ERF_Write1DProfiles_stag.cpp:25
void WriteGenericPlotfileHeaderWithTerrain(std::ostream &HeaderFile, int nlevels, const amrex::Vector< amrex::BoxArray > &bArray, const amrex::Vector< std::string > &varnames, const amrex::Vector< amrex::Geometry > &my_geom, amrex::Real time, const amrex::Vector< int > &level_steps, const amrex::Vector< amrex::IntVect > &my_ref_ratio, const std::string &versionName, const std::string &levelPrefix, const std::string &mfPrefix) const
Definition: ERF_Plotfile.cpp:1770
amrex::Vector< amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > > d_sponge_ptrs
Definition: ERF.H:1232
amrex::Vector< amrex::Vector< amrex::Real > > h_rhoqt_src
Definition: ERF.H:1201
amrex::Vector< long > dt_mri_ratio
Definition: ERF.H:772
amrex::Vector< amrex::Vector< std::unique_ptr< amrex::MultiFab > > > Tau
Definition: ERF.H:862
amrex::Vector< std::unique_ptr< amrex::MultiFab > > vel_t_avg
Definition: ERF.H:784
void sum_energy_quantities(amrex::Real time)
Definition: ERF_WriteScalarProfiles.cpp:318
amrex::Vector< std::unique_ptr< amrex::MultiFab > > SFS_q2fx3_lev
Definition: ERF.H:875
static amrex::Real dt_max
Definition: ERF.H:992
amrex::Vector< amrex::Vector< std::unique_ptr< amrex::iMultiFab > > > lmask_lev
Definition: ERF.H:869
void CreateWeatherDataGeomBoxArrayDistMap(const std::string &filename, amrex::Geometry &geom_weather, amrex::BoxArray &nba, amrex::DistributionMapping &dm)
Definition: ERF_WeatherDataInterpolation.cpp:155
AMREX_FORCE_INLINE int NumSamplePointLogs() noexcept
Definition: ERF.H:1371
amrex::Vector< amrex::Real > h_havg_pressure
Definition: ERF.H:1236
void update_diffusive_arrays(int lev, const amrex::BoxArray &ba, const amrex::DistributionMapping &dm)
Definition: ERF_MakeNewArrays.cpp:428
static int verbose
Definition: ERF.H:1122
amrex::Vector< amrex::Real > h_havg_qc
Definition: ERF.H:1238
amrex::Vector< std::unique_ptr< ERFPhysBCFunct_w > > physbcs_w
Definition: ERF.H:795
static int nghost_eb_basic()
Definition: ERF.H:1550
void Advance(int lev, amrex::Real time, amrex::Real dt_lev, int iteration, int ncycle)
Definition: ERF_Advance.cpp:23
static std::string column_file_name
Definition: ERF.H:1183
amrex::Vector< std::unique_ptr< std::fstream > > samplelinelog
Definition: ERF.H:1523
amrex::Vector< std::unique_ptr< amrex::MultiFab > > terrain_blanking
Definition: ERF.H:902
std::unique_ptr< Microphysics > micro
Definition: ERF.H:815
int m_plot2d_int_2
Definition: ERF.H:1017
amrex::Vector< amrex::MultiFab > base_state
Definition: ERF.H:913
AMREX_FORCE_INLINE amrex::YAFluxRegister * getAdvFluxReg(int lev)
Definition: ERF.H:1328
int m_plot3d_int_1
Definition: ERF.H:1014
amrex::Vector< amrex::Real > h_havg_density
Definition: ERF.H:1234
amrex::Vector< std::unique_ptr< amrex::MultiFab > > z_phys_nd_new
Definition: ERF.H:896
amrex::Vector< amrex::Real > fixed_fast_dt
Definition: ERF.H:996
bool metgrid_retain_sfc
Definition: ERF.H:1164
amrex::Vector< std::unique_ptr< amrex::MultiFab > > Qv_prim
Definition: ERF.H:800
static int sum_interval
Definition: ERF.H:1127
static int pert_interval
Definition: ERF.H:1128
void restart()
Definition: ERF.cpp:1629
int last_plot3d_file_step_2
Definition: ERF.H:956
amrex::Vector< std::unique_ptr< amrex::MultiFab > > SFS_q1fx2_lev
Definition: ERF.H:874
amrex::IntVect & SampleLine(int i)
Definition: ERF.H:1404
amrex::Vector< amrex::MultiFab > rV_new
Definition: ERF.H:807
std::string PlotFileName(int lev) const
void write_1D_profiles(amrex::Real time)
Definition: ERF_Write1DProfiles.cpp:17
void initialize_integrator(int lev, amrex::MultiFab &cons_mf, amrex::MultiFab &vel_mf)
Definition: ERF_MakeNewArrays.cpp:634
amrex::Vector< std::unique_ptr< amrex::MultiFab > > thin_zforce
Definition: ERF.H:953
amrex::Vector< std::string > plot3d_var_names_2
Definition: ERF.H:1038
amrex::Vector< amrex::BCRec > domain_bcs_type
Definition: ERF.H:929
amrex::Array< amrex::LinOpBCType, AMREX_SPACEDIM > get_projection_bc(amrex::Orientation::Side side) const noexcept
Definition: ERF_SolveWithMLMG.cpp:17
amrex::Vector< std::unique_ptr< amrex::MultiFab > > Qr_prim
Definition: ERF.H:801
void post_timestep(int nstep, amrex::Real time, amrex::Real dt_lev)
Definition: ERF.cpp:586
amrex::Real m_subvol_per
Definition: ERF.H:1023
std::unique_ptr< amrex::MultiFab > mf_MUB
Definition: ERF.H:1173
std::string pp_prefix
Definition: ERF.H:495
std::string SampleLineLogName(int i) const noexcept
The filename of the ith samplelinelog file.
Definition: ERF.H:1535
amrex::Vector< amrex::Vector< std::unique_ptr< amrex::MultiFab > > > sst_lev
Definition: ERF.H:867
amrex::Vector< std::string > plot2d_var_names_1
Definition: ERF.H:1039
amrex::Vector< std::unique_ptr< amrex::MultiFab > > thin_xforce
Definition: ERF.H:951
AMREX_FORCE_INLINE int NumSamplePoints() noexcept
Definition: ERF.H:1398
std::unique_ptr< amrex::MultiFab > mf_C2H
Definition: ERF.H:1172
bool metgrid_use_sfc
Definition: ERF.H:1163
amrex::Real m_plot2d_per_2
Definition: ERF.H:1022
static amrex::Real bndry_output_planes_per
Definition: ERF.H:1188
amrex::Real m_check_per
Definition: ERF.H:1035
void init_custom(int lev)
Definition: ERF_InitCustom.cpp:26
amrex::Vector< std::unique_ptr< IRadiation > > rad
Definition: ERF.H:835
std::unique_ptr< ProblemBase > prob
Definition: ERF.H:759
amrex::Vector< int > num_files_at_level
Definition: ERF.H:762
void WriteMultiLevelPlotfileWithTerrain(const std::string &plotfilename, int nlevels, const amrex::Vector< const amrex::MultiFab * > &mf, const amrex::Vector< const amrex::MultiFab * > &mf_nd, const amrex::Vector< std::string > &varnames, const amrex::Vector< amrex::Geometry > &my_geom, amrex::Real time, const amrex::Vector< int > &level_steps, const amrex::Vector< amrex::IntVect > &my_ref_ratio, const std::string &versionName="HyperCLaw-V1.1", const std::string &levelPrefix="Level_", const std::string &mfPrefix="Cell", const amrex::Vector< std::string > &extra_dirs=amrex::Vector< std::string >()) const
Definition: ERF_Plotfile.cpp:1683
void init_bcs()
Definition: ERF_InitBCs.cpp:20
int profile_int
Definition: ERF.H:1029
bool metgrid_debug_quiescent
Definition: ERF.H:1155
amrex::Vector< std::unique_ptr< ERFPhysBCFunct_u > > physbcs_u
Definition: ERF.H:793
amrex::Vector< amrex::Real > t_new
Definition: ERF.H:769
bool destag_profiles
Definition: ERF.H:1030
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > stretched_dz_d
Definition: ERF.H:911
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > th_bc_data
Definition: ERF.H:736
amrex::Vector< amrex::Real > t_avg_cnt
Definition: ERF.H:785
void FillPatchFineLevel(int lev, amrex::Real time, const amrex::Vector< amrex::MultiFab * > &mfs_vel, const amrex::Vector< amrex::MultiFab * > &mfs_mom, const amrex::MultiFab &old_base_state, const amrex::MultiFab &new_base_state, bool fillset=true, bool cons_only=false)
Definition: ERF_FillPatch.cpp:20
std::string DerDataLogName(int i) const noexcept
Definition: ERF.H:1529
void remake_zphys(int lev, amrex::Real time, std::unique_ptr< amrex::MultiFab > &temp_zphys_nd)
Definition: ERF_MakeNewArrays.cpp:589
int m_check_int
Definition: ERF.H:1034
amrex::Vector< std::unique_ptr< amrex::MultiFab > > solar_zenith
Definition: ERF.H:848
amrex::Real estTimeStep(int lev, long &dt_fast_ratio) const
Definition: ERF_ComputeTimestep.cpp:54
static std::string sponge_type
Definition: ERF.H:1140
static amrex::Real startCPUTime
Definition: ERF.H:1416
amrex::Vector< std::unique_ptr< amrex::MultiFab > > SFS_diss_lev
Definition: ERF.H:873
amrex::Vector< amrex::MultiFab > rU_old
Definition: ERF.H:804
void FillIntermediatePatch(int lev, amrex::Real time, const amrex::Vector< amrex::MultiFab * > &mfs_vel, const amrex::Vector< amrex::MultiFab * > &mfs_mom, int ng_cons, int ng_vel, bool cons_only, int icomp_cons, int ncomp_cons)
Definition: ERF_FillIntermediatePatch.cpp:28
amrex::Vector< amrex::Real > t_old
Definition: ERF.H:770
AMREX_FORCE_INLINE int NumSampleLines() noexcept
Definition: ERF.H:1411
amrex::Vector< std::unique_ptr< amrex::MultiFab > > Theta_prim
Definition: ERF.H:799
amrex::Vector< std::unique_ptr< amrex::MultiFab > > z_t_rk
Definition: ERF.H:899
void make_subdomains(const amrex::BoxList &ba, amrex::Vector< amrex::BoxArray > &bins)
Definition: ERF_MakeSubdomains.cpp:6
void compute_divergence(int lev, amrex::MultiFab &rhs, amrex::Array< amrex::MultiFab const *, AMREX_SPACEDIM > rho0_u_const, amrex::Geometry const &geom_at_lev)
Definition: ERF_ComputeDivergence.cpp:10
std::string MakeVTKFilename(int nstep)
Definition: ERF_Write1DProfiles.cpp:574
void fill_rhs(amrex::MultiFab &rhs_mf, const amrex::MultiFab &state_mf, amrex::Real time, const amrex::Geometry &geom)
void appendPlotVariables(const std::string &pp_plot_var_names, amrex::Vector< std::string > &plot_var_names)
Definition: ERF_Plotfile.cpp:184
amrex::Vector< std::unique_ptr< ERFPhysBCFunct_v > > physbcs_v
Definition: ERF.H:794
static amrex::Real column_per
Definition: ERF.H:1180
amrex::Vector< std::string > tot_e_datalogname
Definition: ERF.H:1517
static amrex::Real stop_time
Definition: ERF.H:978
static int output_bndry_planes
Definition: ERF.H:1186
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > d_v_geos
Definition: ERF.H:1211
void update_terrain_arrays(int lev)
Definition: ERF_MakeNewArrays.cpp:623
static std::string nc_bdy_file
Definition: ERF.H:1146
bool metgrid_interp_theta
Definition: ERF.H:1160
static amrex::Vector< amrex::Vector< std::string > > nc_init_file
Definition: ERF.H:1143
void init_only(int lev, amrex::Real time)
Definition: ERF.cpp:1672
static int input_bndry_planes
Definition: ERF.H:1192
static StateInterpType interpolation_type
Definition: ERF.H:1137
amrex::Gpu::DeviceVector< amrex::Real > d_havg_qc
Definition: ERF.H:1244
void AverageDown()
Definition: ERF_AverageDown.cpp:16
amrex::Vector< amrex::Vector< amrex::Real > > h_v_geos
Definition: ERF.H:1210
bool projection_has_dirichlet(amrex::Array< amrex::LinOpBCType, AMREX_SPACEDIM > bcs) const
Definition: ERF_PoissonSolve_tb.cpp:8
bool regrid_level_0_on_restart
Definition: ERF.H:1005
amrex::Vector< std::unique_ptr< amrex::MultiFab > > Lwave_onegrid
Definition: ERF.H:920
InputSoundingData input_sounding_data
Definition: ERF.H:727
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > d_rhoqt_src
Definition: ERF.H:1202
std::unique_ptr< SampleData > data_sampler
Definition: ERF.H:1510
amrex::MultiFab fine_mask
Definition: ERF.H:1264
void Write2DPlotFile(int which, PlotFileType plotfile_type, amrex::Vector< std::string > plot_var_names)
Definition: ERF_Plotfile.cpp:1859
amrex::Gpu::DeviceVector< amrex::Real > d_havg_density
Definition: ERF.H:1240
int last_subvol
Definition: ERF.H:959
void init_from_hse(int lev)
Definition: ERF_InitFromHSE.cpp:32
const std::string datetime_format
Definition: ERF.H:981
amrex::Vector< amrex::YAFluxRegister * > advflux_reg
Definition: ERF.H:924
int metgrid_force_sfc_k
Definition: ERF.H:1167
amrex::Vector< amrex::Vector< amrex::Real > > h_rhotheta_src
Definition: ERF.H:1198
static int ng_pres_hse
Definition: ERF.H:1195
static amrex::Real bndry_output_planes_start_time
Definition: ERF.H:1189
bool real_extrap_w
Definition: ERF.H:1149
static amrex::Real cfl
Definition: ERF.H:987
void project_velocity(int lev, amrex::Real dt)
Definition: ERF_PoissonSolve.cpp:10
void solve_with_mlmg(int lev, amrex::Vector< amrex::MultiFab > &rhs, amrex::Vector< amrex::MultiFab > &p, amrex::Vector< amrex::Array< amrex::MultiFab, AMREX_SPACEDIM >> &fluxes)
Definition: ERF_SolveWithMLMG.cpp:40
amrex::Gpu::DeviceVector< amrex::BCRec > domain_bcs_type_d
Definition: ERF.H:930
amrex::Vector< std::unique_ptr< ForestDrag > > m_forest_drag
Definition: ERF.H:1251
const int datwidth
Definition: ERF.H:965
void derive_diag_profiles_stag(amrex::Real time, amrex::Gpu::HostVector< amrex::Real > &h_avg_u, amrex::Gpu::HostVector< amrex::Real > &h_avg_v, amrex::Gpu::HostVector< amrex::Real > &h_avg_w, amrex::Gpu::HostVector< amrex::Real > &h_avg_rho, amrex::Gpu::HostVector< amrex::Real > &h_avg_th, amrex::Gpu::HostVector< amrex::Real > &h_avg_ksgs, amrex::Gpu::HostVector< amrex::Real > &h_avg_Kmv, amrex::Gpu::HostVector< amrex::Real > &h_avg_Khv, amrex::Gpu::HostVector< amrex::Real > &h_avg_qv, amrex::Gpu::HostVector< amrex::Real > &h_avg_qc, amrex::Gpu::HostVector< amrex::Real > &h_avg_qr, amrex::Gpu::HostVector< amrex::Real > &h_avg_wqv, amrex::Gpu::HostVector< amrex::Real > &h_avg_wqc, amrex::Gpu::HostVector< amrex::Real > &h_avg_wqr, amrex::Gpu::HostVector< amrex::Real > &h_avg_qi, amrex::Gpu::HostVector< amrex::Real > &h_avg_qs, amrex::Gpu::HostVector< amrex::Real > &h_avg_qg, amrex::Gpu::HostVector< amrex::Real > &h_avg_uu, amrex::Gpu::HostVector< amrex::Real > &h_avg_uv, amrex::Gpu::HostVector< amrex::Real > &h_avg_uw, amrex::Gpu::HostVector< amrex::Real > &h_avg_vv, amrex::Gpu::HostVector< amrex::Real > &h_avg_vw, amrex::Gpu::HostVector< amrex::Real > &h_avg_ww, amrex::Gpu::HostVector< amrex::Real > &h_avg_uth, amrex::Gpu::HostVector< amrex::Real > &h_avg_vth, amrex::Gpu::HostVector< amrex::Real > &h_avg_wth, amrex::Gpu::HostVector< amrex::Real > &h_avg_thth, amrex::Gpu::HostVector< amrex::Real > &h_avg_ku, amrex::Gpu::HostVector< amrex::Real > &h_avg_kv, amrex::Gpu::HostVector< amrex::Real > &h_avg_kw, amrex::Gpu::HostVector< amrex::Real > &h_avg_p, amrex::Gpu::HostVector< amrex::Real > &h_avg_pu, amrex::Gpu::HostVector< amrex::Real > &h_avg_pv, amrex::Gpu::HostVector< amrex::Real > &h_avg_pw, amrex::Gpu::HostVector< amrex::Real > &h_avg_wthv)
Definition: ERF_Write1DProfiles_stag.cpp:296
amrex::Real sampler_per
Definition: ERF.H:1509
amrex::Vector< amrex::BoxArray > ba1d
Definition: ERF.H:1169
InputSpongeData input_sponge_data
Definition: ERF.H:730
amrex::Vector< amrex::Vector< amrex::BoxArray > > subdomains
Definition: ERF.H:1258
std::string restart_chkfile
Definition: ERF.H:984
bool metgrid_use_below_sfc
Definition: ERF.H:1162
void init_zphys(int lev, amrex::Real time)
Definition: ERF_MakeNewArrays.cpp:518
amrex::Vector< std::string > sampleptlogname
Definition: ERF.H:1520
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > xvel_bc_data
Definition: ERF.H:733
void InitData_pre()
Definition: ERF.cpp:781
amrex::IntVect & SamplePoint(int i)
Definition: ERF.H:1391
bool use_datetime
Definition: ERF.H:980
amrex::Vector< std::array< amrex::Real, 2 > > hurricane_track_xy
Definition: ERF.H:153
amrex::Vector< amrex::Vector< amrex::Real > > h_u_geos
Definition: ERF.H:1207
void InitializeLevelFromData(int lev, const amrex::MultiFab &initial_data)
std::string subvol_file
Definition: ERF.H:1012
int rad_datalog_int
Definition: ERF.H:851
void sum_derived_quantities(amrex::Real time)
Definition: ERF_WriteScalarProfiles.cpp:187
void sum_integrated_quantities(amrex::Real time)
Definition: ERF_WriteScalarProfiles.cpp:15
amrex::Vector< std::string > datalogname
Definition: ERF.H:1515
void initHSE()
Initialize HSE.
Definition: ERF_Init1D.cpp:142
static amrex::Real column_loc_x
Definition: ERF.H:1181
amrex::Vector< std::unique_ptr< amrex::MultiFab > > ax
Definition: ERF.H:885
amrex::Vector< std::unique_ptr< amrex::MultiFab > > z_phys_nd
Definition: ERF.H:881
void MakeDiagnosticAverage(amrex::Vector< amrex::Real > &h_havg, amrex::MultiFab &S, int n)
Definition: ERF.cpp:2325
amrex::Real metgrid_proximity
Definition: ERF.H:1165
void setRecordDerDataInfo(int i, const std::string &filename)
Definition: ERF.H:1447
amrex::Vector< amrex::Real > h_havg_temperature
Definition: ERF.H:1235
void advance_dycore(int level, amrex::Vector< amrex::MultiFab > &state_old, amrex::Vector< amrex::MultiFab > &state_new, amrex::MultiFab &xvel_old, amrex::MultiFab &yvel_old, amrex::MultiFab &zvel_old, amrex::MultiFab &xvel_new, amrex::MultiFab &yvel_new, amrex::MultiFab &zvel_new, amrex::MultiFab &source, amrex::MultiFab &xmom_src, amrex::MultiFab &ymom_src, amrex::MultiFab &zmom_src, amrex::MultiFab &buoyancy, amrex::Geometry fine_geom, amrex::Real dt, amrex::Real time)
Definition: ERF_AdvanceDycore.cpp:38
amrex::Vector< std::unique_ptr< std::fstream > > sampleptlog
Definition: ERF.H:1519
void FillWeatherDataMultiFab(const std::string &filename, const amrex::Geometry &geom_weather, const amrex::BoxArray &nba, const amrex::DistributionMapping &dm, amrex::Vector< amrex::MultiFab > &weather_forecast_data)
Definition: ERF_WeatherDataInterpolation.cpp:451
void poisson_wall_dist(int lev)
Definition: ERF_PoissonWallDist.cpp:20
amrex::Vector< std::unique_ptr< amrex::MultiFab > > detJ_cc_src
Definition: ERF.H:891
amrex::Gpu::DeviceVector< amrex::Real > d_havg_pressure
Definition: ERF.H:1242
std::string plot2d_file_1
Definition: ERF.H:1010
amrex::Vector< std::string > der_datalogname
Definition: ERF.H:1516
amrex::Vector< std::unique_ptr< amrex::MultiFab > > SmnSmn_lev
Definition: ERF.H:864
const amrex::Vector< std::string > derived_names
Definition: ERF.H:1048
void solve_with_gmres(int lev, const amrex::Box &subdomain, amrex::MultiFab &rhs, amrex::MultiFab &p, amrex::Array< amrex::MultiFab, AMREX_SPACEDIM > &fluxes, amrex::MultiFab &ax_sub, amrex::MultiFab &ay_sub, amrex::MultiFab &znd_sub)
Definition: ERF_SolveWithGMRES.cpp:12
amrex::Vector< std::unique_ptr< amrex::MultiFab > > ay_src
Definition: ERF.H:893
void sample_points(int lev, amrex::Real time, amrex::IntVect cell, amrex::MultiFab &mf)
Definition: ERF_WriteScalarProfiles.cpp:527
AMREX_FORCE_INLINE std::ostream & DataLog(int i)
Definition: ERF.H:1335
void writeJobInfo(const std::string &dir) const
Definition: ERF_WriteJobInfo.cpp:9
void solve_with_EB_mlmg(int lev, amrex::Vector< amrex::MultiFab > &rhs, amrex::Vector< amrex::MultiFab > &p, amrex::Vector< amrex::Array< amrex::MultiFab, AMREX_SPACEDIM >> &fluxes)
Definition: ERF_SolveWithEBMLMG.cpp:19
void ComputeDt(int step=-1)
Definition: ERF_ComputeTimestep.cpp:11
amrex::Vector< int > nsubsteps
Definition: ERF.H:766
amrex::Vector< std::unique_ptr< amrex::iMultiFab > > yflux_imask
Definition: ERF.H:946
amrex::Vector< amrex::MultiFab > rW_new
Definition: ERF.H:809
amrex::Vector< std::unique_ptr< amrex::MultiFab > > lon_m
Definition: ERF.H:722
std::unique_ptr< WriteBndryPlanes > m_w2d
Definition: ERF.H:1248
AMREX_FORCE_INLINE std::ostream & SampleLineLog(int i)
Definition: ERF.H:1378
amrex::Vector< amrex::Vector< amrex::MultiFab * > > lsm_flux
Definition: ERF.H:833
amrex::Vector< std::string > plot3d_var_names_1
Definition: ERF.H:1037
std::string SamplePointLogName(int i) const noexcept
The filename of the ith sampleptlog file.
Definition: ERF.H:1532
void turbPert_update(const int lev, const amrex::Real dt)
Definition: ERF_InitTurbPert.cpp:12
void InitData_post()
Definition: ERF.cpp:855
void refinement_criteria_setup()
Definition: ERF_Tagging.cpp:227
void InterpWeatherDataOntoMesh(const amrex::Geometry &geom_weather, amrex::MultiFab &weather_forecast_interp)
Definition: ERF_WeatherDataInterpolation.cpp:249
static int nghost_eb_volume()
Definition: ERF.H:1554
bool metgrid_debug_dry
Definition: ERF.H:1157
static AMREX_FORCE_INLINE int ComputeGhostCells(const SolverChoice &sc)
Definition: ERF.H:1270
static int bndry_output_planes_interval
Definition: ERF.H:1187
void init_geo_wind_profile(const std::string input_file, amrex::Vector< amrex::Real > &u_geos, amrex::Gpu::DeviceVector< amrex::Real > &u_geos_d, amrex::Vector< amrex::Real > &v_geos, amrex::Gpu::DeviceVector< amrex::Real > &v_geos_d, const amrex::Geometry &lgeom, const amrex::Vector< amrex::Real > &zlev_stag)
Definition: ERF_InitGeowind.cpp:10
void ERF_shared()
Definition: ERF.cpp:119
amrex::Vector< std::string > plot2d_var_names_2
Definition: ERF.H:1040
static int nghost_eb_full()
Definition: ERF.H:1557
void RemakeLevel(int lev, amrex::Real time, const amrex::BoxArray &ba, const amrex::DistributionMapping &dm) override
Definition: ERF_MakeNewLevel.cpp:444
static void GotoNextLine(std::istream &is)
Definition: ERF_Checkpoint.cpp:16
void MakeNewLevelFromScratch(int lev, amrex::Real time, const amrex::BoxArray &ba, const amrex::DistributionMapping &dm) override
Definition: ERF_MakeNewLevel.cpp:25
void init_stuff(int lev, const amrex::BoxArray &ba, const amrex::DistributionMapping &dm, amrex::Vector< amrex::MultiFab > &lev_new, amrex::Vector< amrex::MultiFab > &lev_old, amrex::MultiFab &tmp_base_state, std::unique_ptr< amrex::MultiFab > &tmp_zphys_nd)
Definition: ERF_MakeNewArrays.cpp:24
amrex::Vector< amrex::MultiFab > zmom_crse_rhs
Definition: ERF.H:813
bool metgrid_debug_isothermal
Definition: ERF.H:1156
void initSponge()
Initialize sponge profiles.
Definition: ERF_InitSponge.cpp:35
amrex::Vector< std::unique_ptr< std::fstream > > der_datalog
Definition: ERF.H:1513
amrex::Real m_plot3d_per_2
Definition: ERF.H:1020
void advance_lsm(int lev, amrex::MultiFab &cons_in, amrex::MultiFab &xvel_in, amrex::MultiFab &yvel_in, const amrex::Real &dt_advance)
Definition: ERF_AdvanceLSM.cpp:5
amrex::Vector< std::unique_ptr< std::fstream > > tot_e_datalog
Definition: ERF.H:1514
amrex::Vector< std::unique_ptr< amrex::MultiFab > > sinPhi_m
Definition: ERF.H:724
void derive_diag_profiles(amrex::Real time, amrex::Gpu::HostVector< amrex::Real > &h_avg_u, amrex::Gpu::HostVector< amrex::Real > &h_avg_v, amrex::Gpu::HostVector< amrex::Real > &h_avg_w, amrex::Gpu::HostVector< amrex::Real > &h_avg_rho, amrex::Gpu::HostVector< amrex::Real > &h_avg_th, amrex::Gpu::HostVector< amrex::Real > &h_avg_ksgs, amrex::Gpu::HostVector< amrex::Real > &h_avg_Kmv, amrex::Gpu::HostVector< amrex::Real > &h_avg_Khv, amrex::Gpu::HostVector< amrex::Real > &h_avg_qv, amrex::Gpu::HostVector< amrex::Real > &h_avg_qc, amrex::Gpu::HostVector< amrex::Real > &h_avg_qr, amrex::Gpu::HostVector< amrex::Real > &h_avg_wqv, amrex::Gpu::HostVector< amrex::Real > &h_avg_wqc, amrex::Gpu::HostVector< amrex::Real > &h_avg_wqr, amrex::Gpu::HostVector< amrex::Real > &h_avg_qi, amrex::Gpu::HostVector< amrex::Real > &h_avg_qs, amrex::Gpu::HostVector< amrex::Real > &h_avg_qg, amrex::Gpu::HostVector< amrex::Real > &h_avg_uu, amrex::Gpu::HostVector< amrex::Real > &h_avg_uv, amrex::Gpu::HostVector< amrex::Real > &h_avg_uw, amrex::Gpu::HostVector< amrex::Real > &h_avg_vv, amrex::Gpu::HostVector< amrex::Real > &h_avg_vw, amrex::Gpu::HostVector< amrex::Real > &h_avg_ww, amrex::Gpu::HostVector< amrex::Real > &h_avg_uth, amrex::Gpu::HostVector< amrex::Real > &h_avg_vth, amrex::Gpu::HostVector< amrex::Real > &h_avg_wth, amrex::Gpu::HostVector< amrex::Real > &h_avg_thth, amrex::Gpu::HostVector< amrex::Real > &h_avg_ku, amrex::Gpu::HostVector< amrex::Real > &h_avg_kv, amrex::Gpu::HostVector< amrex::Real > &h_avg_kw, amrex::Gpu::HostVector< amrex::Real > &h_avg_p, amrex::Gpu::HostVector< amrex::Real > &h_avg_pu, amrex::Gpu::HostVector< amrex::Real > &h_avg_pv, amrex::Gpu::HostVector< amrex::Real > &h_avg_pw, amrex::Gpu::HostVector< amrex::Real > &h_avg_wthv)
Definition: ERF_Write1DProfiles.cpp:190
int real_width
Definition: ERF.H:1147
void MakeNewLevelFromCoarse(int lev, amrex::Real time, const amrex::BoxArray &ba, const amrex::DistributionMapping &dm) override
Definition: ERF_MakeNewLevel.cpp:280
amrex::Vector< std::unique_ptr< amrex::MultiFab > > ax_src
Definition: ERF.H:892
amrex::Vector< std::unique_ptr< amrex::MultiFab > > z_phys_cc_src
Definition: ERF.H:890
void input_sponge(int lev)
Definition: ERF_InitSponge.cpp:17
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > d_u_geos
Definition: ERF.H:1208
void Define_ERFFillPatchers(int lev)
Definition: ERF.cpp:2397
AMREX_FORCE_INLINE int NumDataLogs() noexcept
Definition: ERF.H:1349
static PlotFileType plotfile3d_type_2
Definition: ERF.H:1133
void setRecordEnergyDataInfo(int i, const std::string &filename)
Definition: ERF.H:1460
TurbulentPerturbation turbPert
Definition: ERF.H:1090
amrex::Vector< amrex::MultiFab > rW_old
Definition: ERF.H:808
void ImposeBCsOnPhi(int lev, amrex::MultiFab &phi, const amrex::Box &subdomain)
Definition: ERF_ImposeBCsOnPhi.cpp:12
void FillCoarsePatch(int lev, amrex::Real time)
Definition: ERF_FillCoarsePatch.cpp:21
void ClearLevel(int lev) override
Definition: ERF_MakeNewLevel.cpp:666
std::unique_ptr< SurfaceLayer > m_SurfaceLayer
Definition: ERF.H:1250
static PlotFileType plotfile2d_type_2
Definition: ERF.H:1135
static amrex::Vector< amrex::AMRErrorTag > ref_tags
Definition: ERF.H:1256
void make_physbcs(int lev)
Definition: ERF_MakeNewArrays.cpp:656
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > d_w_subsid
Definition: ERF.H:1205
amrex::Vector< ERFFillPatcher > FPr_w
Definition: ERF.H:859
int real_set_width
Definition: ERF.H:1148
amrex::Vector< std::unique_ptr< amrex::MultiFab > > SFS_hfx2_lev
Definition: ERF.H:872
amrex::Vector< amrex::Vector< amrex::Real > > zlevels_stag
Definition: ERF.H:878
AMREX_FORCE_INLINE int NumDerDataLogs() noexcept
Definition: ERF.H:1356
amrex::EBFArrayBoxFactory const & EBFactory(int lev) const noexcept
Definition: ERF.H:1546
amrex::Vector< amrex::Vector< amrex::MultiFab * > > lsm_data
Definition: ERF.H:832
amrex::Vector< amrex::Vector< amrex::Real > > stretched_dz_h
Definition: ERF.H:910
int m_plot2d_int_1
Definition: ERF.H:1016
void WriteCheckpointFile() const
Definition: ERF_Checkpoint.cpp:26
int sampler_interval
Definition: ERF.H:1508
amrex::Real volWgtSumMF(int lev, const amrex::MultiFab &mf, int comp, bool finemask)
Definition: ERF_WriteScalarProfiles.cpp:651
static int output_1d_column
Definition: ERF.H:1178
bool metgrid_debug_psfc
Definition: ERF.H:1158
amrex::Vector< std::unique_ptr< amrex::MultiFab > > az_src
Definition: ERF.H:894
static int fixed_mri_dt_ratio
Definition: ERF.H:997
amrex::Vector< amrex::Real > dt
Definition: ERF.H:771
static amrex::Real init_shrink
Definition: ERF.H:989
void WriteVTKPolyline(const std::string &filename, amrex::Vector< std::array< amrex::Real, 2 >> &points_xy)
Definition: ERF_Write1DProfiles.cpp:593
void Write3DPlotFile(int which, PlotFileType plotfile_type, amrex::Vector< std::string > plot_var_names)
Definition: ERF_Plotfile.cpp:263
void advance_radiation(int lev, amrex::MultiFab &cons_in, const amrex::Real &dt_advance)
Definition: ERF_AdvanceRadiation.cpp:5
static amrex::Real dt_max_initial
Definition: ERF.H:991
static bool is_it_time_for_action(int nstep, amrex::Real time, amrex::Real dt, int action_interval, amrex::Real action_per)
Definition: ERF_WriteScalarProfiles.cpp:747
static std::string nc_low_file
Definition: ERF.H:1152
amrex::Vector< amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > > d_rayleigh_ptrs
Definition: ERF.H:1229
void InitData()
Definition: ERF.cpp:772
void advance_microphysics(int lev, amrex::MultiFab &cons_in, const amrex::Real &dt_advance, const int &iteration, const amrex::Real &time)
Definition: ERF_AdvanceMicrophysics.cpp:5
amrex::Vector< std::unique_ptr< amrex::MultiFab > > Lwave
Definition: ERF.H:918
amrex::Vector< std::unique_ptr< amrex::MultiFab > > cosPhi_m
Definition: ERF.H:724
void ParameterSanityChecks()
Definition: ERF.cpp:2152
void derive_stress_profiles(amrex::Gpu::HostVector< amrex::Real > &h_avg_tau11, amrex::Gpu::HostVector< amrex::Real > &h_avg_tau12, amrex::Gpu::HostVector< amrex::Real > &h_avg_tau13, amrex::Gpu::HostVector< amrex::Real > &h_avg_tau22, amrex::Gpu::HostVector< amrex::Real > &h_avg_tau23, amrex::Gpu::HostVector< amrex::Real > &h_avg_tau33, amrex::Gpu::HostVector< amrex::Real > &h_avg_hfx3, amrex::Gpu::HostVector< amrex::Real > &h_avg_q1fx3, amrex::Gpu::HostVector< amrex::Real > &h_avg_q2fx3, amrex::Gpu::HostVector< amrex::Real > &h_avg_diss)
Definition: ERF_Write1DProfiles.cpp:475
int cf_width
Definition: ERF.H:854
int last_plot2d_file_step_1
Definition: ERF.H:957
amrex::Vector< amrex::Vector< amrex::Vector< amrex::Real > > > h_rayleigh_ptrs
Definition: ERF.H:1223
amrex::Vector< std::unique_ptr< amrex::iMultiFab > > zflux_imask
Definition: ERF.H:947
bool m_expand_plotvars_to_unif_rr
Definition: ERF.H:1013
amrex::Real m_plot3d_per_1
Definition: ERF.H:1019
int plot_file_on_restart
Definition: ERF.H:962
void Construct_ERFFillPatchers(int lev)
Definition: ERF.cpp:2371
void post_update(amrex::MultiFab &state_mf, amrex::Real time, const amrex::Geometry &geom)
amrex::Vector< int > num_boxes_at_level
Definition: ERF.H:761
static amrex::Real sub_cfl
Definition: ERF.H:988
static void print_error(MPI_Comm, const std::string &msg)
Definition: ERF_ConsoleIO.cpp:43
static int ng_dens_hse
Definition: ERF.H:1194
std::unique_ptr< ReadBndryPlanes > m_r2d
Definition: ERF.H:1249
amrex::Vector< std::unique_ptr< amrex::MultiFab > > ay
Definition: ERF.H:886
void project_momenta(int lev, amrex::Real dt, amrex::Vector< amrex::MultiFab > &vars)
Definition: ERF_PoissonSolve.cpp:41
amrex::Vector< amrex::Vector< amrex::MultiFab * > > qmoist
Definition: ERF.H:816
amrex::Vector< std::unique_ptr< amrex::MultiFab > > SFS_q1fx1_lev
Definition: ERF.H:874
static amrex::Real getCPUTime()
Definition: ERF.H:1420
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > zvel_bc_data
Definition: ERF.H:735
void derive_stress_profiles_stag(amrex::Gpu::HostVector< amrex::Real > &h_avg_tau11, amrex::Gpu::HostVector< amrex::Real > &h_avg_tau12, amrex::Gpu::HostVector< amrex::Real > &h_avg_tau13, amrex::Gpu::HostVector< amrex::Real > &h_avg_tau22, amrex::Gpu::HostVector< amrex::Real > &h_avg_tau23, amrex::Gpu::HostVector< amrex::Real > &h_avg_tau33, amrex::Gpu::HostVector< amrex::Real > &h_avg_hfx3, amrex::Gpu::HostVector< amrex::Real > &h_avg_q1fx3, amrex::Gpu::HostVector< amrex::Real > &h_avg_q2fx3, amrex::Gpu::HostVector< amrex::Real > &h_avg_diss)
Definition: ERF_Write1DProfiles_stag.cpp:605
LandSurface lsm
Definition: ERF.H:831
static amrex::Real change_max
Definition: ERF.H:990
void AverageDownTo(int crse_lev, int scomp, int ncomp)
Definition: ERF_AverageDown.cpp:36
void setRecordSampleLineInfo(int i, int lev, amrex::IntVect &cell, const std::string &filename)
Definition: ERF.H:1490
void setSpongeRefFromSounding(bool restarting)
Set sponge mean profiles from input sounding.
Definition: ERF_InitSponge.cpp:65
amrex::Vector< std::unique_ptr< amrex::MultiFab > > detJ_cc_new
Definition: ERF.H:897
void WeatherDataInterpolation(const amrex::Real time)
Definition: ERF_WeatherDataInterpolation.cpp:511
amrex::Array< std::string, 2 *AMREX_SPACEDIM > domain_bc_type
Definition: ERF.H:933
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > yvel_bc_data
Definition: ERF.H:734
void init_uniform(int lev)
Definition: ERF_InitUniform.cpp:17
static void writeBuildInfo(std::ostream &os)
Definition: ERF_WriteJobInfo.cpp:137
amrex::GpuArray< ERF_BC, AMREX_SPACEDIM *2 > phys_bc_type
Definition: ERF.H:942
std::string check_file
Definition: ERF.H:1033
amrex::Vector< amrex::Vector< amrex::MultiFab > > initial_state
Definition: ERF.H:154
amrex::Vector< amrex::IntVect > samplepoint
Definition: ERF.H:1521
amrex::Vector< amrex::Real > h_havg_qv
Definition: ERF.H:1237
static void print_usage(MPI_Comm, std::ostream &)
Definition: ERF_ConsoleIO.cpp:26
amrex::Vector< amrex::MultiFab > rV_old
Definition: ERF.H:806
amrex::Vector< std::unique_ptr< amrex::MultiFab > > lat_m
Definition: ERF.H:722
std::unique_ptr< amrex::MultiFab > mf_C1H
Definition: ERF.H:1171
int last_plot3d_file_step_1
Definition: ERF.H:955
amrex::Vector< std::unique_ptr< amrex::MultiFab > > Hwave
Definition: ERF.H:917
amrex::Vector< int > istep
Definition: ERF.H:765
const int datprecision
Definition: ERF.H:966
amrex::Vector< std::unique_ptr< amrex::iMultiFab > > xflux_imask
Definition: ERF.H:945
void derive_upwp(amrex::Vector< amrex::Real > &h_havg)
static PlotFileType plotfile2d_type_1
Definition: ERF.H:1134
int metgrid_order
Definition: ERF.H:1166
bool finished_wave
Definition: ERF.H:921
void ReadCheckpointFile()
Definition: ERF_Checkpoint.cpp:439
bool writeNow(const amrex::Real cur_time, const amrex::Real dt, const int nstep, const int plot_int, const amrex::Real plot_per)
Definition: ERF.cpp:2448
amrex::Vector< amrex::Vector< amrex::MultiFab > > vars_old
Definition: ERF.H:778
amrex::MultiFab & build_fine_mask(int lev)
Definition: ERF_WriteScalarProfiles.cpp:710
amrex::Vector< std::unique_ptr< ERFPhysBCFunct_base > > physbcs_base
Definition: ERF.H:796
void init_thin_body(int lev, const amrex::BoxArray &ba, const amrex::DistributionMapping &dm)
Definition: ERF_MakeNewLevel.cpp:705
static PlotFileType plotfile3d_type_1
Definition: ERF.H:1132
AMREX_FORCE_INLINE std::ostream & SamplePointLog(int i)
Definition: ERF.H:1364
amrex::Array< amrex::Array< amrex::Real, AMREX_SPACEDIM *2 >, AMREX_SPACEDIM+NBCVAR_max > m_bc_neumann_vals
Definition: ERF.H:939
amrex::Vector< amrex::Vector< std::unique_ptr< amrex::MultiFab > > > tsk_lev
Definition: ERF.H:868
void setRecordSamplePointInfo(int i, int lev, amrex::IntVect &cell, const std::string &filename)
Definition: ERF.H:1473
static int column_interval
Definition: ERF.H:1179
static void print_summary(std::ostream &)
int m_plot3d_int_2
Definition: ERF.H:1015
amrex::Vector< std::unique_ptr< amrex::MultiFab > > qheating_rates
Definition: ERF.H:836
void turbPert_amplitude(const int lev)
Definition: ERF_InitTurbPert.cpp:32
void sample_lines(int lev, amrex::Real time, amrex::IntVect cell, amrex::MultiFab &mf)
Definition: ERF_WriteScalarProfiles.cpp:563
void initializeMicrophysics(const int &)
Definition: ERF.cpp:1594
bool plot_lsm
Definition: ERF.H:1026
const amrex::Vector< std::string > cons_names
Definition: ERF.H:1041
void ReadCheckpointFileSurfaceLayer()
Definition: ERF_Checkpoint.cpp:1003
void timeStep(int lev, amrex::Real time, int iteration)
Definition: ERF_TimeStep.cpp:17
Definition: ERF_LandSurface.H:15
@ pres
Definition: ERF_Kessler.H:25
@ qv
Definition: ERF_Kessler.H:28
@ T
Definition: ERF_IndexDefines.H:110
Definition: ERF_ConsoleIO.cpp:12
real(c_double), parameter, private pi
Definition: ERF_module_mp_morr_two_moment.F90:100
AdvType moistscal_horiz_adv_type
Definition: ERF_AdvStruct.H:399
AdvType dycore_vert_adv_type
Definition: ERF_AdvStruct.H:396
AdvType moistscal_vert_adv_type
Definition: ERF_AdvStruct.H:400
AdvType dryscal_horiz_adv_type
Definition: ERF_AdvStruct.H:397
AdvType dycore_horiz_adv_type
Definition: ERF_AdvStruct.H:395
AdvType dryscal_vert_adv_type
Definition: ERF_AdvStruct.H:398
Definition: ERF_DataStruct.H:123
bool use_num_diff
Definition: ERF_DataStruct.H:866
AdvChoice advChoice
Definition: ERF_DataStruct.H:785
static TerrainType terrain_type
Definition: ERF_DataStruct.H:768
Definition: ERF_TurbPertStruct.H:22