C ======================================================================
      SUBROUTINE DSDYN
C ----------------------------------------------------------------------
C   Title:Trex
C
C   *** FORTRAN 77 ***
C
C   Associated files:
C      TF.dat
C      TF.dso.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     INTGL3   : Integrator
C     APLIN    : nonlinear transfer characteristic
C     PI3      : Proportional/Integral Controller
C     LIMIT    : Hard Limiter
C     RAMP3    : Rate Limiter
C     TSAT21   : Transformer Saturation Subroutine
C     G6P200   : Valve Group Model
C     ESYS65   : 3-Phase RRL Source model w/ Zero sequence
C
      REAL    INTGL3
      REAL    APLIN
      REAL    PI3
      REAL    LIMIT
      REAL    RAMP3
C ======================================================================
C   Variable Declarations
C ----------------------------------------------------------------------
      REAL    TIME,DELT,PRINT,FINTIM,STOR,VAR,CON,PGB
      INTEGER ICH,NEXC
      REAL    AR, KB12, Freq, Vrms, IROUTER, IRINNER
      REAL    KB34, Ptf, ILS, RT_1, RT_2, VT2(6), IT2(6)
      REAL    GI, RT_3, RT_4, VT3(6), IT3(6), RT_5, AM
      REAL    RT_6, RT_7, VT4(6), IT4(6), RT_8, RT_9
      REAL    RT_10, VT1(6), IT1(6), RT_25, RT_26, RT_27
      REAL    RT_28, RT_29, RT_30, RT_31, RT_32, RT_33
      REAL    RT_34, RT_35, RT_36, RT_37, RT_38, RT_39
      REAL    RT_40, RT_41, RT_42, RT_43, RT_44, RT_45
      REAL    RT_46, RT_47, RT_48, RT_49, RT_50, RT_51
      REAL    RT_52, RT_53, RT_54, RT_55, RT_56, RT_57
      REAL    RT_58, RT_59, RT_60, RT_61, RT_62, RT_63
      REAL    RT_64, RT_65, RT_66, RT_67, RT_68, RT_69
      REAL    RT_70, RT_71, RT_72, RT_73, RT_74, RT_75
      REAL    RT_76, RT_77, RT_78, RT_79, RT_80, RT_81
      REAL    RT_82, RT_83, RT_84, RT_85, RT_86, RT_87
      REAL    RT_88, RT_89, RT_90, RT_91, RT_92, RT_93
      REAL    RT_94, RT_95, RT_96, RT_97, RT_98, RT_99
      REAL    RT_100, RT_101, RT_102, RT_103, RT_104
      REAL    RT_105, RT_106, RT_107, RT_108, RT_109
      REAL    RT_110, RT_111, RT_112, RT_113, RT_114
      REAL    RT_115, RT_116, RT_117, RT_118, RT_119
      REAL    RT_120, RT_121, RT_122, RT_123, RT_124
      REAL    RT_125, RT_126, RT_127, RT_128, RT_129
      REAL    RT_130, RT_131, RT_132, RT_133, RT_134
      REAL    RT_135, RT_136, RT_137, X_C(10), Y_C(10)
      REAL    RT_138, RT_139, RT_140, PGAIN, ITAU, RT_141
      REAL    RT_142, RT_143, VDO, RT_144, RT_145, RT_146
      REAL    RT_147, RT_148, RT_149, RT_150, RT_151
      REAL    RT_152, RT_153, RT_154, NPP, RT_155, RT_156
      REAL    IDS, X_KNEE, C_INST, P_MES, Q_MES, E_RMS
      REAL    THETA_D, A_CURR, B_CURR, C_CURR, A_MAG
      REAL    B_MAG, C_MAG, SRC_V, SRC_F, SRC_P
      INTEGER FP2(6), FP3(6), FP4(6), FP1(6), IT_3, IT_4
      INTEGER IT_5, IT_6, IT_7, IT_8, IT_9, IT_10
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   Subsystem: "MG"
C ----------------------------------------------------------------------
      RT_26 = 0.0
      RT_27 = VAR(1)
      RT_30 = 76.360
      RT_36 = 1.576712e+04
      RT_38 = 1.20
      RT_40 = 2.376e+03
      RT_42 = 3.0
      RT_44 = 3.3410
      RT_47 = 2.718e+09
      RT_48 = 5.546995e+05
      RT_49 = 13.80
      Vrms = RT_49
      RT_29 = CON(1)
