C ====================================================================== SUBROUTINE DSOUT C ---------------------------------------------------------------------- C Title:Trex C C *** FORTRAN 77 *** C C Associated files: C TF.dat C TF.dsd.f C C ====================================================================== C Include and Common Block Declarations C ---------------------------------------------------------------------- INCLUDE 'emt.d' INCLUDE 'emt.e' COMMON /S1/TIME,DELT,ICH,PRINT,FINTIM COMMON /S2/STOR(ND10),NEXC COMMON /S4/VAR(ND11),CON(ND12),PGB(ND13) REAL PI_, TWO_PI, PI_BY3, PI2_BY3 DATA PI_/3.141592654/, TWO_PI/6.283185307/ DATA PI_BY3/1.047197551/, PI2_BY3/2.094395102/ C C ====================================================================== C Function Declarations C ---------------------------------------------------------------------- C REALP2 : Real Pole C REAL REALP2 C ====================================================================== C Variable Declarations C ---------------------------------------------------------------------- REAL TIME,DELT,PRINT,FINTIM,STOR,VAR,CON,PGB INTEGER ICH,NEXC REAL Ptf, ILS, Vload, Id, Vd2, Vd3, Vd4, Vd1 REAL RT_11, RT_12, RT_13, RT_14, RT_15, RT_16 REAL RT_17, RT_18, Id1, Id2, Id3, Id4, Iload REAL Vrouter, Vrinner, RT_19, RT_20, RT_21 REAL RT_22, RT_23, RT_24, IDS C C ====================================================================== C CON() definitions C ---------------------------------------------------------------------- C CON( 1) = Reference of "RT_29" before definition. C CON( 2) = Reference of "Ptf" before definition. C CON( 3) = Reference of "RT_31" before definition. C CON( 4) = Reference of "ILS" before definition. C CON( 5) = Reference of "IDS" before definition. C C ====================================================================== C VAR() definitions C ---------------------------------------------------------------------- C VAR(1) Slider Desc="Pmotor" Group="" Min=-1.12e+07 Max=1.12e+07 InitValue=0.0 Units="MW" C VAR(2) Slider Desc="Comm Inductance" Group="GenCtrl" Min=1.0e-05 Max=5.0e-05 InitValue=2.26338e-05 Units="H" C VAR(3) Switch Desc="Control Mode" Group="GenCtrl" InitValue=0 P0=0.0 P1=1.0 OnText="ON_V" OffText="OFF_I" C VAR(4) Slider Desc="min Gamma" Group="GenCtrl" Min=0.0 Max=90.0 InitValue=15.0 Units="deg" C VAR(5) Slider Desc="Intgrl Gain" Group="GenCtrl" Min=0.0 Max=1000.0 InitValue=0.0 Units="farads" C VAR(6) Slider Desc="Prop Gain" Group="GenCtrl" Min=0.0 Max=100.0 InitValue=0.0 Units="ohms" C VAR(7) Switch Desc="PS12" Group="" InitValue=1 P0=0.0 P1=1.0 OnText="ON" OffText="OFF" C VAR(8) Switch Desc="PS34" Group="" InitValue=1 P0=0.0 P1=1.0 OnText="ON" OffText="OFF" C VAR(9) Slider Desc="GI_PLO" Group="GenCtrl" Min=0.0 Max=200.0 InitValue=100.0 Units="" C C ====================================================================== C PGB() definitions C ---------------------------------------------------------------------- C PGB(1) Output Desc="Pnet" Group="" Max=1.0e+09 Min=-1.0e+09 Units="Watt" C PGB(2) Output Desc="MG Freq" Group="" Max=100.0 Min=0.0 Units="Hz" C PGB(3) Output Desc="Wnet" Group="" Max=5.0e+09 Min=-5.0e+09 Units="Joule" C PGB(4) Output Desc="Tinner" Group="" Max=200.0 Min=0.0 Units="Deg C" C PGB(5) Output Desc="Ratinner" Group="" Max=10.0 Min=1.0 Units="p.u." C PGB(6) Output Desc="Action" Group="" Max=1.0e+10 Min=0.0 Units="A^2-sec" C PGB(7) Output Desc="Touter" Group="" Max=200.0 Min=0.0 Units="Deg C" C PGB(8) Output Desc="Ratouter" Group="" Max=10.0 Min=1.0 Units="p.u." C PGB(9) Output Desc="Reference" Group="" Max=30.0 Min=-30.0 Units="kA" C PGB(10) Output Desc="KB" Group="" Max=2.0 Min=-2.0 Units="" C PGB(11) Output Desc="Comp" Group="" Max=2.0 Min=-2.0 Units="" C PGB(12) Output Desc="Q" Group="" Max=2.0 Min=-2.0 Units="" C PGB(13) Output Desc="Ierror" Group="" Max=30.0 Min=0.0 Units="kA" C PGB(14) Output Desc="Vrequest" Group="" Max=10.0 Min=-10.0 Units="kV" C PGB(15) Output Desc="Areq" Group="" Max=180.0 Min=0.0 Units="deg" C PGB(16) Output Desc="Alpha Command" Group="" Max=180.0 Min=0.0 Units="deg" C PGB(17) Output Desc="Amax" Group="" Max=180.0 Min=0.0 Units="deg" C PGB(18) Output Desc="Measured Alpha" Group="" Max=120.0 Min=0.0 Units="deg" C PGB(19) Output Desc="Measured Gamma" Group="" Max=180.0 Min=0.0 Units="deg" C PGB(20) Output Desc="Firing Pulse (3,1)" Group="" Max=2.0 Min=-2.0 Units="" C PGB(21) Output Desc="Vthy(3,1)" Group="" Max=2.0 Min=-2.0 Units="kV" C PGB(22) Output Desc="Ithy(3,1)" Group="" Max=2.0 Min=-2.0 Units="kA" C PGB(23) Output Desc="Firing Pulse (4,1)" Group="" Max=2.0 Min=-2.0 Units="" C PGB(24) Output Desc="Vthy(4,1)" Group="" Max=2.0 Min=-2.0 Units="kV" C PGB(25) Output Desc="Ithy(4,1)" Group="" Max=2.0 Min=-2.0 Units="kA" C PGB(26) Output Desc="Firing Pulse (2,1)" Group="" Max=2.0 Min=-2.0 Units="" C PGB(27) Output Desc="Vthy(2,1)" Group="" Max=2.0 Min=-2.0 Units="kV" C PGB(28) Output Desc="Ithy(2,1)" Group="" Max=2.0 Min=-2.0 Units="kA" C PGB(29) Output Desc="Ithy(1,1)" Group="" Max=2.0 Min=-2.0 Units="kA" C PGB(30) Output Desc="Ithy(1,2)" Group="" Max=2.0 Min=-2.0 Units="kA" C PGB(31) Output Desc="Ithy(1,3)" Group="" Max=2.0 Min=-2.0 Units="kA" C PGB(32) Output Desc="Ithy(1,4)" Group="" Max=2.0 Min=-2.0 Units="kA" C PGB(33) Output Desc="Ithy(1,5)" Group="" Max=2.0 Min=-2.0 Units="kA" C PGB(34) Output Desc="Ithy(1,6)" Group="" Max=2.0 Min=-2.0 Units="kA" C PGB(35) Output Desc="Firing Pulse (1,1)" Group="" Max=2.0 Min=-2.0 Units="" C PGB(36) Output Desc="Vthy(1,1)" Group="" Max=2.0 Min=-2.0 Units="kV" C PGB(37) Output Desc="Load Voltage" Group="Main/DC" Max=5.0 Min=-5.0 Units="kV" C PGB(38) Output Desc="Load Voltage (smooth)" Group="Main/DC" Max=5.0 Min=-5.0 Units="kV" C PGB(39) Output Desc="Iload" Group="Main/DC" Max=2.0 Min=-2.0 Units="kA" C PGB(40) Output Desc="Id" Group="Main/DC" Max=2.0 Min=-2.0 Units="kA" C PGB(41) Output Desc="Id1" Group="" Max=20.0 Min=0.0 Units="kA" C PGB(42) Output Desc="Id2" Group="" Max=20.0 Min=0.0 Units="kA" C PGB(43) Output Desc="Id3" Group="" Max=20.0 Min=0.0 Units="kA" C PGB(44) Output Desc="Id4" Group="" Max=20.0 Min=0.0 Units="kA" C PGB(45) Output Desc="OH Power" Group="Main/DC" Max=2.0e+08 Min=-2.0e+08 Units="Watt" C PGB(46) Output Desc="Vd1" Group="" Max=2.0 Min=-2.0 Units="kV" C PGB(47) Output Desc="Vd2" Group="" Max=2.0 Min=-2.0 Units="kV" C PGB(48) Output Desc="Vd3" Group="" Max=2.0 Min=-2.0 Units="kV" C PGB(49) Output Desc="Vd4" Group="" Max=2.0 Min=-2.0 Units="kV" C C ====================================================================== C C ====================================================================== C Subsystem: "MG" C ---------------------------------------------------------------------- C ====================================================================== C Subsystem: "TF Inner Leg R(I^2T)" C ---------------------------------------------------------------------- C ====================================================================== C Subsystem: "TF Outer Leg R(I^2T)" C ---------------------------------------------------------------------- C ====================================================================== C Subsystem: "Mode Control" C ---------------------------------------------------------------------- C ====================================================================== C Subsystem: "TF Power Supply" C ---------------------------------------------------------------------- Id = CDC(4,7,1) Id2 = CDCNG(7,1) Id3 = CDCNG(8,1) Id4 = CDCNG(9,1) Iload = CDCNG(10,1) Id1 = CDCNG(11,1) Vload = VDC(4,1)-VDC(31,1) Vd2 = VDC(29,1)-VDC(31,1) Vd3 = VDC(30,1)-VDC(31,1) Vd4 = VDC(43,1)-VDC(31,1) Vd1 = VDC(28,1)-VDC(31,1) Vrouter = VDC(4,1)-VDC(6,1) Vrinner = VDC(5,1)-VDC(31,1) C Gain RT_11 = 0.10 * VLOAD C C Real Pole RT_12 = REALP2( 1.0 , 0.0010, RT_11) C PGB(37) = RT_11 PGB(38) = RT_12 C Gain RT_13 = 10.0 * ILOAD C PGB(39) = RT_13 ILS = RT_13 C Gain RT_14 = 10.0 * ID1 C PGB(40) = RT_14 IDS = RT_14 C Gain RT_15 = 10.0 * ID2 C C Gain RT_16 = 10.0 * ID1 C C Gain RT_17 = 10.0 * ID4 C C Gain RT_18 = 10.0 * ID3 C PGB(41) = RT_16 PGB(42) = RT_15 PGB(43) = RT_18 PGB(44) = RT_17 C RT_19 = VLOAD * ILOAD C C Gain RT_20 = 1.0e+06 * RT_19 C Ptf = RT_20 PGB(45) = RT_20 C Gain RT_21 = 0.10 * VD2 C C Gain RT_22 = 0.10 * VD1 C C Gain RT_23 = 0.10 * VD4 C C Gain RT_24 = 0.10 * VD3 C PGB(46) = RT_22 PGB(47) = RT_21 PGB(48) = RT_24 PGB(49) = RT_23 CON(2) = Ptf CON(4) = ILS CON(5) = IDS C ====================================================================== RETURN END