diff --git a/src/spicelib/devices/bsim4/b4ld.c b/src/spicelib/devices/bsim4/b4ld.c
index c849d6393..ef08bc422 100644
--- a/src/spicelib/devices/bsim4/b4ld.c
+++ b/src/spicelib/devices/bsim4/b4ld.c
@@ -97,8 +97,8 @@ BSIM4model *model;
 BSIM4model *model = (BSIM4model*)inModel;
 BSIM4instance *here;
 #endif
-double ceqgstot, dgstot_dvd, dgstot_dvg, dgstot_dvs, dgstot_dvb;
-double ceqgdtot, dgdtot_dvd, dgdtot_dvg, dgdtot_dvs, dgdtot_dvb;
+double ceqgstot, ceqgstot_fvk_sNodePrime, ceqgstot_fvk_sNode, dgstot_dvd, dgstot_dvg, dgstot_dvs, dgstot_dvb;
+double ceqgdtot, ceqgdtot_fvk_dNodePrime, ceqgdtot_fvk_dNode, dgdtot_dvd, dgdtot_dvg, dgdtot_dvs, dgdtot_dvb;
 double gstot, gstotd, gstotg, gstots, gstotb, gspr, Rs, Rd;
 double gdtot, gdtotd, gdtotg, gdtots, gdtotb, gdpr;
 double vgs_eff, vgd_eff, dvgs_eff_dvg, dvgd_eff_dvg;
@@ -127,7 +127,6 @@ double delvbd_jct, delvbs_jct, vbs_jct, vbd_jct;
 double SourceSatCurrent, DrainSatCurrent;
 double ag0, qgb, von, cbhat, VgstNVt, ExpVgst;
 double ceqqb, ceqqb_fvk, ceqqd, ceqqd_fvk, ceqqg, ceqqg_fvk, ceqqjd=0.0, ceqqjd_fvk = 0.0, ceqqjs=0.0, ceqqjs_fvk = 0.0, ceq, geq;
-double ceqqb_SnodePrime_fvk = 0.0, ceqqd_SnodePrime_fvk = 0.0, ceqqg_SnodePrime_fvk = 0.0, ceqqjd_SnodePrime_fvk = 0.0, ceqqjs_SnodePrime_fvk = 0.0 ;
 double cdrain, cdhat, ceqdrn, ceqdrn_fvk, ceqbd, ceqbd_fvk, ceqbs, ceqbs_fvk, ceqjd, ceqjd_fvk, ceqjs, ceqjs_fvk, gjbd, gjbs;
 double czbd, czbdsw, czbdswg, czbs, czbssw, czbsswg, evbd, evbs, arg, sarg;
 double delvbd, delvbs, delvds, delvgd, delvgs;
@@ -4614,15 +4613,12 @@ line860:
 
           ceqqg = cqgate - gcggb * vgb + gcgdb * vbd + gcgsb * vbs;
           ceqqg_fvk = cqgate ;
-          ceqqg_SnodePrime_fvk = ceqqg ;
           ceqqd = cqdrn - gcdgb * vgb - gcdgmb * vgmb + (gcddb + gcdbdb)
                 * vbd - gcdbdb * vbd_jct + gcdsb * vbs;
           ceqqd_fvk = cqdrn ;
-          ceqqd_SnodePrime_fvk = ceqqd ;
           ceqqb = cqbody - gcbgb * vgb - gcbgmb * vgmb
                 + gcbdb * vbd + gcbsb * vbs;
           ceqqb_fvk = cqbody ;
