/* acload.f -- translated by f2c (version of 3 February 1990  3:36:42).
   You must link the resulting object file with the libraries:
	-lF77 -lI77 -lm -lc   (in that order)
*/

#include "f2c.h"

/* Common Block Declarations */

struct {
    integer ielmnt, isbckt, nsbckt, iunsat, nunsat, itemps, numtem, isens, 
	    nsens, ifour, nfour, ifield, icode, idelim, icolum, insize, 
	    junode, lsbkpt, numbkp, iorder, jmnode, iur, iuc, ilc, ilr, 
	    numoff, isr, nmoffc, iseq, iseq1, neqn, nodevs, ndiag, iswap, 
	    iequa, macins, lvnim1, lx0, lvn, lynl, lyu, lyl, lx1, lx2, lx3, 
	    lx4, lx5, lx6, lx7, ld0, ld1, ltd, imynl, imvn, lcvn, nsnod, 
	    nsmat, nsval, icnod, icmat, icval, loutpt, lpol, lzer, irswpf, 
	    irswpr, icswpf, icswpr, irpt, jcpt, irowno, jcolno, nttbr, nttar, 
	    lvntmp;
} tabinf_;

#define tabinf_1 tabinf_

struct {
    integer locate[50], jelcnt[50], nunods, ncnods, numnod, nstop, nut, nlt, 
	    nxtrm, ndist, ntlin, ibr, numvs, numalt, numcyc;
} cirdat_;

#define cirdat_1 cirdat_

struct {
    doublereal omega, time, delta, delold[7], ag[7], vt, xni, egfet, xmu, 
	    sfactr;
    integer mode, modedc, icalc, initf, method, iord, maxord, noncon, iterno, 
	    itemno, nosolv, modac, ipiv, ivmflg, ipostp, iscrch, iofile;
} status_;

#define status_1 status_

struct {
    doublereal fstart, fstop, fincr, skw2, refprl, spw2;
    integer jacflg, idfreq, inoise, nosprt, nosout, nosin, idist, idprt;
} ac_;

#define ac_1 ac_

struct {
    doublereal value[200000];
} blank_;

#define blank_1 blank_

/*<       subroutine acload >*/
/* Subroutine */ int acload_()
{
    /* System generated locals */
    integer i_1;

    /* Builtin functions */
    double sin(), cos();

    /* Local variables */
    static doublereal cggb, area, cbsb, cgsb, xcgb, xcgd;
    static integer ndim;
    static doublereal gmbs;
    static integer locm;
    static doublereal gcpr, gdpr;
    static integer lmat, locv, loct, locy;
    static doublereal xceq, gspr, gepr;
    static integer iptr;
    static doublereal xcpi, xcmu, xcbx, xccs, xnrm, xrev, xqco, xcgs, rval, 
	    xval;
    static integer node1, node2;
    extern /* Subroutine */ int zero8_();
    static integer i;
    static doublereal capbd, xcddb, xcbgb, xcdgb, capbs, xcmcb, xcggb, xcgdb, 
	    xcdsb, xcsdb, xcsgb, xcgsb, xcbdb, xcdxd, xcbsb, xcgxd;
    static integer locyi, itype;
    static doublereal xcgxs, xcdxs, xcsxd, xcssb, xcsxs, z0, y0, gm, go, td, 
	    gx, xl;
#define nodplc ((integer *)&blank_1)
#define cvalue ((complex *)&blank_1)
    static doublereal devmod, xw, covlgs, covlgd, covlgb, gbd, gcd, gcg, ggd, 
	    arg, gbs;
    static integer loc;
    static doublereal geq, gpi, gds, val, ggs, xgd, gcs, gmu, xbd, xgm, xqc, 
	    xgs, xbs, cbdb, cbgb, cgdb;

/*<       implicit double precision (a-h,o-z) >*/

/*     this routine zeroes-out and then loads the complex coefficient */
/*     matrix */

/* spice version 2g.6  sccsid=tabinf 3/15/83 */
/*<       common /tabinf/ ielmnt,isbckt,nsbckt,iunsat,nunsat,itemps,numtem, >*/
/*<      1   isens,nsens,ifour,nfour,ifield,icode,idelim,icolum,insize, >*/
/*<      2   junode,lsbkpt,numbkp,iorder,jmnode,iur,iuc,ilc,ilr,numoff,isr, >*/
/*<      3   nmoffc,iseq,iseq1,neqn,nodevs,ndiag,iswap,iequa,macins,lvnim1, >*/
/*<      4   lx0,lvn,lynl,lyu,lyl,lx1,lx2,lx3,lx4,lx5,lx6,lx7,ld0,ld1,ltd, >*/
/*<      5   imynl,imvn,lcvn,nsnod,nsmat,nsval,icnod,icmat,icval, >*/
/*<      6   loutpt,lpol,lzer,irswpf,irswpr,icswpf,icswpr,irpt,jcpt, >*/
/*<      7   irowno,jcolno,nttbr,nttar,lvntmp >*/
/* spice version 2g.6  sccsid=cirdat 3/15/83 */
/*<       common /cirdat/ locate(50),jelcnt(50),nunods,ncnods,numnod,nstop, >*/
/*<      1   nut,nlt,nxtrm,ndist,ntlin,ibr,numvs,numalt,numcyc >*/
/* spice version 2g.6  sccsid=status 3/15/83 */
/*<       common /status/ omega,time,delta,delold(7),ag(7),vt,xni,egfet, >*/
/*<      1   xmu,sfactr,mode,modedc,icalc,initf,method,iord,maxord,noncon, >*/
/*<      2   iterno,itemno,nosolv,modac,ipiv,ivmflg,ipostp,iscrch,iofile >*/
/* spice version 2g.6  sccsid=ac 3/15/83 */
/*<       common /ac/ fstart,fstop,fincr,skw2,refprl,spw2,jacflg,idfreq, >*/
/*<      1   inoise,nosprt,nosout,nosin,idist,idprt >*/
/* spice version 2g.6  sccsid=blank 3/15/83 */
/*<       common /blank/ value(200000) >*/
/*<       integer nodplc(64) >*/
/*<       complex cvalue(32) >*/
/*<       equivalence (value(1),nodplc(1),cvalue(1)) >*/


/*<       complex cval >*/

/*  zero y matrix and current vector */

/*<       call zero8(value(lvn+1),nstop+nttbr) >*/
    i_1 = cirdat_1.nstop + tabinf_1.nttbr;
    zero8_(&blank_1.value[tabinf_1.lvn], &i_1);
/*<       call zero8(value(imvn+1),nstop+nttbr) >*/
    i_1 = cirdat_1.nstop + tabinf_1.nttbr;
    zero8_(&blank_1.value[tabinf_1.imvn], &i_1);

/*  resistors */

/*<       loc=locate(1) >*/
    loc = cirdat_1.locate[0];
/*<    20 if ((loc.eq.0).or.(nodplc(loc+8).ne.0)) go to 30 >*/
L20:
    if (loc == 0 || nodplc[loc + 7] != 0) {
	goto L30;
    }
/*<       locv=nodplc(loc+1) >*/
    locv = nodplc[loc];
/*<       val=value(locv+1) >*/
    val = blank_1.value[locv];
/*<       locy=lynl+nodplc(loc+6) >*/
    locy = tabinf_1.lynl + nodplc[loc + 5];
/*<       value(locy)=value(locy)+val >*/
    blank_1.value[locy - 1] += val;
/*<       locy=lynl+nodplc(loc+7) >*/
    locy = tabinf_1.lynl + nodplc[loc + 6];
/*<       value(locy)=value(locy)+val >*/
    blank_1.value[locy - 1] += val;
/*<       locy=lynl+nodplc(loc+4) >*/
    locy = tabinf_1.lynl + nodplc[loc + 3];
/*<       value(locy)=value(locy)-val >*/
    blank_1.value[locy - 1] -= val;
/*<       locy=lynl+nodplc(loc+5) >*/
    locy = tabinf_1.lynl + nodplc[loc + 4];
/*<       value(locy)=value(locy)-val >*/
    blank_1.value[locy - 1] -= val;
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 20 >*/
    goto L20;

/*  capacitors */

/*<    30 loc=locate(2) >*/
L30:
    loc = cirdat_1.locate[1];
/*<    40 if ((loc.eq.0).or.(nodplc(loc+12).ne.0)) go to 50 >*/
L40:
    if (loc == 0 || nodplc[loc + 11] != 0) {
	goto L50;
    }
/*<       locv=nodplc(loc+1) >*/
    locv = nodplc[loc];
/*<       val=omega*value(locv+1) >*/
    val = status_1.omega * blank_1.value[locv];
/*<       locyi=imynl+nodplc(loc+10) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 9];
/*<       value(locyi)=value(locyi)+val >*/
    blank_1.value[locyi - 1] += val;
/*<       locyi=imynl+nodplc(loc+11) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 10];
/*<       value(locyi)=value(locyi)+val >*/
    blank_1.value[locyi - 1] += val;
/*<       locyi=imynl+nodplc(loc+5) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 4];
/*<       value(locyi)=value(locyi)-val >*/
    blank_1.value[locyi - 1] -= val;
/*<       locyi=imynl+nodplc(loc+6) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 5];
/*<       value(locyi)=value(locyi)-val >*/
    blank_1.value[locyi - 1] -= val;
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 40 >*/
    goto L40;

/*  inductors */

/*<    50 loc=locate(3) >*/
L50:
    loc = cirdat_1.locate[2];
/*<    60 if ((loc.eq.0).or.(nodplc(loc+14).ne.0)) go to 70 >*/
L60:
    if (loc == 0 || nodplc[loc + 13] != 0) {
	goto L70;
    }
/*<       locv=nodplc(loc+1) >*/
    locv = nodplc[loc];
/*<       val=omega*value(locv+1) >*/
    val = status_1.omega * blank_1.value[locv];
/*<       locyi=imynl+nodplc(loc+13) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 12];
/*<       locy=lynl+nodplc(loc+13) >*/
    locy = tabinf_1.lynl + nodplc[loc + 12];
/*<       value(locy)=0.0d0 >*/
    blank_1.value[locy - 1] = 0.;
/*<       value(locyi)=-val >*/
    blank_1.value[locyi - 1] = -val;
/*<       locy=lynl+nodplc(loc+6) >*/
    locy = tabinf_1.lynl + nodplc[loc + 5];
/*<       locyi=imynl+nodplc(loc+6) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 5];
/*<       value(locy)=1.0d0 >*/
    blank_1.value[locy - 1] = 1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+7) >*/
    locy = tabinf_1.lynl + nodplc[loc + 6];
