Dengda98 · Dengda98 · Apr 10, 2025 · Feb 23, 2025 · Feb 28, 2025 · Mar 3, 2025
diff --git a/.github/workflows/build.yml b/.github/workflows/build.yml
@@ -315,6 +315,18 @@ jobs:
           python plot_compare_pygrt.py
         continue-on-error: true  # 即使失败，仍然标记为成功
 
+      - name: Test strain stress (dynamic)
+        working-directory: ${{ env.EXAMPLE_COPY_PATH }}/compute_strain_stress/dynamic
+        run: |
+          chmod +x *.sh
+          ./run_grt.sh
+
+      - name: Test static displacement (static)
+        working-directory: ${{ env.EXAMPLE_COPY_PATH }}/compute_strain_stress/static
+        run: |
+          chmod +x *.sh
+          ./run_stgrt.sh
+
       - name: Remove the test files
         run: |
           rm -rf ${{ env.EXAMPLE_COPY_PATH }}

diff --git a/.github/workflows/testbuild.yml b/.github/workflows/testbuild.yml
@@ -281,6 +281,18 @@ jobs:
           ./run_grt.sh
           python run_pygrt.py
 
+      - name: Test strain stress (dynamic)
+        working-directory: ${{ env.EXAMPLE_COPY_PATH }}/compute_strain_stress/dynamic
+        run: |
+          chmod +x *.sh
+          ./run_grt.sh
+
+      - name: Test static displacement (static)
+        working-directory: ${{ env.EXAMPLE_COPY_PATH }}/compute_strain_stress/static
+        run: |
+          chmod +x *.sh
+          ./run_stgrt.sh
+
       - name: Remove the test files
         run: |
           rm -rf ${{ env.EXAMPLE_COPY_PATH }}

