home *** CD-ROM | disk | FTP | other *** search
/ HAM Radio 1 / HamRadio.cdr / tech / design5 / optosc.lst < prev    next >
Encoding:
File List  |  1979-12-31  |  19.0 KB  |  468 lines
  1.                                                                          Page   1
  2.                                                                         01-01-80
  3.                                                                         00:05:05
  4. D Line# 1     7                            Microsoft FORTRAN77 V3.31 August 1985
  5.       1 C
  6.       2 C     OPTOSC PROGRAM (OPTIMUM OSCILLATOR)
  7.       3 C
  8.       4 C     COMPUTATION OF THE OPTIMAL NETWORK FOR A
  9.       5 C     TRANSISTORIZED WITH OR WITHOUT DIELECTRIC RESONATOR.
  10.       6 C
  11.       7 C
  12.       8 C     PROGRAM BY: PATRICK CHAMPAGNE
  13.       9 C                 DEFENSE COMMUNICATIONS DIVISION
  14.      10 C                 CANADIAN MARCONI COMPANY
  15.      11 C              2442 TRENTON, MONTREAL, QUEBEC, CANADA
  16.      12 C
  17.      13 C     2 POSSIBLE OSCILLATION CONDITIONS ARE ANALYZED:
  18.      14 C     A) TRANSISTOR POTENTIALLY UNSTABLE
  19.      15 C     B) TRANSISTOR UNCONDITIONNALLY STABLE (THE COMPUTATION OF
  20.      16 C        THE OPTIMUM NETWORK FOR NEGATIVE CONDUCTANCE IS BASED ON
  21.      17 C        R.SPENCE'S "LINEAR ACTIVE NETWORKS")
  22.      18 C
  23.      19 C
  24.      20 C     THE FOLLOWING PARAMETERS ARE USED IN THE PROGRAM:
  25.      21 C               Sij    :S PARAMETERS OF THE TRANSISTOR
  26.      22 C               Yij    :Y PARAMETERS OF THE TRANSISTOR
  27.      23 C               EEFF   :EFFECTIVE DIELECTRIC CONSTANT OF MICROSTRIP
  28.      24 C               WL     :WAVELENGTH
  29.      25 C               BGi    :FEEDBACK SUSCEPTANCE FOR MAX. NEG. CONDUCTANCE
  30.      26 C               GNi    :MAXIMUM NEGATIVE CONDUCTANCES
  31.      27 C               ZIN    :RESONATOR IMPEDANCE AS SEEN BY TRANSISTOR
  32.      28 C               XIN    :RESONATOR ADMITTANCE AS SEEN BY TRANSISTOR
  33.      29 C               B0     :TRANSISTOR OUTPUT SUSCEPTANCE
  34.      30 C               XM(i,j):INDEFINITE ADMITTANCE MATRIX
  35.      31 C               YE,Y0,
  36.      32 C               YC,YB  :TRANSISTOR ADMITTANCES
  37.      33 C
  38.      34 $DEBUG
  39.      35       COMPLEX S11,S21,S12,S22,Y11,Y12,Y21,Y22,Y1,Y2,Y3,Y4,DEL
  40.      36       COMPLEX XM(3,3),D,YE,YB,YC,Y0,GM,ZIN,XIN
  41.      37       CHARACTER*10 AA
  42.      38       WRITE(*,3)
  43.      39       OPEN(6,FILE='PRN',STATUS='OLD')
  44.      40       WRITE(6,700)
  45.      41 700   FORMAT(' *****************************************************')
  46.      42       WRITE(6,710)
  47.      43 710   FORMAT(' *                                                   *')
  48.      44       WRITE(6,3)
  49.      45 3     FORMAT(' *     COMPUTATION OF TRANSISTORIZED OSCILLATORS     * ')
  50.      46       WRITE(6,710)
  51.      47       WRITE(6,700)
  52.      48       WRITE(*,7)
  53.      49 7     FORMAT(1H0,'TRANSISTOR TYPE?')
  54.      50       READ(*,9)AA
  55.      51 9     FORMAT(A)
  56.      52       WRITE(*,11)AA
  57.      53       WRITE(6,11)AA
  58.      54 11    FORMAT(1H0,'THE TRANSISTOR USED IS ',A)
  59.      55       WRITE(*,13)
  60.      56 13    FORMAT(1H0,'BIAS CONDITIONS (VCE,IC) OR (VDS,ID)?')
  61.                                                                          Page   2
  62.                                                                         01-01-80
  63.                                                                         00:05:05
  64. D Line# 1     7                            Microsoft FORTRAN77 V3.31 August 1985
  65.      57       READ(*,17)VCE,AIC
  66.      58 17    FORMAT(2F6.3)
  67.      59       WRITE(6,19)
  68.      60 19    FORMAT(1H0,'BIAS CONDITIONS')
  69.      61       WRITE(6,21)VCE,AIC
  70.      62 21    FORMAT(1H0,'VCE(VDS)=',F6.3,'V',3X,'IC(ID)=',F6.3,'MA')
  71.      63       WRITE(*,61)
  72.      64 61    FORMAT(1H0,'OPERATING FREQUENCY (IN GHZ)')
  73.      65       READ (*,62) FREQ
  74.      66 62    FORMAT(F6.3)
  75.      67       WRITE(*,63) FREQ
  76.      68       WRITE(6,63) FREQ
  77.      69 63    FORMAT(1H0,'THE OPERATING FREQUENCY IS ',F6.3,' GHZ')
  78.      70 C
  79.      71 C     THE COMPUTATION OF THE EFFECTIVE DIELECTRIC CONSTANT IS BASED
  80.      72 C     ON EDWARDS' 'FOUNDATIONS FOR MICROSTRIP DESIGN'
  81.      73 C
  82.      74       T1=0.43*FREQ**2-0.009*FREQ**3
  83.      75       T2=(.635/50)**1.33
  84.      76       T1=1+T1*T2
  85.      77 C
  86.      78 C     E EFF. =6.53 IS USED FOR F=0
  87.      79 C
  88.      80       T1=(9.8-6.53)/T1
  89.      81       EEFF=9.8-T1
  90.      82       WL=298/(FREQ*EEFF**.5)
  91.      83  
  92.      84       WRITE(6,96)WL
  93.      85 96    FORMAT(1H0,'THE WAVELENGTH IS ',F6.3,' MM')
  94.      86 C
  95.      87 C     THE TRANSISTOR CONFIGURATION CORRESPOND TO THE CONDITION
  96.      88 C     FOR WHICH THE S-PARAMETERS ARE GIVEN
  97.      89 C
  98.      90       WRITE(6,10)
  99.      91       WRITE(*,10)
  100.      92 10    FORMAT(1H0,'TRANSISTOR CONFIGURATION')
  101.      93       WRITE(*,12)
  102.      94       WRITE(6,12)
  103.      95 12    FORMAT(1H0,'1=CB,CG; 2=CE,CS; 3=CC,CD')
  104.      96       READ(*,23) ICONF
  105.      97 23    FORMAT(I2)
  106.      98       WRITE(6,14) ICONF
  107.      99 14    FORMAT(1H0,10X,I2)
  108.     100 C
  109.     101 C     S PARAMETERS READING
  110.     102 C
  111.     103       WRITE(*,15)
  112.     104       WRITE(6,15)
  113.     105 15    FORMAT(1H0,'S11:AMPLITUDE AND PHASE (DEGREES)')
  114.     106       READ(*,5) P11
  115.     107       READ(*,5) ARG11
  116.     108       WRITE(6,27)P11,ARG11
  117.     109 27    FORMAT(1H0,5X,2(F8.3,5X))
  118.     110       WRITE(*,20)
  119.     111       WRITE(6,20)
  120.     112 20    FORMAT(1H0,'S21:AMPLITUDE AND PHASE')
  121.                                                                          Page   3
  122.                                                                         01-01-80
  123.                                                                         00:05:05
  124. D Line# 1     7                            Microsoft FORTRAN77 V3.31 August 1985
  125.     113       READ(*,5) P21
  126.     114       READ(*,5) ARG21
  127.     115       WRITE(6,27) P21,ARG21
  128.     116       WRITE(*,25)
  129.     117       WRITE(6,25)
  130.     118 25    FORMAT(1H0,'S12:AMPLITUDE AND PHASE')
  131.     119       READ(*,5) P12
  132.     120       READ(*,5) ARG12
  133.     121       WRITE(6,27) P12,ARG12
  134.     122       WRITE(*,30)
  135.     123       WRITE(6,30)
  136.     124 30    FORMAT(1H0,'S22:AMPLITUDE AND PHASE')
  137.     125       READ(*,5) P22
  138.     126       READ(*,5) ARG22
  139.     127       WRITE(6,27) P22,ARG22
  140.     128 5     FORMAT (F8.3)
  141.     129 C
  142.     130 C     TRANSFORMATION OF ANGLES FROM DEGREES TO RADIANS
  143.     131 C
  144.     132       XPI=3.1416/180
  145.     133       ARG11=ARG11*XPI
  146.     134       ARG21=ARG21*XPI
  147.     135       ARG12=ARG12*XPI
  148.     136       ARG22=ARG22*XPI
  149.     137 C
  150.     138 C     COMPLEX S PARAMETERS
  151.     139 C
  152.     140       S11=CMPLX(P11*COS(ARG11),P11*SIN(ARG11))
  153.     141       S21=CMPLX(P21*COS(ARG21),P21*SIN(ARG21))
  154.     142       S12=CMPLX(P12*COS(ARG12),P12*SIN(ARG12))
  155.     143       S22=CMPLX(P22*COS(ARG22),P22*SIN(ARG22))
  156.     144 C
  157.     145 C     ROLLETT STABILITY FACTOR
  158.     146 C
  159.     147       DEL=S11*S22-S12*S21
  160.     148       FS=1+CABS(DEL)**2-CABS(S11)**2-CABS(S22)**2
  161.     149       FS=FS/(2*CABS(S21)*CABS(S12))
  162.     150       WRITE(*,32) FS
  163.     151       WRITE(6,32) FS
  164.     152 32    FORMAT(1H0,'STABILITY FACTOR K=',F8.3,)
  165.     153 C
  166.     154 C     K >= TO 1 CORRSPONDS TO UNCONDITIONNALLY STABLE
  167.     155 C
  168.     156 C     Y PARAMETERS
  169.     157 C
  170.     158       D=(1+S11)*(1+S22)-S12*S21
  171.     159       Y11=((1+S22)*(1-S11)+S12*S21)/D
  172.     160       Y12=-(2*S12)/D
  173.     161       Y21=-(2*S21)/D
  174.     162       Y22=((1+S11)*(1-S22)+S12*S21)/D
  175.     163       WRITE(*,35) Y11
  176.     164       WRITE(6,35) Y11
  177.     165 35    FORMAT(///,1H0,'Y11=',2(F8.3,3X))
  178.     166       WRITE(*,40) Y12
  179.     167       WRITE(6,40) Y12
  180.     168 40    FORMAT(1H0,'Y12=',2(F8.3,3X))
  181.                                                                          Page   4
  182.                                                                         01-01-80
  183.                                                                         00:05:05
  184. D Line# 1     7                            Microsoft FORTRAN77 V3.31 August 1985
  185.     169       WRITE(*,45) Y21
  186.     170       WRITE(6,45) Y21
  187.     171 45    FORMAT(1H0,'Y21=',2(F8.3,3X))
  188.     172       WRITE(*,50) Y22
  189.     173       WRITE(6,50) Y22
  190.     174 50    FORMAT(1H0,'Y22=',2(F8.3,3X))
  191.     175 C
  192.     176 C     IF THE TRANSISTOR IS POTENTIALLY UNSTABLE, THE SERIES FEEDBACK
  193.     177 C     NETWORK (CONTAINING ONLY ONE SUSCEPTANCE), WHICH MAXIMIZES THE
  194.     178 C     NEGATIVE CONDUCTANCE, IS CALCULATED.
  195.     179 C
  196.     180       IF (FS.GE.1)GOTO 200
  197.     181       WRITE(*,52)
  198.     182       WRITE(6,52)
  199.     183 52    FORMAT(1H0,'ONE SUSCEPTANCE NETWORK CALCULATION')
  200.     184       G11=REAL(Y11)
  201.     185       G22=REAL(Y22)
  202.     186       B11=AIMAG(Y11)
  203.     187       B22=AIMAG(Y22)
  204.     188       YM=REAL(Y12*Y21)
  205.     189       YL=CABS(Y12*Y21)
  206.     190       YI=AIMAG(Y12*Y21)
  207.     191 C
  208.     192 C     2 SOLUTIONS ARE POSSIBLE FOR SERIES FEEDBACK SUSCEPTANCE.
  209.     193 C     FOR EACH OF THE SOLUTIONS, THE NEGATIVE CONDACTANCE IS
  210.     194 C     COMPUTED AND THE LARGEST VALUE WILL BE USED.
  211.     195 C
  212.     196 C     THE COMPUTATION IS DONE FOR BOTH POSSIBLE OUTPUT PORT.
  213.     197 C     THE USER MUST DECIDE, BASED ONT THESE RESULTS, WHICH WILL
  214.     198 C     THE CHOSEN CONFIGURATION.
  215.     199 C
  216.     200       BG1=-B11+(G11*YI)/(YM+YL)
  217.     201       Y1=CMPLX(0.,BG1)
  218.     202       GN1=G22-(YM+YL)/(2*G11)
  219.     203       BG2=-B11+(G11*YI)/(YM-YL)
  220.     204       Y2=CMPLX(0.,BG2)
  221.     205       GN2=G22-(YM-YL)/(2*G11)
  222.     206       WRITE(*,60)GN1,GN2
  223.     207       WRITE(6,60)GN1,GN2
  224.     208 60    FORMAT(1H0,'NEG. CONDUCTANCES :GN1= ',F8.3,'  GN2=',F8.3)
  225.     209       IF(GN1.LE.GN2)GOTO 65
  226.     210       WRITE(*,70)
  227.     211       WRITE(6,70)
  228.     212 70    FORMAT(1H0,'GN2 MORE NEG. THAN GN1, GNMAX=GN2')
  229.     213 80    GNMAX=GN2
  230.     214       BG=BG2
  231.     215       GOTO 85
  232.     216 65    WRITE(*,90)
  233.     217       WRITE(6,90)
  234.     218 90    FORMAT(1H0,'GN1 MORE NEG. THAN  GN2, GNMAX=GN1')
  235.     219 75    GNMAX=GN1
  236.     220       BG=BG1
  237.     221 85    WRITE(*,95)
  238.     222       WRITE(6,95)
  239.     223 95    FORMAT(1H0,'CIRCUIT PARAMETERS FOR AN OUTPUT ON PORT 2')
  240.     224       B0=B22-AIMAG((Y12*Y21)/(Y11+CMPLX(0.,BG)))
  241.                                                                          Page   5
  242.                                                                         01-01-80
  243.                                                                         00:05:05
  244. D Line# 1     7                            Microsoft FORTRAN77 V3.31 August 1985
  245.     225       WRITE(*,100) GNMAX,BG,B0
  246.     226       WRITE(6,100) GNMAX,BG,B0
  247.     227 C
  248.     228 C     COMPUTATION OF THE DISTANCE BETWEEN DR AND TRANSISTOR.
  249.     229 C     AN ALUMINA WITH 50OHMS LINES IS ASSUMED.
  250.     230 C
  251.     231 100   FORMAT(1H0,5X,'GNMAX=',F8.3,'  BG=',F8.3,'  B0=',F8.3)
  252.     232       DIST=ATAN(BG)
  253.     233       IF (DIST.GE.0.0)GOTO 102
  254.     234       DIST=DIST+3.1416
  255.     235 102   DIST=DIST*WL/6.2832
  256.     236       WRITE(*,97)DIST
  257.     237       WRITE(6,97)DIST
  258.     238 97    FORMAT(1H0,'THE DISTANCE BETWEEN DR AND TRANS.=',F6.3,' MM')
  259.     239 C
  260.     240 C     COMPUTATION OF OUTPUT ADMITTANCES FOR REAL DR
  261.     241 C
  262.     242       WRITE(6,620)
  263.     243       WRITE(*,620)
  264.     244 620   FORMAT(/,1H0,'ADMITTANCES WITH REAL DR')
  265.     245       WRITE(*,800)
  266.     246 800   FORMAT(1H0,' ENTER S11 MAX, S11 MIN AND # OF COMPUTATIONS')
  267.     247       READ (*,810) CLIM1,CLIM2,M
  268.     248 810   FORMAT(2F6.3,I2)
  269.     249       DRHO=(CLIM1-CLIM2)/M
  270.     250       M=M+1
  271.     251       DO 820 I=1,M
  272. 1   252       RHO=CLIM2+DRHO*(I-1)
  273. 1   253       RES=(((100*RHO)/(1-RHO))+50.)/50.
  274. 1   254       ZIN=RES+CMPLX(0.,1.)*TAN(DIST*6.2832/WL)
  275. 1   255       ZIN=ZIN/(1+CMPLX(0.,RES)*TAN(DIST*6.2832/WL))
  276. 1   256       XIN=1/ZIN
  277. 1   257       BOUT=B22-AIMAG((Y12*Y21)/(Y11+XIN))
  278. 1   258       GNEG=G22-REAL((Y12*Y21)/(Y11+XIN))
  279. 1   259       WRITE(*,630)RHO,GNEG,BOUT
  280. 1   260       WRITE(6,630)RHO,GNEG,BOUT
  281. 1   261 630   FORMAT(/,' RHO=',F6.3,3X,'GNEG=',F6.3,3X,'BOUT=',F6.3)
  282. 1   262       WRITE(6,615)ZIN,XIN
  283. 1   263 615   FORMAT(' ZIN=',2F6.3,5X,'XIN=',2F6.3)
  284. 1   264 820   CONTINUE
  285.     265 880   WRITE(*,850)
  286.     266 850   FORMAT(1H0,'ENTER S11 AND DIST(MM),(S11<0 ==>FIN)')
  287.     267       READ (*,825) RHO,DIST
  288.     268 825   FORMAT(2F6.3)
  289.     269       IF (RHO.LT.0)GOTO 900
  290.     270       RES=(((100*RHO)/(1-RHO))+50.)/50.
  291.     271       ZIN=RES+CMPLX(0.,1.)*TAN(DIST*6.2832/WL)
  292.     272       ZIN=ZIN/(1+CMPLX(0.,RES)*TAN(DIST*6.2832/WL))
  293.     273       XIN=1/ZIN
  294.     274       BOUT=B22-AIMAG((Y12*Y21)/(Y11+XIN))
  295.     275       GNEG=G22-REAL((Y12*Y21)/(Y11+XIN))
  296.     276       WRITE(*,630)RHO,GNEG,BOUT
  297.     277       WRITE(6,630)RHO,GNEG,BOUT
  298.     278       WRITE(6,890)ZIN,XIN,DIST
  299.     279 890   FORMAT(' ZIN=',2F6.3,5X,'XIN=',2F6.3,3X,'DIST=',F6.3)
  300.     280       GOTO 880
  301.                                                                          Page   6
  302.                                                                         01-01-80
  303.                                                                         00:05:05
  304. D Line# 1     7                            Microsoft FORTRAN77 V3.31 August 1985
  305.     281 900   BG3=-B22+(G22*YI)/(YM+YL)
  306.     282       Y3=CMPLX(0.,BG3)
  307.     283       GN3=G11-(YM+YL)/(2*G22)
  308.     284       BG4=-B22+(G22*YI)/(YM-YL)
  309.     285       Y4=CMPLX(0.,BG4)
  310.     286       GN4=G11-(YM-YL)/(2*G22)
  311.     287       WRITE(*,105)GN3,GN4
  312.     288       WRITE(6,105)GN3,GN4
  313.     289 105   FORMAT(///,1H0,'NEGATIVE CONDUCTANCES: GN3=',F8.3,5X,'GN4=',F8.3)
  314.     290       IF(GN3.LE.GN4)GOTO 120
  315.     291       WRITE(*,110)
  316.     292       WRITE(6,110)
  317.     293 110   FORMAT(1H0,'GN4 MORE NEG. THAN GN3, GNMAX=GN4')
  318.     294 115   GNMAXA=GN4
  319.     295       BGA=BG4
  320.     296       GOTO 135
  321.     297 120   WRITE(*,125)
  322.     298       WRITE(6,125)
  323.     299 125   FORMAT(1H0,'GN3 MORE NEG. THAN GN4, GNMAX=GN3')
  324.     300 130   GNMAXA=GN3
  325.     301       BGA=BG3
  326.     302 135   WRITE(*,140)
  327.     303       WRITE(6,140)
  328.     304 140   FORMAT(1H0,'CIRCUIT PARAMETERS FOR AN OUTPUT ON PORT 1')
  329.     305       B0A=B11-AIMAG((Y12*Y21)/(Y22+CMPLX(0.,BGA)))
  330.     306       WRITE(*,100) GNMAXA,BGA,B0A
  331.     307       WRITE(6,100) GNMAXA,BGA,B0A
  332.     308       DIST=ATAN(BGA)
  333.     309       IF (DIST.GE.0.0) GOTO 103
  334.     310       DIST=DIST+3.1416
  335.     311 103   DIST=DIST*WL/6.2832
  336.     312       WRITE(*,142)DIST
  337.     313       WRITE(6,142)DIST
  338.     314 142   FORMAT(1H0,'THE DISTANCE BETWEEN DR AND TRANS.=',F6.3,' MM')
  339.     315       WRITE(6,620)
  340.     316       WRITE(*,620)
  341.     317       WRITE(*,800)
  342.     318       READ (*,810) CLIM1,CLIM2,M
  343.     319       DRHO=(CLIM1-CLIM2)/M
  344.     320       M=M+1
  345.     321       DO 920 I=1,M
  346. 1   322       RHO=CLIM2+DRHO*(I-1)
  347. 1   323       RES=(((100*RHO)/(1-RHO))+50.)/50.
  348. 1   324       ZIN=RES+CMPLX(0.,1.)*TAN(DIST*6.2832/WL)
  349. 1   325       ZIN=ZIN/(1+CMPLX(0.,RES)*TAN(DIST*6.2832/WL))
  350. 1   326       XIN=1/ZIN
  351. 1   327       BOUT=B11-AIMAG((Y12*Y21)/(Y22+XIN))
  352. 1   328       GNEG=G11-REAL((Y12*Y21)/(Y22+XIN))
  353. 1   329       WRITE(*,630)RHO,GNEG,BOUT
  354. 1   330       WRITE(6,630)RHO,GNEG,BOUT
  355. 1   331       WRITE(6,615)ZIN,XIN
  356. 1   332 920   CONTINUE
  357.     333 1010  WRITE(*,850)
  358.     334       READ (*,825) RHO,DIST
  359.     335       IF(RHO.LT.0)GOTO 200
  360.     336       RES=(((100*RHO)/(1-RHO))+50.)/50.
  361.                                                                          Page   7
  362.                                                                         01-01-80
  363.                                                                         00:05:05
  364. D Line# 1     7                            Microsoft FORTRAN77 V3.31 August 1985
  365.     337       ZIN=RES+CMPLX(0.,1.)*TAN(DIST*6.2832/WL)
  366.     338       ZIN=ZIN/(1+CMPLX(0.,RES)*TAN(DIST*6.2832/WL))
  367.     339       XIN=1/ZIN
  368.     340       BOUT=B11-AIMAG((Y12*Y21)/(Y22+XIN))
  369.     341       GNEG=G11-REAL((Y12*Y21)/(Y22+XIN))
  370.     342       WRITE(*,630)RHO,GNEG,BOUT
  371.     343       WRITE(6,630)RHO,GNEG,BOUT
  372.     344       WRITE(6,890)ZIN,XIN,DIST
  373.     345       GOTO 1010
  374.     346 200   WRITE(*,150)
  375.     347       IF(FS.LT.1) GOTO 500
  376.     348       WRITE(6,150)
  377.     349 150   FORMAT(///,1H0,'THREE ADMITTANCES NETWORK CALCULATION')
  378.     350       IF(ICONF-2)205,215,225
  379.     351 C
  380.     352 C     THE ANALYSIS IS BASED ON SMALL-SIGNAL S-PARAMETERS. THE NETWORK
  381.     353 C     CONSISTING OF THREE ADMITTANCES WHICH LEADS TO A MAXIMUM
  382.     354 C     LOADING OSCILLATOR IS COMPUTED. EACH OF THE THREE POSSIBLE OUTPUT
  383.     355 C     PORTS IS CONSIDERED.
  384.     356 C
  385.     357 205   XM(1,1)=Y11
  386.     358       XM(3,1)=Y21
  387.     359       XM(2,1)=-(Y11+Y21)
  388.     360       XM(1,3)=Y12
  389.     361       XM(1,2)=-(Y12+Y11)
  390.     362       XM(3,3)=Y22
  391.     363       XM(2,3)=-(Y22+Y12)
  392.     364       XM(3,2)=-(Y22+Y21)
  393.     365       XM(2,2)=-(XM(3,2)+XM(1,2))
  394.     366       WRITE(*,210)
  395.     367       WRITE(6,210)
  396.     368 210   FORMAT(1H0,15X,'Y(I,J) MATRIX FROM CB OR CG')
  397.     369       GOTO 235
  398.     370 215   XM(2,2)=Y11
  399.     371       XM(2,3)=Y12
  400.     372       XM(3,2)=Y21
  401.     373       XM(3,3)=Y22
  402.     374       XM(2,1)=-(Y11+Y12)
  403.     375       XM(3,1)=-(Y21+Y22)
  404.     376       XM(1,2)=-(Y11+Y21)
  405.     377       XM(1,3)=-(Y12+Y22)
  406.     378       XM(1,1)=-(XM(1,2)+XM(1,3))
  407.     379       WRITE(*,220)
  408.     380       WRITE(6,220)
  409.     381 220   FORMAT(1H0,15X,'Y(I,J) MATRIX FROM CE OR CS')
  410.     382       GOTO 235
  411.     383 225   XM(1,1)=Y22
  412.     384       XM(1,2)=Y21
  413.     385       XM(2,1)=Y12
  414.     386       XM(2,2)=Y11
  415.     387       XM(1,3)=-(Y22+Y21)
  416.     388       XM(2,3)=-(Y12+Y11)
  417.     389       XM(3,1)=-(Y22+Y12)
  418.     390       XM(3,2)=-(Y21+Y11)
  419.     391       XM(3,3)=-(XM(3,1)+XM(3,2))
  420.     392       WRITE(*,230)
  421.                                                                          Page   8
  422.                                                                         01-01-80
  423.                                                                         00:05:05
  424. D Line# 1     7                            Microsoft FORTRAN77 V3.31 August 1985
  425.     393       WRITE(6,230)
  426.     394 230   FORMAT(1H0,15X,'Y(I,J) MATRIX FROM CC OR CD')
  427.     395 235   DO 240 I=1,3
  428. 1   396       WRITE(*,245)(XM(I,J),J=1,3)
  429. 1   397       WRITE(6,245)(XM(I,J),J=1,3)
  430. 1   398 245   FORMAT(1H0,3(4X,2(F8.3)))
  431. 1   399 240   CONTINUE
  432.     400       GM=(1/3)*(REAL(XM(1,1))+REAL(XM(2,2))+REAL(XM(3,3)))
  433.     401       IF (CABS(GM).GE.0) GOTO 250
  434.     402       WRITE(*,255)
  435.     403       WRITE(6,255)
  436.     404 255   FORMAT(1H0,2X,'AVERAGE AUTO-CONDUCTANCE <0, INACTIVE DEVICE')
  437.     405       GOTO 500
  438.     406 250   YE=0.5*(XM(2,2)+XM(3,3)-XM(1,1))
  439.     407       YB=0.5*(XM(3,3)+XM(1,1)-XM(2,2))
  440.     408       YC=0.5*(XM(1,1)+XM(2,2)-XM(3,3))
  441.     409       Y0=0.5*(XM(1,2)-XM(2,1))
  442.     410       T=AIMAG(Y0)**2
  443.     411       T=T-(REAL(YE)*REAL(YB)+REAL(YB)*REAL(YC)+REAL(YC)*REAL(YE))
  444.     412 C
  445.     413 C     BASE NETWORK
  446.     414 C
  447.     415       XF=REAL(Y0)/AIMAG(Y0)
  448.     416       WRITE(*,265)
  449.     417       WRITE(6,265)
  450.     418 265   FORMAT(//,1H0,'BASE SUSCEPTANCES OF THE NETWORK')
  451.     419       BBC=-XF*REAL(YE)-AIMAG(YE)
  452.     420       BEC=-XF*REAL(YB)-AIMAG(YB)
  453.     421       BEB=-XF*REAL(YC)-AIMAG(YC)
  454.     422       WRITE(*,270)BBC
  455.     423       WRITE(6,270)BBC
  456.     424 270   FORMAT(//,1H0,'BBC=',2X,F8.3)
  457.     425       WRITE(*,280)BEC
  458.     426       WRITE(6,280)BEC
  459.     427 280   FORMAT(1H0,'BEC=',2X,F8.3)
  460.     428       WRITE(*,290)BEB
  461.     429       WRITE(6,290)BEB
  462.     430 290   FORMAT(1H0,'BEB=',2X,F8.3)
  463.     431 C
  464.     432 C     COMPUTATION OF SUSCEPTANCE TO BE ADDED IN
  465.     433 C     PARALL WITH MAXIMUM LOAD
  466.     434 C
  467.     435       WRITE(6,260)
  468.     436