/*<       locyi=imynl+nodplc(loc+7) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 6];
/*<       value(locy)=-1.0d0 >*/
    blank_1.value[locy - 1] = -1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+8) >*/
    locy = tabinf_1.lynl + nodplc[loc + 7];
/*<       locyi=imynl+nodplc(loc+8) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 7];
/*<       value(locy)=1.0d0 >*/
    blank_1.value[locy - 1] = 1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+9) >*/
    locy = tabinf_1.lynl + nodplc[loc + 8];
/*<       locyi=imynl+nodplc(loc+9) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 8];
/*<       value(locy)=-1.0d0 >*/
    blank_1.value[locy - 1] = -1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 60 >*/
    goto L60;

/*  mutual inductors */

/*<    70 loc=locate(4) >*/
L70:
    loc = cirdat_1.locate[3];
/*<    80 if ((loc.eq.0).or.(nodplc(loc+6).ne.0)) go to 90 >*/
L80:
    if (loc == 0 || nodplc[loc + 5] != 0) {
	goto L90;
    }
/*<       locv=nodplc(loc+1) >*/
    locv = nodplc[loc];
/*<       val=omega*value(locv+1) >*/
    val = status_1.omega * blank_1.value[locv];
/*<       locy=lynl+nodplc(loc+4) >*/
    locy = tabinf_1.lynl + nodplc[loc + 3];
/*<       locyi=imynl+nodplc(loc+4) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 3];
/*<       value(locy)=0.0d0 >*/
    blank_1.value[locy - 1] = 0.;
/*<       value(locyi)=-val >*/
    blank_1.value[locyi - 1] = -val;
/*<       locy=lynl+nodplc(loc+5) >*/
    locy = tabinf_1.lynl + nodplc[loc + 4];
/*<       locyi=imynl+nodplc(loc+5) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 4];
/*<       value(locy)=0.0d0 >*/
    blank_1.value[locy - 1] = 0.;
/*<       value(locyi)=-val >*/
    blank_1.value[locyi - 1] = -val;
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 80 >*/
    goto L80;

/*  nonlinear voltage controlled current sources */

/*<    90 loc=locate(5) >*/
L90:
    loc = cirdat_1.locate[4];
/*<    95 if ((loc.eq.0).or.(nodplc(loc+13).ne.0)) go to 100 >*/
L95:
    if (loc == 0 || nodplc[loc + 12] != 0) {
	goto L100;
    }
/*<       ndim=nodplc(loc+4) >*/
    ndim = nodplc[loc + 3];
/*<       lmat=nodplc(loc+7) >*/
    lmat = nodplc[loc + 6];
/*<       loct=lx0+nodplc(loc+12)+2 >*/
    loct = tabinf_1.lx0 + nodplc[loc + 11] + 2;
/*<       do 97 i=1,ndim >*/
    i_1 = ndim;
    for (i = 1; i <= i_1; ++i) {
/*<       val=value(loct) >*/
	val = blank_1.value[loct - 1];
/*<       loct=loct+2 >*/
	loct += 2;
/*<       locy=lynl+nodplc(lmat+1) >*/
	locy = tabinf_1.lynl + nodplc[lmat];
/*<       value(locy)=value(locy)+val >*/
	blank_1.value[locy - 1] += val;
/*<       locy=lynl+nodplc(lmat+2) >*/
	locy = tabinf_1.lynl + nodplc[lmat + 1];
/*<       value(locy)=value(locy)-val >*/
	blank_1.value[locy - 1] -= val;
/*<       locy=lynl+nodplc(lmat+3) >*/
	locy = tabinf_1.lynl + nodplc[lmat + 2];
/*<       value(locy)=value(locy)-val >*/
	blank_1.value[locy - 1] -= val;
/*<       locy=lynl+nodplc(lmat+4) >*/
	locy = tabinf_1.lynl + nodplc[lmat + 3];
/*<       value(locy)=value(locy)+val >*/
	blank_1.value[locy - 1] += val;
/*<       lmat=lmat+4 >*/
	lmat += 4;
/*<    97 continue >*/
/* L97: */
    }
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 95 >*/
    goto L95;

/*  nonlinear voltage controlled voltage sources */

/*<   100 loc=locate(6) >*/
L100:
    loc = cirdat_1.locate[5];
/*<   105 if ((loc.eq.0).or.(nodplc(loc+14).ne.0)) go to 110 >*/
L105:
    if (loc == 0 || nodplc[loc + 13] != 0) {
	goto L110;
    }
/*<       ndim=nodplc(loc+4) >*/
    ndim = nodplc[loc + 3];
/*<       lmat=nodplc(loc+8) >*/
    lmat = nodplc[loc + 7];
/*<       loct=lx0+nodplc(loc+13)+3 >*/
    loct = tabinf_1.lx0 + nodplc[loc + 12] + 3;
/*<       locy=lynl+nodplc(lmat+1) >*/
    locy = tabinf_1.lynl + nodplc[lmat];
/*<       locyi=imynl+nodplc(lmat+1) >*/
    locyi = tabinf_1.imynl + nodplc[lmat];
/*<       value(locy)=+1.0d0 >*/
    blank_1.value[locy - 1] = 1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(lmat+2) >*/
    locy = tabinf_1.lynl + nodplc[lmat + 1];
/*<       locyi=imynl+nodplc(lmat+2) >*/
    locyi = tabinf_1.imynl + nodplc[lmat + 1];
/*<       value(locy)=-1.0d0 >*/
    blank_1.value[locy - 1] = -1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(lmat+3) >*/
    locy = tabinf_1.lynl + nodplc[lmat + 2];
/*<       locyi=imynl+nodplc(lmat+3) >*/
    locyi = tabinf_1.imynl + nodplc[lmat + 2];
/*<       value(locy)=+1.0d0 >*/
    blank_1.value[locy - 1] = 1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(lmat+4) >*/
    locy = tabinf_1.lynl + nodplc[lmat + 3];
/*<       locyi=imynl+nodplc(lmat+4) >*/
    locyi = tabinf_1.imynl + nodplc[lmat + 3];
/*<       value(locy)=-1.0d0 >*/
    blank_1.value[locy - 1] = -1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       lmat=lmat+4 >*/
    lmat += 4;
/*<       do 107 i=1,ndim >*/
    i_1 = ndim;
    for (i = 1; i <= i_1; ++i) {
/*<       val=value(loct) >*/
	val = blank_1.value[loct - 1];
/*<       loct=loct+2 >*/
	loct += 2;
/*<       locy=lynl+nodplc(lmat+1) >*/
	locy = tabinf_1.lynl + nodplc[lmat];
/*<       value(locy)=value(locy)-val >*/
	blank_1.value[locy - 1] -= val;
/*<       locy=lynl+nodplc(lmat+2) >*/
	locy = tabinf_1.lynl + nodplc[lmat + 1];
/*<       value(locy)=value(locy)+val >*/
	blank_1.value[locy - 1] += val;
/*<       lmat=lmat+2 >*/
	lmat += 2;
/*<   107 continue >*/
/* L107: */
    }
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 105 >*/
    goto L105;

