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])