C
      RT_28 = RT_27 * RT_29
C
      Ptf = CON(2)
      RT_31 = CON(3)
      RT_25 = - RT_26 - Ptf + RT_28 - RT_31
C
      PGB(1) = RT_25
C  Integrator
      RT_45 = INTGL3( 2.718e+09, 0.0, 4.25e+09, 1.0, RT_25)
C
C
      IF (ABS(RT_48) .LT. 1.0E-10) THEN
         IF (RT_48 .LT. 0.0)  THEN
            RT_33 = -1.0E10 * RT_45
         ELSE
            RT_33 =  1.0E10 * RT_45
         ENDIF
      ELSE
         RT_33 = RT_45 / RT_48
      ENDIF
C
C
      IF (RT_33 .LE. 0.0) THEN
         RT_34 = 0.0
      ELSE
         RT_34 = SQRT(RT_33)
      ENDIF
C
C
      IF (ABS(RT_30) .LT. 1.0E-10) THEN
         IF (RT_30 .LT. 0.0)  THEN
            RT_29 = -1.0E10 * RT_34
         ELSE
            RT_29 =  1.0E10 * RT_34
         ENDIF
      ELSE
         RT_29 = RT_34 / RT_30
      ENDIF
C
C
      RT_35 = RT_34 * RT_36
C
C
      RT_37 = RT_34 ** RT_38
C
C
      RT_39 = RT_37 * RT_40
C
C
      RT_41 = RT_34 ** RT_42
C
C
      RT_43 = RT_41 * RT_44
C
      RT_32 = + RT_39 + RT_43
C
      RT_31 = + RT_35 + RT_32
C
      Freq = RT_34
      PGB(2) = RT_34
      RT_46 = + RT_47 - RT_45
C
      PGB(3) = RT_46
      CON(1) = RT_29
      CON(3) = RT_31
C ======================================================================
C   Subsystem: "TF Inner Leg R(I^2T)"
C ----------------------------------------------------------------------
      RT_51 = 4.5496e-15
      RT_54 = 0.136080
      RT_57 = 2.0934e+14
      RT_60 = 2.8709e+11
      RT_61 = 10.0
      RT_62 = 6.2949e+13
      RT_64 = 5.3309e-32
      RT_66 = 6.742e-04
      RT_69 = 0.00410
      RT_72 = 1.0
C
      RT_56 = RT_57 * RT_61
C
C
      RT_58 = RT_61 * RT_61
C
C
      RT_59 = RT_58 * RT_60
C
      RT_50 = + RT_62 + RT_56 - RT_59
C
      ILS = CON(4)
C
      IF (ABS(RT_66) .LT. 1.0E-10) THEN
         IF (RT_66 .LT. 0.0)  THEN
            RT_77 = -1.0E10 * ILS
         ELSE
            RT_77 =  1.0E10 * ILS
         ENDIF
      ELSE
         RT_77 = ILS / RT_66
      ENDIF
C
C
      RT_75 = RT_77 * RT_77
C
C  Gain
      RT_76 = 1000.0 * RT_75
C
C  Integrator
      RT_78 = INTGL3( 0.0, 0.0, 1.0e+20, 0.0010, RT_76)
C
      RT_63 = + RT_78 + RT_50
C
C
      RT_52 = RT_63 * RT_63
C
C
      RT_53 = RT_52 * RT_64
C
C
      RT_74 = RT_51 * RT_63
C
      RT_55 = + RT_54 + RT_74 + RT_53