-          ceqqb_SnodePrime_fvk = ceqqb ;
 
           if (here->BSIM4rgateMod == 3)
           {
@@ -4636,10 +4632,8 @@ line860:
           if (here->BSIM4rbodyMod)
           {   ceqqjs = *(ckt->CKTstate0 + here->BSIM4cqbs) + gcsbsb * vbs_jct;
               ceqqjs_fvk = *(ckt->CKTstate0 + here->BSIM4cqbs) ;
-              ceqqjs_SnodePrime_fvk = ceqqjs ;
               ceqqjd = *(ckt->CKTstate0 + here->BSIM4cqbd) + gcdbdb * vbd_jct;
               ceqqjd_fvk = *(ckt->CKTstate0 + here->BSIM4cqbd) ;
-              ceqqjd_SnodePrime_fvk = ceqqjd ;
           }
 
           if (here->BSIM4trnqsMod)
@@ -5350,23 +5344,42 @@ line900:
 ///////////////////////////////////////////
 
         /* KCL - Non-Linear and Dynamic Linear Parts */
-        *(ckt->CKTfvk+here->BSIM4dNodePrime) -= m * (ceqjd_fvk - ceqbd_fvk - ceqdrn_fvk - ceqqd_fvk + Idtoteq_fvk) ; /* NO ceqgdtot in FVK */
-        *(here->KCLcurrentdNodePrimeRHS_1) = -(m * ceqjd_fvk) ;
-        *(here->KCLcurrentdNodePrimeRHS_2) = m * ceqbd_fvk ;
-        *(here->KCLcurrentdNodePrimeRHS_3) = m * ceqdrn_fvk ;
-        *(here->KCLcurrentdNodePrimeRHS_4) = m * ceqqd_fvk ;
-        *(here->KCLcurrentdNodePrimeRHS_5) = -(m * Idtoteq_fvk) ;
+        ceqgdtot_fvk_dNodePrime = (- gdtotd * *(ckt->CKTrhsOld+here->BSIM4dNodePrime)
+                                   - gdtot * *(ckt->CKTrhsOld+here->BSIM4dNode)
+                                   - gdtotg * *(ckt->CKTrhsOld+here->BSIM4gNodePrime)
+                                   - gdtots * *(ckt->CKTrhsOld+here->BSIM4sNodePrime)
+                                   - gdtotb * *(ckt->CKTrhsOld+here->BSIM4bNodePrime))
+                                  - ceqgdtot ; //A*V-b
+        *(ckt->CKTfvk+here->BSIM4dNodePrime) += m * ceqgdtot_fvk_dNodePrime ;
+        *(here->KCLcurrentdNodePrimeRHS_1) = m * ceqgdtot_fvk_dNodePrime ;
+        *(ckt->CKTfvk+here->BSIM4dNodePrime) -= m * (ceqjd_fvk - ceqbd_fvk - ceqdrn_fvk - ceqqd_fvk + Idtoteq_fvk) ;
+        *(here->KCLcurrentdNodePrimeRHS_2) = -(m * ceqjd_fvk) ;
+        *(here->KCLcurrentdNodePrimeRHS_3) = m * ceqbd_fvk ;
+        *(here->KCLcurrentdNodePrimeRHS_4) = m * ceqdrn_fvk ;
+        *(here->KCLcurrentdNodePrimeRHS_5) = m * ceqqd_fvk ;
+        *(here->KCLcurrentdNodePrimeRHS_6) = -(m * Idtoteq_fvk) ;
 
-        *(ckt->CKTfvk+here->BSIM4gNodePrime) += m * (ceqqg_fvk + Igtoteq_fvk) ; /* NO ceqgcrg in FVK */
+        *(ckt->CKTfvk+here->BSIM4gNodePrime) += m * (ceqqg_fvk - ceqgcrg + Igtoteq_fvk) ;
         *(here->KCLcurrentgNodePrimeRHS_1) = m * ceqqg_fvk ;
-        *(here->KCLcurrentgNodePrimeRHS_2) = m * Igtoteq_fvk ;
+        *(here->KCLcurrentgNodePrimeRHS_2) = -(m * ceqgcrg) ;
+        *(here->KCLcurrentgNodePrimeRHS_3) = m * Igtoteq_fvk ;
 
-        if (here->BSIM4rgateMod == 3)
+        if (here->BSIM4rgateMod == 2)
         {
-            *(ckt->CKTfvk+here->BSIM4gNodeMid) += m * ceqqgmid_fvk ; /* NO ceqgcrg in FVK */
-            *(here->KCLcurrentgNodeMidRHS) = m * ceqqgmid_fvk ;
+            *(ckt->CKTfvk+here->BSIM4gNodeExt) += m * ceqgcrg ;
+            *(here->KCLcurrentgNodeExtRHS) = m * ceqgcrg ;
+        } else if (here->BSIM4rgateMod == 3) {
+            *(ckt->CKTfvk+here->BSIM4gNodeMid) += m * (ceqqgmid_fvk + ceqgcrg) ;
+            *(here->KCLcurrentgNodeMidRHS_1) = m * ceqqgmid_fvk ;
+            *(here->KCLcurrentgNodeMidRHS_2) = m * ceqgcrg ;
         }
 
+        ceqgstot_fvk_sNodePrime = (- gstotd * *(ckt->CKTrhsOld+here->BSIM4dNodePrime)
+                                   - gstotg * *(ckt->CKTrhsOld+here->BSIM4gNodePrime)
+                                   - gstots * *(ckt->CKTrhsOld+here->BSIM4sNodePrime)
+                                   - gstot * *(ckt->CKTrhsOld+here->BSIM4sNode)
+                                   - gstotb * *(ckt->CKTrhsOld+here->BSIM4bNodePrime))
+                                  - ceqgstot ; //A*V-b
         if (!here->BSIM4rbodyMod)
         {
             *(ckt->CKTfvk+here->BSIM4bNodePrime) -= m * (ceqbd_fvk + ceqbs_fvk - ceqjd_fvk - ceqjs_fvk - ceqqb_fvk + Ibtoteq_fvk) ;
@@ -5377,17 +5390,19 @@ line900:
             *(here->KCLcurrentbNodePrimeRHS_5) = m * ceqqb_fvk ;
             *(here->KCLcurrentbNodePrimeRHS_6) = -(m * Ibtoteq_fvk) ;
 
+            *(ckt->CKTfvk+here->BSIM4sNodePrime) += m * ceqgstot_fvk_sNodePrime ;
+            *(here->KCLcurrentsNodePrimeRHS_1) = m * ceqgstot_fvk_sNodePrime ;
             *(ckt->CKTfvk+here->BSIM4sNodePrime) -= m * (ceqdrn_fvk - ceqbs_fvk + ceqjs_fvk
                                                          + ceqqg_fvk + ceqqb_fvk + ceqqd_fvk
-                                                         + ceqqgmid_fvk + Istoteq_fvk) ; /* NO ceqgstot in FVK */
-            *(here->KCLcurrentsNodePrimeRHS_1) = -(m * ceqdrn_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_2) = m * ceqbs_fvk ;
-            *(here->KCLcurrentsNodePrimeRHS_3) = -(m * ceqjs_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_4) = -(m * ceqqg_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_5) = -(m * ceqqb_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_6) = -(m * ceqqd_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_7) = -(m * ceqqgmid_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_8) = -(m * Istoteq_fvk) ;
+                                                         + ceqqgmid_fvk + Istoteq_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_2) = -(m * ceqdrn_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_3) = m * ceqbs_fvk ;
+            *(here->KCLcurrentsNodePrimeRHS_4) = -(m * ceqjs_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_5) = -(m * ceqqg_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_6) = -(m * ceqqb_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_7) = -(m * ceqqd_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_8) = -(m * ceqqgmid_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_9) = -(m * Istoteq_fvk) ;
         } else {
             *(ckt->CKTfvk+here->BSIM4dbNode) += m * (ceqjd_fvk + ceqqjd_fvk) ;
             *(here->KCLcurrentdbNodeRHS_1) = m * ceqjd_fvk ;
@@ -5403,19 +5418,42 @@ line900:
             *(here->KCLcurrentsbNodeRHS_1) = m * ceqjs_fvk ;
             *(here->KCLcurrentsbNodeRHS_2) = m * ceqqjs_fvk ;
 
+            *(ckt->CKTfvk+here->BSIM4sNodePrime) += m * ceqgstot_fvk_sNodePrime ;
+            *(here->KCLcurrentsNodePrimeRHS_1) = m * ceqgstot_fvk_sNodePrime ;
             *(ckt->CKTfvk+here->BSIM4sNodePrime) -= m * (ceqdrn_fvk - ceqbs_fvk + ceqjs_fvk + ceqqd_fvk
                                                          + ceqqg_fvk + ceqqb_fvk + ceqqjd_fvk + ceqqjs_fvk
-                                                         + ceqqgmid_fvk + Istoteq_fvk) ; /* NO ceqgstot in FVK */
-            *(here->KCLcurrentsNodePrimeRHS_1) = -(m * ceqdrn_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_2) = m * ceqbs_fvk ;
-            *(here->KCLcurrentsNodePrimeRHS_3) = -(m * ceqjs_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_4) = -(m * ceqqd_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_5) = -(m * ceqqg_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_6) = -(m * ceqqb_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_7) = -(m * ceqqjd_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_8) = -(m * ceqqjs_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_9) = -(m * ceqqgmid_fvk) ;
-            *(here->KCLcurrentsNodePrimeRHS_10) = -(m * Istoteq_fvk) ;
+                                                         + ceqqgmid_fvk + Istoteq_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_2) = -(m * ceqdrn_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_3) = m * ceqbs_fvk ;
+            *(here->KCLcurrentsNodePrimeRHS_4) = -(m * ceqjs_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_5) = -(m * ceqqd_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_6) = -(m * ceqqg_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_7) = -(m * ceqqb_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_8) = -(m * ceqqjd_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_9) = -(m * ceqqjs_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_10) = -(m * ceqqgmid_fvk) ;
+            *(here->KCLcurrentsNodePrimeRHS_11) = -(m * Istoteq_fvk) ;
+        }
+
+        if (model->BSIM4rdsMod)
+        {
+            ceqgdtot_fvk_dNode = (+ gdtotg * *(ckt->CKTrhsOld+here->BSIM4gNodePrime)
+                                  + gdtots * *(ckt->CKTrhsOld+here->BSIM4sNodePrime)
+                                  + gdtotb * *(ckt->CKTrhsOld+here->BSIM4bNodePrime)
+                                  + gdtotd * *(ckt->CKTrhsOld+here->BSIM4dNodePrime)
+                                  + gdtot * *(ckt->CKTrhsOld+here->BSIM4dNode))
+                                 - ceqgdtot ; //A*V-b
+            *(ckt->CKTfvk+here->BSIM4dNode) += m * ceqgdtot_fvk_dNode ;
+            *(here->KCLcurrentdNodeRHS) = m * ceqgdtot_fvk_dNode ;
+
+            ceqgstot_fvk_sNode = (+ gstotd * *(ckt->CKTrhsOld+here->BSIM4dNodePrime)
+                                  + gstotg * *(ckt->CKTrhsOld+here->BSIM4gNodePrime)
+                                  + gstotb * *(ckt->CKTrhsOld+here->BSIM4bNodePrime)
+                                  + gstots * *(ckt->CKTrhsOld+here->BSIM4sNodePrime)
+                                  + gstot * *(ckt->CKTrhsOld+here->BSIM4sNode))
+                                 - ceqgstot ; //A*V-b
+            *(ckt->CKTfvk+here->BSIM4sNode) += m * ceqgstot_fvk_sNode ;
+            *(here->KCLcurrentsNodeRHS) = m * ceqgstot_fvk_sNode ;
         }
 
         if (here->BSIM4trnqsMod)
