Skip to content
GitLab
Projects
Groups
Snippets
Help
Loading...
Help
What's new
10
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Open sidebar
GPIAG-Software
IFOS3D
Commits
5bb8c6d7
Commit
5bb8c6d7
authored
Apr 18, 2016
by
Tilman Steinweg
Browse files
Options
Browse Files
Download
Plain Diff
Merge branch 'develop'
parents
f42c903a
29d95e46
Changes
28
Hide whitespace changes
Inline
Side-by-side
Showing
28 changed files
with
2536 additions
and
308 deletions
+2536
-308
.gitignore
.gitignore
+2
-0
par/Makefile
par/Makefile
+38
-0
par/in_and_out/ifos3d.inp
par/in_and_out/ifos3d.inp
+11
-6
par/in_and_out/ifos3d_inv_all_parameters.json
par/in_and_out/ifos3d_inv_all_parameters.json
+132
-0
par/in_and_out/ifos3d_inv_all_parameters_commented.json
par/in_and_out/ifos3d_inv_all_parameters_commented.json
+240
-0
par/in_and_out/ifos3d_toy.inp
par/in_and_out/ifos3d_toy.inp
+12
-7
par/in_and_out/ifos3d_toy.json
par/in_and_out/ifos3d_toy.json
+100
-0
par/in_and_out/ifos3d_toy_FW.json
par/in_and_out/ifos3d_toy_FW.json
+73
-0
par/receiver/receiver.dat
par/receiver/receiver.dat
+20
-0
par/run_toy_example.sh
par/run_toy_example.sh
+9
-0
src/Makefile
src/Makefile
+15
-9
src/checkfd_ssg.c
src/checkfd_ssg.c
+1
-1
src/exchange_par.c
src/exchange_par.c
+5
-5
src/fd.h
src/fd.h
+30
-1
src/globvar.h
src/globvar.h
+4
-3
src/hh_toy_start.c
src/hh_toy_start.c
+218
-0
src/hh_toy_true.c
src/hh_toy_true.c
+0
-0
src/ifos3d.c
src/ifos3d.c
+33
-21
src/json_parser.c
src/json_parser.c
+463
-0
src/part_model.c
src/part_model.c
+20
-1
src/read_par.c
src/read_par.c
+11
-21
src/read_par_json.c
src/read_par_json.c
+797
-0
src/receiver.c
src/receiver.c
+168
-134
src/saveseis.c
src/saveseis.c
+14
-13
src/snapmerge.c
src/snapmerge.c
+80
-53
src/sources.c
src/sources.c
+4
-4
src/wavelet.c
src/wavelet.c
+6
-6
src/writepar.c
src/writepar.c
+30
-23
No files found.
.gitignore
0 → 100644
View file @
5bb8c6d7
*.kdev4
*.swp
par/Makefile
0 → 100644
View file @
5bb8c6d7
# This is par/Makefile
#--------------------------------------------------
# the command "make" or "make install" installs the
# cseife library and the main code of ifos3d
#--------------------------------------------------
default
:
../ ifos3d
.PHONY
:
libcseife
libcseife
:
$(MAKE)
-C
../libcseife
install
.PHONY
:
ifos3d
ifos3d
:
libcseife
ifdef
MODEL
$(MAKE)
-C
../src/ ifos3d
MODEL_SCR
=
$(MODEL)
else
$(MAKE)
-C
../src/ ifos3d
endif
.PHONY
:
install
install
:
clean-bin ifos3d
.PHONY
:
clean
clean
:
rm
-rf
../bin/
*
$(MAKE)
-C
../libcseife clean
$(MAKE)
-C
../src/ clean
.PHONY
:
clean-bin
clean-bin
:
rm
-rf
../bin/
*
.PHONY
:
reinstall
reinstall
:
clean ifos3d
par/in_and_out/ifos3d.inp
View file @
5bb8c6d7
...
@@ -34,16 +34,16 @@ timestep_(in_seconds)_(DT) = 5.0e-5
...
@@ -34,16 +34,16 @@ timestep_(in_seconds)_(DT) = 5.0e-5
#
#
#--------------------Source---------------------------------------
#--------------------Source---------------------------------------
# Shape_of_source-signal:
# Shape_of_source-signal:
(ricker=1;fumue=2;from_SIGNAL_FILE=3;SIN**3=4;deltapulse=5)_(
QUELLART
) = 4
(ricker=1;fumue=2;from_SIGNAL_FILE=3;SIN**3=4;deltapulse=5)_(
SOURCE_SHAPE
) = 4
point_source_(explosive=1;force_in_x=2;in_y=3;in_z=4;custom=5)_(
QUELL
TYP) = 4
point_source_(explosive=1;force_in_x=2;in_y=3;in_z=4;custom=5)_(
SOURCE_
TYP
E
) = 4
# If
QUELL
TYP <5 the following two lines are ignored
# If
SOURCE_
TYP
E
<5 the following two lines are ignored
force_angle_between_x_y_(in_degree)_(ALPHA) = 45.0
force_angle_between_x_y_(in_degree)_(ALPHA) = 45.0
force_angle_between_x_z_(in_degree)_(BETA) = 45.0
force_angle_between_x_z_(in_degree)_(BETA) = 45.0
# Plane wave (PW) excitation,if PLANE_WAVE_DEPTH>0, SRCREC is treated as 0
# Plane wave (PW) excitation,if PLANE_WAVE_DEPTH>0, SRCREC is treated as 0
depth_of_PW_excitation_(no<=0)_(in_meter)_(PLANE_WAVE_DEPTH) = 0.0
depth_of_PW_excitation_(no<=0)_(in_meter)_(PLANE_WAVE_DEPTH) = 0.0
dip_of_PW_from_vertical_(in_degrees)_(PHI) =0.0
dip_of_PW_from_vertical_(in_degrees)_(PHI) =0.0
duration_of_source-signal_PW_(in_seconds)_(TS) = 0.0033
duration_of_source-signal_PW_(in_seconds)_(TS) = 0.0033
# External signal input instead of
QUELLART
# External signal input instead of
SOURCE_SHAPE
SIGNAL_FILE = ?
SIGNAL_FILE = ?
read_source_positions_from_SOURCE_FILE_(yes=1)_(SRCREC) = 1
read_source_positions_from_SOURCE_FILE_(yes=1)_(SRCREC) = 1
SOURCE_FILE = ./sources/sources_toy.dat
SOURCE_FILE = ./sources/sources_toy.dat
...
@@ -103,11 +103,16 @@ output_of_seismograms_(SEISMO) = 1
...
@@ -103,11 +103,16 @@ output_of_seismograms_(SEISMO) = 1
# SEISMO=3: curl and div
# SEISMO=3: curl and div
# SEISMO=4: everything
# SEISMO=4: everything
# Warning: "curl" is not really curl in 3D
# Warning: "curl" is not really curl in 3D
read_receiver_positions_from_file_(yes=1)_(READREC) = 0
read_receiver_positions_(READREC) = 2
# receiver line READREC=0
# from_file READREC=1
# rec array READREC=2
REC_FILE = ./receiver/receiver.dat
REC_FILE = ./receiver/receiver.dat
reference_point_for_receiver_coordinate_system_(REFREC) = 0.0 , 0.0 , 0.0
reference_point_for_receiver_coordinate_system_(REFREC) = 0.0 , 0.0 , 0.0
# if READREC=1 the following
three line
s are ignored
# if READREC=1 the following
receiver option
s are ignored
# Note that z denotes the vertical direction !
# Note that z denotes the vertical direction !
#
# --------------------Receiver line -------------------------------
position_of_first_receiver_(in_m)_(XREC1,YREC1,ZREC1) = 90.0 , 90.0 , 90.0
position_of_first_receiver_(in_m)_(XREC1,YREC1,ZREC1) = 90.0 , 90.0 , 90.0
position_of_last_receiver_(in_m)_(XREC2,YREC2,ZREC2) = 90.0 , 90.0 , 90.0
position_of_last_receiver_(in_m)_(XREC2,YREC2,ZREC2) = 90.0 , 90.0 , 90.0
distance_between_two_adjacent_receivers_(in_gridpoints)_(NGEOPH) = 1
distance_between_two_adjacent_receivers_(in_gridpoints)_(NGEOPH) = 1
...
...
par/in_and_out/ifos3d_inv_all_parameters.json
0 → 100644
View file @
5bb8c6d7
{
"MODELING PARAMETERS"
:
"comment"
,
"Domain Decomposition"
:
"comment"
,
"NPROCX"
:
"2"
,
"NPROCY"
:
"2"
,
"NPROCZ"
:
"2"
,
"3-D Grid"
:
"comment"
,
"NX"
:
"160"
,
"NY"
:
"160"
,
"NZ"
:
"184"
,
"DX"
:
"0.8"
,
"DY"
:
"0.8"
,
"DZ"
:
"0.8"
,
"FD order"
:
"comment"
,
"FDORDER"
:
"4"
,
"FDCOEFF"
:
"2"
,
"Time Stepping"
:
"comment"
,
"TIME"
:
"0.06"
,
"DT"
:
"5.0e-05"
,
"Source"
:
"comment"
,
"SOURCE_SHAPE"
:
"4"
,
"SOURCE_TYPE"
:
"4"
,
"ALPHA"
:
"45.0"
,
"BETA"
:
"45.0"
,
"PLANE_WAVE_DEPTH"
:
"0.0"
,
"PHI"
:
"0.0"
,
"TS"
:
"0.0033"
,
"SIGNAL_FILE"
:
"./STF/stf.su"
,
"SRCREC"
:
"1"
,
"SOURCE_FILE"
:
"./sources/sources_toy.dat"
,
"RUN_MULTIPLE_SHOTS"
:
"1"
,
"Model"
:
"comment"
,
"READMOD"
:
"0"
,
"MFILE"
:
"model/toy"
,
"Q-approximation"
:
"comment"
,
"L"
:
"0"
,
"FL1"
:
"1000.0"
,
"TAU"
:
"0.000001"
,
"Free Surface"
:
"comment"
,
"FREE_SURF"
:
"0"
,
"Absorbing Boundary"
:
"comment"
,
"ABS_TYPE"
:
"1"
,
"FW"
:
"10"
,
"DAMPING"
:
"8.0"
,
"VPPML"
:
"6200.0"
,
"FPML"
:
"200.00"
,
"BOUNDARY"
:
"0"
,
"Snapshots"
:
"comment"
,
"SNAP"
:
"0"
,
"TSNAP1"
:
"0.01"
,
"TSNAP2"
:
"0.24"
,
"TSNAPINC"
:
"0.0075"
,
"IDX, IDY, IDZ"
:
"1,1,1"
,
"SNAP_FORMAT"
:
"4"
,
"SNAP_FILE"
:
"./snap/back"
,
"SNAP_PLANE"
:
"1"
,
"Receiver"
:
"comment"
,
"SEISMO"
:
"1"
,
"READREC"
:
"0"
,
"REC_FILE"
:
"./receiver/receiver.dat"
,
"REFRECX, REFRECY, REFRECZ"
:
"0.0 , 0.0, 0.0"
,
"XREC1, YREC1, ZREC1"
:
"90.0 , 90.0, 90.0"
,
"XREC2, YREC2, ZREC2"
:
"90.0 , 90.0, 90.0"
,
"NGEOPH"
:
"1"
,
"Receiver array"
:
"comment"
,
"REC_ARRAY"
:
"1"
,
"REC_ARRAY_DEPTH"
:
"24.0"
,
"REC_ARRAY_DIST"
:
"30.0"
,
"DRX"
:
"10"
,
"DRY"
:
"10"
,
"Seismograms"
:
"comment"
,
"NDT"
:
"1"
,
"NDTSHIFT"
:
"0"
,
"SEIS_FORMAT"
:
"1"
,
"SEIS_FILE"
:
"su/IFOS"
,
"Method"
:
"comment"
,
"METHOD"
:
"0"
,
"INVERSION PARAMETERS"
:
"comment"
,
"In- and Output Files"
:
"comment"
,
"GRAD_FILE"
:
"./grad/toy_grad"
,
"MOD_OUT_FILE"
:
"./model/toy"
,
"SEIS_OBS_FILE"
:
"./su_obs/obs_toy"
,
"EXTOBS"
:
"0"
,
"INV_FILE"
:
"./in_and_out/workflow_toy.dat"
,
"HESS_FILE"
:
"./hess/toy_hess"
,
"General"
:
"comment"
,
"ITMIN, ITMAX"
:
"1 , 80"
,
"FILT"
:
"1"
,
"NFMAX"
:
"5"
,
"TAST"
:
"100"
,
"VP0, VS0, RHO0"
:
"6200.0, 3600.0, 2800.0"
,
"WEIGHT_VP,WEIGHT_VS,WEIGHT_RHO"
:
"1.0, 1.0, 0.0"
,
"Steplength estimation"
:
"comment"
,
"NSHOTS_STEP"
:
"4"
,
"TESTSTEP"
:
"0.02"
,
"Gradient preconditioning"
:
"comment"
,
"DAMPTYPE"
:
"2"
,
"Hessian"
:
"comment"
,
"HESS"
:
"0"
,
"READ_HESS"
:
"0"
,
"REC_HESS"
:
"1"
,
"WATER_HESS_VP, WATER_HESS_VS, WATER_HESS_RHO"
:
"0.0192, 0.0192, 1.0e-14"
,
"L-BFGS"
:
"comment"
,
"LBFGS"
:
"0"
,
"NUMPAR"
:
"2"
,
"BFGSNUM"
:
"5"
}
par/in_and_out/ifos3d_inv_all_parameters_commented.json
0 → 100644
View file @
5bb8c6d7
{
"MODELING PARAMETERS"
:
"comment"
,
"Note that z denotes the vertical direction !"
:
"comment"
,
"Domain Decomposition"
:
"comment"
,
"NPROCX"
:
"2"
,
"NPROCY"
:
"2"
,
"NPROCZ"
:
"2"
,
"3-D Grid"
:
"comment"
,
"NX"
:
"160"
,
"NY"
:
"160"
,
"NZ"
:
"184"
,
"DX"
:
"0.8"
,
"DY"
:
"0.8"
,
"DZ"
:
"0.8"
,
"Caution: first gridpoint is not at {0.0,0.0,0.0} but at {dx,dy,dz}!"
:
"comment"
,
"FD order"
:
"comment"
,
"FDORDER"
:
"4"
,
"possible values are 2, 4, 6, 8, 10, 12"
:
"comment"
,
"fd_coefficients_(Taylor=1;Holberg=2)"
:
"comment"
,
"FDCOEFF"
:
"2"
,
"Time Stepping"
:
"comment"
,
"time_of_wave_propagation_(in_sec)"
:
"comment"
,
"TIME"
:
"0.06"
,
"timestep_(in_seconds)"
:
"comment"
,
"DT"
:
"5.0e-05"
,
"Source"
:
"comment"
,
"Shape_of_source-signal:"
:
"comment"
,
"SOURCE_SHAPE"
:
"4"
,
"SOURCE_SHAPE values: ricker=1;fumue=2;from_SIGNAL_FILE=3;SIN**3=4;deltapulse=5"
:
"comment"
,
"External signal input instead of SOURCE_SHAPE"
:
"comment"
,
"SIGNAL_FILE"
:
"./STF/stf.su"
,
"SOURCE_TYPE"
:
"4"
,
"SOURCE_TYPE values (point_source): explosive=1;force_in_x=2;force_in_y=3;force_in_z=4;custom=5"
:
"comment"
,
"If SOURCE_TYPE <5 the following two lines are ignored"
:
"comment"
,
"ALPHA"
:
"45.0"
,
"BETA"
:
"45.0"
,
"ALPHA=force_angle_between_x_y_(in_degree) (BETA between x_z)"
:
"comment"
,
"read_source_positions_from_SOURCE_FILE_(yes=1;Plane wave=2) "
:
"comment"
,
"SRCREC"
:
"1"
,
"SOURCE_FILE"
:
"./sources/sources_toy.dat"
,
"run_multiple_shots_defined_in_SOURCE_FILE_(yes=1)"
:
"comment"
,
"RUN_MULTIPLE_SHOTS"
:
"1"
,
"Following lines can be ignored if SRCREC!=2"
:
"comment"
,
"Plane wave (PW) excitation,if PLANE_WAVE_DEPTH>0, SRCREC is treated as 0"
:
"comment"
,
"depth_of_PW_excitation_(no<=0)_(in_meter)"
:
"comment"
,
"PLANE_WAVE_DEPTH"
:
"0.0"
,
"dip_of_PW_from_vertical_(in_degrees)"
:
"comment"
,
"PHI"
:
"0.0"
,
"duration_of_source-signal_PW_(in_seconds)"
:
"comment"
,
"TS"
:
"0.0033"
,
"Model"
:
"comment"
,
"read_model_from_MFILE(yes=1)"
:
"comment"
,
"READMOD"
:
"0"
,
"MFILE"
:
"model/toy"
,
"Q-approximation"
:
"comment"
,
"Number_of_relaxation_mechanisms"
:
"comment"
,
"L"
:
"0"
,
"L_Relaxation_frequencies"
:
"comment"
,
"FL1"
:
"1000.0"
,
"Tau_value_for_entire_model"
:
"comment"
,
"TAU"
:
"0.000001"
,
"Free Surface"
:
"comment"
,
"free_surface_(yes=1)"
:
"comment"
,
"FREE_SURF"
:
"0"
,
"Absorbing Boundary"
:
"comment"
,
"type_of_boundary_condition_(ABS_TYPE)_(PML=1/ABS=2)"
:
"comment"
,
"ABS_TYPE"
:
"1"
,
"width_of_absorbing_frame_(in_grid_points)_(No<=0)"
:
"comment"
,
"FW"
:
"10"
,
"percentage_of_amplitude_decay_at_outer_edge"
:
"comment"
,
"DAMPING"
:
"8.0"
,
"parameter for PML:"
:
"comment"
,
"P_wave_velocity_near_the_grid_boundary(in_m/s)"
:
"comment"
,
"VPPML"
:
"6200.0"
,
"Dominant_frequency_(in_Hz)"
:
"comment"
,
"FPML"
:
"200.00"
,
"apply_periodic_boundary_condition_at_edges"
:
"comment"
,
"BOUNDARY"
:
"0"
,
"Snapshots"
:
"comment"
,
"output_of_snapshots(yes>0)"
:
"comment"
,
"SNAP"
:
"0"
,
"output of particle velocities: SNAP=1"
:
"comment"
,
"output of pressure field: SNAP=2"
:
"comment"
,
"output of curl and divergence energy: SNAP=3"
:
"comment"
,
"output of both particle velocities and energy : SNAP=4"
:
"comment"
,
"first_snapshot_(in_sec)"
:
"comment"
,
"TSNAP1"
:
"0.01"
,
"last_snapshot_(in_sec)"
:
"comment"
,
"TSNAP2"
:
"0.24"
,
"increment_(in_sec)"
:
"comment"
,
"TSNAPINC"
:
"0.0075"
,
"increment_x,y,z-direction"
:
"comment"
,
"IDX, IDY, IDZ"
:
"1,1,1"
,
"data-format(ASCII(2);BINARY(3))"
:
"comment"
,
"SNAP_FORMAT"
:
"4"
,
"basic_filename"
:
"comment"
,
"SNAP_FILE"
:
"./snap/back"
,
"if SNAP !=3 the following line is ignored"
:
"comment"
,
"SNAP_PLANE"
:
"1"
,
"output of snapshots as energy wihout sign SNAP_PLANE=1"
:
"comment"
,
"energy with sign true for x-z-plane SNAP_PLANE=2"
:
"comment"
,
"energy with sign true for x-y-plane SNAP_PLANE=3"
:
"comment"
,
"energy with sign true for y-z-plane SNAP_PLANE=4"
:
"comment"
,
"Receiver"
:
"comment"
,
"output_of_seismograms"
:
"comment"
,
"SEISMO"
:
"1"
,
"SEISMO=0: no seismograms"
:
"comment"
,
"SEISMO=1: particle-velocities"
:
"comment"
,
"SEISMO=2: pressure (hydrophones)"
:
"comment"
,
"SEISMO=3: curl and div"
:
"comment"
,
"SEISMO=4: everything"
:
"comment"
,
"Warning: "
curl
" is not really curl in 3D"
:
"comment"
,
"read_receiver_positions_from_file_(yes=1)"
:
"comment"
,
"READREC"
:
"0"
,
"REC_FILE"
:
"./receiver/receiver.dat"
,
"reference_point_for_receiver_coordinate_system"
:
"comment"
,
"REFRECX, REFRECY, REFRECZ"
:
"0.0 , 0.0, 0.0"
,
"if READREC=1 the following three lines are ignored"
:
"comment"
,
"Note that z denotes the vertical direction !"
:
"comment"
,
"position_of_first_receiver_(in_m)"
:
"comment"
,
"XREC1, YREC1, ZREC1"
:
"90.0 , 90.0, 90.0"
,
"position_of_last_receiver_(in_m)"
:
"comment"
,
"XREC2, YREC2, ZREC2"
:
"90.0 , 90.0, 90.0"
,
"distance_between_two_adjacent_receivers"
:
"comment"
,
"NGEOPH"
:
"1"
,
"Receiver array"
:
"comment"
,
"parameters for horizontal plane of receivers"
:
"comment"
,
"number_of_planes_(no<=0)"
:
"comment"
,
"REC_ARRAY"
:
"1"
,
"depth_of_first_(upper)_plane_(in_m)"
:
"comment"
,
"REC_ARRAY_DEPTH"
:
"24.0"
,
"vertical_distance_between_planes_(in_m)"
:
"comment"
,
"REC_ARRAY_DIST"
:
"30.0"
,
"distance_between_receivers_in_x/y-direction"
:
"comment"
,
"DRX"
:
"10"
,
"DRY"
:
"10"
,
"Seismograms"
:
"comment"
,
"samplingrate_and_timelag_(in_timesteps!)"
:
"comment"
,
"NDT"
:
"1"
,
"NDTSHIFT"
:
"0"
,
"write every (ndt)th sample, leaving ndtshift samples at the beginning out"
:
"comment"
,
"default: ndt=1 and ndtshift=0. (ndt=0 is set to ndt=1, ndt<0 are set to -ndt.)"
:
"comment"
,
"data-format"
:
"comment"
,
"SEIS_FORMAT"
:
"1"
,
"0: SEG-Y (ASCII-text/native 4-byte-floats (IEEE on PC)/little endian on PC)"
:
"comment"
,
"1: SU (native 4-byte-floats (IEEE on PC)/little endian on PC)"
:
"comment"
,
"2: TEXTUAL (native ASCII)"
:
"comment"
,
"3: BINARY (IEEE-4-byte-floats on PC/little endian on PC)"
:
"comment"
,
"4: SEG-Y (ASCII-text/native 4-byte-floats (IEEE on PC)/little endian on PC)"
:
"comment"
,
"5: SEG-Y (ASCII-text/IBM-4-byte-floats on PC/big endian on PC) "
:
"comment"
,
"basic_filename"
:
"comment"
,
"SEIS_FILE"
:
"su/IFOS"
,
"Method"
:
"comment"
,
"METHOD"
:
"0"
,
"0: only forward simulation"
:
"comment"
,
"1: conjugate_gradient_FWI"
:
"comment"
,
"INVERSION PARAMETERS"
:
"comment"
,
"In- and Output Files"
:
"comment"
,
"gradient_filename"
:
"comment"
,
"GRAD_FILE"
:
"./grad/toy_grad"
,
"model_output_filename"
:
"comment"
,
"MOD_OUT_FILE"
:
"./model/toy"
,
"observed_data_fileneame"
:
"comment"
,
"SEIS_OBS_FILE"
:
"./su_obs/obs_toy"
,
"external_or_internal_observed_data"
:
"comment"
,
"EXTOBS"
:
"0"
,
"inversion_parameter_file"
:
"comment"
,
"INV_FILE"
:
"./in_and_out/workflow_toy.dat"
,
"hessian_file"
:
"comment"
,
"HESS_FILE"
:
"./hess/toy_hess"
,
"General"
:
"comment"
,
"minimum/maximum_iteration_number >0"
:
"comment"
,
"ITMIN, ITMAX"
:
"1 , 80"
,
"filtering_(yes=1/no=0)"
:
"comment"
,
"FILT"
:
"1"
,
"maximum_number_frequencies_per_iteration"
:
"comment"
,
"NFMAX"
:
"5"
,
"number_of_timestep_per_period_used_for_inversion"
:
"comment"
,
"TAST"
:
"100"
,
"average_model_parameter"
:
"comment"
,
"VP0, VS0, RHO0"
:
"6200.0, 3600.0, 2800.0"
,
"velocities in m/s, density in kg/m³"
:
"comment"
,
"parameter_class_weighting_factors_for_vp/vs/rho_"
:
"comment"
,
"WEIGHT_VP,WEIGHT_VS,WEIGHT_RHO"
:
"1.0, 1.0, 0.0"
,
"choose from 1.0 (full update) to 0.0 (no update)"
:
"comment"
,
"Steplength estimation"
:
"comment"
,
"number_of_shots_used_for_steplength_estimation"
:
"comment"
,
"NSHOTS_STEP"
:
"4"
,
"initial_test_steplength"
:
"comment"
,
"TESTSTEP"
:
"0.02"
,
"Gradient preconditioning"
:
"comment"
,
"Type_of_preconditioning"
:
"comment"
,
"DAMPTYPE"
:
"2"
,
"0: no damping; 1: Circular taper around receivers;"
:
"comment"
,
"2: Cicular taper around sources and receivers; 3: Tapering of source and receiver planes"
:
"comment"
,
"Hessian"
:
"comment"
,
"Apply_Hessian_(yes=1/no=0)"
:
"comment"
,
"HESS"
:
"0"
,
"Read_Hessian_from_file"
:
"comment"
,
"READ_HESS"
:
"0"
,
"Part_of_receivers_used_for_Hessian"
:
"comment"
,
"REC_HESS"
:
"1"
,
"(Each REC_HESS Receiver is used to calculate the Hessian, not yet implemented)"
:
"comment"
,
"Water_level_Hessian_for_vp/vs/rho"
:
"comment"
,
"WATER_HESS_VP, WATER_HESS_VS, WATER_HESS_RHO"
:
"0.0192, 0.0192, 1.0e-14"
,
"L-BFGS"
:
"comment"
,
"Apply_L_BFGS"
:
"comment"
,
"LBFGS"
:
"0"
,
"Number_of_inverted_parameters"
:
"comment"
,
"NUMPAR"
:
"2"
,
"Number_iterations_used_for_LBFGS"
:
"comment"
,
"BFGSNUM"
:
"5"
}
par/in_and_out/ifos3d_toy.inp
View file @
5bb8c6d7
...
@@ -34,16 +34,16 @@ timestep_(in_seconds)_(DT) = 5.0e-5
...
@@ -34,16 +34,16 @@ timestep_(in_seconds)_(DT) = 5.0e-5
#
#
#--------------------Source---------------------------------------
#--------------------Source---------------------------------------
# Shape_of_source-signal:
# Shape_of_source-signal:
(ricker=1;fumue=2;from_SIGNAL_FILE=3;SIN**3=4;deltapulse=5)_(
QUELLART
) = 4
(ricker=1;fumue=2;from_SIGNAL_FILE=3;SIN**3=4;deltapulse=5)_(
SOURCE_SHAPE
) = 4
point_source_(explosive=1;force_in_x=2;in_y=3;in_z=4;custom=5)_(
QUELL
TYP) = 4
point_source_(explosive=1;force_in_x=2;in_y=3;in_z=4;custom=5)_(
SOURCE_
TYP
E
) = 4
# If
QUELL
TYP <5 the following two lines are ignored
# If
SOURCE_
TYP
E
<5 the following two lines are ignored
force_angle_between_x_y_(in_degree)_(ALPHA) = 45.0
force_angle_between_x_y_(in_degree)_(ALPHA) = 45.0
force_angle_between_x_z_(in_degree)_(BETA) = 45.0
force_angle_between_x_z_(in_degree)_(BETA) = 45.0
# Plane wave (PW) excitation,if PLANE_WAVE_DEPTH>0, SRCREC is treated as 0
# Plane wave (PW) excitation,if PLANE_WAVE_DEPTH>0, SRCREC is treated as 0
depth_of_PW_excitation_(no<=0)_(in_meter)_(PLANE_WAVE_DEPTH) = 0.0
depth_of_PW_excitation_(no<=0)_(in_meter)_(PLANE_WAVE_DEPTH) = 0.0
dip_of_PW_from_vertical_(in_degrees)_(PHI) =0.0
dip_of_PW_from_vertical_(in_degrees)_(PHI) =0.0
duration_of_source-signal_PW_(in_seconds)_(TS) = 0.0033
duration_of_source-signal_PW_(in_seconds)_(TS) = 0.0033
# External signal input instead of
QUELLART
# External signal input instead of
SOURCE_SHAPE
SIGNAL_FILE = ?
SIGNAL_FILE = ?
read_source_positions_from_SOURCE_FILE_(yes=1)_(SRCREC) = 1
read_source_positions_from_SOURCE_FILE_(yes=1)_(SRCREC) = 1
SOURCE_FILE = ./sources/sources_toy.dat
SOURCE_FILE = ./sources/sources_toy.dat
...
@@ -103,11 +103,16 @@ output_of_seismograms_(SEISMO) = 1
...
@@ -103,11 +103,16 @@ output_of_seismograms_(SEISMO) = 1
# SEISMO=3: curl and div
# SEISMO=3: curl and div
# SEISMO=4: everything
# SEISMO=4: everything
# Warning: "curl" is not really curl in 3D
# Warning: "curl" is not really curl in 3D
read_receiver_positions_from_file_(yes=1)_(READREC) = 0
read_receiver_positions_(READREC) = 2
# receiver line READREC=0
# from_file READREC=1
# rec array READREC=2
REC_FILE = ./receiver/receiver.dat
REC_FILE = ./receiver/receiver.dat
reference_point_for_receiver_coordinate_system_(REFREC) = 0.0 , 0.0 , 0.0
reference_point_for_receiver_coordinate_system_(REFREC) = 0.0 , 0.0 , 0.0
# if READREC=1 the following
three line
s are ignored
# if READREC=1 the following
receiver option
s are ignored
# Note that z denotes the vertical direction !
# Note that z denotes the vertical direction !
#
# --------------------Receiver line -------------------------------
position_of_first_receiver_(in_m)_(XREC1,YREC1,ZREC1) = 90.0 , 90.0 , 90.0
position_of_first_receiver_(in_m)_(XREC1,YREC1,ZREC1) = 90.0 , 90.0 , 90.0
position_of_last_receiver_(in_m)_(XREC2,YREC2,ZREC2) = 90.0 , 90.0 , 90.0
position_of_last_receiver_(in_m)_(XREC2,YREC2,ZREC2) = 90.0 , 90.0 , 90.0
distance_between_two_adjacent_receivers_(in_gridpoints)_(NGEOPH) = 1
distance_between_two_adjacent_receivers_(in_gridpoints)_(NGEOPH) = 1
...
@@ -185,4 +190,4 @@ Water_level_Hessian_for_vp/vs/rho_(WATER_HESS) = 0.0192, 0.0192, 1.0e-14
...
@@ -185,4 +190,4 @@ Water_level_Hessian_for_vp/vs/rho_(WATER_HESS) = 0.0192, 0.0192, 1.0e-14
#------------------------L-BFGS ----------------------------------------
#------------------------L-BFGS ----------------------------------------
Apply_L_BFGS_(LBFGS) = 0
Apply_L_BFGS_(LBFGS) = 0
Number_of_inverted_parameters_(NUMPAR) = 2
Number_of_inverted_parameters_(NUMPAR) = 2
Number_iterations_used_for_LBFGS_(BFGSNUM) = 5
Number_iterations_used_for_LBFGS_(BFGSNUM) = 5
\ No newline at end of file
par/in_and_out/ifos3d_toy.json
0 → 100644
View file @
5bb8c6d7
{
"MODELING PARAMETERS"
:
"comment"
,
"Domain Decomposition"
:
"comment"
,
"NPROCX"
:
"2"
,
"NPROCY"
:
"2"
,
"NPROCZ"
:
"2"
,
"3-D Grid"
:
"comment"
,
"NX"
:
"160"
,
"NY"
:
"160"
,
"NZ"
:
"184"
,
"DX"
:
"0.8"
,
"DY"
:
"0.8"
,
"DZ"
:
"0.8"
,
"FD order"
:
"comment"
,
"FDORDER"
:
"4"
,
"FDCOEFF"
:
"2"
,
"Time Stepping"
:
"comment"
,
"TIME"
:
"0.06"
,
"DT"
:
"5.0e-05"
,
"Source"
:
"comment"
,
"SOURCE_SHAPE"
:
"4"
,
"SOURCE_TYPE"
:
"4"
,
"SRCREC"
:
"1"
,
"SOURCE_FILE"
:
"./sources/sources_toy.dat"
,
"RUN_MULTIPLE_SHOTS"
:
"1"
,
"Model"
:
"comment"
,
"READMOD"
:
"0"
,
"MFILE"
:
"model/toy"
,
"Q-approximation"
:
"comment"
,
"L"
:
"0"
,
"Free Surface"
:
"comment"
,
"FREE_SURF"
:
"0"
,
"Absorbing Boundary"
:
"comment"
,
"ABS_TYPE"
:
"1"
,
"FW"
:
"10"
,
"VPPML"
:
"6200.0"
,
"FPML"
:
"200.00"
,
"Snapshots"
:
"comment"
,
"SNAP"
:
"0"
,
"Receiver"
:
"comment"
,
"SEISMO"
:
"1"
,
"READREC"
:
"2"
,
"Receiver array"
:
"comment"
,
"REC_ARRAY"
:
"1"
,
"REC_ARRAY_DEPTH"
:
"24.0"
,
"REC_ARRAY_DIST"
:
"30.0"
,
"DRX"
:
"10"
,
"DRY"
:
"10"
,
"Seismograms"
:
"comment"
,
"NDT"
: