diff --git a/wetb/hawc2/htc_extensions.py b/wetb/hawc2/htc_extensions.py index b585503ebf7e950f81fe0a92b9c843ec30595f6e..644ac8fbe04f1d0dbbeb9545c06b8c4eb3da3d42 100644 --- a/wetb/hawc2/htc_extensions.py +++ b/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)) diff --git a/wetb/hawc2/tests/test_htc_file.py b/wetb/hawc2/tests/test_htc_file.py index d6b06e2af51804b578a973c8918d5c098c82b1c1..f3fd1afbcc0e1c9a7dd9f4dee8d1ab1987528739 100644 --- a/wetb/hawc2/tests/test_htc_file.py +++ b/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;"""