fix dxyz in filename of add_mann_turbulence

wetb/hawc2/htc_extensions.py

@@ -55,7 +55,7 @@ class HTCDefaults(object):
         if filenames is None:
 
             import numpy as np
-            dxyz = tuple(np.array(box_dimension) / (np.array(no_grid_points) - 1))
+            dxyz = tuple(np.array(box_dimension) / np.array(no_grid_points))
             from wetb.wind.turbulence import mann_turbulence
             filenames = ["./turb/" + mann_turbulence.name_format %
                          ((L, ae23, Gamma, high_frq_compensation) + no_grid_points + dxyz + (seed, uvw))

wetb/hawc2/tests/test_htc_file.py

@@ -159,9 +159,9 @@ class TestHtcFile(unittest.TestCase):
         htcfile.add_mann_turbulence(30.1, 1.1, 3.3, 102, False)
         s = """begin mann;
     create_turb_parameters\t30.1 1.1 3.3 102 0;\tL, alfaeps, gamma, seed, highfrq compensation
-    filename_u\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.23x3.23_s0102u.turb;
-    filename_v\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.23x3.23_s0102v.turb;
-    filename_w\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.23x3.23_s0102w.turb;
+    filename_u\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.12x3.12_s0102u.turb;
+    filename_v\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.12x3.12_s0102v.turb;
+    filename_w\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.12x3.12_s0102w.turb;
     box_dim_u\t16384 0.3662;
     box_dim_v\t32 3.125;
     box_dim_w\t32 3.125;"""