/*  nonlinear current controlled current sources */

/*<   110 loc=locate(7) >*/
L110:
    loc = cirdat_1.locate[6];
/*<   115 if ((loc.eq.0).or.(nodplc(loc+13).ne.0)) go to 120 >*/
L115:
    if (loc == 0 || nodplc[loc + 12] != 0) {
	goto L120;
    }
/*<       ndim=nodplc(loc+4) >*/
    ndim = nodplc[loc + 3];
/*<       lmat=nodplc(loc+7) >*/
    lmat = nodplc[loc + 6];
/*<       loct=lx0+nodplc(loc+12)+2 >*/
    loct = tabinf_1.lx0 + nodplc[loc + 11] + 2;
/*<       do 117 i=1,ndim >*/
    i_1 = ndim;
    for (i = 1; i <= i_1; ++i) {
/*<       val=value(loct) >*/
	val = blank_1.value[loct - 1];
/*<       loct=loct+2 >*/
	loct += 2;
/*<       locy=lynl+nodplc(lmat+1) >*/
	locy = tabinf_1.lynl + nodplc[lmat];
/*<       locyi=imynl+nodplc(lmat+1) >*/
	locyi = tabinf_1.imynl + nodplc[lmat];
/*<       value(locy)=+val >*/
	blank_1.value[locy - 1] = val;
/*<       value(locyi)=0.0d0 >*/
	blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(lmat+2) >*/
	locy = tabinf_1.lynl + nodplc[lmat + 1];
/*<       locyi=imynl+nodplc(lmat+2) >*/
	locyi = tabinf_1.imynl + nodplc[lmat + 1];
/*<       value(locy)=-val >*/
	blank_1.value[locy - 1] = -val;
/*<       value(locyi)=0.0d0 >*/
	blank_1.value[locyi - 1] = 0.;
/*<       lmat=lmat+2 >*/
	lmat += 2;
/*<   117 continue >*/
/* L117: */
    }
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 115 >*/
    goto L115;

/*  nonlinear current controlled voltage sources */

/*<   120 loc=locate(8) >*/
L120:
    loc = cirdat_1.locate[7];
/*<   125 if ((loc.eq.0).or.(nodplc(loc+14).ne.0)) go to 140 >*/
L125:
    if (loc == 0 || nodplc[loc + 13] != 0) {
	goto L140;
    }
/*<       ndim=nodplc(loc+4) >*/
    ndim = nodplc[loc + 3];
/*<       lmat=nodplc(loc+8) >*/
    lmat = nodplc[loc + 7];
/*<       loct=lx0+nodplc(loc+13)+3 >*/
    loct = tabinf_1.lx0 + nodplc[loc + 12] + 3;
/*<       locy=lynl+nodplc(lmat+1) >*/
    locy = tabinf_1.lynl + nodplc[lmat];
/*<       locyi=imynl+nodplc(lmat+1) >*/
    locyi = tabinf_1.imynl + nodplc[lmat];
/*<       value(locy)=+1.0d0 >*/
    blank_1.value[locy - 1] = 1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(lmat+2) >*/
    locy = tabinf_1.lynl + nodplc[lmat + 1];
/*<       locyi=imynl+nodplc(lmat+2) >*/
    locyi = tabinf_1.imynl + nodplc[lmat + 1];
/*<       value(locy)=-1.0d0 >*/
    blank_1.value[locy - 1] = -1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(lmat+3) >*/
    locy = tabinf_1.lynl + nodplc[lmat + 2];
/*<       locyi=imynl+nodplc(lmat+3) >*/
    locyi = tabinf_1.imynl + nodplc[lmat + 2];
/*<       value(locy)=+1.0d0 >*/
    blank_1.value[locy - 1] = 1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(lmat+4) >*/
    locy = tabinf_1.lynl + nodplc[lmat + 3];
/*<       locyi=imynl+nodplc(lmat+4) >*/
    locyi = tabinf_1.imynl + nodplc[lmat + 3];
/*<       value(locy)=-1.0d0 >*/
    blank_1.value[locy - 1] = -1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       lmat=lmat+4 >*/
    lmat += 4;
/*<       do 127 i=1,ndim >*/
    i_1 = ndim;
    for (i = 1; i <= i_1; ++i) {
/*<       val=value(loct) >*/
	val = blank_1.value[loct - 1];
/*<       loct=loct+2 >*/
	loct += 2;
/*<       locy=lynl+nodplc(lmat+i) >*/
	locy = tabinf_1.lynl + nodplc[lmat + i - 1];
/*<       value(locy)=value(locy)-val >*/
	blank_1.value[locy - 1] -= val;
/*<   127 continue >*/
/* L127: */
    }
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 125 >*/
    goto L125;

/*  voltage sources */

/*<   140 loc=locate(9) >*/
L140:
    loc = cirdat_1.locate[8];
/*<   150 if ((loc.eq.0).or.(nodplc(loc+11).ne.0)) go to 160 >*/
L150:
    if (loc == 0 || nodplc[loc + 10] != 0) {
	goto L160;
    }
/*<       locv=nodplc(loc+1) >*/
    locv = nodplc[loc];
/*<       iptr=nodplc(loc+6) >*/
    iptr = nodplc[loc + 5];
/*<       value(lvn+iptr)=value(locv+2) >*/
    blank_1.value[tabinf_1.lvn + iptr - 1] = blank_1.value[locv + 1];
/*<       value(imvn+iptr)=value(locv+3) >*/
    blank_1.value[tabinf_1.imvn + iptr - 1] = blank_1.value[locv + 2];
/*<       locy=lynl+nodplc(loc+7) >*/
    locy = tabinf_1.lynl + nodplc[loc + 6];
/*<       value(locy)=value(locy)+1.0d0 >*/
    blank_1.value[locy - 1] += 1.;
/*<       locy=lynl+nodplc(loc+8) >*/
    locy = tabinf_1.lynl + nodplc[loc + 7];
/*<       value(locy)=value(locy)-1.0d0 >*/
    blank_1.value[locy - 1] += -1.;
/*<       locy=lynl+nodplc(loc+9) >*/
    locy = tabinf_1.lynl + nodplc[loc + 8];
/*<       value(locy)=value(locy)+1.0d0 >*/
    blank_1.value[locy - 1] += 1.;
/*<       locy=lynl+nodplc(loc+10) >*/
    locy = tabinf_1.lynl + nodplc[loc + 9];
/*<       value(locy)=value(locy)-1.0d0 >*/
    blank_1.value[locy - 1] += -1.;
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 150 >*/
    goto L150;

/*  current sources */

/*<   160 loc=locate(10) >*/
L160:
    loc = cirdat_1.locate[9];
/*<   170 if ((loc.eq.0).or.(nodplc(loc+6).ne.0)) go to 200 >*/
L170:
    if (loc == 0 || nodplc[loc + 5] != 0) {
	goto L200;
    }
/*<       locv=nodplc(loc+1) >*/
    locv = nodplc[loc];
/*<       node1=nodplc(loc+2) >*/
    node1 = nodplc[loc + 1];
/*<       node2=nodplc(loc+3) >*/
    node2 = nodplc[loc + 2];
/*<       value(lvn+node1)=value(lvn+node1)-value(locv+2) >*/
    blank_1.value[tabinf_1.lvn + node1 - 1] -= blank_1.value[locv + 1];
/*<       value(imvn+node1)=value(imvn+node1)-value(locv+3) >*/
    blank_1.value[tabinf_1.imvn + node1 - 1] -= blank_1.value[locv + 2];
/*<       value(lvn+node2)=value(lvn+node2)+value(locv+2) >*/
    blank_1.value[tabinf_1.lvn + node2 - 1] += blank_1.value[locv + 1];
/*<       value(imvn+node2)=value(imvn+node2)+value(locv+3) >*/
    blank_1.value[tabinf_1.imvn + node2 - 1] += blank_1.value[locv + 2];
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 170 >*/
    goto L170;

/*  diodes */

/*<   200 loc=locate(11) >*/
L200:
    loc = cirdat_1.locate[10];
/*<   210 if ((loc.eq.0).or.(nodplc(loc+16).ne.0)) go to 250 >*/
L210:
    if (loc == 0 || nodplc[loc + 15] != 0) {
	goto L250;
    }
/*<       locv=nodplc(loc+1) >*/
    locv = nodplc[loc];
/*<       area=value(locv+1) >*/
    area = blank_1.value[locv];
/*<       locm=nodplc(loc+5) >*/
    locm = nodplc[loc + 4];
/*<       locm=nodplc(locm+1) >*/
    locm = nodplc[locm];
/*<       loct=lx0+nodplc(loc+11) >*/
    loct = tabinf_1.lx0 + nodplc[loc + 10];
