forked from williamBarnhart/gCent
-
Notifications
You must be signed in to change notification settings - Fork 0
/
gCent_in.m
34 lines (30 loc) · 1.76 KB
/
gCent_in.m
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
% gCent_in.m
% Input script for driving gCent
%Event parameters
eventID = 'Nevada_20200515';
eventID2 = 'nn00725272'
eventLoc = [-117.875 38.159 0]; %Lon, lat, depth
eventSDR = [73 78 -24]; %Event Strike, dip, rake
eventMag = 7; %Event magnitude
eventDate = [2020 05 15]; % Event date, year, month, day
% Datafiles, split between InSAR, GPS, and optical sensors.
% Include all data files that you wish to include in the inversion, even if they've already been resampled.
% Give full paths to processed directories
%
% InSAR data files: give path to filt_topophase.unw.geo and los.rdr.geo
% files
%
% Optical data: give full path to EW and NS displacements files
% (displacements in meters)
%
% GPS data: give full path to a .mat file
insarDataFiles = {'/Users/wbarnhart/Work_local/EQmonitoring/Nevada_20200515/p144/200510-200516/merged','/Users/wbarnhart/Work_local/EQmonitoring/Nevada_20200515/p64/200511-200517/merged'};
processor = 'MINTPY'; % Can be 'ISCE', 'ISCE_MAI', 'MINTPY', 'ROIPAC', or 'GMT'. GMT inputs should be convertedto .xyz files prior to running using grd2xyz. Compatible only with MINTPY outputs generated from ISCE
iscestack = 'topsStack'; % Specify if ISCE stack processor used is topsStack or alosStack. Leave empty when not using MINTPY processor
opticalDataFilesEW = {};
opticalDataFilesNS = {};
gpsTimeSeriesDir = {};
WORKDIR = ['/Users/wbarnhart/Work_local/EQmonitoring/' eventID];
%Elevation of water in dem.crop file. Used to crop out water areas. To skip
%water masking, set waterElev = [];
waterElev = [];