Added a 'cylheight' parameter whose sole function is to multiply itself by 2 and set 'hcyl'

Author: Martin D. Weinberg

expui/BiorthBasis.cc

@@ -874,6 +874,7 @@ namespace BasisClasses
     "rcylmax",
     "acyl",
     "hcyl",
+    "cylheight",
     "snr",
     "evcut",
     "nmaxfid",
@@ -1130,6 +1131,7 @@ namespace BasisClasses
 
       if (conf["acyl"      ])       acyl  = conf["acyl"      ].as<double>();
       if (conf["hcyl"      ])       hcyl  = conf["hcyl"      ].as<double>();
+      if (conf["cylheight" ])       hcyl  = conf["cylheight" ].as<double>()*2;
       if (conf["lmaxfid"   ])    lmaxfid  = conf["lmaxfid"   ].as<int>();
       if (conf["nmaxfid"   ])    nmaxfid  = conf["nmaxfid"   ].as<int>();
       if (conf["nmax"      ])       nmax  = conf["nmax"      ].as<int>();

src/Cylinder.cc

@@ -29,6 +29,7 @@ Cylinder::valid_keys = {
   "rcylmax",
   "acyl",
   "hcyl",
+  "cylheight",
   "hexp",
   "snr",
   "evcut",
@@ -414,6 +415,7 @@ void Cylinder::initialize()
 
     if (conf["acyl"      ])       acyl  = conf["acyl"      ].as<double>();
     if (conf["hcyl"      ])       hcyl  = conf["hcyl"      ].as<double>();
+    if (conf["cylheight" ])       hcyl  = conf["cylheight" ].as<double>()*2;
     if (conf["hexp"      ])       hexp  = conf["hexp"      ].as<double>();
     if (conf["snr"       ])        snr  = conf["snr"       ].as<double>();
     if (conf["evcut"     ])        rem  = conf["evcut"     ].as<double>();