/*<       gspr=value(locm+2)*area >*/
    gspr = blank_1.value[locm + 1] * area;
/*<       geq=value(loct+2) >*/
    geq = blank_1.value[loct + 1];
/*<       xceq=value(loct+4)*omega >*/
    xceq = blank_1.value[loct + 3] * status_1.omega;
/*<       locy=lynl+nodplc(loc+13) >*/
    locy = tabinf_1.lynl + nodplc[loc + 12];
/*<       value(locy)=value(locy)+gspr >*/
    blank_1.value[locy - 1] += gspr;
/*<       locy=lynl+nodplc(loc+14) >*/
    locy = tabinf_1.lynl + nodplc[loc + 13];
/*<       locyi=imynl+nodplc(loc+14) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 13];
/*<       value(locy)=value(locy)+geq >*/
    blank_1.value[locy - 1] += geq;
/*<       value(locyi)=value(locyi)+xceq >*/
    blank_1.value[locyi - 1] += xceq;
/*<       locy=lynl+nodplc(loc+15) >*/
    locy = tabinf_1.lynl + nodplc[loc + 14];
/*<       locyi=imynl+nodplc(loc+15) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 14];
/*<       value(locy)=value(locy)+geq+gspr >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] + geq + gspr;
/*<       value(locyi)=value(locyi)+xceq >*/
    blank_1.value[locyi - 1] += xceq;
/*<       locy=lynl+nodplc(loc+7) >*/
    locy = tabinf_1.lynl + nodplc[loc + 6];
/*<       value(locy)=value(locy)-gspr >*/
    blank_1.value[locy - 1] -= gspr;
/*<       locy=lynl+nodplc(loc+8) >*/
    locy = tabinf_1.lynl + nodplc[loc + 7];
/*<       locyi=imynl+nodplc(loc+8) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 7];
/*<       value(locy)=value(locy)-geq >*/
    blank_1.value[locy - 1] -= geq;
/*<       value(locyi)=value(locyi)-xceq >*/
    blank_1.value[locyi - 1] -= xceq;
/*<       locy=lynl+nodplc(loc+9) >*/
    locy = tabinf_1.lynl + nodplc[loc + 8];
/*<       value(locy)=value(locy)-gspr >*/
    blank_1.value[locy - 1] -= gspr;
/*<       locy=lynl+nodplc(loc+10) >*/
    locy = tabinf_1.lynl + nodplc[loc + 9];
/*<       locyi=imynl+nodplc(loc+10) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 9];
/*<       value(locy)=value(locy)-geq >*/
    blank_1.value[locy - 1] -= geq;
/*<       value(locyi)=value(locyi)-xceq >*/
    blank_1.value[locyi - 1] -= xceq;
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 210 >*/
    goto L210;

/*  bjts */

/*<   250 loc=locate(12) >*/
L250:
    loc = cirdat_1.locate[11];
/*<   260 if ((loc.eq.0).or.(nodplc(loc+36).ne.0)) go to 300 >*/
L260:
    if (loc == 0 || nodplc[loc + 35] != 0) {
	goto L300;
    }
/*<       locv=nodplc(loc+1) >*/
    locv = nodplc[loc];
/*<       area=value(locv+1) >*/
    area = blank_1.value[locv];
/*<       locm=nodplc(loc+8) >*/
    locm = nodplc[loc + 7];
/*<       locm=nodplc(locm+1) >*/
    locm = nodplc[locm];
/*<       loct=lx0+nodplc(loc+22) >*/
    loct = tabinf_1.lx0 + nodplc[loc + 21];
/*<       gcpr=value(locm+20)*area >*/
    gcpr = blank_1.value[locm + 19] * area;
/*<       gepr=value(locm+19)*area >*/
    gepr = blank_1.value[locm + 18] * area;
/*<       gpi=value(loct+4) >*/
    gpi = blank_1.value[loct + 3];
/*<       gmu=value(loct+5) >*/
    gmu = blank_1.value[loct + 4];
/*<       gm=value(loct+6) >*/
    gm = blank_1.value[loct + 5];
/*<       go=value(loct+7) >*/
    go = blank_1.value[loct + 6];
/*<       xgm=0.0d0 >*/
    xgm = 0.;
/*<       td=value(locm+28) >*/
    td = blank_1.value[locm + 27];
/*<       if(td.eq.0.0d0) go to 270 >*/
    if (td == 0.) {
	goto L270;
    }
/*<       arg=td*omega >*/
    arg = td * status_1.omega;
/*<       gm=gm+go >*/
    gm += go;
/*<       xgm=-gm*dsin(arg) >*/
    xgm = -gm * sin(arg);
/*<       gm=gm*dcos(arg)-go >*/
    gm = gm * cos(arg) - go;
/*<   270 gx=value(loct+16) >*/
L270:
    gx = blank_1.value[loct + 15];
/*<       xcpi=value(loct+9)*omega >*/
    xcpi = blank_1.value[loct + 8] * status_1.omega;
/*<       xcmu=value(loct+11)*omega >*/
    xcmu = blank_1.value[loct + 10] * status_1.omega;
/*<       xcbx=value(loct+15)*omega >*/
    xcbx = blank_1.value[loct + 14] * status_1.omega;
/*<       xccs=value(loct+13)*omega >*/
    xccs = blank_1.value[loct + 12] * status_1.omega;
/*<       xcmcb=value(loct+17)*omega >*/
    xcmcb = blank_1.value[loct + 16] * status_1.omega;
/*<       locy=lynl+nodplc(loc+24) >*/
    locy = tabinf_1.lynl + nodplc[loc + 23];
/*<       value(locy)=value(locy)+gcpr >*/
    blank_1.value[locy - 1] += gcpr;
/*<       locy=lynl+nodplc(loc+25) >*/
    locy = tabinf_1.lynl + nodplc[loc + 24];
/*<       locyi=imynl+nodplc(loc+25) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 24];
/*<       value(locy)=value(locy)+gx >*/
    blank_1.value[locy - 1] += gx;
/*<       value(locyi)=value(locyi)+xcbx >*/
    blank_1.value[locyi - 1] += xcbx;
/*<       locy=lynl+nodplc(loc+26) >*/
    locy = tabinf_1.lynl + nodplc[loc + 25];
/*<       value(locy)=value(locy)+gepr >*/
    blank_1.value[locy - 1] += gepr;
/*<       locy=lynl+nodplc(loc+27) >*/
    locy = tabinf_1.lynl + nodplc[loc + 26];
/*<       locyi=imynl+nodplc(loc+27) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 26];
/*<       value(locy)=value(locy)+gmu+go+gcpr >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] + gmu + go + gcpr;
/*<       value(locyi)=value(locyi)+xcmu+xccs+xcbx >*/
    blank_1.value[locyi - 1] = blank_1.value[locyi - 1] + xcmu + xccs + xcbx;
/*<       locy=lynl+nodplc(loc+28) >*/
    locy = tabinf_1.lynl + nodplc[loc + 27];
/*<       locyi=imynl+nodplc(loc+28) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 27];
/*<       value(locy)=value(locy)+gx+gpi+gmu >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] + gx + gpi + gmu;
/*<       value(locyi)=value(locyi)+xcpi+xcmu+xcmcb >*/
    blank_1.value[locyi - 1] = blank_1.value[locyi - 1] + xcpi + xcmu + xcmcb;

/*<       locy=lynl+nodplc(loc+29) >*/
    locy = tabinf_1.lynl + nodplc[loc + 28];
/*<       locyi=imynl+nodplc(loc+29) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 28];
/*<       value(locy)=value(locy)+gpi+gepr+gm+go >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] + gpi + gepr + gm + go;
/*<       value(locyi)=value(locyi)+xcpi+xgm >*/
    blank_1.value[locyi - 1] = blank_1.value[locyi - 1] + xcpi + xgm;
/*<       locy=lynl+nodplc(loc+10) >*/
    locy = tabinf_1.lynl + nodplc[loc + 9];
/*<       value(locy)=value(locy)-gcpr >*/
    blank_1.value[locy - 1] -= gcpr;
/*<       locy=lynl+nodplc(loc+11) >*/
    locy = tabinf_1.lynl + nodplc[loc + 10];
/*<       value(locy)=value(locy)-gx >*/
    blank_1.value[locy - 1] -= gx;
/*<       locy=lynl+nodplc(loc+12) >*/
    locy = tabinf_1.lynl + nodplc[loc + 11];
/*<       value(locy)=value(locy)-gepr >*/
    blank_1.value[locy - 1] -= gepr;
/*<       locy=lynl+nodplc(loc+13) >*/
    locy = tabinf_1.lynl + nodplc[loc + 12];
/*<       value(locy)=value(locy)-gcpr >*/
    blank_1.value[locy - 1] -= gcpr;
/*<       locy=lynl+nodplc(loc+14) >*/
    locy = tabinf_1.lynl + nodplc[loc + 13];
