MT.py 文件源码

def PlotAppRes3Layers_wrapper(fmin, fmax, nbdata, h1, h2, rhol1, rhol2, rhol3, mul1, mul2, mul3, epsl1, epsl2, epsl3, PlotEnvelope, F_Envelope):

    frangn=frange(np.log10(fmin), np.log10(fmax), nbdata)
    sig3= np.array([0., 0.001, 0.1, 0.001])
    thick3 = np.array([120000., 50., 50.])
    eps3=np.array([1., 1., 1., 1])
    mu3=np.array([1., 1., 1., 1])
    chg3=np.array([0., 0.1, 0., 0.2])
    chg3_0=np.array([0., 0.1, 0., 0.])
    taux3=np.array([0., 0.1, 0., 0.1])
    c3=np.array([1., 1., 1., 1.])

    sig3[1]=1./rhol1
    sig3[2]=1./rhol2
    sig3[3]=1./rhol3
    mu3[1]=mul1
    mu3[2]=mul2
    mu3[3]=mul3
    eps3[1]=epsl1
    eps3[2]=epsl2
    eps3[3]=epsl3
    thick3[1]=h1
    thick3[2]=h2

    PlotAppRes(frangn, thick3, sig3, chg3_0, taux3, c3, mu3, eps3, 3, F_Envelope, PlotEnvelope)