def __init__(self):
self.flat = False
self.tilt = False
self.isomod = modparam.isomod()
self.nlay = 0
self.ngrid = 0
return
import modparam
import numpy as np
import matplotlib.pyplot as plt
#
# para = modparam.para1d()
isomod = modparam.isomod()
# isomod.readmodtxt('./old_code/TEST/Q22A.mod1')
# isomod.init_arr(3)
# isomod.thickness= np.array([0.2, 29.8, 170.])
#
#
# isomod.numbp = np.array([2, 4, 5])
# isomod.mtype = np.array([4, 2, 2])
# isomod.vpvs = np.array([2., 1.75, 1.75])
# isomod.update()
import h5py
dset = h5py.File('./CU_SDT1.0.mod.h5')
data = dset['118.0_-22.0'].value
isomod.parameterize_input(zarr=data[:, 0],
vsarr=data[:, 1],
mohodepth=30.,
seddepth=0.2,
maxdepth=200.)
isomod.update()
# vs = np.append(isomod.vs[:2, 0], isomod.vs[:30, 1])
# vs = np.append(vs, isomod.vs[:60, 2])
# h = np.append(isomod.hArr[:2, 0], isomod.hArr[:30, 1])