/*<       locyi=imynl+nodplc(loc+14) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 13];
/*<       value(locy)=value(locy)-gmu+gm >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] - gmu + gm;
/*<       value(locyi)=value(locyi)-xcmu+xgm >*/
    blank_1.value[locyi - 1] = blank_1.value[locyi - 1] - xcmu + xgm;
/*<       locy=lynl+nodplc(loc+15) >*/
    locy = tabinf_1.lynl + nodplc[loc + 14];
/*<       locyi=imynl+nodplc(loc+15) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 14];
/*<       value(locy)=value(locy)-gm-go >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] - gm - go;
/*<       value(locyi)=value(locyi)-xgm >*/
    blank_1.value[locyi - 1] -= xgm;
/*<       locy=lynl+nodplc(loc+16) >*/
    locy = tabinf_1.lynl + nodplc[loc + 15];
/*<       value(locy)=value(locy)-gx >*/
    blank_1.value[locy - 1] -= gx;
/*<       locy=lynl+nodplc(loc+17) >*/
    locy = tabinf_1.lynl + nodplc[loc + 16];
/*<       locyi=imynl+nodplc(loc+17) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 16];
/*<       value(locy)=value(locy)-gmu >*/
    blank_1.value[locy - 1] -= gmu;
/*<       value(locyi)=value(locyi)-xcmu-xcmcb >*/
    blank_1.value[locyi - 1] = blank_1.value[locyi - 1] - xcmu - xcmcb;
/*<       locy=lynl+nodplc(loc+18) >*/
    locy = tabinf_1.lynl + nodplc[loc + 17];
/*<       locyi=imynl+nodplc(loc+18) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 17];
/*<       value(locy)=value(locy)-gpi >*/
    blank_1.value[locy - 1] -= gpi;
/*<       value(locyi)=value(locyi)-xcpi >*/
    blank_1.value[locyi - 1] -= xcpi;
/*<       locy=lynl+nodplc(loc+19) >*/
    locy = tabinf_1.lynl + nodplc[loc + 18];
/*<       value(locy)=value(locy)-gepr >*/
    blank_1.value[locy - 1] -= gepr;
/*<       locy=lynl+nodplc(loc+20) >*/
    locy = tabinf_1.lynl + nodplc[loc + 19];
/*<       locyi=imynl+nodplc(loc+20) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 19];
/*<       value(locy)=value(locy)-go >*/
    blank_1.value[locy - 1] -= go;
/*<       value(locyi)=value(locyi)+xcmcb >*/
    blank_1.value[locyi - 1] += xcmcb;
/*<       locy=lynl+nodplc(loc+21) >*/
    locy = tabinf_1.lynl + nodplc[loc + 20];
/*<       locyi=imynl+nodplc(loc+21) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 20];
/*<       value(locy)=value(locy)-gpi-gm >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] - gpi - gm;
/*<       value(locyi)=value(locyi)-xcpi-xgm-xcmcb >*/
    blank_1.value[locyi - 1] = blank_1.value[locyi - 1] - xcpi - xgm - xcmcb;
/*<       locyi=imynl+nodplc(loc+31) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 30];
/*<       value(locyi)=value(locyi)+xccs >*/
    blank_1.value[locyi - 1] += xccs;
/*<       locyi=imynl+nodplc(loc+32) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 31];
/*<       value(locyi)=value(locyi)-xccs >*/
    blank_1.value[locyi - 1] -= xccs;
/*<       locyi=imynl+nodplc(loc+33) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 32];
/*<       value(locyi)=value(locyi)-xccs >*/
    blank_1.value[locyi - 1] -= xccs;
/*<       locyi=imynl+nodplc(loc+34) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 33];
/*<       value(locyi)=value(locyi)-xcbx >*/
    blank_1.value[locyi - 1] -= xcbx;
/*<       locyi=imynl+nodplc(loc+35) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 34];
/*<       value(locyi)=value(locyi)-xcbx >*/
    blank_1.value[locyi - 1] -= xcbx;
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 260 >*/
    goto L260;

/*  jfets */

/*<   300 loc=locate(13) >*/
L300:
    loc = cirdat_1.locate[12];
/*<   310 if ((loc.eq.0).or.(nodplc(loc+25).ne.0)) go to 350 >*/
L310:
    if (loc == 0 || nodplc[loc + 24] != 0) {
	goto L350;
    }
/*<       locv=nodplc(loc+1) >*/
    locv = nodplc[loc];
/*<       area=value(locv+1) >*/
    area = blank_1.value[locv];
/*<       locm=nodplc(loc+7) >*/
    locm = nodplc[loc + 6];
/*<       locm=nodplc(locm+1) >*/
    locm = nodplc[locm];
/*<       loct=lx0+nodplc(loc+19) >*/
    loct = tabinf_1.lx0 + nodplc[loc + 18];
/*<       gdpr=value(locm+4)*area >*/
    gdpr = blank_1.value[locm + 3] * area;
/*<       gspr=value(locm+5)*area >*/
    gspr = blank_1.value[locm + 4] * area;
/*<       gm=value(loct+5) >*/
    gm = blank_1.value[loct + 4];
/*<       gds=value(loct+6) >*/
    gds = blank_1.value[loct + 5];
/*<       ggs=value(loct+7) >*/
    ggs = blank_1.value[loct + 6];
/*<       xgs=value(loct+9)*omega >*/
    xgs = blank_1.value[loct + 8] * status_1.omega;
/*<       ggd=value(loct+8) >*/
    ggd = blank_1.value[loct + 7];
/*<       xgd=value(loct+11)*omega >*/
    xgd = blank_1.value[loct + 10] * status_1.omega;
/*<       locy=lynl+nodplc(loc+20) >*/
    locy = tabinf_1.lynl + nodplc[loc + 19];
/*<       value(locy)=value(locy)+gdpr >*/
    blank_1.value[locy - 1] += gdpr;
/*<       locy=lynl+nodplc(loc+21) >*/
    locy = tabinf_1.lynl + nodplc[loc + 20];
/*<       locyi=imynl+nodplc(loc+21) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 20];
/*<       value(locy)=value(locy)+ggd+ggs >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] + ggd + ggs;
/*<       value(locyi)=value(locyi)+xgd+xgs >*/
    blank_1.value[locyi - 1] = blank_1.value[locyi - 1] + xgd + xgs;
/*<       locy=lynl+nodplc(loc+22) >*/
    locy = tabinf_1.lynl + nodplc[loc + 21];
/*<       value(locy)=value(locy)+gspr >*/
    blank_1.value[locy - 1] += gspr;
/*<       locy=lynl+nodplc(loc+23) >*/
    locy = tabinf_1.lynl + nodplc[loc + 22];
/*<       locyi=imynl+nodplc(loc+23) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 22];
/*<       value(locy)=value(locy)+gdpr+gds+ggd >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] + gdpr + gds + ggd;
/*<       value(locyi)=value(locyi)+xgd >*/
    blank_1.value[locyi - 1] += xgd;
/*<       locy=lynl+nodplc(loc+24) >*/
    locy = tabinf_1.lynl + nodplc[loc + 23];
/*<       locyi=imynl+nodplc(loc+24) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 23];
/*<       value(locy)=value(locy)+gspr+gds+gm+ggs >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] + gspr + gds + gm + ggs;

/*<       value(locyi)=value(locyi)+xgs >*/
    blank_1.value[locyi - 1] += xgs;
/*<       locy=lynl+nodplc(loc+9) >*/
    locy = tabinf_1.lynl + nodplc[loc + 8];
/*<       value(locy)=value(locy)-gdpr >*/
    blank_1.value[locy - 1] -= gdpr;
/*<       locy=lynl+nodplc(loc+10) >*/
    locy = tabinf_1.lynl + nodplc[loc + 9];
/*<       locyi=imynl+nodplc(loc+10) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 9];
/*<       value(locy)=value(locy)-ggd >*/
    blank_1.value[locy - 1] -= ggd;
/*<       value(locyi)=value(locyi)-xgd >*/
    blank_1.value[locyi - 1] -= xgd;
/*<       locy=lynl+nodplc(loc+11) >*/
    locy = tabinf_1.lynl + nodplc[loc + 10];
/*<       locyi=imynl+nodplc(loc+11) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 10];
/*<       value(locy)=value(locy)-ggs >*/
    blank_1.value[locy - 1] -= ggs;
/*<       value(locyi)=value(locyi)-xgs >*/
    blank_1.value[locyi - 1] -= xgs;
/*<       locy=lynl+nodplc(loc+12) >*/
    locy = tabinf_1.lynl + nodplc[loc + 11];
/*<       value(locy)=value(locy)-gspr >*/
    blank_1.value[locy - 1] -= gspr;
/*<       locy=lynl+nodplc(loc+13) >*/
    locy = tabinf_1.lynl + nodplc[loc + 12];
