Skip to content

Commit

Permalink
Add initHSE call for init_type==ideal
Browse files Browse the repository at this point in the history
  • Loading branch information
ewquon committed Sep 9, 2023
1 parent 049f84f commit 44e7625
Showing 1 changed file with 3 additions and 2 deletions.
5 changes: 3 additions & 2 deletions Source/ERF.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -793,6 +793,7 @@ ERF::init_only (int lev, Real time)
// The base state is initialized from WRF wrfinput data, output by
// ideal.exe or real.exe
init_from_wrfinput(lev);
if (init_type == "ideal") initHSE();

} else if (init_type == "metgrid") {
// The base state is initialized from data output by WPS metgrid;
Expand All @@ -802,13 +803,13 @@ ERF::init_only (int lev, Real time)
} else if (init_type == "uniform") {
// Initialize a uniform background field and base state based on the
// problem-specified reference density and temperature
initHSE(lev);
init_uniform(lev);
initHSE(lev);
} else {
// No background flow initialization specified, initialize the
// background field to be equal to the base state, calculated from the
// problem-specific erf_init_dens_hse
initHSE(lev);
initHSE(lev); // need to call this first
init_from_hse(lev);
}

Expand Down

0 comments on commit 44e7625

Please sign in to comment.