@@ -5427,17 +5465,69 @@ line900:
 
 
         /* KCL - Static Linear Part */
-        *(ckt->CKTfvk+here->BSIM4dNodePrime) += m * gdpr * (*(ckt->CKTrhsOld+here->BSIM4dNodePrime) - *(ckt->CKTrhsOld+here->BSIM4dNode)) ;
-        *(here->KCLcurrentdNodePrime) = m * gdpr * (*(ckt->CKTrhsOld+here->BSIM4dNodePrime) - *(ckt->CKTrhsOld+here->BSIM4dNode)) ;
+        if (!model->BSIM4rdsMod)
+        {
+            *(ckt->CKTfvk+here->BSIM4dNodePrime) += m * gdpr * (*(ckt->CKTrhsOld+here->BSIM4dNodePrime) - *(ckt->CKTrhsOld+here->BSIM4dNode)) ;
+            *(here->KCLcurrentdNodePrime) = m * gdpr * (*(ckt->CKTrhsOld+here->BSIM4dNodePrime) - *(ckt->CKTrhsOld+here->BSIM4dNode)) ;
 
-        *(ckt->CKTfvk+here->BSIM4dNode) -= m * gdpr * (*(ckt->CKTrhsOld+here->BSIM4dNodePrime) - *(ckt->CKTrhsOld+here->BSIM4dNode)) ;
-        *(here->KCLcurrentdNode) = -(m * gdpr * (*(ckt->CKTrhsOld+here->BSIM4dNodePrime) - *(ckt->CKTrhsOld+here->BSIM4dNode))) ;
+            *(ckt->CKTfvk+here->BSIM4dNode) -= m * gdpr * (*(ckt->CKTrhsOld+here->BSIM4dNodePrime) - *(ckt->CKTrhsOld+here->BSIM4dNode)) ;
+            *(here->KCLcurrentdNode) = -(m * gdpr * (*(ckt->CKTrhsOld+here->BSIM4dNodePrime) - *(ckt->CKTrhsOld+here->BSIM4dNode))) ;
 
-        *(ckt->CKTfvk+here->BSIM4sNodePrime) += m * gspr * (*(ckt->CKTrhsOld+here->BSIM4sNodePrime) - *(ckt->CKTrhsOld+here->BSIM4sNode)) ;
-        *(here->KCLcurrentsNodePrime) = m * gspr * (*(ckt->CKTrhsOld+here->BSIM4sNodePrime) - *(ckt->CKTrhsOld+here->BSIM4sNode)) ;
+            *(ckt->CKTfvk+here->BSIM4sNodePrime) += m * gspr * (*(ckt->CKTrhsOld+here->BSIM4sNodePrime) - *(ckt->CKTrhsOld+here->BSIM4sNode)) ;
+            *(here->KCLcurrentsNodePrime) = m * gspr * (*(ckt->CKTrhsOld+here->BSIM4sNodePrime) - *(ckt->CKTrhsOld+here->BSIM4sNode)) ;
+
+            *(ckt->CKTfvk+here->BSIM4sNode) -= m * gspr * (*(ckt->CKTrhsOld+here->BSIM4sNodePrime) - *(ckt->CKTrhsOld+here->BSIM4sNode)) ;
+            *(here->KCLcurrentsNode) = -(m * gspr * (*(ckt->CKTrhsOld+here->BSIM4sNodePrime) - *(ckt->CKTrhsOld+here->BSIM4sNode))) ;
+        }
+
+        if (here->BSIM4rgateMod == 1)
+        {
+            *(ckt->CKTfvk+here->BSIM4gNodeExt) += m * geltd * (*(ckt->CKTrhsOld+here->BSIM4gNodeExt) - *(ckt->CKTrhsOld+here->BSIM4gNodePrime)) ;
+            *(here->KCLcurrentgNodeExt) = m * geltd * (*(ckt->CKTrhsOld+here->BSIM4gNodeExt) - *(ckt->CKTrhsOld+here->BSIM4gNodePrime)) ;
+
+            *(ckt->CKTfvk+here->BSIM4gNodePrime) -= m * geltd * (*(ckt->CKTrhsOld+here->BSIM4gNodeExt) - *(ckt->CKTrhsOld+here->BSIM4gNodePrime)) ;
+            *(here->KCLcurrentgNodePrime) = -(m * geltd * (*(ckt->CKTrhsOld+here->BSIM4gNodeExt) - *(ckt->CKTrhsOld+here->BSIM4gNodePrime))) ;
+        } else if (here->BSIM4rgateMod == 3) {
+            *(ckt->CKTfvk+here->BSIM4gNodeExt) += m * geltd * (*(ckt->CKTrhsOld+here->BSIM4gNodeExt) - *(ckt->CKTrhsOld+here->BSIM4gNodeMid)) ;
+            *(here->KCLcurrentgNodeExt) = m * geltd * (*(ckt->CKTrhsOld+here->BSIM4gNodeExt) - *(ckt->CKTrhsOld+here->BSIM4gNodeMid)) ;
+
+            *(ckt->CKTfvk+here->BSIM4gNodeMid) -= m * geltd * (*(ckt->CKTrhsOld+here->BSIM4gNodeExt) - *(ckt->CKTrhsOld+here->BSIM4gNodeMid)) ;
+            *(here->KCLcurrentgNodeMid) = -(m * geltd * (*(ckt->CKTrhsOld+here->BSIM4gNodeExt) - *(ckt->CKTrhsOld+here->BSIM4gNodeMid))) ;
+        }
+
+        if (here->BSIM4rbodyMod)
+        {
+            *(ckt->CKTfvk+here->BSIM4dbNode) += m * here->BSIM4grbpd * (*(ckt->CKTrhsOld+here->BSIM4dbNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime)) ;
+            *(here->KCLcurrentdbNode_1) = m * here->BSIM4grbpd * (*(ckt->CKTrhsOld+here->BSIM4dbNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime)) ;
+
+            *(ckt->CKTfvk+here->BSIM4bNodePrime) -= m * here->BSIM4grbpd * (*(ckt->CKTrhsOld+here->BSIM4dbNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime)) ;
+            *(here->KCLcurrentbNodePrime_1) = -(m * here->BSIM4grbpd * (*(ckt->CKTrhsOld+here->BSIM4dbNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime))) ;
+
+            *(ckt->CKTfvk+here->BSIM4dbNode) += m * here->BSIM4grbdb * (*(ckt->CKTrhsOld+here->BSIM4dbNode) - *(ckt->CKTrhsOld+here->BSIM4bNode)) ;
+            *(here->KCLcurrentdbNode_2) = m * here->BSIM4grbdb * (*(ckt->CKTrhsOld+here->BSIM4dbNode) - *(ckt->CKTrhsOld+here->BSIM4bNode)) ;
+
+            *(ckt->CKTfvk+here->BSIM4bNode) -= m * here->BSIM4grbdb * (*(ckt->CKTrhsOld+here->BSIM4dbNode) - *(ckt->CKTrhsOld+here->BSIM4bNode)) ;
+            *(here->KCLcurrentbNode_1) = -(m * here->BSIM4grbdb * (*(ckt->CKTrhsOld+here->BSIM4dbNode) - *(ckt->CKTrhsOld+here->BSIM4bNode))) ;
+
+            *(ckt->CKTfvk+here->BSIM4bNode) += m * here->BSIM4grbpb * (*(ckt->CKTrhsOld+here->BSIM4bNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime)) ;
+            *(here->KCLcurrentbNode_2) = m * here->BSIM4grbpb * (*(ckt->CKTrhsOld+here->BSIM4bNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime)) ;
+
+            *(ckt->CKTfvk+here->BSIM4bNodePrime) -= m * here->BSIM4grbpb * (*(ckt->CKTrhsOld+here->BSIM4bNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime)) ;
+            *(here->KCLcurrentbNodePrime_2) = -(m * here->BSIM4grbpb * (*(ckt->CKTrhsOld+here->BSIM4bNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime))) ;
+
+            *(ckt->CKTfvk+here->BSIM4sbNode) += m * here->BSIM4grbps * (*(ckt->CKTrhsOld+here->BSIM4sbNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime)) ;
+            *(here->KCLcurrentsbNode_1) = m * here->BSIM4grbps * (*(ckt->CKTrhsOld+here->BSIM4sbNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime)) ;
+
+            *(ckt->CKTfvk+here->BSIM4bNodePrime) -= m * here->BSIM4grbps * (*(ckt->CKTrhsOld+here->BSIM4sbNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime)) ;
+            *(here->KCLcurrentbNodePrime_3) = -(m * here->BSIM4grbps * (*(ckt->CKTrhsOld+here->BSIM4sbNode) - *(ckt->CKTrhsOld+here->BSIM4bNodePrime))) ;
+
+            *(ckt->CKTfvk+here->BSIM4sbNode) += m * here->BSIM4grbsb * (*(ckt->CKTrhsOld+here->BSIM4sbNode) - *(ckt->CKTrhsOld+here->BSIM4bNode)) ;
+            *(here->KCLcurrentsbNode_2) = m * here->BSIM4grbsb * (*(ckt->CKTrhsOld+here->BSIM4sbNode) - *(ckt->CKTrhsOld+here->BSIM4bNode)) ;
+
+            *(ckt->CKTfvk+here->BSIM4bNode) -= m * here->BSIM4grbsb * (*(ckt->CKTrhsOld+here->BSIM4sbNode) - *(ckt->CKTrhsOld+here->BSIM4bNode)) ;
+            *(here->KCLcurrentbNode_3) = -(m * here->BSIM4grbsb * (*(ckt->CKTrhsOld+here->BSIM4sbNode) - *(ckt->CKTrhsOld+here->BSIM4bNode))) ;
+        }
 
-        *(ckt->CKTfvk+here->BSIM4sNode) -= m * gspr * (*(ckt->CKTrhsOld+here->BSIM4sNodePrime) - *(ckt->CKTrhsOld+here->BSIM4sNode)) ;
-        *(here->KCLcurrentsNode) = -(m * gspr * (*(ckt->CKTrhsOld+here->BSIM4sNodePrime) - *(ckt->CKTrhsOld+here->BSIM4sNode))) ;
 ///////////////////////////////////////////
 #endif
 
diff --git a/src/spicelib/devices/bsim4/b4node.c b/src/spicelib/devices/bsim4/b4node.c
index 1f4f9355e..2b87aacd3 100644
--- a/src/spicelib/devices/bsim4/b4node.c
+++ b/src/spicelib/devices/bsim4/b4node.c
@@ -34,12 +34,25 @@ BSIM4nodeIsNonLinear (GENmodel *inModel, CKTcircuit *ckt)
 
             ckt->CKTnodeIsLinear [here->BSIM4dNodePrime] = 0 ;
             ckt->CKTnodeIsLinear [here->BSIM4sNodePrime] = 0 ;
-//            ckt->CKTnodeIsLinear [here->BSIM4gNodePrime] = 0 ;
-//            ckt->CKTnodeIsLinear [here->BSIM4gNodeMid] = 0 ;
-            ckt->CKTnodeIsLinear [here->BSIM4dbNode] = 0 ;
+            ckt->CKTnodeIsLinear [here->BSIM4gNodePrime] = 0 ;
             ckt->CKTnodeIsLinear [here->BSIM4bNodePrime] = 0 ;