/*<       value(locy)=value(locy)-gdpr >*/
    blank_1.value[locy - 1] -= gdpr;
/*<       locy=lynl+nodplc(loc+14) >*/
    locy = tabinf_1.lynl + nodplc[loc + 13];
/*<       locyi=imynl+nodplc(loc+14) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 13];
/*<       value(locy)=value(locy)-ggd+gm >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] - ggd + gm;
/*<       value(locyi)=value(locyi)-xgd >*/
    blank_1.value[locyi - 1] -= xgd;
/*<       locy=lynl+nodplc(loc+15) >*/
    locy = tabinf_1.lynl + nodplc[loc + 14];
/*<       value(locy)=value(locy)-gds-gm >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] - gds - gm;
/*<       locy=lynl+nodplc(loc+16) >*/
    locy = tabinf_1.lynl + nodplc[loc + 15];
/*<       locyi=imynl+nodplc(loc+16) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 15];
/*<       value(locy)=value(locy)-ggs-gm >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] - ggs - gm;
/*<       value(locyi)=value(locyi)-xgs >*/
    blank_1.value[locyi - 1] -= xgs;
/*<       locy=lynl+nodplc(loc+17) >*/
    locy = tabinf_1.lynl + nodplc[loc + 16];
/*<       value(locy)=value(locy)-gspr >*/
    blank_1.value[locy - 1] -= gspr;
/*<       locy=lynl+nodplc(loc+18) >*/
    locy = tabinf_1.lynl + nodplc[loc + 17];
/*<       value(locy)=value(locy)-gds >*/
    blank_1.value[locy - 1] -= gds;
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 310 >*/
    goto L310;

/*  mosfets */

/*<   350 loc=locate(14) >*/
L350:
    loc = cirdat_1.locate[13];
/*<   360 if ((loc.eq.0).or.(nodplc(loc+33).ne.0)) go to 400 >*/
L360:
    if (loc == 0 || nodplc[loc + 32] != 0) {
	goto L400;
    }
/*<       locv=nodplc(loc+1) >*/
    locv = nodplc[loc];
/*<       locm=nodplc(loc+8) >*/
    locm = nodplc[loc + 7];
/*<       itype=nodplc(locm+2) >*/
    itype = nodplc[locm + 1];
/*<       locm=nodplc(locm+1) >*/
    locm = nodplc[locm];
/*<       devmod=value(locv+8) >*/
    devmod = blank_1.value[locv + 7];
/*<       xnrm=1.0d0 >*/
    xnrm = 1.;
/*<       xrev=0.0d0 >*/
    xrev = 0.;
/*<       if (devmod.ge.0.0d0) go to 370 >*/
    if (devmod >= 0.) {
	goto L370;
    }
/*<       xnrm=0.0d0 >*/
    xnrm = 0.;
/*<       xrev=1.0d0 >*/
    xrev = 1.;
/*<   370 loct=lx0+nodplc(loc+26) >*/
L370:
    loct = tabinf_1.lx0 + nodplc[loc + 25];
/*<       if (value(locm+7).eq.0.0d0.and. >*/
/*<      1   value(locm+8).eq.0.0d0) go to 375 >*/
    if (blank_1.value[locm + 6] == 0. && blank_1.value[locm + 7] == 0.) {
	goto L375;
    }
/*<       gdpr=value(locm+7) >*/
    gdpr = blank_1.value[locm + 6];
/*<       gspr=value(locm+8) >*/
    gspr = blank_1.value[locm + 7];
/*<       go to 380 >*/
    goto L380;
/*<   375 gdpr=value(locm+16)/value(locv+13) >*/
L375:
    gdpr = blank_1.value[locm + 15] / blank_1.value[locv + 12];
/*<       gspr=value(locm+16)/value(locv+14) >*/
    gspr = blank_1.value[locm + 15] / blank_1.value[locv + 13];
/*<   380 gm=value(loct+7) >*/
L380:
    gm = blank_1.value[loct + 6];
/*<       gds=value(loct+8) >*/
    gds = blank_1.value[loct + 7];
/*<       gmbs=value(loct+9) >*/
    gmbs = blank_1.value[loct + 8];
/*<       gbd=value(loct+10) >*/
    gbd = blank_1.value[loct + 9];
/*<       gbs=value(loct+11) >*/
    gbs = blank_1.value[loct + 10];
/*<       capbd=value(loct+24) >*/
    capbd = blank_1.value[loct + 23];
/*<       capbs=value(loct+26) >*/
    capbs = blank_1.value[loct + 25];
/* c */
/* c    charge oriented model parameters */
/* c */
/*<       xl=value(locv+1)-2.0d0*value(locm+28) >*/
    xl = blank_1.value[locv] - blank_1.value[locm + 27] * 2.;
/*<       xw=value(locv+2) >*/
    xw = blank_1.value[locv + 1];
/*<       xqco=value(locm+35) >*/
    xqco = blank_1.value[locm + 34];
/*<       xqc=value(locv+15) >*/
    xqc = blank_1.value[locv + 14];
/*<       covlgs=value(locm+13)*xw >*/
    covlgs = blank_1.value[locm + 12] * xw;
/*<       covlgd=value(locm+14)*xw >*/
    covlgd = blank_1.value[locm + 13] * xw;
/*<       covlgb=value(locm+15)*xl >*/
    covlgb = blank_1.value[locm + 14] * xl;
/*<       if (xqco.gt.0.5d0) go to 385 >*/
    if (xqco > .5) {
	goto L385;
    }
/*<       cggb=value(loct+18) >*/
    cggb = blank_1.value[loct + 17];
/*<       cgdb=value(loct+19) >*/
    cgdb = blank_1.value[loct + 18];
/*<       cgsb=value(loct+20) >*/
    cgsb = blank_1.value[loct + 19];
/*<       cbgb=value(loct+21) >*/
    cbgb = blank_1.value[loct + 20];
/*<       cbdb=value(loct+22) >*/
    cbdb = blank_1.value[loct + 21];
/*<       cbsb=value(loct+23) >*/
    cbsb = blank_1.value[loct + 22];
/*<       gcg=(cggb+cbgb)*omega >*/
    gcg = (cggb + cbgb) * status_1.omega;
/*<       gcd=(cgdb+cbdb)*omega >*/
    gcd = (cgdb + cbdb) * status_1.omega;
/*<       gcs=(cgsb+cbsb)*omega >*/
    gcs = (cgsb + cbsb) * status_1.omega;
/*<       xcgxd=-xqc*gcg >*/
    xcgxd = -xqc * gcg;
/*<       xcgxs=-(1.0d0-xqc)*gcg >*/
    xcgxs = -(1. - xqc) * gcg;
/*<       xcdxd=-xqc*gcd >*/
    xcdxd = -xqc * gcd;
/*<       xcdxs=-(1.0d0-xqc)*gcd >*/
    xcdxs = -(1. - xqc) * gcd;
/*<       xcsxd=-xqc*gcs >*/
    xcsxd = -xqc * gcs;
/*<       xcsxs=-(1.0d0-xqc)*gcs >*/
    xcsxs = -(1. - xqc) * gcs;
/*<       xcdgb=xcgxd-covlgd*omega >*/
    xcdgb = xcgxd - covlgd * status_1.omega;
/*<       xcddb=xcdxd+(capbd+covlgd)*omega >*/
    xcddb = xcdxd + (capbd + covlgd) * status_1.omega;
/*<       xcdsb=xcsxd >*/
    xcdsb = xcsxd;
/*<       xcsgb=xcgxs-covlgs*omega >*/
    xcsgb = xcgxs - covlgs * status_1.omega;
/*<       xcsdb=xcdxs >*/
    xcsdb = xcdxs;
/*<       xcssb=xcsxs+(capbs+covlgs)*omega >*/
    xcssb = xcsxs + (capbs + covlgs) * status_1.omega;
/*<       xcggb=(cggb+covlgd+covlgs+covlgb)*omega >*/
    xcggb = (cggb + covlgd + covlgs + covlgb) * status_1.omega;
/*<       xcgdb=(cgdb-covlgd)*omega >*/
    xcgdb = (cgdb - covlgd) * status_1.omega;
/*<       xcgsb=(cgsb-covlgs)*omega >*/
    xcgsb = (cgsb - covlgs) * status_1.omega;
/*<       xcbgb=(cbgb-covlgb)*omega >*/
    xcbgb = (cbgb - covlgb) * status_1.omega;
/*<       xcbdb=(cbdb-capbd)*omega >*/
    xcbdb = (cbdb - capbd) * status_1.omega;
/*<       xcbsb=(cbsb-capbs)*omega >*/
    xcbsb = (cbsb - capbs) * status_1.omega;
/*<       go to 390 >*/
    goto L390;

/*     meyer"s model parameters */