diff --git a/example/check_upar/README b/example/check_upar/README
@@ -0,0 +1 @@
+Use difference to check spatial derivatives calculation.
diff --git a/example/check_upar/compare_uir.pdf b/example/check_upar/compare_uir.pdf
diff --git a/example/check_upar/compare_uiz.pdf b/example/check_upar/compare_uiz.pdf
diff --git a/example/check_upar/milrow b/example/check_upar/milrow
@@ -0,0 +1,9 @@
+0.2  3.4  1.7  2.3  9e10  9e10
+0.6  3.7  1.9  2.4  9e10  9e10
+0.5  4.2  2.1  2.4  9e10  9e10
+0.5  4.6  2.3  2.5  9e10  9e10
+0.7  4.9    2.8  2.6  9e10  9e10
+0.5  5.1  2.9  2.7  9e10  9e10
+6.0  5.9  3.3  2.7  9e10  9e10
+28.  6.9  4.0  2.8  9e10  9e10
+-0.1   8.2  4.7  3.2  9e10  9e10
diff --git a/example/check_upar/plot.sh b/example/check_upar/plot.sh
@@ -0,0 +1,69 @@
+#!/bin/bash
+
+SAC_DISPLAY_COPYRIGHT=0
+
+out="compare_pdf"
+rm -rf $out
+mkdir -p $out
+
+diff=0.001
+depsrc=2
+deprcv=3.2
+deprcv2=`echo $deprcv+$diff | bc | awk '{printf("%.3f", $1)}'`
+dist=10
+dist2=`echo $dist+$diff | bc | awk '{printf("%.3f", $1)}'`
+
+GRN="GRN"
+
+# --------------------- ui_z --------------------------------------
+for ch in $(ls $GRN/milrow_${depsrc}_${deprcv2}_${dist}); do
+echo $ch
+sac <<EOF
+r $GRN/milrow_${depsrc}_${deprcv2}_${dist}/$ch
+subf $GRN/milrow_${depsrc}_${deprcv}_${dist}/$ch
+div -$diff
+w uiz_diff.sac
+r uiz_diff.sac
+r more $GRN/milrow_${depsrc}_${deprcv}_${dist}/z${ch:0:3}.sac
+int
+qdp off
+color red inc list blue red 
+p2
+saveimg $out/compare_$ch.pdf
+q
+EOF
+
+done
+
+rm uiz_diff.sac
+
+pdftk $out/* cat output compare_uiz.pdf
+
+rm -rf $out
+mkdir -p $out
+
+# --------------------- ui_z --------------------------------------
+for ch in $(ls $GRN/milrow_${depsrc}_${deprcv}_${dist2}); do
+echo $ch
+sac <<EOF
+r $GRN/milrow_${depsrc}_${deprcv}_${dist2}/$ch
+subf $GRN/milrow_${depsrc}_${deprcv}_${dist}/$ch
+div $diff
+w uir_diff.sac
+r uir_diff.sac
+r more $GRN/milrow_${depsrc}_${deprcv}_${dist}/r${ch:0:3}.sac
+int
+qdp off
+color red inc list blue red 
+p2
+saveimg $out/compare_$ch.pdf
+q
+EOF
+
+done
+
+rm uir_diff.sac
+
+pdftk $out/* cat output compare_uir.pdf
+
+rm -rf $out
diff --git a/example/check_upar/run_grt.sh b/example/check_upar/run_grt.sh
@@ -0,0 +1,19 @@
+#!/bin/bash 
+
+depsrc=2
+nt=1024
+dt=0.01
+
+modname="milrow"
+out="GRN"
+
+
+diff=0.001
+deprcv=3.2
+deprcv2=`echo $deprcv+$diff | bc | awk '{printf("%.3f", $1)}'`
+dist=10
+dist2=`echo $dist+$diff | bc | awk '{printf("%.3f", $1)}'`
+
+grt -M${modname} -O${out} -N${nt}/${dt} -D${depsrc}/$deprcv -R$dist -e
+grt -M${modname} -O${out} -N${nt}/${dt} -D${depsrc}/$deprcv2 -R$dist
+grt -M${modname} -O${out} -N${nt}/${dt} -D${depsrc}/$deprcv -R$dist2
diff --git a/example/compute_strain_stress/dynamic/milrow b/example/compute_strain_stress/dynamic/milrow
@@ -0,0 +1,9 @@
+0.2  3.4  1.7  2.3  9e10  9e10
+0.6  3.7  1.9  2.4  9e10  9e10
+0.5  4.2  2.1  2.4  9e10  9e10
+0.5  4.6  2.3  2.5  9e10  9e10
+0.7  4.9    2.8  2.6  9e10  9e10
+0.5  5.1  2.9  2.7  9e10  9e10
+6.0  5.9  3.3  2.7  9e10  9e10
+28.  6.9  4.0  2.8  9e10  9e10
+-0.1   8.2  4.7  3.2  9e10  9e10
diff --git a/example/compute_strain_stress/dynamic/run_grt.sh b/example/compute_strain_stress/dynamic/run_grt.sh
@@ -0,0 +1,26 @@
+#!/bin/bash
+
+
+dist=10
+depsrc=2
+deprcv=0
+
+nt=1024
+dt=0.01
+
+modname="milrow"
+out="GRN"
+
+grt -M${modname} -O${out} -N${nt}/${dt} -D${depsrc}/${deprcv} -R${dist} -e 
+
+# Fault
+S="u1e8"
+stk=77
+dip=88
+rak=99
+az=30
+grt.syn -G${out}/${modname}_${depsrc}_${deprcv}_${dist} -A$az -S$S -M$stk/$dip/$rak -Dp/0.6 -Osyn_dc -e
+
+
+grt.strain syn_dc/out
+grt.stress syn_dc/out
diff --git a/example/compute_strain_stress/static/milrow b/example/compute_strain_stress/static/milrow
@@ -0,0 +1,9 @@
+0.2  3.4  1.7  2.3  9e10  9e10
+0.6  3.7  1.9  2.4  9e10  9e10
+0.5  4.2  2.1  2.4  9e10  9e10
+0.5  4.6  2.3  2.5  9e10  9e10
+0.7  4.9    2.8  2.6  9e10  9e10
+0.5  5.1  2.9  2.7  9e10  9e10
+6.0  5.9  3.3  2.7  9e10  9e10
+28.  6.9  4.0  2.8  9e10  9e10
+-0.1   8.2  4.7  3.2  9e10  9e10
diff --git a/example/compute_strain_stress/static/run_stgrt.sh b/example/compute_strain_stress/static/run_stgrt.sh
@@ -0,0 +1,24 @@
+#!/bin/bash
+
+depsrc=2
+deprcv=0
+
+x1=-3
+x2=3
+nx=10
+
+y1=-3
+y2=3
+ny=10
+
+stgrt -Mmilrow -D${depsrc}/${deprcv} -X$x1/$x2/$nx -Y$y1/$y2/$ny -e > grn 
+
+# Fault
+S="u1e8"
+stk=77
+dip=88
+rak=99
+stgrt.syn -S$S -M$stk/$dip/$rak -e -N < grn > syn 
+
+stgrt.strain < syn > strain
+stgrt.stress < syn > stress
diff --git a/example/compute_strain_stress/static2/README b/example/compute_strain_stress/static2/README
@@ -0,0 +1,2 @@
+结果对比：
+    郝金来, 姚振兴, 2012. 均匀弹性分层介质模型中的同震位移、应变以及应力[J]. 地球物理学报, 55(5): 1682-1694.
diff --git a/example/compute_strain_stress/static2/halfspace2 b/example/compute_strain_stress/static2/halfspace2
@@ -0,0 +1 @@
+0.0  4.0 2.7 2.5 9e30 9e30
diff --git a/example/compute_strain_stress/static2/mod b/example/compute_strain_stress/static2/mod
@@ -0,0 +1,5 @@
+20 5.8 3.35 2.7 9e30 9e30
+20 6.8 3.90 2.8 9e30 9e30
+20 7.8 4.45 2.9 9e30 9e30
+20 8.8 5.00 3.0 9e30 9e30
+20 9.8 5.55 3.1 9e30 9e30
diff --git a/example/compute_strain_stress/static2/plot.ipynb b/example/compute_strain_stress/static2/plot.ipynb
diff --git a/example/compute_strain_stress/static2/run.sh b/example/compute_strain_stress/static2/run.sh
@@ -0,0 +1,31 @@
+#!/bin/bash 
+
+
+mod="halfspace2"
+stgrt -M$mod -D20/16 -X3/3/1 -Y-10/10/41 -e > grn1
+
+stgrt.syn -M0/70/0 -Su1e6 -e -N < grn1 > syn_1_stkslip
+stgrt.stress < syn_1_stkslip > stress_1_stkslip
+
+stgrt.syn -M0/70/90 -Su1e6 -e -N < grn1 > syn_1_dipslip
+stgrt.stress < syn_1_dipslip > stress_1_dipslip
+
+
+
+stgrt -M$mod -D10/14 -X3/3/1 -Y-10/10/41 -e > grn2
+
+stgrt.syn -M0/70/0 -Su1e6 -e -N < grn2 > syn_2_stkslip
+stgrt.stress < syn_2_stkslip > stress_2_stkslip
+
+stgrt.syn -M0/70/90 -Su1e6 -e -N < grn2 > syn_2_dipslip
+stgrt.stress < syn_2_dipslip > stress_2_dipslip
+
+
+mod="mod"
+stgrt -M$mod -D5/0 -X3/3/1 -Y-10/10/41 -e > grn3
+
+stgrt.syn -M0/70/0 -Su1e6 -e -N < grn3 > syn_3_stkslip
+stgrt.stress < syn_3_stkslip > stress_3_stkslip
+
+stgrt.syn -M0/70/90 -Su1e6 -e -N < grn3 > syn_3_dipslip
+stgrt.stress < syn_3_dipslip > stress_3_dipslip
diff --git a/example/static_disp/disp_dc.pdf b/example/static_disp/disp_dc.pdf
diff --git a/example/static_disp/disp_exp.pdf b/example/static_disp/disp_exp.pdf
diff --git a/example/static_disp/disp_sf.pdf b/example/static_disp/disp_sf.pdf
diff --git a/example/static_disp/halfspace2 b/example/static_disp/halfspace2
@@ -0,0 +1 @@
+0.0  4.0 2.7 2.5 9e30 9e30
diff --git a/example/static_disp/run_stgrt.sh b/example/static_disp/run_stgrt.sh
@@ -0,0 +1,73 @@
+#!/bin/bash 
+
+depsrc=2
+deprcv=0
+
+x1=-3
+x2=3
+nx=20
+
+y1=-3
+y2=3
+ny=20
+
+stgrt -Mhalfspace2 -D${depsrc}/${deprcv} -X$x1/$x2/$nx -Y$y1/$y2/$ny > grn 
+
+# Fault
+S="1e23"
+stk=60
+dip=90
+rak=0
+stgrt.syn -S$S -M$stk/$dip/$rak -N < grn > syn 
+
+gmt set FONT_TITLE 9p
+gmt begin disp_dc pdf
+    gmt xyz2grd syn -GsynZ.nc -R$y1/$y2/$x1/$x2 -I$ny+n/$nx+n -i0,1,2 -:
+    gmt xyz2grd syn -GsynN.nc -R$y1/$y2/$x1/$x2 -I$ny+n/$nx+n -i0,1,3 -:
+    gmt xyz2grd syn -GsynE.nc -R$y1/$y2/$x1/$x2 -I$ny+n/$nx+n -i0,1,4 -:
+
+    gmt basemap -Baf -BWSen+t"Fault, $stk/$dip/$rak, $S" -JX5c/5c -R$y1/$y2/$x1/$x2
+    gmt grdimage synZ.nc
+    gmt grdvector synE.nc synN.nc -Q0.1c+e+jc+h1+gblack -Si0.3c
+
+    gmt meca -Sa0.5c <<EOF
+0 0 $deprcv $stk $dip $rak 5
+EOF
+
+    gmt colorbar -Bx+l"Z (cm)"
+gmt end
+
+
+# Single Force
+S="1e17"
+fn=1
+fe=1
+fz=0
+stgrt.syn -S$S -F$fn/$fe/$fz -N < grn > syn 
+gmt begin disp_sf pdf
+    gmt xyz2grd syn -GsynZ.nc -R$y1/$y2/$x1/$x2 -I$ny+n/$nx+n -i0,1,2 -:
+    gmt xyz2grd syn -GsynN.nc -R$y1/$y2/$x1/$x2 -I$ny+n/$nx+n -i0,1,3 -:
+    gmt xyz2grd syn -GsynE.nc -R$y1/$y2/$x1/$x2 -I$ny+n/$nx+n -i0,1,4 -:
+
+    gmt basemap -Baf -BWSen+t"Force, $fn/$fe/$fz, $S" -JX5c/5c -R$y1/$y2/$x1/$x2
+    gmt grdimage synZ.nc
+    gmt grdvector synE.nc synN.nc -Q0.1c+e+jc+h1+gblack -Si0.5c
+
+    gmt colorbar -Bx+l"Z (cm)"
+gmt end
+
+# Explosion
+S="1e23"
+stgrt.syn -S$S -N < grn > syn 
+gmt begin disp_exp pdf
+    gmt xyz2grd syn -GsynZ.nc -R$y1/$y2/$x1/$x2 -I$ny+n/$nx+n -i0,1,2 -:
+    gmt xyz2grd syn -GsynN.nc -R$y1/$y2/$x1/$x2 -I$ny+n/$nx+n -i0,1,3 -:
+    gmt xyz2grd syn -GsynE.nc -R$y1/$y2/$x1/$x2 -I$ny+n/$nx+n -i0,1,4 -:
+
+    gmt basemap -Baf -BWSen+t"Explosion, $S" -JX5c/5c -R$y1/$y2/$x1/$x2
+    gmt grdimage synZ.nc
+    gmt grdvector synE.nc synN.nc -Q0.1c+e+jc+h1+gblack -Si0.3c
+
+    gmt colorbar -Bx+l"Z (cm)"
+gmt end
+
diff --git a/example/view_integ_stats/README b/example/view_integ_stats/README
@@ -0,0 +1 @@
+View slow convergency when source and receiver are at close or same depth.
diff --git a/example/view_integ_stats/halfspace2 b/example/view_integ_stats/halfspace2
@@ -0,0 +1 @@
+0.0  4.0 2.7 2.5 9e30 9e30
diff --git a/example/view_integ_stats/run_stgrt.sh b/example/view_integ_stats/run_stgrt.sh
@@ -0,0 +1,4 @@
+#!/bin/bash 
+
+
+stgrt -Mhalfspace2 -D0.1/0 -X-3/3/10 -Y-5/5/20 -e -S
diff --git a/example/view_integ_stats/view.ipynb b/example/view_integ_stats/view.ipynb
diff --git a/pygrt/C_extension/Makefile b/pygrt/C_extension/Makefile
@@ -48,6 +48,7 @@ CFLAGS :=  $(LINK_STATIC) $(LFFT_FLAGS) -lm -O3 \
 SUBDIRS = \
     src/common \
 	src/dynamic \
+	src/static \
 	src/travt
 
 
@@ -56,12 +57,12 @@ all: objs progs libs
 
 objs:
 	for dir in $(SUBDIRS); do \
-		$(MAKE) -C $$dir CC="$(CC)" CFLAGS="$(CFLAGS)" objs; \
+		$(MAKE) -C $$dir CC="$(CC)" CFLAGS="$(CFLAGS)" objs || exit 1; \
 	done
 
 progs: objs 
 	for dir in $(SUBDIRS); do \
-		$(MAKE) -C $$dir CC="$(CC)" CFLAGS="$(CFLAGS)" progs; \
+		$(MAKE) -C $$dir CC="$(CC)" CFLAGS="$(CFLAGS)" progs || exit 1; \
 	done
 
 

diff --git a/pygrt/C_extension/include/common/const.h b/pygrt/C_extension/include/common/const.h
@@ -72,7 +72,8 @@ typedef int MYINT;  ///< 整数
     #define FIVEQUARTERPI  3.92699082f  ///< \f$ \frac{5\pi}{4} \f$
     #define SEVENQUARTERPI  5.49778714f   ///< \f$ \frac{7\pi}{4} \f$
     #define INV_SQRT_TWO 0.70710678f   ///< \f$ \frac{1}{\sqrt{2}} \f$ 
-
+    #define DEG1 0.017453293f  ///< \f$ \frac{\pi}{180} \f$ 
+
 #else 
     typedef double _Complex MYCOMPLEX;
     typedef double MYREAL;
@@ -109,7 +110,8 @@ typedef int MYINT;  ///< 整数
     #define FIVEQUARTERPI  3.9269908169872414  ///< \f$ \frac{5\pi}{4} \f$
     #define SEVENQUARTERPI  5.497787143782138   ///< \f$ \frac{7\pi}{4} \f$
     #define INV_SQRT_TWO 0.7071067811865475   ///< \f$ \frac{1}{\sqrt{2}} \f$ 
-
+    #define DEG1 0.017453292519943295  ///< \f$ \frac{\pi}{180} \f$ 
+
 #endif
 
 #define CZERO CMPLX(RZERO, RZERO)  ///< 0.0 + j0.0
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1 @@
		Use difference to check spatial derivatives calculation.
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1,2 @@
		结果对比：
		郝金来, 姚振兴, 2012. 均匀弹性分层介质模型中的同震位移、应变以及应力[J]. 地球物理学报, 55(5): 1682-1694.
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1 @@
		View slow convergency when source and receiver are at close or same depth.
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1,4 @@
		#!/bin/bash


		stgrt -Mhalfspace2 -D0.1/0 -X-3/3/10 -Y-5/5/20 -e -S