dr*1000\dz

pyat/czcore.py

@@ -170,7 +170,7 @@ def call_Bellhop_Rays(freq, sd, dr, bathm, ssp, sed, base,azi):
     run_type = 'R'  # incoherently sum beams, see AT docs for more info
     nbeams = 301
     alpha = np.linspace(-15, 15, 2)  # min and max launch angle -20 degrees to 20 degrees
-    box = Box(Zmax, 100)  # bound the region you let the beams go, depth in meters and range in km
+    box = Box(Zmax, Rmax+10)  # bound the region you let the beams go, depth in meters and range in km
     deltas = 0  # length step of ray trace, 0 means automatically choose
     beam = Beam(RunType=run_type, Nbeams=nbeams, alpha=alpha, box=box, deltas=deltas)  # package
     """删除ssp 已修改 8.26"""

pyat/ram.py

@@ -23,8 +23,8 @@ def call_ram(freq, sd, bathm, ssp, sed, dr=0.1):
     rmax = max(bathm.r)             # Maximum receiving range (unit: m)
     zmax = round(max(bathm.d) * 2)  # Maximum receiving depth (uint: m)
     zmplt = round(max(bathm.d)*1.1)     # Maximum depth of plot (uint: m)
-    dz = 1500 / freq / 15           # Step in depth direction: wavelength / 15
-    # dr = 10 * dz                     # Step in range direction
+    dz = 10          # Step in depth direction: wavelength / 15
+    dr = 1000* dr                     # Step in range direction
     Zmax = calZmax(ssp, sed)
     rd = Zmax
     pos = PosRAM(sd, rd, rmax, zmax, dr, dz, zmplt)