/*<   385 xcgs=(value(loct+12)+covlgs)*omega >*/
L385:
    xcgs = (blank_1.value[loct + 11] + covlgs) * status_1.omega;
/*<       xcgd=(value(loct+14)+covlgd)*omega >*/
    xcgd = (blank_1.value[loct + 13] + covlgd) * status_1.omega;
/*<       xcgb=(value(loct+16)+covlgb)*omega >*/
    xcgb = (blank_1.value[loct + 15] + covlgb) * status_1.omega;
/*<       xbd=capbd*omega >*/
    xbd = capbd * status_1.omega;
/*<       xbs=capbs*omega >*/
    xbs = capbs * status_1.omega;
/* c */
/* c    do the mapping from meyer"s model into charge oriented model */
/* c */
/*<       xcggb=xcgd+xcgs+xcgb >*/
    xcggb = xcgd + xcgs + xcgb;
/*<       xcgdb=-xcgd >*/
    xcgdb = -xcgd;
/*<       xcgsb=-xcgs >*/
    xcgsb = -xcgs;
/*<       xcbgb=-xcgb >*/
    xcbgb = -xcgb;
/*<       xcbdb=-xbd >*/
    xcbdb = -xbd;
/*<       xcbsb=-xbs >*/
    xcbsb = -xbs;
/*<       xcddb=xcgd+xbd >*/
    xcddb = xcgd + xbd;
/*<       xcssb=xcgs+xbs >*/
    xcssb = xcgs + xbs;
/* c    xcgsb=-xcgb */
/*<       xcdgb=-xcgd >*/
    xcdgb = -xcgd;
/*<       xcsgb=-xcgs >*/
    xcsgb = -xcgs;
/*<       xcdsb=0.0d0 >*/
    xcdsb = 0.;
/*<       xcsdb=0.0d0 >*/
    xcsdb = 0.;
/* c */
/*<   390 locyi=imynl+nodplc(loc+28) >*/
L390:
    locyi = tabinf_1.imynl + nodplc[loc + 27];
/*<       value(locyi)=value(locyi)+xcggb >*/
    blank_1.value[locyi - 1] += xcggb;
/*<       locyi=imynl+nodplc(loc+30) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 29];
/*<       value(locyi)=value(locyi)-xcbgb-xcbdb-xcbsb >*/
    blank_1.value[locyi - 1] = blank_1.value[locyi - 1] - xcbgb - xcbdb - 
	    xcbsb;
/*<       locyi=imynl+nodplc(loc+31) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 30];
/*<       value(locyi)=value(locyi)+xcddb >*/
    blank_1.value[locyi - 1] += xcddb;
/*<       locyi=imynl+nodplc(loc+32) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 31];
/*<       value(locyi)=value(locyi)+xcssb >*/
    blank_1.value[locyi - 1] += xcssb;
/*<       locyi=imynl+nodplc(loc+11) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 10];
/*<       value(locyi)=value(locyi)-xcggb-xcgdb-xcgsb >*/
    blank_1.value[locyi - 1] = blank_1.value[locyi - 1] - xcggb - xcgdb - 
	    xcgsb;
/*<       locyi=imynl+nodplc(loc+12) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 11];
/*<       value(locyi)=value(locyi)+xcgdb >*/
    blank_1.value[locyi - 1] += xcgdb;
/*<       locyi=imynl+nodplc(loc+13) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 12];
/*<       value(locyi)=value(locyi)+xcgsb >*/
    blank_1.value[locyi - 1] += xcgsb;
/*<       locyi=imynl+nodplc(loc+15) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 14];
/*<       value(locyi)=value(locyi)+xcbgb >*/
    blank_1.value[locyi - 1] += xcbgb;
/*<       locyi=imynl+nodplc(loc+16) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 15];
/*<       value(locyi)=value(locyi)+xcbdb >*/
    blank_1.value[locyi - 1] += xcbdb;
/*<       locyi=imynl+nodplc(loc+17) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 16];
/*<       value(locyi)=value(locyi)+xcbsb >*/
    blank_1.value[locyi - 1] += xcbsb;
/*<       locyi=imynl+nodplc(loc+19) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 18];
/*<       value(locyi)=value(locyi)+xcdgb >*/
    blank_1.value[locyi - 1] += xcdgb;
/*<       locyi=imynl+nodplc(loc+20) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 19];
/*<       value(locyi)=value(locyi)-xcdgb-xcddb-xcdsb >*/
    blank_1.value[locyi - 1] = blank_1.value[locyi - 1] - xcdgb - xcddb - 
	    xcdsb;
/*<       locyi=imynl+nodplc(loc+21) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 20];
/*<       value(locyi)=value(locyi)+xcdsb >*/
    blank_1.value[locyi - 1] += xcdsb;
/*<       locyi=imynl+nodplc(loc+22) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 21];
/*<       value(locyi)=value(locyi)+xcsgb >*/
    blank_1.value[locyi - 1] += xcsgb;
/*<       locyi=imynl+nodplc(loc+24) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 23];
/*<       value(locyi)=value(locyi)-xcsgb-xcsdb-xcssb >*/
    blank_1.value[locyi - 1] = blank_1.value[locyi - 1] - xcsgb - xcsdb - 
	    xcssb;
/*<       locyi=imynl+nodplc(loc+25) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 24];
/*<       value(locyi)=value(locyi)+xcsdb >*/
    blank_1.value[locyi - 1] += xcsdb;
/*<       locy=lynl+nodplc(loc+27) >*/
    locy = tabinf_1.lynl + nodplc[loc + 26];
/*<       value(locy)=value(locy)+gdpr >*/
    blank_1.value[locy - 1] += gdpr;
/*<       locy=lynl+nodplc(loc+29) >*/
    locy = tabinf_1.lynl + nodplc[loc + 28];
/*<       value(locy)=value(locy)+gspr >*/
    blank_1.value[locy - 1] += gspr;
/*<       locy=lynl+nodplc(loc+30) >*/
    locy = tabinf_1.lynl + nodplc[loc + 29];
/*<       value(locy)=value(locy)+gbd+gbs >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] + gbd + gbs;
/*<       locy=lynl+nodplc(loc+31) >*/
    locy = tabinf_1.lynl + nodplc[loc + 30];
/*<       value(locy)=value(locy)+gdpr+gds+gbd+xrev*(gm+gmbs) >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] + gdpr + gds + gbd + 
	    xrev * (gm + gmbs);
/*<       locy=lynl+nodplc(loc+32) >*/
    locy = tabinf_1.lynl + nodplc[loc + 31];
/*<       value(locy)=value(locy)+gspr+gds+gbs+xnrm*(gm+gmbs) >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] + gspr + gds + gbs + 
	    xnrm * (gm + gmbs);
/*<       locy=lynl+nodplc(loc+10) >*/
    locy = tabinf_1.lynl + nodplc[loc + 9];
/*<       value(locy)=value(locy)-gdpr >*/
    blank_1.value[locy - 1] -= gdpr;
/*<       locy=lynl+nodplc(loc+14) >*/
    locy = tabinf_1.lynl + nodplc[loc + 13];
/*<       value(locy)=value(locy)-gspr >*/
    blank_1.value[locy - 1] -= gspr;
/*<       locy=lynl+nodplc(loc+16) >*/
    locy = tabinf_1.lynl + nodplc[loc + 15];
/*<       value(locy)=value(locy)-gbd >*/
    blank_1.value[locy - 1] -= gbd;
/*<       locy=lynl+nodplc(loc+17) >*/
    locy = tabinf_1.lynl + nodplc[loc + 16];
/*<       value(locy)=value(locy)-gbs >*/
    blank_1.value[locy - 1] -= gbs;
/*<       locy=lynl+nodplc(loc+18) >*/
    locy = tabinf_1.lynl + nodplc[loc + 17];
/*<       value(locy)=value(locy)-gdpr >*/
    blank_1.value[locy - 1] -= gdpr;
/*<       locy=lynl+nodplc(loc+19) >*/
    locy = tabinf_1.lynl + nodplc[loc + 18];
/*<       value(locy)=value(locy)+(xnrm-xrev)*gm >*/
    blank_1.value[locy - 1] += (xnrm - xrev) * gm;
/*<       locy=lynl+nodplc(loc+20) >*/
    locy = tabinf_1.lynl + nodplc[loc + 19];
/*<       value(locy)=value(locy)-gbd+(xnrm-xrev)*gmbs >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] - gbd + (xnrm - xrev) * 
	    gmbs;
/*<       locy=lynl+nodplc(loc+21) >*/
    locy = tabinf_1.lynl + nodplc[loc + 20];
/*<       value(locy)=value(locy)-gds-xnrm*(gm+gmbs) >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] - gds - xnrm * (gm + 
	    gmbs);
/*<       locy=lynl+nodplc(loc+22) >*/
    locy = tabinf_1.lynl + nodplc[loc + 21];