C
      RT_65 = AMAX1(RT_55,RT_61)
C
      PGB(4) = RT_65
      RT_67 = + RT_65 - RT_61
C
C
      RT_68 = RT_67 * RT_69
C
C
      RT_70 = ILS * RT_68
C
      RT_71 = + RT_68 + RT_72
C
      PGB(5) = RT_71
C  Gain
      RT_73 = 0.10 * RT_70
C
      IRINNER = RT_73
C
      RT_79 = ILS * ILS
C
C  Gain
      RT_80 = 1000.0 * RT_79
C
C  Integrator
      RT_81 = INTGL3( 0.0, 0.0, 1.0e+20, 0.0010, RT_80)
C
      PGB(6) = RT_81
C ======================================================================
C   Subsystem: "TF Outer Leg R(I^2T)"
C ----------------------------------------------------------------------
      RT_83 = 4.5496e-15
      RT_86 = 0.136080
      RT_89 = 2.0934e+14
      RT_92 = 2.8709e+11
      RT_93 = 10.0
      RT_94 = 6.2949e+13
      RT_96 = 5.3309e-32
      RT_98 = 0.0037356
      RT_101 = 0.00410
      RT_104 = 1.0
C
      RT_88 = RT_89 * RT_93
C
C
      RT_90 = RT_93 * RT_93
C
C
      RT_91 = RT_90 * RT_92
C
      RT_82 = + RT_94 + RT_88 - RT_91
C
C
      IF (ABS(RT_98) .LT. 1.0E-10) THEN
         IF (RT_98 .LT. 0.0)  THEN
            RT_109 = -1.0E10 * ILS
         ELSE
            RT_109 =  1.0E10 * ILS
         ENDIF
      ELSE
         RT_109 = ILS / RT_98
      ENDIF
C
C
      RT_107 = RT_109 * RT_109
C
C  Gain
      RT_108 = 1000.0 * RT_107
C
C  Integrator
      RT_110 = INTGL3( 0.0, 0.0, 1.0e+20, 0.0010, RT_108)
C
      RT_95 = + RT_110 + RT_82
C
C
      RT_84 = RT_95 * RT_95
C
C
      RT_85 = RT_84 * RT_96
C
C
      RT_106 = RT_83 * RT_95
C
      RT_87 = + RT_86 + RT_106 + RT_85
C
      RT_97 = AMAX1(RT_87,RT_93)
C
      PGB(7) = RT_97
      RT_99 = + RT_97 - RT_93
C
C
      RT_100 = RT_99 * RT_101
C
C
      RT_102 = ILS * RT_100
C
      RT_103 = + RT_100 + RT_104
C
      PGB(8) = RT_103
C  Gain
      RT_105 = 0.10 * RT_102
C
      IROUTER = RT_105
C ======================================================================
C   Subsystem: "Mode Control"
C ----------------------------------------------------------------------
      RT_124 = 13.80
      RT_125 = 0.750
      RT_128 = VAR(2)
      RT_129 = VAR(3)
      RT_131 = VAR(4)
      RT_132 = TIME
      RT_133 = TIME
      IT_3 = 0
      IT_4 = 1
      IT_5 = -7
      IT_6 = 0
      IT_10 = 33
      RT_144 = VAR(5)
      RT_145 = 1.0
      PGAIN = VAR(6)
      RT_147 = 1.0
      RT_149 = VAR(7)
      RT_150 = TIME
      RT_154 = VAR(8)
      RT_155 = 2.0
C  Cosine Function
      RT_119 = COS(RT_131 * 0.0174533)
C
C
      IF (ABS(RT_125) .LT. 1.0E-10) THEN
         IF (RT_125 .LT. 0.0)  THEN
            RT_123 = -1.0E10 * RT_124
         ELSE
            RT_123 =  1.0E10 * RT_124
         ENDIF
      ELSE
         RT_123 = RT_124 / RT_125
      ENDIF
