Comments o PETOOL.docVIP

Comments on PETOOL Chengyou yin Hefei Electronic Engineering Institute cyouyin@ Recently I have read your two papers. They are as follows. O. Ozgun, “Recursive two-way parabolic equation approach for modeling terrain effects in tropospheric propagation,” IEEE Trans. Antennas Propag., vol. 57, no. 9, pp. 2706–2713, 2009. Ozlem Ozgun, Gokhan Apaydin , Mustafa Kuzuoglu, Levent Sevgi, “PETOOL: MATLAB-based one-way and two-way split-step parabolic equation tool for radiowave propagation over variable terrain”, Computer Physics Communications 182 (2011) 2638–2654 I have acquired a lot of benefits from your papers. Then I also read your software PETOOL that are included in [2]. Now I think there are some questions in the software. Now I present in the following to discuss with author or readers that have read PETOOL. 1. about refractivity profile function In the SSPE-function programming, 962 m = n.*n - 1 + 2.*inputs.z / inputs.ae; This sentence is only for narrow-angle PE, not suitable for wide-angle PE. Wide-angle PE is (1) And the modified refractivity is (2) The SSPE-function programming’s sentences 254 255 is for wide-angle PE. %wide-angle 254 prop = exp(1j*beta*delx*(sqrt(1-p.^2./beta^2)-1)); 255 envpr = exp(1j * beta * delx * (mn-1)); % ENVIRONMENT TERM Therefore the programming’s sentence 962 may be incorrect. 2. about iteration computation of C1x、C2x In the SSPE-function programming, 331 dmft.c1x = exp(1j*delx*sqrt(beta^2+ad)); 333 dmft.c2x = exp(1j*delx*sqrt(beta^2+ad)); 335 dmft.c3x = exp(1j*delx*sqrt(beta^2+ad)); These sentences are used to compute the coefficients of C1x、C2x iteration for DMFT. However in wide-angle situation, the propagation factor should be (3) (4) So I think the correct sentences should be dmft.c1x = exp(1j*delx*(sqrt(beta^2+ad)-beta)); dmft.c2x = exp(1j*delx*(sqrt(beta^2+ad)-beta)); dmft.c3x = exp(1j*delx*(sqrt(beta^2+ad)-beta)); 3. about p