Kessler microphysics updated (erf-model#1326)

* Changes to EOS.H with moisture * Change all functions in EOS.H to work with moisture * Adding qp for buoyancy type 1 * Updating docs for buoyancy term * Updating docs for buoyancy term * Update docs for buoyancy * Avoiding extra divide * Update buoyancy docs * Some corrections to buoyancy docs * Some corrections to Docs in buoyancy * Correcting docs issue * Correcting docs * adding Kessler microphysics docs * adding Kessler microphysics docs * adding Kessler microphysics docs * adding Kessler microphysics docs * Removing unused vars * Updating prob for squall line * Updating prob for squall line * Minor change to makefile * Minor change * Finalizing Kessler * Commenting out IC file writing * Updating prob.H with host device macros * Correcting argument list * Correcting argument list * Correcting CUDA compilation errors * Correcting CUDA compilation errors * Correcting CUDA compilation errors --------- Co-authored-by: Mahesh Natarajan <[email protected]> Co-authored-by: Ann Almgren <[email protected]>
atmyers · Nov 30, 2023 · 49b69fa · 49b69fa
1 parent 0e5e440
commit 49b69fa
Show file tree

Hide file tree

Showing 9 changed files with 117 additions and 195 deletions.
diff --git a/Exec/SquallLine_2D/inputs_moisture b/Exec/SquallLine_2D/inputs_moisture
@@ -45,7 +45,7 @@ amr.max_level       = 0       # maximum level number allowed
 # CHECKPOINT FILES
 amr.check_file      = chk        # root name of checkpoint file
 amr.check_int       = 1000       # number of timesteps between checkpoints
-#amr.restart         = chk12000
+#amr.restart         = chk09000
 
 # PLOTFILES
 erf.plot_file_1         = plt        # root name of plotfile
@@ -70,7 +70,7 @@ erf.molec_diff_type = "ConstantAlpha"
 erf.rho0_trans = 1.0 # [kg/m^3], used to convert input diffusivities
 erf.dynamicViscosity = 200.0 # [kg/(m-s)] ==> nu = 75.0 m^2/s
 erf.alpha_T = 00.0 # [m^2/s]
-erf.alpha_C = 50.0
+erf.alpha_C = 25.0
 
 erf.moisture_model = "Kessler"
 erf.use_moist_background = true
@@ -79,6 +79,13 @@ erf.moistscal_horiz_adv_string = "Centered_2nd"
 erf.moistscal_vert_adv_string = "Centered_2nd"
 
 # PROBLEM PARAMETERS (optional)
-prob.T_0 = 300.0
-prob.U_0 = 0
-prob.T_pert = 3
+prob.z_tr = 12000.0
+prob.height = 1200.0
+prob.theta_0 = 300.0
+prob.theta_tr = 343.0
+prob.T_tr = 213.0
+prob.x_c = 0.0
+prob.z_c = 1500.0
+prob.x_r = 4000.0
+prob.z_r = 1500.0
+prob.theta_c = 3.0
diff --git a/Exec/SquallLine_2D/inputs_moisture_outflow b/Exec/SquallLine_2D/inputs_moisture_outflow
diff --git a/Exec/SquallLine_2D/prob.H b/Exec/SquallLine_2D/prob.H
@@ -13,16 +13,20 @@
 #include "TileNoZ.H"
 
 struct ProbParm : ProbParmDefaults {
-  amrex::Real T_0 = 300.0; // surface temperature == mean potential temperature
-  amrex::Real U_0 = 10.0;
-  amrex::Real V_0 = 0.0;
-  amrex::Real x_c = 0.0; // center of thermal perturbation
-  amrex::Real z_c = 3200.0;
-  amrex::Real x_r = 1000.0;
-  amrex::Real z_r = 1000.0;
-  amrex::Real T_pert = 5.0; // perturbation temperature
-  // overridden physical constants
+  amrex::Real z_tr = 12000.0;
+  amrex::Real height = 1200.0;
+  amrex::Real theta_0 = 300.0;
+  amrex::Real theta_tr = 343.0;
+  amrex::Real T_tr = 213.0;
+
+  amrex::Real x_c = 0.0;
+  amrex::Real z_c = 1.5e3;
+  amrex::Real x_r = 4.0e3;
+  amrex::Real z_r = 1.5e3;
+  amrex::Real theta_c = 3.0;
+
   amrex::Real C_p = 1004.0;
+  amrex::Real Theta_0 = 300.0;
 }; // namespace ProbParm
 
 
@@ -63,7 +67,6 @@ public:
 
     void erf_init_dens_hse_moist (amrex::MultiFab& rho_hse,
                                   std::unique_ptr<amrex::MultiFab>& z_phys_nd,
-                                  std::unique_ptr<amrex::MultiFab>& z_phys_cc,
                                   amrex::Geometry const& geom) override;
 
     void init_custom_terrain (
@@ -79,6 +82,46 @@ public:
         amrex::Vector<amrex::Real>& thetabar,
         amrex::Geometry      const& geom) override;
 
+    amrex::Real compute_theta (amrex::Real z);
+
+    amrex::Real compute_p_k (amrex::Real& p_k,
+                             const amrex::Real p_k_minus_1,
+                             amrex::Real& theta_k,
+                             amrex::Real& rho_k,
+                             amrex::Real& q_v_k,
+                             amrex::Real& T_dp,
+                             amrex::Real& T_b,
+                             const amrex::Real dz,
+                             const amrex::Real z,
+                             const amrex::Real rho_k_minus_1);
+
+    amrex::Real compute_F (const amrex::Real& p_k,
+                           const amrex::Real& p_k_minus_1,
+                           amrex::Real &theta_k,
+                           amrex::Real& rho_k,
+                           amrex::Real& q_v_k,
+                           amrex::Real& T_dp,
+                           amrex::Real& T_b,
+                           const amrex::Real& dz,
+                           const amrex::Real& z,
+                           const amrex::Real& rho_k_minus_1);
+
+    void compute_rho (const amrex::Real& z,
+                      const amrex::Real& pressure,
+                      amrex::Real &theta,
+                      amrex::Real& rho,
+                      amrex::Real& q_v,
+                      amrex::Real& T_dp,
+                      amrex::Real& T_b);
+
+    void init_isentropic_hse_no_terrain(amrex::Real *theta,
+                                        amrex::Real* r,
+                                        amrex::Real* p,
+                                        amrex::Real *q_v,
+                                        const amrex::Real& dz,
+                                        const amrex::Real&  prob_lo_z,
+                                        const int& khi);
+
 protected:
     std::string name () override { return "Supercell"; }