C
C
      IF (ABS(RT_123) .LT. 1.0E-10) THEN
         IF (RT_123 .LT. 0.0)  THEN
            RT_122 = -1.0E10 * Vrms
         ELSE
            RT_122 =  1.0E10 * Vrms
         ENDIF
      ELSE
         RT_122 = Vrms / RT_123
      ENDIF
C
C  Gain
      RT_127 = 6.28320 * Freq
C
C
      RT_126 = RT_128 * RT_127
C
C  Gain
      RT_130 = 1.41420 * RT_126
C
C
      IF (ABS(RT_122) .LT. 1.0E-10) THEN
         IF (RT_122 .LT. 0.0)  THEN
            RT_121 = -1.0E10 * RT_130
         ELSE
            RT_121 =  1.0E10 * RT_130
         ENDIF
      ELSE
         RT_121 = RT_130 / RT_122
      ENDIF
C
C  Step_Xfer_Function
      IF (RT_133 .LT. 0.050) THEN
         RT_134 = 0.0
      ELSE
         RT_134 = 1.0
      ENDIF
C
      X_C(1) = 0.0
      X_C(2) = 0.050
      X_C(3) = 1.250
      X_C(4) = 1.2510
      X_C(5) = 2.0
      X_C(6) = 3.0
      X_C(7) = 0.2110
      X_C(8) = 0.60
      X_C(9) = 1.10
      X_C(10) = 1.50
      Y_C(1) = 0.0
      Y_C(2) = 71.1720
      Y_C(3) = 71.1720
      Y_C(4) = 0.0
      Y_C(5) = 0.0
      Y_C(6) = 0.0
      Y_C(7) = 19.750
      Y_C(8) = 0.0
      Y_C(9) = 22.60
      Y_C(10) = 0.0
C Transfer Function
      RT_137 = APLIN(6,X_C,Y_C,RT_132)
C
      RT_111 = + RT_137 - ILS
C
      PGB(9) = RT_137
      IDS = CON(5)
C  COMPARATOR
      IF ( IDS .GT. FLOAT(IT_10)) THEN
         IT_9 = 1.0
      ELSE
         IT_9 = 0.0
      ENDIF
C
C  RSFlipFlop
       IF(TIME.LT.DELT) THEN
          STOR(NEXC+1)=0
       ELSE
          IF ((IT_9.NE.0).AND.(IT_6.EQ.0)) THEN
             STOR(NEXC+1)=1
          ELSE IF ((IT_6.NE.0).AND.(IT_9.EQ.0)) THEN
             STOR(NEXC+1)=0
          END IF
       END IF
       IT_7=STOR(NEXC+1)
       IT_8=1-IT_7
       NEXC=NEXC+1
C
C  Selector
      IF (IT_7 .EQ. 1) THEN
         RT_136 = FLOAT(IT_5)
      ELSE
         RT_136 = FLOAT(IT_4)
      ENDIF
C
C  Selector
      IF (NINT(RT_134) .EQ. 1) THEN
         RT_135 = RT_136
      ELSE
         RT_135 = FLOAT(IT_3)
      ENDIF
C
      PGB(10) = RT_135
      PGB(11) = FLOAT(IT_9)
      PGB(12) = FLOAT(IT_7)
      PGB(13) = RT_111
C  Gain
      VDO = 1.350 * RT_122
C
C
      IF (ABS(RT_144) .LT. 1.0E-10) THEN
         IF (RT_144 .LT. 0.0)  THEN
            ITAU = -1.0E10 * RT_145
         ELSE
            ITAU =  1.0E10 * RT_145
         ENDIF
      ELSE
         ITAU = RT_145 / RT_144
      ENDIF
C
C  PI_Controller
      RT_140 = PI3( PGAIN , ITAU , -1.10 , 1.10, 1.10 , RT_111)
C
C  Selector
      IF (NINT(RT_129) .EQ. 1) THEN
         RT_112 = RT_137
      ELSE
         RT_112 = RT_140
      ENDIF
C
      PGB(14) = RT_140
C
      RT_146 = RT_147 * VDO
