你已经派生过 gprMax
镜像自地址
https://gitee.com/sunhf/gprMax.git
已同步 2025-08-07 23:14:03 +08:00
Updated syntax
这个提交包含在:
@@ -6,7 +6,7 @@
|
|||||||
from pathlib import Path
|
from pathlib import Path
|
||||||
|
|
||||||
import gprMax
|
import gprMax
|
||||||
from user_libs.GPRAntennaModels.GSSI import antenna_like_GSSI_400
|
from toolboxes.GPRAntennaModels.GSSI import antenna_like_GSSI_400
|
||||||
|
|
||||||
import numpy as np
|
import numpy as np
|
||||||
|
|
||||||
@@ -105,24 +105,27 @@ for i in range(1, 51):
|
|||||||
snap = gprMax.Snapshot(p1=(0, y / 2, 0), p2=(x, y / 2 + dl, z), dl=(dl, dl, dl),
|
snap = gprMax.Snapshot(p1=(0, y / 2, 0), p2=(x, y / 2 + dl, z), dl=(dl, dl, dl),
|
||||||
filename=Path(*parts[:-1], parts[-1] + '_' + str(i)).name,
|
filename=Path(*parts[:-1], parts[-1] + '_' + str(i)).name,
|
||||||
time=i * tw / 50)
|
time=i * tw / 50)
|
||||||
# scene.add(snap)
|
scene.add(snap)
|
||||||
|
|
||||||
# create a geometry view of the main grid and the sub grid stitched together
|
# create a geometry view of the main grid and the sub grid stitched together
|
||||||
gv = gprMax.GeometryView(p1=(0, 0, 0),
|
gv1 = gprMax.GeometryView(p1=(sg_x0, sg_y0, sg_z0), p2=(sg_x1, sg_y1, sg_z1),
|
||||||
p2=(1, 1, 1),
|
dl=(dl_s, dl_s, dl_s),
|
||||||
dl=dl,
|
filename=fn.with_suffix('').parts[-1] + '_sg',
|
||||||
filename=fn.with_suffix('').parts[-1],
|
output_type='n')
|
||||||
output_type='f',
|
#sg.add(gv1)
|
||||||
multi_grid=True)
|
|
||||||
|
|
||||||
# create a geometry view of the main grid and the sub grid stitched together
|
gv3 = gprMax.GeometryView(p1=(sg_x0, sg_y0, 1.512), p2=(sg_x1, sg_y1, 1.513),
|
||||||
gv_normal = gprMax.GeometryView(p1=(0, 0, 0),
|
dl=(dl_s, dl_s, dl_s),
|
||||||
|
filename=fn.with_suffix('').parts[-1] + '_sg',
|
||||||
|
output_type='f')
|
||||||
|
sg.add(gv3)
|
||||||
|
|
||||||
|
gv2 = gprMax.GeometryView(p1=(0, 0, 0),
|
||||||
p2=domain.props.p1,
|
p2=domain.props.p1,
|
||||||
dl=dl,
|
dl=dl,
|
||||||
filename=fn.with_suffix('').parts[-1] + '_voxels',
|
filename=fn.with_suffix('').parts[-1],
|
||||||
output_type='n')
|
output_type='n')
|
||||||
# scene.add(gv)
|
scene.add(gv2)
|
||||||
scene.add(gv_normal)
|
|
||||||
|
|
||||||
# half space material
|
# half space material
|
||||||
layer_m = gprMax.Material(er=9, se=0.397e-3, mr=1, sm=0, id='soil_2')
|
layer_m = gprMax.Material(er=9, se=0.397e-3, mr=1, sm=0, id='soil_2')
|
||||||
@@ -134,4 +137,4 @@ scene.add(fb)
|
|||||||
sr = gprMax.AddSurfaceRoughness(p1=(0, 0, 1), p2=(3, 1, 1), frac_dim=1.5, weighting=(1, 1), limits=(0.4, 1.2), fractal_box_id='fbox', seed=1)
|
sr = gprMax.AddSurfaceRoughness(p1=(0, 0, 1), p2=(3, 1, 1), frac_dim=1.5, weighting=(1, 1), limits=(0.4, 1.2), fractal_box_id='fbox', seed=1)
|
||||||
scene.add(sr)
|
scene.add(sr)
|
||||||
|
|
||||||
gprMax.run(scenes=[scene], n=1, geometry_only=False, outputfile=fn, subgrid=True, autotranslate=True)
|
gprMax.run(scenes=[scene], n=1, geometry_only=True, outputfile=fn, subgrid=True, autotranslate=True)
|
||||||
|
在新工单中引用
屏蔽一个用户