def test_FUV(): x = loadbunch(fname, masked=False) # Switch epoch from Matlab to Python x.t -= 366 x.t0 -= 366 for i, flag in enumerate(x.flags): F, U, V = FUV(x.t, x.t0, x.lind - 1, x.lat, flag) print(f"i: {i} ngflags: {flag}") # We use broadcasting instead of replication, so # we need to test only against the first row of # the octave output in such cases. if F.shape[0] == 1: sub = (i, slice(0, 1)) else: sub = (i, ) assert_array_almost_equal(F, x.Fo[sub]) assert_array_almost_equal(U, x.Uo[sub]) if V.shape[0] == 1: sub = (i, slice(0, 1)) else: sub = (i, ) assert_array_almost_equal(V, x.Vo[sub])
def test_FUV(): x = loadbunch(fname, masked=False) # Switch epoch from Matlab to Python x.t -= 366 x.t0 -= 366 for i, flag in enumerate(x.flags): F, U, V = FUV(x.t, x.t0, x.lind-1, x.lat, flag) print('i: {} ngflags: {}'.format(i, flag)) # We use broadcasting instead of replication, so # we need to test only against the first row of # the octave output in such cases. if F.shape[0] == 1: sub = (i, slice(0, 1)) else: sub = (i,) assert_array_almost_equal(F, x.Fo[sub]) assert_array_almost_equal(U, x.Uo[sub]) if V.shape[0] == 1: sub = (i, slice(0, 1)) else: sub = (i,) assert_array_almost_equal(V, x.Vo[sub])
def test_FUV(): x = loadbunch(fname, masked=False) for i, flag in enumerate(x.flags): F, U, V = FUV(x.t, x.t0, x.lind - 1, x.lat, flag) print('i:', i, "ngflgs:", flag) # We use broadcasting instead of replication, so # we need to test only against the first row of # the octave output in such cases. if F.shape[0] == 1: sub = (i, slice(0, 1)) else: sub = (i, ) assert_array_almost_equal(F, x.Fo[sub]) assert_array_almost_equal(U, x.Uo[sub]) if V.shape[0] == 1: sub = (i, slice(0, 1)) else: sub = (i, ) assert_array_almost_equal(V, x.Vo[sub])
from __future__ import (absolute_import, division, print_function) import os from utide.utilities import loadbunch __version__ = '0.1.0' _base_dir = os.path.join(os.path.dirname(__file__), 'data') _ut_constants_fname = os.path.join(_base_dir, 'ut_constants.npz') # At least for now, use NaNs rather than masked arrays. ut_constants = loadbunch(_ut_constants_fname, masked=False) # A list of strings is much easier to work with. constit_names = [n.strip() for n in ut_constants.const.name] # Make a dictionary for index lookups. constit_index_dict = dict([(name, i) for (i, name) in enumerate(constit_names)]) from ._solve import solve # noqa from ._reconstruct import reconstruct # noqa __all__ = ['solve', 'reconstruct']
print('Size of boundary: ', len(bnd)) seg_types = [0] * len(bnd) adcirc.order_boundary(bnd,seg_types) #adcirc.update(os.path.join(data_files_dir,'vdatum','vdatum_fl_sab_adcirc54.nc')) def parse_string(name): lista = [e.decode().strip() for e in name.tolist()] return ''.join(lista) names = [] const = adcirc.Dataset.variables['tidenames'][:] for name in const: names.append(parse_string(name.data)) con_info = loadbunch(_ut_constants_fname)['const'] from utide import _ut_constants_fname from utide.utilities import loadbunch con_info = loadbunch(_ut_constants_fname)['const'] k = 0 ind_nc, ind_ttide = [], [] const_name = [e.strip() for e in con_info['name'].tolist()] consts = ['STEADY', 'M2', 'S2', 'N2', 'K1', 'O1', 'P1', 'M4', 'M6'] for name in consts: try: if name == 'STEADY': indx = const_name.index('Z0')
print('Size of boundary: ', len(bnd)) seg_types = [0] * len(bnd) adcirc.order_boundary(bnd, seg_types) def parse_string(name): lista = [e.decode().strip() for e in name.tolist()] return ''.join(lista) names = [] const = adcirc.Dataset.variables['tidenames'][:] for name in const: names.append(parse_string(name)) #changed from name.data con_info = loadbunch(_ut_constants_fname)['const'] k = 0 ind_nc, ind_ttide = [], [] const_name = [e.strip() for e in con_info['name'].tolist()] consts = ['STEADY', 'M2', 'S2', 'N2', 'K1', 'O1', 'P1', 'M4', 'M6'] for name in consts: try: if name == 'STEADY': indx = const_name.index('Z0') else: indx = const_name.index(name) k += 1 ind_ttide.append(indx)
from __future__ import (absolute_import, division, print_function) import os import numpy as np from utide.utilities import loadbunch datadir = os.path.join(os.path.dirname(__file__)) fname = 'ut_constants' data = loadbunch( os.path.join(datadir, '{}.mat'.format(fname)), masked=False) np.savez('{}.npz'.format(fname), **data) fname = 'FUV0' data = loadbunch( os.path.join(datadir, '{}.mat'.format(fname)), masked=False) np.savez('{}.npz'.format(fname), **data)
def test_name_A_g(coef, matfile): coef_mat = loadbunch(os.path.join(datapath, matfile)).coef coef_mat = convert_unicode_arrays(coef_mat) assert_array_equal(coef.name, coef_mat.name) assert_array_almost_equal(coef.A, coef_mat.A) assert_array_almost_equal(coef.g, coef_mat.g)
import os import numpy as np from utide.utilities import loadbunch datadir = os.path.join(os.path.dirname(__file__)) fname = "ut_constants" data = loadbunch(os.path.join(datadir, f"{fname}.mat"), masked=False) np.savez(f"{fname}.npz", **data) fname = "FUV0" data = loadbunch(os.path.join(datadir, f"{fname}.mat"), masked=False) np.savez(f"{fname}.npz", **data)
def test_name_A_g(coef, matfile): coef_mat = loadbunch(os.path.join(datapath, matfile)).coef coef_mat = convert_unicode_arrays(coef_mat) assert_array_equal(coef.name, coef_mat.name) assert_array_almost_equal(coef.A, coef_mat.A) assert_array_almost_equal(coef.g, coef_mat.g)