Change order of compute_Background_PE optional arguments

The dimensional rho option should not be after the mapping options
since the mapping options will need to be specified for an
unmapped dimensional density case. This now makes the optional
argument ordering more logical.

src/Science.cpp

@@ -296,7 +296,7 @@ bool compare_pairs( pair<double, double> a, pair<double, double> b ) {
 // Compute Background Potential Energy (BPE)
 void compute_Background_PE(double & BPE_tot, TArrayn::DTArray & rho,
         int Nx, int Ny, int Nz, double Lx, double Ly, double Lz, double g,
-        double rho_0, int iter, bool mapped, Array<double,1> hill, bool dimensional_rho) {
+        double rho_0, int iter, bool dimensional_rho, bool mapped, Array<double,1> hill) {
     // Tensor variables for indexing
     blitz::firstIndex ii;
     blitz::secondIndex jj;

src/Science.hpp

@@ -42,7 +42,7 @@ void dissipation(TArrayn::DTArray & diss, TArrayn::DTArray & u, TArrayn::DTArray
 // Background Potential Energy (BPE)
 void compute_Background_PE(double & BPE_tot, TArrayn::DTArray & rho, int Nx, int Ny, int Nz,
         double Lx, double Ly, double Lz, double g, double rho_0, int iter,
-        bool mapped = false, Array<double,1> hill = Array<double,1>(), bool dimensional_rho = false);
+        bool dimensional_rho = false, bool mapped = false, Array<double,1> hill = Array<double,1>());
 // Internal energy converted to BPE
 void compute_BPE_from_internal(double & phi_i, TArrayn::DTArray & rho,
         double kappa_rho, double rho_0, double g, int Nz,