-            ckt->CKTnodeIsLinear [here->BSIM4sbNode] = 0 ;
-            ckt->CKTnodeIsLinear [here->BSIM4qNode] = 0 ;
+
+            if (!here->BSIM4rbodyMod)
+            {
+                ckt->CKTnodeIsLinear [here->BSIM4dbNode] = 0 ;
+                ckt->CKTnodeIsLinear [here->BSIM4sbNode] = 0 ;
+            }
+
+            if (model->BSIM4rdsMod)
+            {
+                ckt->CKTnodeIsLinear [here->BSIM4dNode] = 0 ;
+                ckt->CKTnodeIsLinear [here->BSIM4sNode] = 0 ;
+            }
+
+            if (here->BSIM4trnqsMod)
+            {
+                ckt->CKTnodeIsLinear [here->BSIM4qNode] = 0 ;
+            }
 
 
             /* KCL - Non-Linear and Dynamic Linear Parts */
@@ -48,11 +61,19 @@ BSIM4nodeIsNonLinear (GENmodel *inModel, CKTcircuit *ckt)
             error = CKTmkCurKCL (ckt, here->BSIM4dNodePrime, &(here->KCLcurrentdNodePrimeRHS_3)) ;
             error = CKTmkCurKCL (ckt, here->BSIM4dNodePrime, &(here->KCLcurrentdNodePrimeRHS_4)) ;
             error = CKTmkCurKCL (ckt, here->BSIM4dNodePrime, &(here->KCLcurrentdNodePrimeRHS_5)) ;
+            error = CKTmkCurKCL (ckt, here->BSIM4dNodePrime, &(here->KCLcurrentdNodePrimeRHS_6)) ;
+
             error = CKTmkCurKCL (ckt, here->BSIM4gNodePrime, &(here->KCLcurrentgNodePrimeRHS_1)) ;
             error = CKTmkCurKCL (ckt, here->BSIM4gNodePrime, &(here->KCLcurrentgNodePrimeRHS_2)) ;
+            error = CKTmkCurKCL (ckt, here->BSIM4gNodePrime, &(here->KCLcurrentgNodePrimeRHS_3)) ;
 
-            if (here->BSIM4rgateMod == 3)
-                error = CKTmkCurKCL (ckt, here->BSIM4gNodeMid, &(here->KCLcurrentgNodeMidRHS)) ;
+            if (here->BSIM4rgateMod == 2)
+            {
+                error = CKTmkCurKCL (ckt, here->BSIM4gNodeExt, &(here->KCLcurrentgNodeExtRHS)) ;
+            } else if (here->BSIM4rgateMod == 3) {
+                error = CKTmkCurKCL (ckt, here->BSIM4gNodeMid, &(here->KCLcurrentgNodeMidRHS_1)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4gNodeMid, &(here->KCLcurrentgNodeMidRHS_2)) ;
+            }
 
             if (!here->BSIM4rbodyMod)
             {
@@ -62,6 +83,7 @@ BSIM4nodeIsNonLinear (GENmodel *inModel, CKTcircuit *ckt)
                 error = CKTmkCurKCL (ckt, here->BSIM4bNodePrime, &(here->KCLcurrentbNodePrimeRHS_4)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4bNodePrime, &(here->KCLcurrentbNodePrimeRHS_5)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4bNodePrime, &(here->KCLcurrentbNodePrimeRHS_6)) ;
+
                 error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrimeRHS_1)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrimeRHS_2)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrimeRHS_3)) ;
