Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- ~/src3/nec2/nec-2/nec $ cat nec2_vs_nec2small.diff
- --- nec2.f 2019-03-23 23:32:40.036835977 -0400
- +++ nec2small.f 2019-03-23 22:40:40.559978030 -0400
- @@ -29,17 +29,14 @@
- C
- C ***
- IMPLICIT REAL (A-H,O-Z)
- -
- -C Change all these params
- +C Change all these params and LD
- C .. N is number segments
- C .. M is number patches
- -C .. N2M = N+2*M, N3M = N+3*M
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- CHARACTER AIN*2, ATST*2, INFILE*256, OTFILE*256
- C***
- C INTEGER AIN,ATST,PNET
- -
- C..Command Line parsing..
- INTEGER NARGS
- LOGICAL FISVALID
- @@ -54,7 +51,7 @@
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
- - COMMON /CMB/ CM(1000000)
- + COMMON /CMB/ CM(90000)
- COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
- &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
- COMMON /SAVE/ IP( N2M), KCOM, COM(20,5), EPSR, SIG, SCRWLT,
- @@ -103,7 +100,6 @@
- C***
- DATA LOADMX, NSMAX, NETMX/200,150,150/, NORMF/200/
- 706 CONTINUE
- -
- C..Fortran 77 Extension to get number args..
- NARGS=IARGC()
- C..Defaults for terminal IO
- @@ -147,14 +143,13 @@
- & OTFILE)
- STOP
- +C***
- 705 CONTINUE
- CALL SECNDS(EXTIM)
- FJ=(0.,1.)
- -C** LD is max permissible N+M (see NM param)
- -C** param from original NEC code
- LD=NM
- NXA(1)=0
- - IRESRV=1000000
- + IRESRV=90000
- C***
- 1 KCOM=0
- IFRTMW=0
- @@ -1104,7 +1099,7 @@
- C ARC GENERATES SEGMENT GEOMETRY DATA FOR AN ARC OF NS SEGMENTS
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
- @@ -1201,7 +1196,7 @@
- C CURRENT VECTOR CUR.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CUR, CURX, VQDS, CURD, CCJ, VSANT, VQD, CS1, CS2
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- @@ -1297,7 +1292,7 @@
- C ***
- IMPLICIT REAL (A-H,O-Z)
- C CMNGF FILLS INTERACTION MATRICIES B, C, AND D FOR N.G.F. SOLUTION
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CB, CC, CD, ZARRAY, EXK, EYK, EZK, EXS, EYS, EZS, EXC
- &, EYC, EZC
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- @@ -1579,7 +1574,7 @@
- C
- C CMSET SETS UP THE COMPLEX STRUCTURE MATRIX IN THE ARRAY CM
- C
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CM, ZARRAY, ZAJ, EXK, EYK, EZK, EXS,
- &EYS, EZS, EXC, EYC, EZC, SSX, D, DETER
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- @@ -1694,7 +1689,7 @@
- C ***
- IMPLICIT REAL (A-H,O-Z)
- C CMSS COMPUTES MATRIX ELEMENTS FOR SURFACE-SURFACE INTERACTIONS.
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX G11, G12, G21, G22, CM, EXK, EYK, EZK, EXS, EYS, EZS,
- & EXC, EYC, EZC
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- @@ -1781,7 +1776,7 @@
- C ***
- IMPLICIT REAL (A-H,O-Z)
- C COMPUTES MATRIX ELEMENTS FOR E ALONG WIRES DUE TO PATCH CURRENT
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CM, ZRATI, ZRATI2, T1, EXK, EYK, EZK, EXS, EYS, EZS,
- &EXC, EYC, EZC, EMEL, CW, FRATI
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- @@ -1934,7 +1929,7 @@
- C CMWS COMPUTES MATRIX ELEMENTS FOR WIRE-SURFACE INTERACTIONS
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CM, CW, ETK, ETS, ETC, EXK, EYK, EZK, EXS, EYS, EZS,
- &EXC, EYC, EZC
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- @@ -2025,7 +2020,7 @@
- C
- C CMWW COMPUTES MATRIX ELEMENTS FOR WIRE-WIRE INTERACTIONS
- C
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CM, CW, ETK, ETS, ETC, EXK, EYK, EZK, EXS, EYS, EZS,
- &EXC, EYC, EZC
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- @@ -2145,7 +2140,7 @@
- C CONNECT SETS UP SEGMENT CONNECTION DATA IN ARRAYS ICON1 AND ICON2
- C BY SEARCHING FOR SEGMENT ENDS THAT ARE IN CONTACT.
- C
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
- @@ -2455,7 +2450,7 @@
- C
- C COUPLE COMPUTES THE MAXIMUM COUPLING BETWEEN PAIRS OF SEGMENTS.
- C
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX Y11A, Y12A, CUR, Y11, Y12, Y22, YL, YIN, ZL, ZIN, RHO
- &, VQD, VSANT, VQDS
- COMMON /YPARM/ NCOUP, ICOUP, NCTAG(5), NCSEG(5), Y11A(5), Y12A(
- @@ -2533,7 +2528,7 @@
- C DATAGN IS THE MAIN ROUTINE FOR INPUT OF GEOMETRY DATA.
- C
- C***
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- C***
- CHARACTER *2 GM, ATST
- CHARACTER *1 IFX,IFY,IFZ,IPT
- @@ -2924,7 +2919,7 @@
- C COMPUTE NEAR E FIELDS OF A SEGMENT WITH SINE, COSINE, AND
- C CONSTANT CURRENTS. GROUND EFFECT INCLUDED.
- C
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX TXK, TYK, TZK, TXS, TYS, TZS, TXC, TYC, TZC, EXK, EYK
- &, EZK, EXS, EYS, EZS, EXC, EYC, EZC, EPX, EPY, ZRATI, REFS, REFPS
- &, ZRSIN, ZRATX, T1, ZSCRN, ZRATI2, TEZS, TERS, TEZC, TERC, TEZK,
- @@ -3160,7 +3155,7 @@
- IMPLICIT REAL (A-H,O-Z)
- C COMPUTE E FIELD OF SINE, COSINE, AND CONSTANT CURRENT FILAMENTS BY
- C THIN WIRE APPROXIMATION.
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CON, GZ1, GZ2, GP1, GP2, GZP1, GZP2, EZS, ERS, EZC,
- &ERC, EZK, ERK
- COMMON /TMI/ ZPK, RKB2, IJX
- @@ -3207,7 +3202,7 @@
- C COMPUTE E FIELD OF SINE, COSINE, AND CONSTANT CURRENT FILAMENTS BY
- C EXTENDED THIN WIRE APPROXIMATION.
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CON, GZ1, GZ2, GZP1, GZP2, GR1, GR2, GRP1, GRP2, EZS,
- & EZC, ERS, ERC, GRK1, GRK2, EZK, ERK, GZZ1, GZZ2
- COMMON /TMI/ ZPK, RKB2, IJX
- @@ -3284,7 +3279,7 @@
- C DOUBLE PRECISION 6/4/85
- C
- SUBROUTINE ETMNS( P1, P2, P3, P4, P5, P6, IPR, E)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- C ***
- C
- C ETMNS FILLS THE ARRAY E WITH THE NEGATIVE OF THE ELECTRIC FIELD
- @@ -3682,7 +3677,7 @@
- C TEXT. (MATRIX TRANSPOSED.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX A, D, ARJ
- DIMENSION A( NDIM, NDIM), IP( NDIM)
- COMMON /SCRATM/ D( N2M)
- @@ -3892,7 +3887,7 @@
- C FBLOCK SETS PARAMETERS FOR OUT-OF-CORE SOLUTION FOR THE PRIMARY
- C MATRIX (A)
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX SSX, DETER
- COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
- &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
- @@ -4000,7 +3995,7 @@
- C FBNGF SETS THE BLOCKING PARAMETERS FOR THE B, C, AND D ARRAYS FOR
- C OUT-OF-CORE STORAGE.
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
- &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
- IRESX= IRESRV- IMAT
- @@ -4077,7 +4072,7 @@
- C THE FACTOR EXP(J*K*R)/(R/LAMDA) NOT INCLUDED
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CIX, CIY, CIZ, EXA, ETH, EPH, CONST, CCX, CCY, CCZ,
- &CDP, CUR
- COMPLEX ZRATI, ZRSIN, RRV, RRH, RRV1, RRH1, RRV2, RRH2,
- @@ -4292,7 +4287,7 @@
- C CALCULATES THE XYZ COMPONENTS OF THE ELECTRIC FIELD DUE TO
- C SURFACE CURRENTS
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CT, CONS, SCUR, EX, EY, EZ
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- @@ -4328,7 +4323,7 @@
- C GF COMPUTES THE INTEGRAND EXP(JKR)/(KR) FOR NUMERICAL INTEGRATION.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /TMI/ ZPK, RKB2, IJ
- ZDK= ZK- ZPK
- RK= SQRT( RKB2+ ZDK* ZDK)
- @@ -4352,7 +4347,7 @@
- C GFIL READS THE N.G.F. FILE
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- integer*4 COM
- COMPLEX CM, SSX, ZRATI, ZRATI2, T1, ZARRAY, AR1, AR2, AR3,
- &EPSCF, FRATI
- @@ -4509,7 +4504,7 @@
- C GFLD COMPUTES THE RADIATED FIELD INCLUDING GROUND WAVE.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CUR, EPI, CIX, CIY, CIZ, EXA, XX1, XX2, U, U2, ERV,
- &EZV, ERH, EPH
- COMPLEX EZH, EX, EY, ETH, UX, ERD
- @@ -4663,7 +4658,7 @@
- C WRITE N.G.F. FILE
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- integer*4 COM
- COMPLEX CM, SSX, ZRATI, ZRATI2, T1, ZARRAY, AR1, AR2, AR3,
- &EPSCF, FRATI
- @@ -4759,7 +4754,7 @@
- C ***
- C INTEGRAND FOR H FIELD OF A WIRE
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /TMH/ ZPK, RHKS
- RS= ZK- ZPK
- RS= RHKS+ RS* RS
- @@ -4783,7 +4778,7 @@
- C (PROC. IRE, SEPT., 1937, PP.1203,1236.)
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX FJ, TPJ, U2, U, RK1, RK2, T1, T2, T3, T4, P1, RV, OMR
- &, W, F, Q1, RH, V, G, XR1, XR2, X1, X2, X3, X4, X5, X6, X7, EZV,
- &ERV, EZH, ERH, EPH, XX1, XX2, ECON, FBAR
- @@ -4860,7 +4855,7 @@
- SUBROUTINE GX( ZZ, RH, XK, GZ, GZP)
- C ***
- C SEGMENT END CONTRIBUTIONS FOR THIN WIRE APPROX.
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX GZ, GZP
- R2= ZZ* ZZ+ RH* RH
- R= SQRT( R2)
- @@ -4876,7 +4871,7 @@
- &)
- C ***
- C SEGMENT END CONTRIBUTIONS FOR EXT. THIN WIRE APPROX.
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX GZ, C1, C2, C3, G1, G1P, G2, G2P, G3, GZP
- R2= ZZ* ZZ+ RH* RH
- R= SQRT( R2)
- @@ -4923,7 +4918,7 @@
- C SUBROUTINE HELIX GENERATES SEGMENT GEOMETRY DATA FOR A HELIX OF NS
- C SEGMENTS
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
- @@ -5086,7 +5081,7 @@
- C ***
- C HINTG COMPUTES THE H FIELD OF A PATCH CURRENT
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX EXK, EYK, EZK, EXS, EYS, EZS, EXC, EYC, EZC, ZRATI,
- &ZRATI2, GAM, F1X, F1Y, F1Z, F2X, F2Y, F2Z, RRV, RRH, T1, FRATI
- COMMON /DATAJ/ S, B, XJ, YJ, ZJ, CABJ, SABJ, SALPJ, EXK, EYK,
- @@ -5176,7 +5171,7 @@
- C HSFLD COMPUTES THE H FIELD FOR CONSTANT, SINE, AND COSINE CURRENT
- C ON A SEGMENT INCLUDING GROUND EFFECTS.
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX EXK, EYK, EZK, EXS, EYS, EZS, EXC, EYC, EZC, ZRATI,
- &ZRATI2, T1, HPK, HPS, HPC, QX, QY, QZ, RRV, RRH, ZRATX, FRATI
- COMMON /DATAJ/ S, B, XJ, YJ, ZJ, CABJ, SABJ, SALPJ, EXK, EYK,
- @@ -5290,7 +5285,7 @@
- C ***
- C CALCULATES H FIELD OF SINE COSINE, AND CONSTANT CURRENT OF SEGMENT
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX FJ, FJK, EKR1, EKR2, T1, T2, CONS, HPS, HPC, HPK
- DIMENSION FJX(2), FJKX(2)
- EQUIVALENCE(FJ,FJX),(FJK,FJKX)
- @@ -5351,7 +5346,7 @@
- C 4 FUNCTIONS AT THE POINT (X,Y).
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX F1, F2, F3, F4, A, B, C, D, FX1, FX2, FX3, FX4, P1,
- &P2, P3, P4, A11, A12, A13, A14, A21, A22, A23, A24, A31, A32, A33
- &, A34, A41, A42, A43, A44, B11, B12, B13, B14, B21, B22, B23, B24
- @@ -5620,7 +5615,7 @@
- C TAG NUMBER ITAGI. IF ITAGI=0 SEGMENT NUMBER M IS RETURNED.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
- @@ -5659,7 +5654,7 @@
- C RALSTONS TEXT.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX A, D, AJR
- INTEGER R, R1, R2, PJ, PR
- LOGICAL L1, L2, L3
- @@ -5767,7 +5762,7 @@
- C TYPES OF LOADING
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX ZARRAY, ZT, TPCJ, ZINT
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- @@ -5914,7 +5909,7 @@
- C BLOCKS OF DESCENDING ORDER.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX A, B, Y, SUM
- COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
- &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
- @@ -5984,7 +5979,7 @@
- C S/R WHICH UNSCRAMBLES, SCRAMBLED FACTORED MATRIX
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX A, TEMP
- COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
- &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
- @@ -6053,7 +6048,7 @@
- C RESPECTIVELY, THEN SHIFTED BY XS,YS,ZS
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
- @@ -6158,7 +6153,7 @@
- C THE STRUCTURE CURRENTS HAVE BEEN COMPUTED.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX EX, EY, EZ, CUR, ACX, BCX, CCX, EXK, EYK, EZK, EXS,
- &EYS, EZS, EXC, EYC, EZC, ZRATI, ZRATI2, T1, FRATI
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- @@ -6285,7 +6280,7 @@
- C PRESENT.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX CMN, RHNT, YMIT, RHS, ZPED, EINC, VSANT, VLT, CUR,
- &VSRC, RHNX, VQD, VQDS, CUX, CM, CMB, CMC, CMD
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- @@ -6629,7 +6624,7 @@
- C ***
- C COMPUTE NEAR E OR H FIELDS OVER A RANGE OF POINTS
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX EX, EY, EZ
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- @@ -6732,7 +6727,7 @@
- C THE STRUCTURE CURRENTS HAVE BEEN COMPUTED.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEXHX,HY,HZ,CUR,ACX, BCX, CCX, EXK, EYK, EZK, EXS, EYS,
- &EZS, EXC, EYC, EZC
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- @@ -6817,7 +6812,7 @@
- C ***
- C PATCH GENERATES AND MODIFIES PATCH GEOMETRY DATA
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
- @@ -7018,7 +7013,7 @@
- C ***
- C INTEGRATE OVER PATCHES AT WIRE CONNECTION POINT
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX EXK, EYK, EZK, EXS, EYS, EZS, EXC, EYC, EZC, E, E1,
- &E2, E3, E4, E5, E6, E7, E8, E9
- COMMON /DATAJ/ S, B, XJ, YJ, ZJ, CABJ, SABJ, SALPJ, EXK, EYK,
- @@ -7111,7 +7106,7 @@
- C PRNT SETS UP THE PRINT FORMATS FOR IMPEDANCE LOADING
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- CHARACTER*6 IFORM, IVAR
- CHARACTER *(*) IA
- DIMENSION IVAR(13), IA(1), IFORM(8), IN(3), INT(3), FL(6), FLT(6
- @@ -7184,7 +7179,7 @@
- C ***
- C FILL INCIDENT FIELD ARRAY FOR CHARGE DISCONTINUITY VOLTAGE SOURCE
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX VQDS, CURD, CCJ, V, EXK, EYK, EZK, EXS, EYS, EZS, EXC
- &, EYC, EZC, ETK, ETS, ETC, VSANT, VQD, E, ZARRAY
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- @@ -7314,7 +7309,7 @@
- C ***
- C COMPUTE RADIATION PATTERN, GAIN, NORMALIZED GAIN
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- C INTEGER HBLK,HCIR,HCLIF
- CHARACTER*6 IGNTP, IGAX, IGTP, HPOL, HCIR, HCLIF, HBLK
- CHARACTER*6 ISENS
- @@ -7592,7 +7587,7 @@
- SUBROUTINE READGM( GM, I1, I2, X1, Y1, Z1, X2, Y2, Z2, RAD)
- C ***
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- INTEGER*4 NTOT
- INTEGER*4 NINT
- INTEGER*4 NFLT
- @@ -7682,7 +7677,7 @@
- SUBROUTINE READMN( GM, I1, I2, I3, I4, F1, F2, F3, F4, F5, F6)
- C ***
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- INTEGER*4 NTOT
- INTEGER*4 NINT
- INTEGER*4 NFLT
- @@ -7772,7 +7767,7 @@
- C REBLOCK ARRAY B IN N.G.F. SOLUTION FROM BLOCKS OF ROWS ON TAPE14
- C TO BLOCKS OF COLUMNS ON TAPE16
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX B, BX
- COMMON /MATPAR/ ICASE, NBLOKS, NPBLK, NLAST, NBLSYM, NPSYM,
- &NLSYM, IMAT, ICASX, NBBX, NPBX, NLBX, NBBL, NPBL, NLBL
- @@ -7809,7 +7804,7 @@
- C STRUCTURE TO COMPLETE A SYMMETRIC STRUCTURE.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
- @@ -8137,7 +8132,7 @@
- C ***
- C COMPUTE COMPONENT OF BASIS FUNCTION I ON SEGMENT IS.
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
- @@ -8275,7 +8270,7 @@
- C THE SOURCE SEGMENT AT T RELATIVE TO THE SEGMENT CENTER.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX E, ERV, EZV, ERH, EZH, EPH, T1, EXK, EYK, EZK, EXS,
- &EYS, EZS, EXC, EYC, EZC, XX1, XX2, U, U2, ZRATI, ZRATI2, FRATI,
- &ER, ET, HRV, HZV, HRH
- @@ -8395,7 +8390,7 @@
- C ***
- C SOLVE FOR CURRENT IN N.G.F. PROCEDURE
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX A, B, C, D, SUM, XY, Y
- COMMON /SCRATM/ Y( N2M)
- COMMON /SEGJ/ AX(30), BX(30), CX(30), JCO(30), JSNO, ISCON(50),
- @@ -8521,7 +8516,7 @@
- C SOLUTION IS RETURNED THROUGH VECTOR B. (MATRIX TRANSPOSED.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX A, B, Y, SUM
- INTEGER PI
- COMMON /SCRATM/ Y( N2M)
- @@ -8567,7 +8562,7 @@
- C MATRIX EQ.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX A, B, Y, SUM, SSX
- COMMON /SMAT/ SSX(16,16)
- COMMON /SCRATM/ Y( N2M)
- @@ -8685,7 +8680,7 @@
- C ***
- C COMPUTE BASIS FUNCTION I
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
- @@ -8827,7 +8822,7 @@
- C TEST FOR CONVERGENCE IN NUMERICAL INTEGRATION
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- DEN= ABS( F2R)
- TR= ABS( F2I)
- IF( DEN.LT. TR) DEN= TR
- @@ -8847,7 +8842,7 @@
- C ***
- C COMPUTE THE COMPONENTS OF ALL BASIS FUNCTIONS ON SEGMENT J
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
- @@ -8896,7 +8891,7 @@
- C CALCULATES THE ELECTRIC FIELD DUE TO UNIT CURRENT IN THE T1 AND T2
- C DIRECTIONS ON A PATCH
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMPLEX EXK, EYK, EZK, EXS, EYS, EZS, EXC, EYC, EZC, ZRATI,
- &ZRATI2, T1, ER, Q1, Q2, RRV, RRH, EDP, FRATI
- COMMON /DATAJ/ S, B, XJ, YJ, ZJ, CABJ, SABJ, SALPJ, EXK, EYK,
- @@ -8987,7 +8982,7 @@
- C WIRE OF NS SEGMENTS.
- C
- IMPLICIT REAL (A-H,O-Z)
- - PARAMETER ( NM=10000, N2M=15000, N3M=20000)
- + PARAMETER ( NM=600, N2M=800, N3M=1000)
- COMMON /DATA/ LD, N1, N2, N, NP, M1, M2, M, MP, X( NM), Y( NM),
- &Z( NM), SI( NM), BI( NM), ALP( NM), BET( NM), ICON1( N2M), ICON2(
- & N2M), ITAG( N2M), ICONX( NM), WLAM, IPSYM
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement