def check_data():
for mkey in variables.keys():
for mmkey in list(variables[mkey].keys()):
for month in range(1, 13):
NCDFS = np.array(
NetCDFSWAN.load(
os.path.join(
tmpFolder,
f'2000/{str(month)}/results/{mkey}.mat'))[mmkey])
start = int(NCDFS[0][0] //
10) # actual value is similar to the index
end = start + NCDFS.shape[0]
v = variables[mkey][mmkey][start:end] # 0 - 745
np.testing.assert_array_equal(NCDFS, v)
print(f"{mkey}.mat ok")
for i, station in enumerate(stations):
n = stations[station]["nsnodes"]
sts = 0
for month in range(1, 13):
NCDFS = NetCDFSWAN.load(
os.path.join(
tmpFolder,
f'2000/{str(month)}/results/{station}.spc'))["spectra"]
ts = NCDFS.shape[0]
spg = spcgroup["spectra"][i, :n]
for node in range(n):
spg_n = spg[node][sts:sts + ts]
NCDFS_n = NCDFS[:, node]
np.testing.assert_array_equal(spg_n, NCDFS_n)
#print(f"station {station} node {node}:", NCDFS_n.shape, spg_n.shape, " start index:", sts)
sts += ts - 1
print(f"{station}.spc ok")
def test_load(): print(NetCDFSWAN.load(os.path.join(tmpFolder,"2000/1/results/HS.mat"))['Hsig'].shape) print(NetCDFSWAN.load(os.path.join(tmpFolder,"2000/1/results/WIND.mat"))) print(NetCDFSWAN.load(os.path.join(tmpFolder,"Mesh/dummy.bot")).shape) print(NetCDFSWAN.load(os.path.join(tmpFolder,"Mesh/dummy.ele")).shape)