@@ -73,12 +95,15 @@ BSIM4nodeIsNonLinear (GENmodel *inModel, CKTcircuit *ckt)
             } else {
                 error = CKTmkCurKCL (ckt, here->BSIM4dbNode, &(here->KCLcurrentdbNodeRHS_1)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4dbNode, &(here->KCLcurrentdbNodeRHS_2)) ;
+
                 error = CKTmkCurKCL (ckt, here->BSIM4bNodePrime, &(here->KCLcurrentbNodePrimeRHS_1)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4bNodePrime, &(here->KCLcurrentbNodePrimeRHS_2)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4bNodePrime, &(here->KCLcurrentbNodePrimeRHS_3)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4bNodePrime, &(here->KCLcurrentbNodePrimeRHS_4)) ;
+
                 error = CKTmkCurKCL (ckt, here->BSIM4sbNode, &(here->KCLcurrentsbNodeRHS_1)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4sbNode, &(here->KCLcurrentsbNodeRHS_2)) ;
+
                 error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrimeRHS_1)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrimeRHS_2)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrimeRHS_3)) ;
@@ -89,18 +114,54 @@ BSIM4nodeIsNonLinear (GENmodel *inModel, CKTcircuit *ckt)
                 error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrimeRHS_8)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrimeRHS_9)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrimeRHS_10)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrimeRHS_11)) ;
+            }
+
+            if (model->BSIM4rdsMod)
+            {
+                error = CKTmkCurKCL (ckt, here->BSIM4dNode, &(here->KCLcurrentdNodeRHS)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4sNode, &(here->KCLcurrentsNodeRHS)) ;
             }
 
             if (here->BSIM4trnqsMod)
+            {
                 error = CKTmkCurKCL (ckt, here->BSIM4qNode, &(here->KCLcurrentqNodeRHS_1)) ;
                 error = CKTmkCurKCL (ckt, here->BSIM4qNode, &(here->KCLcurrentqNodeRHS_2)) ;
+            }
 
 
             /* KCL - Static Linear Part */
-            error = CKTmkCurKCL (ckt, here->BSIM4dNodePrime, &(here->KCLcurrentdNodePrime)) ;
-            error = CKTmkCurKCL (ckt, here->BSIM4dNode, &(here->KCLcurrentdNode)) ;
-            error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrime)) ;
-            error = CKTmkCurKCL (ckt, here->BSIM4sNode, &(here->KCLcurrentsNode)) ;
+            if (!model->BSIM4rdsMod)
+            {
+                error = CKTmkCurKCL (ckt, here->BSIM4dNodePrime, &(here->KCLcurrentdNodePrime)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4dNode, &(here->KCLcurrentdNode)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4sNodePrime, &(here->KCLcurrentsNodePrime)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4sNode, &(here->KCLcurrentsNode)) ;
+            } else {
+            }
+
+            if (here->BSIM4rgateMod == 1)
+            {
+                error = CKTmkCurKCL (ckt, here->BSIM4gNodeExt, &(here->KCLcurrentgNodeExt)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4gNodePrime, &(here->KCLcurrentgNodePrime)) ;
+            } else if (here->BSIM4rgateMod == 3) {
+                error = CKTmkCurKCL (ckt, here->BSIM4gNodeExt, &(here->KCLcurrentgNodeExt)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4gNodeMid, &(here->KCLcurrentgNodeMid)) ;
+            }
+
+            if (here->BSIM4rbodyMod)
+            {
+                error = CKTmkCurKCL (ckt, here->BSIM4dbNode, &(here->KCLcurrentdbNode_1)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4dbNode, &(here->KCLcurrentdbNode_2)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4bNodePrime, &(here->KCLcurrentbNodePrime_1)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4bNodePrime, &(here->KCLcurrentbNodePrime_2)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4bNodePrime, &(here->KCLcurrentbNodePrime_3)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4bNode, &(here->KCLcurrentbNode_1)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4bNode, &(here->KCLcurrentbNode_2)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4bNode, &(here->KCLcurrentbNode_3)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4sbNode, &(here->KCLcurrentsbNode_1)) ;
+                error = CKTmkCurKCL (ckt, here->BSIM4sbNode, &(here->KCLcurrentsbNode_2)) ;
+            }
         }
     }
 
