CB_013_G_K_inputs.m

contents

function VBR = CB_013_G_K_inputs()
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    % CB_013_G_K_inputs.m
    %
    %  Specify unrelaxed shear and bulk moduli.
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    % specify state variables as usual
    VBR = struct();
    VBR.in.SV.T_K = linspace(800, 1000, 4);
    sz_T = size(VBR.in.SV.T_K);
    VBR.in.SV.P_GPa = linspace(2, 3, 4);
    VBR.in.SV.rho = 3300 * ones(sz_T);
    VBR.in.SV.phi = 0.01 * ones(sz_T);
    VBR.in.SV.sig_MPa = 1 * ones(sz_T);
    VBR.in.SV.dg_um = 1e4 * ones(sz_T);
    VBR.in.SV.f = [0.01, 0.1];

    % specify methods as usual
    VBR.in.elastic.methods_list={'anharmonic';'anh_poro';};
    VBR.in.anelastic.methods_list={'eburgers_psp';'andrade_psp';'xfit_mxw'};
    VBR.in.viscous.methods_list={'HZK2011'};

    % also specify the unrelaxed moduli at elevated temperature, pressure.
    % you could, instead, load these from a file!
    VBR.in.elastic.Gu_TP = linspace(50, 60, 4) * 1e9; % shear modulus
    VBR.in.elastic.Ku_TP = VBR.in.elastic.Gu_TP * 1.5; % bulk modulus

    % call the VBRc
    VBR = VBR_spine(VBR);
end