Advertisement
Guest User

Untitled

a guest
Mar 24th, 2019
198
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 28.47 KB | None | 0 0
  1. ~/src3/nec2/nec-2/nec $ cat nec2_vs_nec2small.diff
  2. --- nec2.f 2019-03-23 23:32:40.036835977 -0400
  3. +++ nec2small.f 2019-03-23 22:40:40.559978030 -0400
  4. @@ -29,17 +29,14 @@
  5. C
  6. C ***
  7. IMPLICIT REAL (A-H,O-Z)
  8. -
  9. -C Change all these params
  10. +C Change all these params and LD
  11. C .. N is number segments
  12. C .. M is number patches
  13. -C .. N2M = N+2*M, N3M = N+3*M
  14. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  15. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  16.  
  17. CHARACTER AIN*2, ATST*2, INFILE*256, OTFILE*256
  18. C***
  19. C INTEGER AIN,ATST,PNET
  20. -
  21. C..Command Line parsing..
  22. INTEGER NARGS
  23. LOGICAL FISVALID
  24. @@ -54,7 +51,7 @@
  25. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  26. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  27. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
  28. - COMMON /CMB/ CM(1000000)
  29. + COMMON /CMB/ CM(90000)
  30. COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
  31. &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
  32. COMMON /SAVE/ IP( N2M), KCOM, COM(20,5), EPSR, SIG, SCRWLT,
  33. @@ -103,7 +100,6 @@
  34. C***
  35. DATA LOADMX, NSMAX, NETMX/200,150,150/, NORMF/200/
  36. 706 CONTINUE
  37. -
  38. C..Fortran 77 Extension to get number args..
  39. NARGS=IARGC()
  40. C..Defaults for terminal IO
  41. @@ -147,14 +143,13 @@
  42. & OTFILE)
  43. STOP
  44.  
  45. +C***
  46. 705 CONTINUE
  47. CALL SECNDS(EXTIM)
  48. FJ=(0.,1.)
  49. -C** LD is max permissible N+M (see NM param)
  50. -C** param from original NEC code
  51. LD=NM
  52. NXA(1)=0
  53. - IRESRV=1000000
  54. + IRESRV=90000
  55. C***
  56. 1 KCOM=0
  57. IFRTMW=0
  58. @@ -1104,7 +1099,7 @@
  59. C ARC GENERATES SEGMENT GEOMETRY DATA FOR AN ARC OF NS SEGMENTS
  60. C
  61. IMPLICIT REAL (A-H,O-Z)
  62. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  63. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  64. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  65. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  66. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
  67. @@ -1201,7 +1196,7 @@
  68. C CURRENT VECTOR CUR.
  69. C
  70. IMPLICIT REAL (A-H,O-Z)
  71. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  72. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  73. COMPLEX CUR, CURX, VQDS, CURD, CCJ, VSANT, VQD, CS1, CS2
  74. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  75. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  76. @@ -1297,7 +1292,7 @@
  77. C ***
  78. IMPLICIT REAL (A-H,O-Z)
  79. C CMNGF FILLS INTERACTION MATRICIES B, C, AND D FOR N.G.F. SOLUTION
  80. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  81. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  82. COMPLEX CB, CC, CD, ZARRAY, EXK, EYK, EZK, EXS, EYS, EZS, EXC
  83. &, EYC, EZC
  84. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  85. @@ -1579,7 +1574,7 @@
  86. C
  87. C CMSET SETS UP THE COMPLEX STRUCTURE MATRIX IN THE ARRAY CM
  88. C
  89. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  90. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  91. COMPLEX CM, ZARRAY, ZAJ, EXK, EYK, EZK, EXS,
  92. &EYS, EZS, EXC, EYC, EZC, SSX, D, DETER
  93. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  94. @@ -1694,7 +1689,7 @@
  95. C ***
  96. IMPLICIT REAL (A-H,O-Z)
  97. C CMSS COMPUTES MATRIX ELEMENTS FOR SURFACE-SURFACE INTERACTIONS.
  98. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  99. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  100. COMPLEX G11, G12, G21, G22, CM, EXK, EYK, EZK, EXS, EYS, EZS,
  101. & EXC, EYC, EZC
  102. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  103. @@ -1781,7 +1776,7 @@
  104. C ***
  105. IMPLICIT REAL (A-H,O-Z)
  106. C COMPUTES MATRIX ELEMENTS FOR E ALONG WIRES DUE TO PATCH CURRENT
  107. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  108. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  109. COMPLEX CM, ZRATI, ZRATI2, T1, EXK, EYK, EZK, EXS, EYS, EZS,
  110. &EXC, EYC, EZC, EMEL, CW, FRATI
  111. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  112. @@ -1934,7 +1929,7 @@
  113. C CMWS COMPUTES MATRIX ELEMENTS FOR WIRE-SURFACE INTERACTIONS
  114. C
  115. IMPLICIT REAL (A-H,O-Z)
  116. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  117. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  118. COMPLEX CM, CW, ETK, ETS, ETC, EXK, EYK, EZK, EXS, EYS, EZS,
  119. &EXC, EYC, EZC
  120. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  121. @@ -2025,7 +2020,7 @@
  122. C
  123. C CMWW COMPUTES MATRIX ELEMENTS FOR WIRE-WIRE INTERACTIONS
  124. C
  125. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  126. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  127. COMPLEX CM, CW, ETK, ETS, ETC, EXK, EYK, EZK, EXS, EYS, EZS,
  128. &EXC, EYC, EZC
  129. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  130. @@ -2145,7 +2140,7 @@
  131. C CONNECT SETS UP SEGMENT CONNECTION DATA IN ARRAYS ICON1 AND ICON2
  132. C BY SEARCHING FOR SEGMENT ENDS THAT ARE IN CONTACT.
  133. C
  134. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  135. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  136. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  137. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  138. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
  139. @@ -2455,7 +2450,7 @@
  140. C
  141. C COUPLE COMPUTES THE MAXIMUM COUPLING BETWEEN PAIRS OF SEGMENTS.
  142. C
  143. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  144. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  145. COMPLEX Y11A, Y12A, CUR, Y11, Y12, Y22, YL, YIN, ZL, ZIN, RHO
  146. &, VQD, VSANT, VQDS
  147. COMMON /YPARM/ NCOUP, ICOUP, NCTAG(5), NCSEG(5), Y11A(5), Y12A(
  148. @@ -2533,7 +2528,7 @@
  149. C DATAGN IS THE MAIN ROUTINE FOR INPUT OF GEOMETRY DATA.
  150. C
  151. C***
  152. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  153. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  154. C***
  155. CHARACTER *2 GM, ATST
  156. CHARACTER *1 IFX,IFY,IFZ,IPT
  157. @@ -2924,7 +2919,7 @@
  158. C COMPUTE NEAR E FIELDS OF A SEGMENT WITH SINE, COSINE, AND
  159. C CONSTANT CURRENTS. GROUND EFFECT INCLUDED.
  160. C
  161. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  162. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  163. COMPLEX TXK, TYK, TZK, TXS, TYS, TZS, TXC, TYC, TZC, EXK, EYK
  164. &, EZK, EXS, EYS, EZS, EXC, EYC, EZC, EPX, EPY, ZRATI, REFS, REFPS
  165. &, ZRSIN, ZRATX, T1, ZSCRN, ZRATI2, TEZS, TERS, TEZC, TERC, TEZK,
  166. @@ -3160,7 +3155,7 @@
  167. IMPLICIT REAL (A-H,O-Z)
  168. C COMPUTE E FIELD OF SINE, COSINE, AND CONSTANT CURRENT FILAMENTS BY
  169. C THIN WIRE APPROXIMATION.
  170. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  171. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  172. COMPLEX CON, GZ1, GZ2, GP1, GP2, GZP1, GZP2, EZS, ERS, EZC,
  173. &ERC, EZK, ERK
  174. COMMON /TMI/ ZPK, RKB2, IJX
  175. @@ -3207,7 +3202,7 @@
  176. C COMPUTE E FIELD OF SINE, COSINE, AND CONSTANT CURRENT FILAMENTS BY
  177. C EXTENDED THIN WIRE APPROXIMATION.
  178. IMPLICIT REAL (A-H,O-Z)
  179. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  180. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  181. COMPLEX CON, GZ1, GZ2, GZP1, GZP2, GR1, GR2, GRP1, GRP2, EZS,
  182. & EZC, ERS, ERC, GRK1, GRK2, EZK, ERK, GZZ1, GZZ2
  183. COMMON /TMI/ ZPK, RKB2, IJX
  184. @@ -3284,7 +3279,7 @@
  185. C DOUBLE PRECISION 6/4/85
  186. C
  187. SUBROUTINE ETMNS( P1, P2, P3, P4, P5, P6, IPR, E)
  188. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  189. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  190. C ***
  191. C
  192. C ETMNS FILLS THE ARRAY E WITH THE NEGATIVE OF THE ELECTRIC FIELD
  193. @@ -3682,7 +3677,7 @@
  194. C TEXT. (MATRIX TRANSPOSED.
  195. C
  196. IMPLICIT REAL (A-H,O-Z)
  197. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  198. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  199. COMPLEX A, D, ARJ
  200. DIMENSION A( NDIM, NDIM), IP( NDIM)
  201. COMMON /SCRATM/ D( N2M)
  202. @@ -3892,7 +3887,7 @@
  203. C FBLOCK SETS PARAMETERS FOR OUT-OF-CORE SOLUTION FOR THE PRIMARY
  204. C MATRIX (A)
  205. IMPLICIT REAL (A-H,O-Z)
  206. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  207. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  208. COMPLEX SSX, DETER
  209. COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
  210. &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
  211. @@ -4000,7 +3995,7 @@
  212. C FBNGF SETS THE BLOCKING PARAMETERS FOR THE B, C, AND D ARRAYS FOR
  213. C OUT-OF-CORE STORAGE.
  214. IMPLICIT REAL (A-H,O-Z)
  215. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  216. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  217. COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
  218. &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
  219. IRESX= IRESRV- IMAT
  220. @@ -4077,7 +4072,7 @@
  221. C THE FACTOR EXP(J*K*R)/(R/LAMDA) NOT INCLUDED
  222. C
  223. IMPLICIT REAL (A-H,O-Z)
  224. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  225. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  226. COMPLEX CIX, CIY, CIZ, EXA, ETH, EPH, CONST, CCX, CCY, CCZ,
  227. &CDP, CUR
  228. COMPLEX ZRATI, ZRSIN, RRV, RRH, RRV1, RRH1, RRV2, RRH2,
  229. @@ -4292,7 +4287,7 @@
  230. C CALCULATES THE XYZ COMPONENTS OF THE ELECTRIC FIELD DUE TO
  231. C SURFACE CURRENTS
  232. IMPLICIT REAL (A-H,O-Z)
  233. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  234. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  235. COMPLEX CT, CONS, SCUR, EX, EY, EZ
  236. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  237. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  238. @@ -4328,7 +4323,7 @@
  239. C GF COMPUTES THE INTEGRAND EXP(JKR)/(KR) FOR NUMERICAL INTEGRATION.
  240. C
  241. IMPLICIT REAL (A-H,O-Z)
  242. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  243. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  244. COMMON /TMI/ ZPK, RKB2, IJ
  245. ZDK= ZK- ZPK
  246. RK= SQRT( RKB2+ ZDK* ZDK)
  247. @@ -4352,7 +4347,7 @@
  248. C GFIL READS THE N.G.F. FILE
  249. C
  250. IMPLICIT REAL (A-H,O-Z)
  251. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  252. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  253. integer*4 COM
  254. COMPLEX CM, SSX, ZRATI, ZRATI2, T1, ZARRAY, AR1, AR2, AR3,
  255. &EPSCF, FRATI
  256. @@ -4509,7 +4504,7 @@
  257. C GFLD COMPUTES THE RADIATED FIELD INCLUDING GROUND WAVE.
  258. C
  259. IMPLICIT REAL (A-H,O-Z)
  260. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  261. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  262. COMPLEX CUR, EPI, CIX, CIY, CIZ, EXA, XX1, XX2, U, U2, ERV,
  263. &EZV, ERH, EPH
  264. COMPLEX EZH, EX, EY, ETH, UX, ERD
  265. @@ -4663,7 +4658,7 @@
  266. C WRITE N.G.F. FILE
  267. C
  268. IMPLICIT REAL (A-H,O-Z)
  269. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  270. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  271. integer*4 COM
  272. COMPLEX CM, SSX, ZRATI, ZRATI2, T1, ZARRAY, AR1, AR2, AR3,
  273. &EPSCF, FRATI
  274. @@ -4759,7 +4754,7 @@
  275. C ***
  276. C INTEGRAND FOR H FIELD OF A WIRE
  277. IMPLICIT REAL (A-H,O-Z)
  278. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  279. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  280. COMMON /TMH/ ZPK, RHKS
  281. RS= ZK- ZPK
  282. RS= RHKS+ RS* RS
  283. @@ -4783,7 +4778,7 @@
  284. C (PROC. IRE, SEPT., 1937, PP.1203,1236.)
  285. C
  286. IMPLICIT REAL (A-H,O-Z)
  287. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  288. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  289. COMPLEX FJ, TPJ, U2, U, RK1, RK2, T1, T2, T3, T4, P1, RV, OMR
  290. &, W, F, Q1, RH, V, G, XR1, XR2, X1, X2, X3, X4, X5, X6, X7, EZV,
  291. &ERV, EZH, ERH, EPH, XX1, XX2, ECON, FBAR
  292. @@ -4860,7 +4855,7 @@
  293. SUBROUTINE GX( ZZ, RH, XK, GZ, GZP)
  294. C ***
  295. C SEGMENT END CONTRIBUTIONS FOR THIN WIRE APPROX.
  296. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  297. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  298. COMPLEX GZ, GZP
  299. R2= ZZ* ZZ+ RH* RH
  300. R= SQRT( R2)
  301. @@ -4876,7 +4871,7 @@
  302. &)
  303. C ***
  304. C SEGMENT END CONTRIBUTIONS FOR EXT. THIN WIRE APPROX.
  305. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  306. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  307. COMPLEX GZ, C1, C2, C3, G1, G1P, G2, G2P, G3, GZP
  308. R2= ZZ* ZZ+ RH* RH
  309. R= SQRT( R2)
  310. @@ -4923,7 +4918,7 @@
  311. C SUBROUTINE HELIX GENERATES SEGMENT GEOMETRY DATA FOR A HELIX OF NS
  312. C SEGMENTS
  313. IMPLICIT REAL (A-H,O-Z)
  314. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  315. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  316. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  317. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  318. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
  319. @@ -5086,7 +5081,7 @@
  320. C ***
  321. C HINTG COMPUTES THE H FIELD OF A PATCH CURRENT
  322. IMPLICIT REAL (A-H,O-Z)
  323. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  324. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  325. COMPLEX EXK, EYK, EZK, EXS, EYS, EZS, EXC, EYC, EZC, ZRATI,
  326. &ZRATI2, GAM, F1X, F1Y, F1Z, F2X, F2Y, F2Z, RRV, RRH, T1, FRATI
  327. COMMON /DATAJ/ S, B, XJ, YJ, ZJ, CABJ, SABJ, SALPJ, EXK, EYK,
  328. @@ -5176,7 +5171,7 @@
  329. C HSFLD COMPUTES THE H FIELD FOR CONSTANT, SINE, AND COSINE CURRENT
  330. C ON A SEGMENT INCLUDING GROUND EFFECTS.
  331. IMPLICIT REAL (A-H,O-Z)
  332. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  333. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  334. COMPLEX EXK, EYK, EZK, EXS, EYS, EZS, EXC, EYC, EZC, ZRATI,
  335. &ZRATI2, T1, HPK, HPS, HPC, QX, QY, QZ, RRV, RRH, ZRATX, FRATI
  336. COMMON /DATAJ/ S, B, XJ, YJ, ZJ, CABJ, SABJ, SALPJ, EXK, EYK,
  337. @@ -5290,7 +5285,7 @@
  338. C ***
  339. C CALCULATES H FIELD OF SINE COSINE, AND CONSTANT CURRENT OF SEGMENT
  340. IMPLICIT REAL (A-H,O-Z)
  341. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  342. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  343. COMPLEX FJ, FJK, EKR1, EKR2, T1, T2, CONS, HPS, HPC, HPK
  344. DIMENSION FJX(2), FJKX(2)
  345. EQUIVALENCE(FJ,FJX),(FJK,FJKX)
  346. @@ -5351,7 +5346,7 @@
  347. C 4 FUNCTIONS AT THE POINT (X,Y).
  348. C
  349. IMPLICIT REAL (A-H,O-Z)
  350. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  351. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  352. COMPLEX F1, F2, F3, F4, A, B, C, D, FX1, FX2, FX3, FX4, P1,
  353. &P2, P3, P4, A11, A12, A13, A14, A21, A22, A23, A24, A31, A32, A33
  354. &, A34, A41, A42, A43, A44, B11, B12, B13, B14, B21, B22, B23, B24
  355. @@ -5620,7 +5615,7 @@
  356. C TAG NUMBER ITAGI. IF ITAGI=0 SEGMENT NUMBER M IS RETURNED.
  357. C
  358. IMPLICIT REAL (A-H,O-Z)
  359. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  360. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  361. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  362. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  363. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
  364. @@ -5659,7 +5654,7 @@
  365. C RALSTONS TEXT.
  366. C
  367. IMPLICIT REAL (A-H,O-Z)
  368. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  369. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  370. COMPLEX A, D, AJR
  371. INTEGER R, R1, R2, PJ, PR
  372. LOGICAL L1, L2, L3
  373. @@ -5767,7 +5762,7 @@
  374. C TYPES OF LOADING
  375. C
  376. IMPLICIT REAL (A-H,O-Z)
  377. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  378. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  379. COMPLEX ZARRAY, ZT, TPCJ, ZINT
  380. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  381. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  382. @@ -5914,7 +5909,7 @@
  383. C BLOCKS OF DESCENDING ORDER.
  384. C
  385. IMPLICIT REAL (A-H,O-Z)
  386. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  387. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  388. COMPLEX A, B, Y, SUM
  389. COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
  390. &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
  391. @@ -5984,7 +5979,7 @@
  392. C S/R WHICH UNSCRAMBLES, SCRAMBLED FACTORED MATRIX
  393. C
  394. IMPLICIT REAL (A-H,O-Z)
  395. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  396. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  397. COMPLEX A, TEMP
  398. COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
  399. &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
  400. @@ -6053,7 +6048,7 @@
  401. C RESPECTIVELY, THEN SHIFTED BY XS,YS,ZS
  402. C
  403. IMPLICIT REAL (A-H,O-Z)
  404. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  405. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  406. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  407. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  408. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
  409. @@ -6158,7 +6153,7 @@
  410. C THE STRUCTURE CURRENTS HAVE BEEN COMPUTED.
  411. C
  412. IMPLICIT REAL (A-H,O-Z)
  413. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  414. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  415. COMPLEX EX, EY, EZ, CUR, ACX, BCX, CCX, EXK, EYK, EZK, EXS,
  416. &EYS, EZS, EXC, EYC, EZC, ZRATI, ZRATI2, T1, FRATI
  417. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  418. @@ -6285,7 +6280,7 @@
  419. C PRESENT.
  420. C
  421. IMPLICIT REAL (A-H,O-Z)
  422. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  423. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  424. COMPLEX CMN, RHNT, YMIT, RHS, ZPED, EINC, VSANT, VLT, CUR,
  425. &VSRC, RHNX, VQD, VQDS, CUX, CM, CMB, CMC, CMD
  426. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  427. @@ -6629,7 +6624,7 @@
  428. C ***
  429. C COMPUTE NEAR E OR H FIELDS OVER A RANGE OF POINTS
  430. IMPLICIT REAL (A-H,O-Z)
  431. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  432. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  433. COMPLEX EX, EY, EZ
  434. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  435. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  436. @@ -6732,7 +6727,7 @@
  437. C THE STRUCTURE CURRENTS HAVE BEEN COMPUTED.
  438. C
  439. IMPLICIT REAL (A-H,O-Z)
  440. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  441. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  442. COMPLEXHX,HY,HZ,CUR,ACX, BCX, CCX, EXK, EYK, EZK, EXS, EYS,
  443. &EZS, EXC, EYC, EZC
  444. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  445. @@ -6817,7 +6812,7 @@
  446. C ***
  447. C PATCH GENERATES AND MODIFIES PATCH GEOMETRY DATA
  448. IMPLICIT REAL (A-H,O-Z)
  449. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  450. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  451. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  452. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  453. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
  454. @@ -7018,7 +7013,7 @@
  455. C ***
  456. C INTEGRATE OVER PATCHES AT WIRE CONNECTION POINT
  457. IMPLICIT REAL (A-H,O-Z)
  458. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  459. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  460. COMPLEX EXK, EYK, EZK, EXS, EYS, EZS, EXC, EYC, EZC, E, E1,
  461. &E2, E3, E4, E5, E6, E7, E8, E9
  462. COMMON /DATAJ/ S, B, XJ, YJ, ZJ, CABJ, SABJ, SALPJ, EXK, EYK,
  463. @@ -7111,7 +7106,7 @@
  464. C PRNT SETS UP THE PRINT FORMATS FOR IMPEDANCE LOADING
  465. C
  466. IMPLICIT REAL (A-H,O-Z)
  467. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  468. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  469. CHARACTER*6 IFORM, IVAR
  470. CHARACTER *(*) IA
  471. DIMENSION IVAR(13), IA(1), IFORM(8), IN(3), INT(3), FL(6), FLT(6
  472. @@ -7184,7 +7179,7 @@
  473. C ***
  474. C FILL INCIDENT FIELD ARRAY FOR CHARGE DISCONTINUITY VOLTAGE SOURCE
  475. IMPLICIT REAL (A-H,O-Z)
  476. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  477. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  478. COMPLEX VQDS, CURD, CCJ, V, EXK, EYK, EZK, EXS, EYS, EZS, EXC
  479. &, EYC, EZC, ETK, ETS, ETC, VSANT, VQD, E, ZARRAY
  480. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  481. @@ -7314,7 +7309,7 @@
  482. C ***
  483. C COMPUTE RADIATION PATTERN, GAIN, NORMALIZED GAIN
  484. IMPLICIT REAL (A-H,O-Z)
  485. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  486. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  487. C INTEGER HBLK,HCIR,HCLIF
  488. CHARACTER*6 IGNTP, IGAX, IGTP, HPOL, HCIR, HCLIF, HBLK
  489. CHARACTER*6 ISENS
  490. @@ -7592,7 +7587,7 @@
  491. SUBROUTINE READGM( GM, I1, I2, X1, Y1, Z1, X2, Y2, Z2, RAD)
  492. C ***
  493. IMPLICIT REAL (A-H,O-Z)
  494. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  495. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  496. INTEGER*4 NTOT
  497. INTEGER*4 NINT
  498. INTEGER*4 NFLT
  499. @@ -7682,7 +7677,7 @@
  500. SUBROUTINE READMN( GM, I1, I2, I3, I4, F1, F2, F3, F4, F5, F6)
  501. C ***
  502. IMPLICIT REAL (A-H,O-Z)
  503. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  504. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  505. INTEGER*4 NTOT
  506. INTEGER*4 NINT
  507. INTEGER*4 NFLT
  508. @@ -7772,7 +7767,7 @@
  509. C REBLOCK ARRAY B IN N.G.F. SOLUTION FROM BLOCKS OF ROWS ON TAPE14
  510. C TO BLOCKS OF COLUMNS ON TAPE16
  511. IMPLICIT REAL (A-H,O-Z)
  512. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  513. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  514. COMPLEX B, BX
  515. COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
  516. &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
  517. @@ -7809,7 +7804,7 @@
  518. C STRUCTURE TO COMPLETE A SYMMETRIC STRUCTURE.
  519. C
  520. IMPLICIT REAL (A-H,O-Z)
  521. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  522. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  523. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  524. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  525. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
  526. @@ -8137,7 +8132,7 @@
  527. C ***
  528. C COMPUTE COMPONENT OF BASIS FUNCTION I ON SEGMENT IS.
  529. IMPLICIT REAL (A-H,O-Z)
  530. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  531. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  532. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  533. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  534. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
  535. @@ -8275,7 +8270,7 @@
  536. C THE SOURCE SEGMENT AT T RELATIVE TO THE SEGMENT CENTER.
  537. C
  538. IMPLICIT REAL (A-H,O-Z)
  539. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  540. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  541. COMPLEX E, ERV, EZV, ERH, EZH, EPH, T1, EXK, EYK, EZK, EXS,
  542. &EYS, EZS, EXC, EYC, EZC, XX1, XX2, U, U2, ZRATI, ZRATI2, FRATI,
  543. &ER, ET, HRV, HZV, HRH
  544. @@ -8395,7 +8390,7 @@
  545. C ***
  546. C SOLVE FOR CURRENT IN N.G.F. PROCEDURE
  547. IMPLICIT REAL (A-H,O-Z)
  548. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  549. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  550. COMPLEX A, B, C, D, SUM, XY, Y
  551. COMMON /SCRATM/ Y( N2M)
  552. COMMON /SEGJ/ AX(30), BX(30), CX(30), JCO(30), JSNO, ISCON(50),
  553. @@ -8521,7 +8516,7 @@
  554. C SOLUTION IS RETURNED THROUGH VECTOR B. (MATRIX TRANSPOSED.
  555. C
  556. IMPLICIT REAL (A-H,O-Z)
  557. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  558. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  559. COMPLEX A, B, Y, SUM
  560. INTEGER PI
  561. COMMON /SCRATM/ Y( N2M)
  562. @@ -8567,7 +8562,7 @@
  563. C MATRIX EQ.
  564. C
  565. IMPLICIT REAL (A-H,O-Z)
  566. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  567. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  568. COMPLEX A, B, Y, SUM, SSX
  569. COMMON /SMAT/ SSX(16,16)
  570. COMMON /SCRATM/ Y( N2M)
  571. @@ -8685,7 +8680,7 @@
  572. C ***
  573. C COMPUTE BASIS FUNCTION I
  574. IMPLICIT REAL (A-H,O-Z)
  575. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  576. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  577. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  578. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  579. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
  580. @@ -8827,7 +8822,7 @@
  581. C TEST FOR CONVERGENCE IN NUMERICAL INTEGRATION
  582. C
  583. IMPLICIT REAL (A-H,O-Z)
  584. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  585. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  586. DEN= ABS( F2R)
  587. TR= ABS( F2I)
  588. IF( DEN.LT. TR) DEN= TR
  589. @@ -8847,7 +8842,7 @@
  590. C ***
  591. C COMPUTE THE COMPONENTS OF ALL BASIS FUNCTIONS ON SEGMENT J
  592. IMPLICIT REAL (A-H,O-Z)
  593. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  594. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  595. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  596. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  597. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
  598. @@ -8896,7 +8891,7 @@
  599. C CALCULATES THE ELECTRIC FIELD DUE TO UNIT CURRENT IN THE T1 AND T2
  600. C DIRECTIONS ON A PATCH
  601. IMPLICIT REAL (A-H,O-Z)
  602. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  603. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  604. COMPLEX EXK, EYK, EZK, EXS, EYS, EZS, EXC, EYC, EZC, ZRATI,
  605. &ZRATI2, T1, ER, Q1, Q2, RRV, RRH, EDP, FRATI
  606. COMMON /DATAJ/ S, B, XJ, YJ, ZJ, CABJ, SABJ, SALPJ, EXK, EYK,
  607. @@ -8987,7 +8982,7 @@
  608. C WIRE OF NS SEGMENTS.
  609. C
  610. IMPLICIT REAL (A-H,O-Z)
  611. - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
  612. + PARAMETER ( NM=600, N2M=800, N3M=1000)
  613. COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
  614. &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
  615. & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement