RFDFWIM/inp_GPRmodel.m at main · mklayek/RFDFWIM · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
function [ model ] = inp_GPRmodel( )
%inp_GPRmodel - Mainly for sig_model- Must required
% defining input parameters for a Ground-Penetrating Radar
% (GPR) modeling and inversion process.
% --------------------------------------------
%   Copyright @ Mrinal Kanti Layek
%   Senior Research Fellow (Geophysics),
%   Department of Geology and Geophysics,
%   Indian Institute of Technology Kharagpur,
%   Kharagpur - 721302, India
%   https://www.researchgate.net/profile/Mrinal_Layek
%   14.01.2018 (DD.MM.YYYY)
% --------------------------------------------
% On 14 th jan, 2018
path(path, eval('pwd'));
path(path, [eval('pwd') '/forward']);

path(path, [eval('pwd') '/input']);
path(path, [eval('pwd') '/inv']);

%==========================================================================
%==================== Forward Model parameters ============================
%==========================================================================

%
        model.nx =180;
        model.ny =180;
    % spatial sampling interval [m]
    model.dh = 0.05;
    model.n = model.nx * model.ny;

    % define x- and y-axis [m]
    model.x = model.dh.*(1:model.nx);
    model.y = model.dh.*(1:model.ny);
% ------------------------------------:-------------------------------------
    % define modeling frequency [Hz]
%     model.f = 50e6;
    model.s = 1.0;
%     model.omega = 2.0*pi*(model.f+1i*model.s);
    model.mu0 = 4.0*pi*1e-7;

    % define thickness of PML boundary frame [gridpoints]
    model.npml = 10;
    model.a0_pml = 9e8;

    % define free surface
    model.free = 0;

    % basename of input model

% %         model.file_true2 = 'input/CHPDX24.sig';
%     model.file_init = 'input/2_cross_TE_init.sig';
%
%     model.file_true2 = 'input/2_cross_TE_true.sig';
%
%     % file names for source and receiver file
%     model.file_rec = 'input/receiver_2_cross_TE';
%     model.file_src = 'input/source_2_cross_TE';


    % define Caxis clip - Set colormap limits
    model.caxis1 = 10e-3;
    model.caxis2 = 1e-6;
    % define fonts
    model.FSize = 18;
    model.Fweight = 'normal';

    % % define figure size;
    % model.screenx = 1000;
    % model.screeny = 1000;


end