diff --git a/src/spicelib/devices/bsim4/bsim4def.h b/src/spicelib/devices/bsim4/bsim4def.h
index 08a4b0255..d54f49506 100644
--- a/src/spicelib/devices/bsim4/bsim4def.h
+++ b/src/spicelib/devices/bsim4/bsim4def.h
@@ -570,20 +570,30 @@ typedef struct sBSIM4instance
 #endif /* NONOISE */
 
 #ifdef KIRCHHOFF
+//RHS
     double *KCLcurrentdNodePrimeRHS_1 ;
     double *KCLcurrentdNodePrimeRHS_2 ;
     double *KCLcurrentdNodePrimeRHS_3 ;
     double *KCLcurrentdNodePrimeRHS_4 ;
     double *KCLcurrentdNodePrimeRHS_5 ;
+    double *KCLcurrentdNodePrimeRHS_6 ;
+
     double *KCLcurrentgNodePrimeRHS_1 ;
     double *KCLcurrentgNodePrimeRHS_2 ;
-    double *KCLcurrentgNodeMidRHS ;
+    double *KCLcurrentgNodePrimeRHS_3 ;
+
+    double *KCLcurrentgNodeExtRHS ;
+
+    double *KCLcurrentgNodeMidRHS_1 ;
+    double *KCLcurrentgNodeMidRHS_2 ;
+
     double *KCLcurrentbNodePrimeRHS_1 ;
     double *KCLcurrentbNodePrimeRHS_2 ;
     double *KCLcurrentbNodePrimeRHS_3 ;
     double *KCLcurrentbNodePrimeRHS_4 ;
     double *KCLcurrentbNodePrimeRHS_5 ;
     double *KCLcurrentbNodePrimeRHS_6 ;
+
     double *KCLcurrentsNodePrimeRHS_1 ;
     double *KCLcurrentsNodePrimeRHS_2 ;
     double *KCLcurrentsNodePrimeRHS_3 ;
@@ -594,17 +604,39 @@ typedef struct sBSIM4instance
     double *KCLcurrentsNodePrimeRHS_8 ;
     double *KCLcurrentsNodePrimeRHS_9 ;
     double *KCLcurrentsNodePrimeRHS_10 ;
+    double *KCLcurrentsNodePrimeRHS_11 ;
+
     double *KCLcurrentdbNodeRHS_1 ;
     double *KCLcurrentdbNodeRHS_2 ;
+
     double *KCLcurrentsbNodeRHS_1 ;
     double *KCLcurrentsbNodeRHS_2 ;
+
+    double *KCLcurrentdNodeRHS ;
+
+    double *KCLcurrentsNodeRHS ;
+
     double *KCLcurrentqNodeRHS_1 ;
     double *KCLcurrentqNodeRHS_2 ;
 
+//Matrix
     double *KCLcurrentdNodePrime ;
     double *KCLcurrentdNode ;
     double *KCLcurrentsNodePrime ;
     double *KCLcurrentsNode ;
+    double *KCLcurrentgNodePrime ;
+    double *KCLcurrentgNodeMid ;
+    double *KCLcurrentgNodeExt ;
+    double *KCLcurrentdbNode_1 ;
+    double *KCLcurrentdbNode_2 ;
+    double *KCLcurrentbNodePrime_1 ;
+    double *KCLcurrentbNodePrime_2 ;
+    double *KCLcurrentbNodePrime_3 ;
+    double *KCLcurrentbNode_1 ;
+    double *KCLcurrentbNode_2 ;
+    double *KCLcurrentbNode_3 ;
+    double *KCLcurrentsbNode_1 ;
+    double *KCLcurrentsbNode_2 ;
 #endif
 
 } BSIM4instance ;