C
C
      IF (ABS(RT_146) .LT. 1.0E-10) THEN
         IF (RT_146 .LT. 0.0)  THEN
            RT_113 = -1.0E10 * RT_112
         ELSE
            RT_113 =  1.0E10 * RT_112
         ENDIF
      ELSE
         RT_113 = RT_112 / RT_146
      ENDIF
C
C  Hard Limit
      RT_115 = LIMIT(-1.0, 1.0, RT_113)
C
C  ArcCosine Function
      RT_114 = ACOS(RT_115)
C
C  Gain
      RT_141 = 57.30 * RT_114
C
      PGB(15) = RT_141
C  Band_Compare
      IF (RT_150 .LT. 0.050) THEN
         RT_151 = 0.0
      ELSE IF (RT_150 .LE. 1.250) THEN
         RT_151 = 1.0
      ELSE
         RT_151 = -7.0
      ENDIF
C
C
      RT_148 = RT_151 * RT_149
C
      KB12 = RT_148
C
      RT_153 = RT_151 * RT_154
C
      RT_152 = + RT_148 + RT_153
C
      KB34 = RT_153
C
      NPP = RT_152 * RT_155
C
C
      IF (ABS(NPP) .LT. 1.0E-10) THEN
         IF (NPP .LT. 0.0)  THEN
            RT_156 = -1.0E10 * ILS
         ELSE
            RT_156 =  1.0E10 * ILS
         ENDIF
      ELSE
         RT_156 = ILS / NPP
      ENDIF
C
C
      RT_120 = RT_121 * RT_156
C
      RT_118 = + RT_120 - RT_119
C
C  Hard Limit
      RT_117 = LIMIT(-1.0, 1.0, RT_118)
C
C  ArcCosine Function
      RT_116 = ACOS(RT_117)
C
      RT_138 = AMIN1(RT_114,RT_116)
C
C  Rate Limiter
       RT_139 = RAMP3(2.0e+03, 2.0e+03, RT_138)
C
C  Gain
      RT_142 = 57.30 * RT_139
C
      PGB(16) = RT_142
      AR = RT_139
C  Gain
      RT_143 = 57.30 * RT_116
C
      PGB(17) = RT_143
C ======================================================================
C   Subsystem: "TF Power Supply"
C ----------------------------------------------------------------------
C  TRANSFORMER SATURATION SUBROUTINE
      X_KNEE = 1.250
      CALL TSAT21(24,22, 22,11, 11,24, 1,
     +     8.3333333,7.50,0.20,X_KNEE,60.0,1.0,0.010, 0.0, 0)
C
C  TRANSFORMER SATURATION SUBROUTINE
      X_KNEE = 1.250
      CALL TSAT21(18,20, 20,19, 19,18, 1,
     +     8.3333333,7.50,0.20,X_KNEE,60.0,1.0,0.010, 0.0, 0)
C
      GI = VAR(9)
      RT_8 = 0.0
C  G6P200 Valve Group Model
      CALL G6P200( 1 , 12 , 13 , 14 ,
     +           29 , 34 , 36 , 31 , 29 , 35 ,
     +           34 , 31 , 29 , 36 , 35 , 31 ,
     +           AR , 1.0e+04 , GI , 70.0 , NINT(KB12) , 1 , -1 ,
     +           RT_1 , RT_2 , 666.60 , 0.030 ,
     +           0.005070 , 1.0e+08 , 0.01880 , 1.0e+06 , 0.0 )