/*<       value(locy)=value(locy)-(xnrm-xrev)*gm >*/
    blank_1.value[locy - 1] -= (xnrm - xrev) * gm;
/*<       locy=lynl+nodplc(loc+23) >*/
    locy = tabinf_1.lynl + nodplc[loc + 22];
/*<       value(locy)=value(locy)-gspr >*/
    blank_1.value[locy - 1] -= gspr;
/*<       locy=lynl+nodplc(loc+24) >*/
    locy = tabinf_1.lynl + nodplc[loc + 23];
/*<       value(locy)=value(locy)-gbs-(xnrm-xrev)*gmbs >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] - gbs - (xnrm - xrev) * 
	    gmbs;
/*<       locy=lynl+nodplc(loc+25) >*/
    locy = tabinf_1.lynl + nodplc[loc + 24];
/*<       value(locy)=value(locy)-gds-xrev*(gm+gmbs) >*/
    blank_1.value[locy - 1] = blank_1.value[locy - 1] - gds - xrev * (gm + 
	    gmbs);
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 360 >*/
    goto L360;

/*  transmission lines */

/*<   400 loc=locate(17) >*/
L400:
    loc = cirdat_1.locate[16];
/*<   410 if ((loc.eq.0).or.(nodplc(loc+33).ne.0)) go to 1000 >*/
L410:
    if (loc == 0 || nodplc[loc + 32] != 0) {
	goto L1000;
    }
/*<       locv=nodplc(loc+1) >*/
    locv = nodplc[loc];
/*<       z0=value(locv+1) >*/
    z0 = blank_1.value[locv];
/*<       y0=1.0d0/z0 >*/
    y0 = 1. / z0;
/*<       td=value(locv+2) >*/
    td = blank_1.value[locv + 1];
/*<       arg=-omega*td >*/
    arg = -status_1.omega * td;
/*<       rval=dcos(arg) >*/
    rval = cos(arg);
/*<       xval=dsin(arg) >*/
    xval = sin(arg);
/*<       locy=lynl+nodplc(loc+10) >*/
    locy = tabinf_1.lynl + nodplc[loc + 9];
/*<       value(locy)=value(locy)+y0 >*/
    blank_1.value[locy - 1] += y0;
/*<       locy=lynl+nodplc(loc+11) >*/
    locy = tabinf_1.lynl + nodplc[loc + 10];
/*<       locyi=imynl+nodplc(loc+11) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 10];
/*<       value(locy)=-y0 >*/
    blank_1.value[locy - 1] = -y0;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+12) >*/
    locy = tabinf_1.lynl + nodplc[loc + 11];
/*<       locyi=imynl+nodplc(loc+12) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 11];
/*<       value(locy)=-1.0d0 >*/
    blank_1.value[locy - 1] = -1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+13) >*/
    locy = tabinf_1.lynl + nodplc[loc + 12];
/*<       value(locy)=value(locy)+y0 >*/
    blank_1.value[locy - 1] += y0;
/*<       locy=lynl+nodplc(loc+14) >*/
    locy = tabinf_1.lynl + nodplc[loc + 13];
/*<       locyi=imynl+nodplc(loc+14) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 13];
/*<       value(locy)=-1.0d0 >*/
    blank_1.value[locy - 1] = -1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+15) >*/
    locy = tabinf_1.lynl + nodplc[loc + 14];
/*<       locyi=imynl+nodplc(loc+15) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 14];
/*<       value(locy)=-y0 >*/
    blank_1.value[locy - 1] = -y0;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+16) >*/
    locy = tabinf_1.lynl + nodplc[loc + 15];
/*<       locyi=imynl+nodplc(loc+16) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 15];
/*<       value(locy)=+y0 >*/
    blank_1.value[locy - 1] = y0;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+17) >*/
    locy = tabinf_1.lynl + nodplc[loc + 16];
/*<       locyi=imynl+nodplc(loc+17) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 16];
/*<       value(locy)=+1.0d0 >*/
    blank_1.value[locy - 1] = 1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+18) >*/
    locy = tabinf_1.lynl + nodplc[loc + 17];
/*<       locyi=imynl+nodplc(loc+18) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 17];
/*<       value(locy)=+y0 >*/
    blank_1.value[locy - 1] = y0;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+19) >*/
    locy = tabinf_1.lynl + nodplc[loc + 18];
/*<       locyi=imynl+nodplc(loc+19) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 18];
/*<       value(locy)=+1.0d0 >*/
    blank_1.value[locy - 1] = 1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+20) >*/
    locy = tabinf_1.lynl + nodplc[loc + 19];
/*<       locyi=imynl+nodplc(loc+20) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 19];
/*<       value(locy)=-1.0d0 >*/
    blank_1.value[locy - 1] = -1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+21) >*/
    locy = tabinf_1.lynl + nodplc[loc + 20];
/*<       locyi=imynl+nodplc(loc+21) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 20];
/*<       value(locy)=-rval >*/
    blank_1.value[locy - 1] = -rval;
/*<       value(locyi)=-xval >*/
    blank_1.value[locyi - 1] = -xval;
/*<       locy=lynl+nodplc(loc+22) >*/
    locy = tabinf_1.lynl + nodplc[loc + 21];
/*<       locyi=imynl+nodplc(loc+22) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 21];
/*<       value(locy)=+rval >*/
    blank_1.value[locy - 1] = rval;
/*<       value(locyi)=+xval >*/
    blank_1.value[locyi - 1] = xval;
/*<       locy=lynl+nodplc(loc+23) >*/
    locy = tabinf_1.lynl + nodplc[loc + 22];
/*<       locyi=imynl+nodplc(loc+23) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 22];
/*<       value(locy)=+1.0d0 >*/
    blank_1.value[locy - 1] = 1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+24) >*/
    locy = tabinf_1.lynl + nodplc[loc + 23];
/*<       locyi=imynl+nodplc(loc+24) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 23];
/*<       value(locy)=-rval*z0 >*/
    blank_1.value[locy - 1] = -rval * z0;
/*<       value(locyi)=-xval*z0 >*/
    blank_1.value[locyi - 1] = -xval * z0;
/*<       locy=lynl+nodplc(loc+25) >*/
    locy = tabinf_1.lynl + nodplc[loc + 24];
/*<       locyi=imynl+nodplc(loc+25) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 24];
/*<       value(locy)=-rval >*/
    blank_1.value[locy - 1] = -rval;
/*<       value(locyi)=-xval >*/
    blank_1.value[locyi - 1] = -xval;
/*<       locy=lynl+nodplc(loc+26) >*/
    locy = tabinf_1.lynl + nodplc[loc + 25];
/*<       locyi=imynl+nodplc(loc+26) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 25];
/*<       value(locy)=+rval >*/
    blank_1.value[locy - 1] = rval;
/*<       value(locyi)=+xval >*/
    blank_1.value[locyi - 1] = xval;
/*<       locy=lynl+nodplc(loc+27) >*/
    locy = tabinf_1.lynl + nodplc[loc + 26];
/*<       locyi=imynl+nodplc(loc+27) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 26];
/*<       value(locy)=-1.0d0 >*/
    blank_1.value[locy - 1] = -1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+28) >*/
    locy = tabinf_1.lynl + nodplc[loc + 27];
/*<       locyi=imynl+nodplc(loc+28) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 27];
/*<       value(locy)=+1.0d0 >*/
    blank_1.value[locy - 1] = 1.;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+29) >*/
    locy = tabinf_1.lynl + nodplc[loc + 28];
/*<       locyi=imynl+nodplc(loc+29) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 28];
/*<       value(locy)=-rval*z0 >*/
    blank_1.value[locy - 1] = -rval * z0;
/*<       value(locyi)=-xval*z0 >*/
    blank_1.value[locyi - 1] = -xval * z0;
/*<       locy=lynl+nodplc(loc+31) >*/
    locy = tabinf_1.lynl + nodplc[loc + 30];
/*<       locyi=imynl+nodplc(loc+31) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 30];
/*<       value(locy)=-y0 >*/
    blank_1.value[locy - 1] = -y0;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       locy=lynl+nodplc(loc+32) >*/
    locy = tabinf_1.lynl + nodplc[loc + 31];
/*<       locyi=imynl+nodplc(loc+32) >*/
    locyi = tabinf_1.imynl + nodplc[loc + 31];
/*<       value(locy)=-y0 >*/
    blank_1.value[locy - 1] = -y0;
/*<       value(locyi)=0.0d0 >*/
    blank_1.value[locyi - 1] = 0.;
/*<       loc=nodplc(loc) >*/
    loc = nodplc[loc - 1];
/*<       go to 410 >*/
    goto L410;

/*  finished */

/*<  1000 return >*/
L1000:
    return 0;
/*<       end >*/
} /* acload_ */

#undef cvalue
#undef nodplc


