percolating some changes through -

EXP-code · Oct 26, 2024 · 283f9d8 · 283f9d8
1 parent 085f662
commit 283f9d8
Show file tree

Hide file tree

Showing 8 changed files with 9 additions and 9 deletions.
diff --git a/expui/BiorthBasis.cc b/expui/BiorthBasis.cc
@@ -1031,7 +1031,7 @@ namespace BasisClasses
     nmax        = 18;
     mmax        = 6;
     mlim        = std::numeric_limits<int>::max();
-    lmaxfid     = 128;
+    lmaxfid     = 72;
     nmaxfid     = 64;
     ncylnx      = 256;
     ncylny      = 128;

diff --git a/src/expand.H b/src/expand.H
@@ -26,7 +26,7 @@
    KEYWORD        EXPLANATION                            DEFAULT
    ----------     ----------------------------           --------
    lmax =         maximum harmonic order                 4
-   nmax =         maxinum radial order                   10
+   nmax =         maximum radial order                   18
    nlog =         interval between log output            10
    nlist =        interval between p-s output            50
    nsteps =       total number of steps to run           500

diff --git a/utils/ICs/cylcache.cc b/utils/ICs/cylcache.cc
@@ -327,7 +327,7 @@ main(int ac, char **av)
     ("ctype", "DiskHalo radial coordinate scaling type (one of: Linear, Log, Rat)",
      cxxopts::value<string>(ctype)->default_value("Log"))
     ("LMAXFID", "Maximum angular order for spherical basis in adaptive construction of the cylindrical basis",
-     cxxopts::value<int>(LMAXFID)->default_value("32"))
+     cxxopts::value<int>(LMAXFID)->default_value("72"))
     ("NMAXFID", "Maximum radial order for the spherical basis in adapative construction of the cylindrical basis",
      cxxopts::value<int>(NMAXFID)->default_value("32"))
     ("MMAX", "Maximum azimuthal order for the cylindrical basis",
@@ -341,7 +341,7 @@ main(int ac, char **av)
     ("NCYLODD", "Number of vertically odd basis functions per harmonic order",
      cxxopts::value<int>(NCYLODD)->default_value("6"))
     ("NMAX", "Total number of basis functions per harmonic order",
-     cxxopts::value<int>(NMAX)->default_value("20"))
+     cxxopts::value<int>(NMAX)->default_value("18"))
     ("VFLAG", "Diagnostic flag for EmpCylSL",
      cxxopts::value<int>(VFLAG)->default_value("31"))
     ("expcond", "Use analytic target density rather than particle distribution",

diff --git a/utils/ICs/initial.cc b/utils/ICs/initial.cc
@@ -398,7 +398,7 @@ main(int ac, char **av)
     ("LMAX", "Harmonic order for halo expansion",
      cxxopts::value<int>(LMAX)->default_value("18"))
     ("LMAXFID", "Harmonic order for EOF spherical model",
-     cxxopts::value<int>(LMAXFID)->default_value("128"))
+     cxxopts::value<int>(LMAXFID)->default_value("72"))
     ("NMAXFID", "Radial order for EOF spherical model",
      cxxopts::value<int>(NMAXFID)->default_value("64"))
     ("MMAX", "Aximuthal order for Cylindrical expansion",

diff --git a/utils/SL/diskpot.cc b/utils/SL/diskpot.cc
@@ -57,7 +57,7 @@ main(int argc, char** argv)
   int cmap     = 0;
   double scale = 1.0;
 
-  int Lmax=16, Nmax=10;
+  int Lmax=16, Nmax=18;
   int numr=10000;
   double rmin=0.0001, rmax=1.0;
   double delr=0.01, xmax=1.0, zmax=0.1;

diff --git a/utils/SL/oftest.cc b/utils/SL/oftest.cc
@@ -73,7 +73,7 @@ int main(int argc, char** argv)
     ("N,mc",       "number of particles for Monte Carlo realization",
      cxxopts::value<int>(N)->default_value("10000"))
     ("nmax",       "maximum number of radial harmonics in the expansion",
-     cxxopts::value<int>(nmax)->default_value("10"))
+     cxxopts::value<int>(nmax)->default_value("18"))
     ("r,rmin",     "minimum radius for the SL grid",
      cxxopts::value<double>(rmin)->default_value("0.0001"))
     ("R,rmax",     "maximum radius for the SL grid",

diff --git a/utils/SL/qtest.cc b/utils/SL/qtest.cc
@@ -35,7 +35,7 @@ int main(int argc, char** argv)
    ("Lmax", "maximum number of angular harmonics in the expansion",
      cxxopts::value<int>(Lmax)->default_value("2"))
    ("nmax", "maximum number of radial harmonics in the expansion",
-     cxxopts::value<int>(nmax)->default_value("10"))
+     cxxopts::value<int>(nmax)->default_value("18"))
    ("numr", "radial knots for the SL grid",
      cxxopts::value<int>(numr)->default_value("1000"))
    ("rmin", "minimum radius for the SL grid",

diff --git a/utils/SL/slabchk.cc b/utils/SL/slabchk.cc
@@ -29,7 +29,7 @@ main(int argc, char** argv)
     ("K,kmax", "maximum order of in-plane harmonics",
      cxxopts::value<int>(kmax)->default_value("4"))
     ("N,nmax", "maximum number of vertical harmonics",
-     cxxopts::value<int>(nmax)->default_value("10"))
+     cxxopts::value<int>(nmax)->default_value("18"))
     ("n,numz", "size of vertical grid",
      cxxopts::value<int>(numz)->default_value("1000"))
     ("Z,zmax", "maximum extent of vertical grid",