C
      FP2(1) = STOR(NEXC-168)
      FP2(2) = STOR(NEXC-167)
      FP2(3) = STOR(NEXC-166)
      FP2(4) = STOR(NEXC-165)
      FP2(5) = STOR(NEXC-164)
      FP2(6) = STOR(NEXC-163)
      VT2(1) = VDC(34,1)-VDC(29,1)
      VT2(2) = VDC(31,1)-VDC(36,1)
      VT2(3) = VDC(35,1)-VDC(29,1)
      VT2(4) = VDC(31,1)-VDC(34,1)
      VT2(5) = VDC(36,1)-VDC(29,1)
      VT2(6) = VDC(31,1)-VDC(35,1)
      IT2(1) = STOR(NEXC-174)
      IT2(2) = STOR(NEXC-173)
      IT2(3) = STOR(NEXC-172)
      IT2(4) = STOR(NEXC-171)
      IT2(5) = STOR(NEXC-170)
      IT2(6) = STOR(NEXC-169)
C
C  G6P200 Valve Group Model
      CALL G6P200( 1 , 27 , 25 , 26 ,
     +           30 , 40 , 42 , 31 , 30 , 41 ,
     +           40 , 31 , 30 , 42 , 41 , 31 ,
     +           AR , 1.0e+04 , GI , 70.0 , NINT(KB34) , 1 , 0 ,
     +           RT_3 , RT_4 , 666.60 , 0.030 ,
     +           0.005070 , 1.0e+08 , 0.01880 , 1.0e+06 , 0.0 )
C
      FP3(1) = STOR(NEXC-168)
      FP3(2) = STOR(NEXC-167)
      FP3(3) = STOR(NEXC-166)
      FP3(4) = STOR(NEXC-165)
      FP3(5) = STOR(NEXC-164)
      FP3(6) = STOR(NEXC-163)
      VT3(1) = VDC(40,1)-VDC(30,1)
      VT3(2) = VDC(31,1)-VDC(42,1)
      VT3(3) = VDC(41,1)-VDC(30,1)
      VT3(4) = VDC(31,1)-VDC(40,1)
      VT3(5) = VDC(42,1)-VDC(30,1)
      VT3(6) = VDC(31,1)-VDC(41,1)
      IT3(1) = STOR(NEXC-174)
      IT3(2) = STOR(NEXC-173)
      IT3(3) = STOR(NEXC-172)
      IT3(4) = STOR(NEXC-171)
      IT3(5) = STOR(NEXC-170)
      IT3(6) = STOR(NEXC-169)
C
C  G6P200 Valve Group Model
      CALL G6P200( 1 , 27 , 25 , 26 ,
     +           43 , 39 , 37 , 31 , 43 , 38 ,
     +           39 , 31 , 43 , 37 , 38 , 31 ,
     +           AR , 1.0e+04 , GI , 70.0 , NINT(KB34) , 1 , -1 ,
     +           RT_6 , RT_7 , 666.60 , 0.030 ,
     +           0.005070 , 1.0e+08 , 0.01880 , 1.0e+06 , 0.0 )
C
      FP4(1) = STOR(NEXC-168)
      FP4(2) = STOR(NEXC-167)
      FP4(3) = STOR(NEXC-166)
      FP4(4) = STOR(NEXC-165)
      FP4(5) = STOR(NEXC-164)
      FP4(6) = STOR(NEXC-163)
      VT4(1) = VDC(39,1)-VDC(43,1)
      VT4(2) = VDC(31,1)-VDC(37,1)
      VT4(3) = VDC(38,1)-VDC(43,1)
      VT4(4) = VDC(31,1)-VDC(39,1)
      VT4(5) = VDC(37,1)-VDC(43,1)
      VT4(6) = VDC(31,1)-VDC(38,1)
      IT4(1) = STOR(NEXC-174)
      IT4(2) = STOR(NEXC-173)
      IT4(3) = STOR(NEXC-172)
      IT4(4) = STOR(NEXC-171)
      IT4(5) = STOR(NEXC-170)
      IT4(6) = STOR(NEXC-169)
C
C  G6P200 Valve Group Model
      CALL G6P200( 1 , 12 , 13 , 14 ,
     +           28 , 32 , 44 , 31 , 28 , 33 ,
     +           32 , 31 , 28 , 44 , 33 , 31 ,
     +           AR , 1.0e+04 , GI , 70.0 , NINT(KB12) , 1 , 0 ,
     +           RT_9 , RT_10 , 666.660 , 0.030 ,
     +           0.005070 , 1.0e+08 , 0.01880 , 1.0e+06 , 0.0 )
