Cosmetic changes and ready for marchenko paper R1

marchenko/demo/oneD/conv.gnp

+set style data linespoints
+set mytics 10
+set xlabel 'number of iterations'
+set ylabel 'convergence rate'
+set size 2.0,2.0
+set size ratio 0.6
+set grid
+
+set log y
+set nolog x
+
+set term postscript eps font 'Helvetica,12' linewidth 4 fontscale 3
+set output 'convergence.eps'
+plot 'conv.txt' using  1:($2) lw 3 notitle
+

marchenko/demo/oneD/conv.txt

+0  1.000000e+00
+1  8.104102e-01
+2  2.776407e-01
+3  1.775258e-01
+4  1.278046e-01
+5  8.376110e-02
+6  6.221900e-02
+7  4.089906e-02
+8  3.275844e-02
+9  2.070254e-02
+10 1.920658e-02
+11 1.091778e-02
+12 1.282995e-02
+13 6.060715e-03
+14 9.706275e-03
+15 3.603180e-03

marchenko/demo/twoD/initialFocus.scr

+#!/bin/bash
+
+export PATH=$HOME/src/OpenSource/bin:$PATH:
+
+dx=2.5
+dt=0.0005
+
+makemod sizex=6000 sizez=2000 dx=$dx dz=$dx cp0=1900  ro0=1200 \
+        orig=-3000,0 file_base=synclDown.su verbose=2 \
+        intt=def x=-3000,500,3000 z=195,195,195 poly=1 cp=1950 ro=3700 \
+        intt=def x=-3000,3000 z=600,600 poly=0 cp=2050 ro=1750 \
+        intt=def x=-3000,3000 z=680,680 poly=0 cp=2150 ro=2220 \
+        intt=def x=-3000,3000 z=780,780 poly=0 cp=2230 ro=1700 \
+        intt=def x=-3000,-2200,-1500,0,1300,2100,3000 z=520,580,680,840,680,600,500 poly=2 cp=2400 ro=2800 \
+
+makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1 scale=1
+
+dxshot=10
+ishot=300
+nshots=301
+
+export OMP_NUM_THREADS=1
+mkdir -p shots
+mkdir -p jobs
+
+while (( ishot < nshots ))
+do
+
+		(( xsrc = -3000 + ${ishot}*${dxshot} ))
+		echo xsrc=$xsrc
+		file_rcv=shots/iniFocus_z1100_x${xsrc}.su
+
+  cat << EOF > jobs/pbs_$ishot.job 
+#!/bin/bash
+#
+#PBS -q medium
+#PBS -N mod_${xsrc}
+#PBS -j eo 
+#PBS -m n 
+#PBS -l nodes=1
+#PBS -V
+
+
+export PATH=\$HOME/src/OpenSource/bin:\$PATH:
+cd \$PBS_O_WORKDIR
+
+export OMP_NUM_THREADS=4
+	
+		fdelmodc \
+			file_cp=synclDown_cp.su ischeme=1 iorder=4 \
+   			file_den=synclDown_ro.su \
+   			file_src=wave.su \
+   			file_rcv=$file_rcv \
+   			src_type=1 \
+			src_orient=1 \
+			src_injectionrate=1 \
+   			rec_type_vz=0 \
+   			rec_type_p=1 \
+   			rec_int_vz=2 \
+			rec_delay=0.1 \
+   			dtrcv=0.004 \
+   			verbose=2 \
+   			tmod=2.100 \
+   			dxrcv=10.0 \
+   			xrcv1=-3000 xrcv2=3000 \
+   			zrcv1=0 zrcv2=0 \
+   			xsrc=$xsrc zsrc=1100 \
+   			ntaper=200 \
+   			left=2 right=2 top=2 bottom=2
+EOF
+
+        qsub jobs/pbs_$ishot.job 
+
+			(( ishot = $ishot + 1))
+done
+
+
+

utils/demo/syn2d.scr

+#!/bin/bash
+
+export OMP_NUM_THREADS=8
+nt=256
+
+../makewave nt=$nt w=g2 file_out=wave.su w=g2 dt=0.004 
+
+../green c=2000 nt=$nt dt=0.004 zsrc1=500 xsrc1=-60 xsrc2=60 dxsrc=15 file_src=wave.su  file_out=green.su
+
+../green c=2000 nt=$nt dt=0.004 zsrc1=500 xsrc1=-1500 xsrc2=1500 dxsrc=15 file_src=wave.su  file_out=shots.su
+
+
+../syn2d file_syn=green.su file_shot=shots.su  verbose=1 nshots=201 file_cfp=syn.su
+
+