/***  B4SOI 12/16/2010 Released by Tanvir Morshed  ***/


/**********
 * Copyright 2010 Regents of the University of California.  All rights reserved.
 * Authors: 1998 Samuel Fung, Dennis Sinitsky and Stephen Tang
 * Authors: 1999-2004 Pin Su, Hui Wan, Wei Jin, b3soiacld.c
 * Authors: 2005- Hui Wan, Xuemei Xi, Ali Niknejad, Chenming Hu.
 * Authors: 2009- Wenwei Yang, Chung-Hsun Lin, Ali Niknejad, Chenming Hu.
 * Authors: 2010- Tanvir Morshed, Ali Niknejad, Chenming Hu.
 * File: b4soiacld.c
 * Modified by Hui Wan, Xuemei Xi 11/30/2005
 * Modified by Wenwei Yang, Chung-Hsun Lin, Darsen Lu 03/06/2009
 * Modified by Tanvir Morshed 09/22/2009
 **********/

#include "ngspice/ngspice.h"

#include "ngspice/cktdefs.h"
#include "b4soidef.h"
#include "ngspice/sperror.h"
#include "ngspice/suffix.h"


int
B4SOIacLoad(
GENmodel *inModel,
CKTcircuit *ckt)
{
register B4SOImodel *model = (B4SOImodel*)inModel;
register B4SOIinstance *here;
register int selfheat;
double xcggb, xcgdb, xcgsb, xcgeb, xcgT;
double xcdgb, xcddb, xcdsb, xcdeb, xcdT;
double xcsgb, xcsdb, xcssb, xcseb, xcsT;
double xcbgb, xcbdb, xcbsb, xcbeb, xcbT;
double xcegb, xceeb, xceT;
double gdpr, gspr, gds;
double cggb, cgdb, cgsb, cgT;
double cdgb, cddb, cdsb, cdeb, cdT;
double cbgb, cbdb, cbsb, cbeb, cbT;
double ceeb, ceT;
double GSoverlapCap, GDoverlapCap, GEoverlapCap, FwdSum, RevSum, Gm, Gmbs, GmT;
double omega;
double dxpart, sxpart;
double gbbg, gbbdp, gbbb, gbbp, gbbsp, gbbT;
double gddpg, gddpdp, gddpsp, gddpb, gddpT;
double gsspg, gsspdp, gsspsp, gsspb, gsspT;
double gppb, gppp;
double xcTt, cTt, gTtt, gTtg, gTtb, gTtdp, gTtsp;
double EDextrinsicCap, ESextrinsicCap;
double xcedb, xcesb;

/* v3.0 */
double Gme, gddpe, gsspe, gbbe, gTte;

/* v3.1 added variables for RF */
double T0=0.0;
double gcrgd, gcrgg, gcrgs, gcrgb, gcrg;
double xcgmgmb, xcgmdb, xcgmsb, xcgmeb, xcdgmb, xcsgmb, xcegmb;
double geltd;
double gigg, gigd, gigs, gigb, gige, gigT;
double gigpg, gigpp;

/* v3.1.1 bug fix */
double gIstotg, gIstotd, gIstotb, gIstots;
double gIdtotg, gIdtotd, gIdtotb, gIdtots;
double gIgtotg, gIgtotd, gIgtotb, gIgtots;

/* v4.0 */
double xcdbb, xcsbb, xcdbdb, xcsbsb, xcjdbdp, xcjsbsp;
double gstot, gstotd, gstotg, gstots, gstotb;
double gdtot, gdtotd, gdtotg, gdtots, gdtotb;

double m;

    omega = ckt->CKTomega;
    for (; model != NULL; model = model->B4SOInextModel) 
    {    

      for (here = model->B4SOIinstances; here!= NULL;
              here = here->B4SOInextInstance) 
         {
              selfheat = (model->B4SOIshMod == 1) && (here->B4SOIrth0 != 0.0);
              if (here->B4SOImode >= 0) 
              {   Gm = here->B4SOIgm;
                  Gmbs = here->B4SOIgmbs;

/* v3.0 */
                  Gme = here->B4SOIgme;

                  GmT = model->B4SOItype * here->B4SOIgmT;
                  FwdSum = Gm + Gmbs + Gme; /* v3.0 */
                  RevSum = 0.0;

                  cbgb = here->B4SOIcbgb;
                  cbsb = here->B4SOIcbsb;
                  cbdb = here->B4SOIcbdb;
                  cbeb = here->B4SOIcbeb;
                  cbT  = model->B4SOItype * here->B4SOIcbT;

                  ceeb = here->B4SOIceeb;
                  ceT  = model->B4SOItype * here->B4SOIceT;

                  cggb = here->B4SOIcggb;
                  cgsb = here->B4SOIcgsb;
                  cgdb = here->B4SOIcgdb;
                  cgT  = model->B4SOItype * here->B4SOIcgT;

                  cdgb = here->B4SOIcdgb;
                  cdsb = here->B4SOIcdsb;
                  cddb = here->B4SOIcddb;
                  cdeb = here->B4SOIcdeb;
                  cdT  = model->B4SOItype * here->B4SOIcdT;

                  cTt = here->pParam->B4SOIcth;


/* v3.1 bug fix */
                  gigg = here->B4SOIgigg;
                  gigb = here->B4SOIgigb;
                  gige = here->B4SOIgige; 
                  gigs = here->B4SOIgigs;
                  gigd = here->B4SOIgigd;
                  gigT = model->B4SOItype * here->B4SOIgigT;

                  /* v4.1 */
                  gigpg = here->B4SOIgigpg;
                  gigpp = here->B4SOIgigpp;

                  gbbg  = -here->B4SOIgbgs;
                  gbbdp = -here->B4SOIgbds;
                  gbbb  = -here->B4SOIgbbs;
                  gbbp  = -here->B4SOIgbps;
                  gbbT  = -model->B4SOItype * here->B4SOIgbT;
                  gbbe  = -here->B4SOIgbes;

                  if (here->B4SOIrbodyMod) { /* v4.0 */
                          gbbdp = -here->B4SOIgiigidld;
                        gbbb = -here->B4SOIgbgiigbpb;
                  }

                  gbbsp = - ( gbbg + gbbdp + gbbb + gbbp + gbbe);

                  gddpg  = -here->B4SOIgjdg;
                  gddpdp = -here->B4SOIgjdd;
                  if (!here->B4SOIrbodyMod) /* v4.0 */
                          gddpb  = -here->B4SOIgjdb;
                  else
                        gddpb = here->B4SOIgiigidlb;

                  gddpT  = -model->B4SOItype * here->B4SOIgjdT;

/* v3.0 */
                  gddpe  = -here->B4SOIgjde;
                  gddpsp = - ( gddpg + gddpdp + gddpb + gddpe);

                  gsspg  = -here->B4SOIgjsg;
                  gsspdp = -here->B4SOIgjsd;
                  if (!here->B4SOIrbodyMod) /* v4.0 */
                          gsspb  = -here->B4SOIgjsb;
                  else
                        gsspb = 0.0;
                  gsspT  = -model->B4SOItype * here->B4SOIgjsT;

                  gsspe  = 0.0;
                  gsspsp = - (gsspg + gsspdp + gsspb + gsspe);

                       gppb = -here->B4SOIgbpbs;
                  gppp = -here->B4SOIgbpps;

                  gTtg  = here->B4SOIgtempg;
                  gTtb  = here->B4SOIgtempb;
                  gTtdp = here->B4SOIgtempd;
                  gTtt  = here->B4SOIgtempT;

/* v3.0 */
                  gTte  = here->B4SOIgtempe;
                  gTtsp = - (gTtg + gTtb + gTtdp + gTte);


/* v3.1.1 bug fix */
                  if (model->B4SOIigcMod)
                  {  
                     gIstotg = here->B4SOIgIgsg + here->B4SOIgIgcsg;
                     gIstotd = here->B4SOIgIgcsd;
                     gIstots = here->B4SOIgIgss + here->B4SOIgIgcss;
                     gIstotb = here->B4SOIgIgcsb;

                     gIdtotg = here->B4SOIgIgdg + here->B4SOIgIgcdg;
                     gIdtotd = here->B4SOIgIgdd + here->B4SOIgIgcdd;
                     gIdtots = here->B4SOIgIgcds;
                     gIdtotb = here->B4SOIgIgcdb;

                     gIgtotg = gIstotg + gIdtotg;
                     gIgtotd = gIstotd + gIdtotd;
                     gIgtots = gIstots + gIdtots;
                     gIgtotb = gIstotb + gIdtotb;
                  }
                  else
                  {   
                     gIstotg = gIstotd = gIstots = gIstotb = 0.0;
                     gIdtotg = gIdtotd = gIdtots = gIdtotb = 0.0;
                     gIgtotg = gIgtotd = gIgtots = gIgtotb = 0.0;
                  }


                  sxpart = 0.6;
                  dxpart = 0.4;

                  /* v3.1 for RF */
                  if (here->B4SOIrgateMod == 2)
                      T0 = *(ckt->CKTstates[0] + here->B4SOIvges)
                         - *(ckt->CKTstates[0] + here->B4SOIvgs);
                  else if (here->B4SOIrgateMod == 3)
                      T0 = *(ckt->CKTstates[0] + here->B4SOIvgms)
                         - *(ckt->CKTstates[0] + here->B4SOIvgs);
                  if (here->B4SOIrgateMod > 1)
                  {   gcrgd = here->B4SOIgcrgd * T0;
                      gcrgg = here->B4SOIgcrgg * T0;
                      gcrgs = here->B4SOIgcrgs * T0;
                      gcrgb = here->B4SOIgcrgb * T0;
                      gcrgg -= here->B4SOIgcrg;
                      gcrg = here->B4SOIgcrg;
                  }
                  else
                  gcrg = gcrgd = gcrgg = gcrgs = gcrgb = 0.0;
                  /* v3.1 for RF end*/

              } 
              else
              {   Gm = -here->B4SOIgm;
                  Gmbs = -here->B4SOIgmbs;

/* v3.0 */
                  Gme = -here->B4SOIgme;

                  GmT = -model->B4SOItype * here->B4SOIgmT;
                  FwdSum = 0.0;
                  RevSum = -Gm - Gmbs - Gme; /* v3.0 */

                  cdgb = - (here->B4SOIcdgb + here->B4SOIcggb + here->B4SOIcbgb);
                  cdsb = - (here->B4SOIcddb + here->B4SOIcgdb + here->B4SOIcbdb);
                  cddb = - (here->B4SOIcdsb + here->B4SOIcgsb + here->B4SOIcbsb);
                  cdeb = - (here->B4SOIcdeb + here->B4SOIcbeb + here->B4SOIceeb);
                  cdT  = - model->B4SOItype * (here->B4SOIcgT + here->B4SOIcbT
                         + here->B4SOIcdT + here->B4SOIceT);

                  ceeb = here->B4SOIceeb;
                  ceT  = model->B4SOItype * here->B4SOIceT;

                  cggb = here->B4SOIcggb;
                  cgsb = here->B4SOIcgdb;
                  cgdb = here->B4SOIcgsb;
                  cgT  = model->B4SOItype * here->B4SOIcgT;

                  cbgb = here->B4SOIcbgb;
                  cbsb = here->B4SOIcbdb;
                  cbdb = here->B4SOIcbsb;
                  cbeb = here->B4SOIcbeb;
                  cbT  = model->B4SOItype * here->B4SOIcbT;

                  cTt = here->pParam->B4SOIcth;


/* v3.1 bug fix */
                  gigg = here->B4SOIgigg;
                  gigb = here->B4SOIgigb;
                  gige = here->B4SOIgige; 
                  gigs = here->B4SOIgigd; /* v3.1.1 bug fix */
                  gigd = here->B4SOIgigs; /* v3.1.1 bug fix */
                  gigT = model->B4SOItype * here->B4SOIgigT;

                  gigpg = here->B4SOIgigpg;/* bugfix_snps for setting gigpg gigpp*/
                  gigpp = here->B4SOIgigpp; 
                  
                  gbbg  = -here->B4SOIgbgs;
                  gbbb  = -here->B4SOIgbbs;
                  gbbp  = -here->B4SOIgbps;
                  gbbsp = -here->B4SOIgbds;
                  gbbT  = -model->B4SOItype * here->B4SOIgbT;
/* v3.0 */
                  gbbe  = -here->B4SOIgbes;

                  if (here->B4SOIrbodyMod) { /* v4.0 */
                        gbbsp = -here->B4SOIgiigidld;
                        gbbb = -here->B4SOIgbgiigbpb;
                  }

                  gbbdp = - ( gbbg + gbbsp + gbbb + gbbp + gbbe);

                  gddpg  = -here->B4SOIgjsg;
                  gddpsp = -here->B4SOIgjsd;
                  if (!here->B4SOIrbodyMod) /* v4.0 */
                          gddpb  = -here->B4SOIgjsb;
                  else 
                        gddpb = 0.0;

                  gddpT  = -model->B4SOItype * here->B4SOIgjsT;

/* v3.0 */
                  gddpe  = 0.0;
                  gddpdp = - (gddpg + gddpsp + gddpb + gddpe );

                  gsspg  = -here->B4SOIgjdg;
                  gsspsp = -here->B4SOIgjdd;
                  if (!here->B4SOIrbodyMod) /* v4.0 */
                          gsspb  = -here->B4SOIgjdb;
                  else
                        gsspb = here->B4SOIgiigidlb;

                  gsspT  = -model->B4SOItype * here->B4SOIgjdT;

/* v3.0 */
                  gsspe  = -here->B4SOIgjde;
                  gsspdp = - ( gsspg + gsspsp + gsspb + gsspe );

                  gppb = -here->B4SOIgbpbs;
                  gppp = -here->B4SOIgbpps;

                  gTtg  = here->B4SOIgtempg;
                  gTtb  = here->B4SOIgtempb;
                  gTtsp = here->B4SOIgtempd;
                  gTtt  = here->B4SOIgtempT;

/* v3.0 */
                  gTte = here->B4SOIgtempe;
                  gTtdp = - (gTtg + gTtb + gTtsp + gTte);


/* v3.1.1 bug fix */
                  if (model->B4SOIigcMod)
                  {   
                      gIstotg = here->B4SOIgIgsg + here->B4SOIgIgcdg;
                      gIstotd = here->B4SOIgIgcds;
                      gIstots = here->B4SOIgIgss + here->B4SOIgIgcdd;
                      gIstotb = here->B4SOIgIgcdb;

                      gIdtotg = here->B4SOIgIgdg + here->B4SOIgIgcsg;
                      gIdtotd = here->B4SOIgIgdd + here->B4SOIgIgcss;
                      gIdtots = here->B4SOIgIgcsd;
                      gIdtotb = here->B4SOIgIgcsb;

                      gIgtotg = gIstotg + gIdtotg;
                      gIgtotd = gIstotd + gIdtotd;
                      gIgtots = gIstots + gIdtots;
                      gIgtotb = gIstotb + gIdtotb;
                  }
                  else
                  {   
                      gIstotg = gIstotd = gIstots = gIstotb = 0.0;
                      gIdtotg = gIdtotd = gIdtots = gIdtotb = 0.0;
                      gIgtotg = gIgtotd = gIgtots = gIgtotb = 0.0;
                  }


                  sxpart = 0.4;
                  dxpart = 0.6;

                  /* v3.1 for RF */
                  if (here->B4SOIrgateMod == 2)
                      T0 = *(ckt->CKTstates[0] + here->B4SOIvges)
                         - *(ckt->CKTstates[0] + here->B4SOIvgs);
                  else if (here->B4SOIrgateMod == 3)
                      T0 = *(ckt->CKTstates[0] + here->B4SOIvgms)
                         - *(ckt->CKTstates[0] + here->B4SOIvgs);
                  if (here->B4SOIrgateMod > 1)
                  {   gcrgd = here->B4SOIgcrgs * T0;
                      gcrgg = here->B4SOIgcrgg * T0;
                      gcrgs = here->B4SOIgcrgd * T0;
                      gcrgb = here->B4SOIgcrgb * T0;
                      gcrgg -= here->B4SOIgcrg;
                      gcrg = here->B4SOIgcrg;
                  }
                  else
                  gcrg = gcrgd = gcrgg = gcrgs = gcrgb = 0.0;

                  /* v3.1 RF end*/

              }

              if (!model->B4SOIrdsMod)        {
                      gdpr=here->B4SOIdrainConductance;
                      gspr=here->B4SOIsourceConductance;
              }
              else gdpr = gspr = 0.0;


              gds= here->B4SOIgds;

              GSoverlapCap = here->B4SOIcgso;
              GDoverlapCap = here->B4SOIcgdo;
              GEoverlapCap = here->pParam->B4SOIcgeo;

              EDextrinsicCap = here->B4SOIgcde;
              ESextrinsicCap = here->B4SOIgcse;


              /* v3.1 added for RF */
              if (here->B4SOIrgateMod == 3)
              {
                xcgmgmb = (GDoverlapCap + GSoverlapCap + GEoverlapCap ) * omega;
                xcgmdb = -GDoverlapCap * omega;
                xcgmsb = -GSoverlapCap * omega;
                xcgmeb = -GEoverlapCap * omega;

                xcdgmb = xcgmdb;
                xcsgmb = xcgmsb;
                xcegmb = xcgmeb;

                xcedb = -EDextrinsicCap * omega;
                xcdeb = (cdeb - EDextrinsicCap) * omega; 
                xcddb = (cddb + GDoverlapCap + EDextrinsicCap) * omega;
                xceeb = (ceeb + GEoverlapCap + EDextrinsicCap + ESextrinsicCap)
                      * omega;
                xcesb = -ESextrinsicCap * omega;
                xcssb = (GSoverlapCap + ESextrinsicCap - (cgsb + cbsb + cdsb))
                      * omega;

                xcseb = -(cbeb + cdeb + ceeb + ESextrinsicCap) * omega;

                xcegb = 0; /* v3.1 change for RF */
                xceT  =  ceT * omega;
                xcggb = here->B4SOIcggb * omega; 
                xcgdb = cgdb * omega;     
                xcgsb = cgsb * omega;     
                xcgeb = 0; 
                xcgT  = cgT * omega;

                xcdgb = cdgb * omega;     
                xcdsb = cdsb * omega;
                xcdT  = cdT * omega;

                xcsgb = -(cggb + cbgb + cdgb) * omega; 
                xcsdb = -(cgdb + cbdb + cddb) * omega;
                xcsT  = -(cgT + cbT + cdT + ceT) * omega;

                xcbgb = cbgb * omega;
                xcbdb = cbdb * omega;
                xcbsb = cbsb * omega;
                xcbeb = cbeb * omega;
                xcbT  = cbT * omega;

                xcTt = cTt * omega;
              }

              else 
              {
                xcedb = -EDextrinsicCap * omega;
                xcdeb = (cdeb - EDextrinsicCap) * omega;
                xcddb = (cddb + GDoverlapCap + EDextrinsicCap) * omega;
                xceeb = (ceeb + GEoverlapCap + EDextrinsicCap + ESextrinsicCap)
                      * omega;
                xcesb = -ESextrinsicCap * omega;
                xcssb = (GSoverlapCap + ESextrinsicCap - (cgsb + cbsb + cdsb))
                      * omega;

                xcseb = -(cbeb + cdeb + ceeb + ESextrinsicCap) * omega;

                xcegb = (- GEoverlapCap) * omega;
                xceT  =  ceT * omega;
                xcggb = (cggb + GDoverlapCap + GSoverlapCap + GEoverlapCap)
                      * omega;
                xcgdb = (cgdb - GDoverlapCap ) * omega;
                xcgsb = (cgsb - GSoverlapCap) * omega;
                xcgeb = (- GEoverlapCap) * omega;
                xcgT  = cgT * omega;

                xcdgb = (cdgb - GDoverlapCap) * omega;
                xcdsb = cdsb * omega;
                xcdT  = cdT * omega;

                xcsgb = -(cggb + cbgb + cdgb + GSoverlapCap) * omega;
                xcsdb = -(cgdb + cbdb + cddb) * omega;
                xcsT  = -(cgT + cbT + cdT + ceT) * omega;

                xcbgb = cbgb * omega;
                xcbdb = cbdb * omega;
                xcbsb = cbsb * omega;
                xcbeb = cbeb * omega;
                xcbT  = cbT * omega;

                xcTt = cTt * omega;
        
                 /* v3.1 */
                xcdgmb = xcsgmb = xcegmb = 0.0;
                xcgmgmb = xcgmdb = xcgmsb = xcgmeb =0.0;
              }

              if (here->B4SOImode >= 0) { /* v4.0 */
                  if (!here->B4SOIrbodyMod) {
                        xcjdbdp = xcjsbsp = 0.0;
                        xcdbb = -(xcdgb + xcddb + xcdsb + xcdgmb + xcdeb);
                        xcsbb = -(xcsgb + xcsdb + xcssb + xcsgmb + xcseb);
                        xcdbdb = 0.0; xcsbsb = 0.0;
                        xcbdb = here->B4SOIcbdb * omega;
                        xcbsb = here->B4SOIcbsb * omega;
                  }
                  else {
                        xcjdbdp = here->B4SOIcjdb * omega;
                        xcjsbsp = here->B4SOIcjsb * omega;
                        xcdbb = -(xcdgb + xcddb + xcdsb + xcdgmb + xcdeb)
                                 + xcjdbdp;
                        xcsbb = -(xcsgb + xcsdb + xcssb + xcsgmb + xcseb)
                                 + xcjsbsp;
                        xcdbdb = -here->B4SOIcjdb * omega;
                        xcsbsb = -here->B4SOIcjsb * omega;
                        xcbdb = here->B4SOIcbdb * omega - xcdbdb;
                        xcbsb = here->B4SOIcbsb * omega - xcsbsb;
                  }
              }
              else {
                  if (!here->B4SOIrbodyMod) {
                        xcjdbdp = xcjsbsp = 0.0;
                        xcdbb = -(xcdgb + xcddb + xcdsb + xcdgmb + xcdeb);
                        xcsbb = -(xcsgb + xcsdb + xcssb + xcsgmb + xcseb);
                        xcdbdb = 0.0; xcsbsb = 0.0;
                        xcbdb = here->B4SOIcbsb * omega;
                        xcbsb = here->B4SOIcbdb * omega;
                  }
                  else {
                        xcjdbdp = here->B4SOIcjsb * omega;
                        xcjsbsp = here->B4SOIcjdb * omega;
                        xcdbb = -(xcdgb + xcddb + xcdsb + xcdgmb + xcdeb)
                                 + xcjdbdp;
                        xcsbb = -(xcsgb + xcsdb + xcssb + xcsgmb + xcseb)
                                 + xcjsbsp;
                        xcdbdb = -here->B4SOIcjsb * omega;
                        xcsbsb = -here->B4SOIcjdb * omega;
                        xcbdb = here->B4SOIcbsb * omega - xcdbdb;
                        xcbsb = here->B4SOIcbdb * omega - xcsbsb;
                  }

              }

              if (model->B4SOIrdsMod == 1)
              {   gstot = here->B4SOIgstot;
                  gstotd = here->B4SOIgstotd;
                  gstotg = here->B4SOIgstotg;
                  gstots = here->B4SOIgstots - gstot;
                  gstotb = here->B4SOIgstotb;

                  gdtot = here->B4SOIgdtot;
                  gdtotd = here->B4SOIgdtotd - gdtot;
                  gdtotg = here->B4SOIgdtotg;
                  gdtots = here->B4SOIgdtots;
                  gdtotb = here->B4SOIgdtotb;
              }
              else
              {   gstot = gstotd = gstotg = gstots = gstotb = 0.0;
                  gdtot = gdtotd = gdtotg = gdtots = gdtotb = 0.0;
              }

        m = here->B4SOIm;

              /* v3.1 for RF */
              geltd = here->B4SOIgrgeltd;
              if (here->B4SOIrgateMod == 1)
              {  
                *(here->B4SOIGEgePtr) += m * geltd;
                *(here->B4SOIGgePtr) -= m * geltd;
                *(here->B4SOIGEgPtr) -= m * geltd;
                *(here->B4SOIGgPtr) += m * (geltd + gigg + gIgtotg); /* v3.1.1 bug fix */
                *(here->B4SOIGdpPtr) += m * (gigd + gIgtotd); /* v3.1.1 bug fix */
                *(here->B4SOIGspPtr) += m * (gigs + gIgtots); /* v3.1.1 bug fix */
                if (here->B4SOIsoiMod != 2) /* v3.2 */
                    *(here->B4SOIGbPtr) -= m * (-gigb - gIgtotb); /* v3.1.1 bug fix */
              }

              else if (here->B4SOIrgateMod == 2)
              {   
                *(here->B4SOIGEgePtr) += m * gcrg;
                *(here->B4SOIGEgPtr) += m * gcrgg;
                *(here->B4SOIGEdpPtr) += m * gcrgd;
                *(here->B4SOIGEspPtr) += m * gcrgs;
                if (here->B4SOIsoiMod != 2) /* v3.2 */
                *(here->B4SOIGEbPtr) += m * gcrgb; 

                *(here->B4SOIGgePtr) -= m * gcrg;
                *(here->B4SOIGgPtr) -= m * (gcrgg - gigg - gIgtotg); /* v3.1.1 bug fix */
                *(here->B4SOIGdpPtr) -= m * (gcrgd - gigd - gIgtotd); /* v3.1.1 bug fix */
                *(here->B4SOIGspPtr) -= m * (gcrgs - gigs - gIgtots); /* v3.1.1 bug fix */
                if (here->B4SOIsoiMod != 2) /* v3.2 */
                *(here->B4SOIGbPtr) -= m * (gcrgb - gigb - gIgtotb); /* v3.1.1 bug fix */
              } 

              else if (here->B4SOIrgateMod == 3)
              {   
                *(here->B4SOIGEgePtr) += m * geltd;
                *(here->B4SOIGEgmPtr) -= m * geltd;
                *(here->B4SOIGMgePtr) -= m * geltd;
                *(here->B4SOIGMgmPtr) += m * (geltd + gcrg);
                *(here->B4SOIGMgmPtr +1) += m * xcgmgmb;
 
                *(here->B4SOIGMdpPtr) += m * gcrgd;
                *(here->B4SOIGMdpPtr +1) += m * xcgmdb;
                *(here->B4SOIGMgPtr) += m * gcrgg;
                *(here->B4SOIGMspPtr) += m * gcrgs;
                *(here->B4SOIGMspPtr +1) += m * xcgmsb;
                if (here->B4SOIsoiMod != 2) /* v3.2 */
                *(here->B4SOIGMbPtr) += m * gcrgb;
                *(here->B4SOIGMePtr +1) += m * xcgmeb;
 
                *(here->B4SOIDPgmPtr +1) += m * xcdgmb;
                *(here->B4SOIGgmPtr) -= m * gcrg;
                *(here->B4SOISPgmPtr +1) += m * xcsgmb;
                *(here->B4SOIEgmPtr +1) += m * xcegmb;
 
                *(here->B4SOIGgPtr) -= m * (gcrgg - gigg - gIgtotg); /* v3.1.1 bug fix */
                *(here->B4SOIGdpPtr) -= m * (gcrgd - gigd - gIgtotd); /* v3.1.1 bug fix */
                *(here->B4SOIGspPtr) -= m * (gcrgs - gigs - gIgtots); /* v3.1.1 bug fix */
                if (here->B4SOIsoiMod != 2) /* v3.2 */
                *(here->B4SOIGbPtr) -= m * (gcrgb - gigb - gIgtotb); /* v3.1.1 bug fix */
              }
              else
              {   
                *(here->B4SOIGgPtr)  += m * (gigg + gIgtotg); /* v3.1.1 bug fix */
                *(here->B4SOIGdpPtr) += m * (gigd + gIgtotd); /* v3.1.1 bug fix */
                *(here->B4SOIGspPtr) += m * (gigs + gIgtots); /* v3.1.1 bug fix */
                if (here->B4SOIsoiMod != 2) /* v3.2 */
                *(here->B4SOIGbPtr)  += m * (gigb + gIgtotb); /* v3.1.1 bug fix */
              }
              /* v3.1 for RF end*/

              if (model->B4SOIrdsMod)
              {   (*(here->B4SOIDgPtr) += m * gdtotg);
                  (*(here->B4SOIDspPtr) += m * gdtots);
                  (*(here->B4SOISdpPtr) += m * gstotd);
                  (*(here->B4SOISgPtr) += m * gstotg);
                  if (here->B4SOIsoiMod != 2) {
                          (*(here->B4SOIDbPtr) += m * gdtotb);
                          (*(here->B4SOISbPtr) += m * gstotb);
                  }
              }

              *(here->B4SOIEdpPtr +1) += m * xcedb;
              *(here->B4SOIEspPtr +1) += m * xcesb;
              *(here->B4SOIDPePtr +1) += m * xcdeb;
              *(here->B4SOISPePtr +1) += m * xcseb;
              *(here->B4SOIEgPtr  +1) += m * xcegb;
              *(here->B4SOIGePtr  +1) += m * xcgeb;

              *(here->B4SOIEePtr  +1) += m * xceeb;

              *(here->B4SOIGgPtr  +1) += m * xcggb;
              *(here->B4SOIGdpPtr +1) += m * xcgdb;
              *(here->B4SOIGspPtr +1) += m * xcgsb;

              *(here->B4SOIDPgPtr +1) += m * xcdgb;
              *(here->B4SOIDPdpPtr +1) += m * xcddb;
              *(here->B4SOIDPspPtr +1) += m * xcdsb;

              *(here->B4SOISPgPtr +1) += m * xcsgb;
              *(here->B4SOISPdpPtr +1) += m * xcsdb;
              *(here->B4SOISPspPtr +1) += m * xcssb;

/* v3.1 */
              if (here->B4SOIsoiMod != 2) /* v3.2 */
              {
                 *(here->B4SOIBePtr +1) += m * xcbeb;
                 *(here->B4SOIBgPtr +1) += m * xcbgb;
                 *(here->B4SOIBdpPtr +1) += m * xcbdb;
                 *(here->B4SOIBspPtr +1) += m * xcbsb;

                 *(here->B4SOIEbPtr  +1) -= m * (xcegb + xceeb + xcedb + xcesb);
                 *(here->B4SOIGbPtr +1) -= m * (xcggb + xcgdb + xcgsb + xcgeb);
/*                 *(here->B4SOIDPbPtr +1) -= xcdgb + xcddb + xcdsb + xcdeb; 
                                          + xcdgmb; */

                 *(here->B4SOIDPbPtr +1) -= m * -xcdbb; /* v4.0 */

/*                 *(here->B4SOISPbPtr +1) -= xcsgb + xcsdb + xcssb + xcseb
                                            + xcsgmb; */
                 *(here->B4SOISPbPtr +1) -= m * -xcsbb; /* v4.0 */
                 *(here->B4SOIBbPtr +1) -= m * (xcbgb + xcbdb + xcbsb + xcbeb);
              }
/* v3.1 */


              if (selfheat)
              {
                 *(here->B4SOITemptempPtr + 1) += m * xcTt;
                 *(here->B4SOIDPtempPtr + 1) += m * xcdT;
                 *(here->B4SOISPtempPtr + 1) += m * xcsT;
                 *(here->B4SOIBtempPtr + 1) += m * xcbT;
                 *(here->B4SOIEtempPtr + 1) += m * xceT;
                 *(here->B4SOIGtempPtr + 1) += m * xcgT;
              }
                                                               
 
/* v3.0 */
              if (here->B4SOIsoiMod != 0) /* v3.2 */
              {
                 *(here->B4SOIDPePtr) += m * (Gme + gddpe);
                 *(here->B4SOISPePtr) += m * (gsspe - Gme);

                 if (here->B4SOIsoiMod != 2) /* v3.2 */
                 {
                    *(here->B4SOIGePtr) += m * gige;
                    *(here->B4SOIBePtr) -= m * gige;
                 }
              }

              *(here->B4SOIEePtr) += 0.0;

              *(here->B4SOIDPgPtr) += m * (Gm + gddpg - gIdtotg -gdtotg); /* v4.0 */
              *(here->B4SOIDPdpPtr) += m * (gdpr + gds + gddpdp + RevSum - gIdtotd
                                     - gdtotd); /* v4.0 */

              *(here->B4SOIDPspPtr) -= m * (gds + FwdSum - gddpsp + gIdtots 
                                     + gdtots); /* v4.0 */
              *(here->B4SOIDPdPtr) -= m * (gdpr + gdtot);

              *(here->B4SOISPgPtr) -= m * (Gm - gsspg + gIstotg + gstotg); /* v4.0 */
              *(here->B4SOISPdpPtr) -= m * (gds + RevSum - gsspdp + gIstotd
                                     + gstotd); /* v4.0 */
              *(here->B4SOISPspPtr) += m * (gspr + gds + FwdSum + gsspsp - gIstots
                                     - gstots); /* v4.0 */
              *(here->B4SOISPsPtr) -= m * (gspr + gstot);


/* v3.1 */
              if (here->B4SOIsoiMod != 2) /* v3.2 */
              {
                 *(here->B4SOIBePtr) += m * gbbe; /* v3.0 */
                 *(here->B4SOIBgPtr)  += m * (gbbg - gigg); /* v3.1 bug fix */
                 *(here->B4SOIBdpPtr) += m * (gbbdp - gigd); /* v3.1 bug fix */
                 *(here->B4SOIBspPtr) += m * (gbbsp - gigs); /* v3.1 bug fix */
                 *(here->B4SOIBbPtr) += m * (gbbb - gigb); /* v3.1 bug fix */
                 *(here->B4SOISPbPtr) -= m * (Gmbs - gsspb + gIstotb + gstotb); 
                                        /* v4.0 */
                 *(here->B4SOIDPbPtr) -= m * ((-gddpb - Gmbs) + gIdtotb + gdtotb); 
                                        /* v4.0 */
              }
/* v3.1 */



              if (selfheat)
              {
                 *(here->B4SOIDPtempPtr) += m * (GmT + gddpT);
                 *(here->B4SOISPtempPtr) += m * (-GmT + gsspT);
                 *(here->B4SOIBtempPtr) += m * (gbbT - gigT); /* v3.1 bug fix */
                 *(here->B4SOIGtempPtr) += m * gigT; /* v3.1 bug fix */
                 *(here->B4SOITemptempPtr) += m * (gTtt + 1/here->pParam->B4SOIrth);
                 *(here->B4SOITempgPtr) += m * gTtg;
                 *(here->B4SOITempbPtr) += m * gTtb;
                 *(here->B4SOITempdpPtr) += m * gTtdp;
                 *(here->B4SOITempspPtr) += m * gTtsp;

/* v3.0 */       
                 if (here->B4SOIsoiMod != 0) /* v3.2 */
                    *(here->B4SOITempePtr) += m * gTte;
              }


              *(here->B4SOIDdPtr) += m * (gdpr + gdtot);
              *(here->B4SOIDdpPtr) -= m * (gdpr - gdtotd);
              *(here->B4SOISsPtr) += m * (gspr + gstot);
              *(here->B4SOISspPtr) -= m * (gspr -gstots);


              if (here->B4SOIbodyMod == 1) {
                 (*(here->B4SOIBpPtr) -= m * gppp);
                 (*(here->B4SOIPbPtr) += m * gppb);
                 (*(here->B4SOIPpPtr) += m * gppp);
              }

              /* v4.1  Ig_agbcp2 stamping */
              (*(here->B4SOIGgPtr) += gigpg);
              if (here->B4SOIbodyMod == 1)  {
                  (*(here->B4SOIPpPtr) -= m * gigpp);
                  (*(here->B4SOIPgPtr) -= m * gigpg);
                  (*(here->B4SOIGpPtr) += m * gigpp);
              }
              else if(here->B4SOIbodyMod == 2)
              {
                  (*(here->B4SOIBbPtr) -= m * gigpp);
                  (*(here->B4SOIBgPtr) -= m * gigpg);
                  (*(here->B4SOIGbPtr) += m * gigpp);
              }



              /* v4.0 */
              if (here->B4SOIrbodyMod)
              {
                      (*(here->B4SOIDPdbPtr + 1) -= m * xcjdbdp);
                      (*(here->B4SOIDPdbPtr) -= m * here->B4SOIGGjdb);
                (*(here->B4SOISPsbPtr + 1) -= m * xcjsbsp);
                (*(here->B4SOISPsbPtr) -= m * here->B4SOIGGjsb);

                (*(here->B4SOIDBdpPtr + 1) -= m * xcjdbdp);
                (*(here->B4SOIDBdpPtr) -= m * here->B4SOIGGjdb);
                (*(here->B4SOIDBdbPtr + 1) += m * xcjdbdp);
                (*(here->B4SOIDBdbPtr) += m * (here->B4SOIGGjdb
                                        + here->B4SOIgrbdb));

                (*(here->B4SOIDBbPtr) -= m * here->B4SOIgrbdb);
                (*(here->B4SOISBbPtr) -= m * here->B4SOIgrbsb);

                (*(here->B4SOISBspPtr + 1) -= m * xcjsbsp);
                (*(here->B4SOISBspPtr) -= m * here->B4SOIGGjsb);
                (*(here->B4SOISBsbPtr + 1) += m * xcjsbsp);
                (*(here->B4SOISBsbPtr) += m * (here->B4SOIGGjsb 
                                        + here->B4SOIgrbsb));

                (*(here->B4SOIBdbPtr) -= m * here->B4SOIgrbdb);
                (*(here->B4SOIBsbPtr) -= m * here->B4SOIgrbsb);
                (*(here->B4SOIBbPtr) += m * (here->B4SOIgrbsb
                                        + here->B4SOIgrbdb));

              } 

              if (here->B4SOIdebugMod != 0)
              {
                      *(here->B4SOIVbsPtr) += 1;
                      *(here->B4SOIIdsPtr) += 1;
                      *(here->B4SOIIcPtr) += 1;
                      *(here->B4SOIIbsPtr) += 1;
                      *(here->B4SOIIbdPtr) += 1;
                      *(here->B4SOIIiiPtr) += 1;
                      *(here->B4SOIIgidlPtr) += 1;
                      *(here->B4SOIItunPtr) += 1;
                      *(here->B4SOIIbpPtr) += 1;
                      *(here->B4SOICbgPtr) += 1;
                      *(here->B4SOICbbPtr) += 1;
                      *(here->B4SOICbdPtr) += 1;
                      *(here->B4SOIQbfPtr) += 1;
                      *(here->B4SOIQjsPtr) += 1;
                      *(here->B4SOIQjdPtr) += 1;

              }

        }
    }
    return(OK);
}