C
      FP1(1) = STOR(NEXC-168)
      FP1(2) = STOR(NEXC-167)
      FP1(3) = STOR(NEXC-166)
      FP1(4) = STOR(NEXC-165)
      FP1(5) = STOR(NEXC-164)
      FP1(6) = STOR(NEXC-163)
      VT1(1) = VDC(32,1)-VDC(28,1)
      VT1(2) = VDC(31,1)-VDC(44,1)
      VT1(3) = VDC(33,1)-VDC(28,1)
      VT1(4) = VDC(31,1)-VDC(32,1)
      VT1(5) = VDC(44,1)-VDC(28,1)
      VT1(6) = VDC(31,1)-VDC(33,1)
      IT1(1) = STOR(NEXC-174)
      IT1(2) = STOR(NEXC-173)
      IT1(3) = STOR(NEXC-172)
      IT1(4) = STOR(NEXC-171)
      IT1(5) = STOR(NEXC-170)
      IT1(6) = STOR(NEXC-169)
C
C  Gain
      RT_5 = 57.2957795 * RT_10
C
C  Gain
      AM = 57.2957795 * RT_9
C
      PGB(18) = AM
      PGB(19) = RT_5
      PGB(20) = FLOAT(FP3(1))
      PGB(21) = VT3(1)
      PGB(22) = IT3(1)
      PGB(23) = FLOAT(FP4(1))
      PGB(24) = VT4(1)
      PGB(25) = IT4(1)
      PGB(26) = FLOAT(FP2(1))
      PGB(27) = VT2(1)
      PGB(28) = IT2(1)
      PGB(29) = IT1(1)
      PGB(30) = IT1(2)
      PGB(31) = IT1(3)
      PGB(32) = IT1(4)
      PGB(33) = IT1(5)
      PGB(34) = IT1(6)
      PGB(35) = FLOAT(FP1(1))
      PGB(36) = VT1(1)
C  Single Phase Current Source.
      C_INST = IROUTER
      CCIN(4, 1) = CCIN(4, 1) + C_INST
      CCIN(6, 1) = CCIN(6, 1) - C_INST
C
      IF (TIME .LT. DELT) THEN
         IF (ABS(CCIN(4, 1)) .LE. 1.0E-6) CCIN(4, 1) = 1.0E-6
         IF (ABS(CCIN(6, 1)) .LE. 1.0E-6) CCIN(6, 1) = 1.0E-6
      ENDIF
C
C  Single Phase Current Source.
      C_INST = IRINNER
      CCIN(5, 1) = CCIN(5, 1) + C_INST
      CCIN(31, 1) = CCIN(31, 1) - C_INST
C
      IF (TIME .LT. DELT) THEN
         IF (ABS(CCIN(5, 1)) .LE. 1.0E-6) CCIN(5, 1) = 1.0E-6
         IF (ABS(CCIN(31, 1)) .LE. 1.0E-6) CCIN(31, 1) = 1.0E-6
      ENDIF
C
C  3-Phase source: MG
      A_MAG = 1.0
      B_MAG = 1.0
      C_MAG = 1.0
      SRC_V = Vrms
      SRC_F = Freq
      SRC_P = RT_8
      CALL ESYS65(1, 3,2,1,0, 1, 3,2,1,
     +            -1.0, 0.00140, 2.14e-04, 1.0, 0.10,
     +            70.0, SRC_F, 0, SRC_V, 0.0010,
     +            A_MAG, B_MAG, C_MAG,
     +            SRC_P, 13.80, 475.0,
     +            1.0, 230.0, 0.020, 0.050,
     +            1.0, 0.020, 0.050, 0,
     +            P_MES, Q_MES, E_RMS, THETA_D,
     +            A_CURR , B_CURR , C_CURR )
C
C ======================================================================
      RETURN
      END