Skip to content
Snippets Groups Projects
Select Git revision
  • master default protected
  • plane-wave
  • 2.1.2
  • 2.1.1
  • 2.1.0
  • 2.0.0
  • 1.4
  • 1.4.0
  • 1.3.3
  • 1.3.2
  • 1.3.1
  • 1.3.0
  • 1.2.3
  • 1.2.2
  • 1.2.1
  • 1.2.0
  • 1.1.4
  • 1.1.3
  • 1.1.2
  • 1.1.1
  • 1.1.0
  • 1.0.0
22 results
Blame Permalink
decomposition.scr 2.08 KiB
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 






#!/bin/bash
#PBS -N fdelmod
#PBS -q verylong
#PBS -l nodes=1
#PBS -k eo
#PBS -j eo

which makewave
which makemod

cd /Users/jan/src/OpenSource/fdelmodc/demo

#makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1
#makewave w=fw fmin=0 flef=5 frig=80 fmax=100  dt=$dt file_out=wavefw.su nt=4096 t0=0.3

cp=2000
rho=1000
dx=2.5
dt=0.0005

makemod sizex=6000 sizez=4000 dx=$dx dz=$dx cp0=$cp cs0=$cs ro0=$rho \
	orig=-3000,-1000 file_base=syncl.su \
	intt=def x=-3000,0,3000 z=250,250,250 poly=0 cp=2300 ro=5000 \
	intt=def x=-3000,-2000,-1000,-800,0,800,3000 z=650,650,700,750,900,750,600 poly=2 cp=2600 ro=1000 \
	intt=def x=-3000,0,3000 z=1390,1390,1390 poly=0 cp=2000 ro=5000 \
	intt=def x=-3000,0,3000 z=1490,1490,1490 poly=0 cp=2400 ro=2000 

makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1


../fdelmodc \
    file_cp=syncl_cp.su ischeme=1 \
    file_den=syncl_ro.su \
    file_src=wave.su \
    file_rcv=shot_fd_Fz_zsrc1150.su \
	nshot=2 dxshot=500 \
    src_type=7 \
	src_orient=1 \
	src_injectionrate=1 \
    rec_type_ud=1 \
    dtrcv=0.0040 \
	rec_delay=0.1 \
    verbose=2 \
    tmod=3.10 \
    dxrcv=10 \
    xrcv1=-2250 xrcv2=2250 \
    zrcv1=1150 zrcv2=1150 \
    xsrc=0 zsrc=0 \
    ntaper=200 \
    left=4 right=4 top=4 bottom=4


supsgraph < anglerp.su hbox=2 wbox=4 style=normal \
	labelsize=10 label2='energy' label1='angle in degrees' \
	titlesize=-1 f1=0 d1=1 d1num=10.0 x1beg=1 > anglerp.eps
	
supsimage < shot_fd_Fz_zsrc1150_rvz.su \
    wbox=2 hbox=4 titlesize=-1 labelsize=10 verbose=1 \
    label1="time [s]" label2="lateral position [m]" > shot_fd_Fz_zsrc1150_rvz.eps

supsimage < shot_fd_Fz_zsrc1150_rp.su \
    wbox=2 hbox=4 titlesize=-1 labelsize=10 verbose=1 \
    label1="time [s]" label2="lateral position [m]" > shot_fd_Fz_zsrc1150_rp.eps

supsimage < shot_fd_Fz_zsrc1150_ru.su \
    wbox=2 hbox=4 titlesize=-1 labelsize=10 verbose=1 \
    label1="time [s]" label2="lateral position [m]" > shot_fd_Fz_zsrc1150_ru.eps

supsimage < shot_fd_Fz_zsrc1150_rd.su \
    wbox=2 hbox=4 titlesize=-1 labelsize=10 verbose=1 \
    label1="time [s]" label2="lateral position [m]" > shot_fd_Fz_zsrc1150_rd.eps