b1dset.c

spice中支持多层次元件模型仿真的可单独运行的插件源码

	    MultDeriv(&d_qb,&d_qb,&d_Arg3);
	    d_qb.value += Vfb + Phi;
	    TimesDeriv(&d_dummy,&d_Vth0,-1.0);
	    PlusDeriv(&d_qb,&d_qb,&d_dummy);
	    TimesDeriv(&d_qb,&d_qb,WLCox);
            qd =  - WLCox * (0.5 * (VgsVth -  Arg1) +
                 Arg1 * Argl6);
	    MultDeriv(&d_dummy,&d_Arg1,&d_Argl6);
	    TimesDeriv(&d_qd,&d_Arg1,-1.0);
	    PlusDeriv(&d_qd,&d_qd,&d_VgsVth);
	    TimesDeriv(&d_qd,&d_qd,0.5);
	    PlusDeriv(&d_qd,&d_qd,&d_dummy);
	    TimesDeriv(&d_qd,&d_qd,-WLCox);
            goto finished;
        } else if( vds >= VdsPinchoff ) {    /* saturation region   */
            Args1 = 1.0 /  (3*A);
	    TimesDeriv(&d_Args1,&d_A,3.0);
	    InvDeriv(&d_Args1,&d_Args1);

            qg = WLCox * (vgs - Vfb - Phi - VgsVth * Args1);
	    MultDeriv(&d_qg,&d_VgsVth,&d_Args1);
	    d_qg.value += Vfb + Phi - vgs;
	    d_qg.d1_p -= 1.0;
	    TimesDeriv(&d_qg,&d_qg,-WLCox);
            qb = WLCox * (Vfb + Phi - Vth0 + (1.0 - A) * VgsVth * Args1);
	    TimesDeriv(&d_dummy,&d_A,-1.0);
	    d_dummy.value += 1.0;
	    MultDeriv(&d_qb,&d_VgsVth, &d_dummy);
	    MultDeriv(&d_qb,&d_qb,&d_Args1);
	    d_qb.value += Vfb + Phi;
	    TimesDeriv(&d_dummy,&d_Vth0,-1.0);
	    PlusDeriv(&d_qb,&d_qb,&d_dummy);
	    TimesDeriv(&d_qb,&d_qb,WLCox);
            qd = -co4v15*WLCox*VgsVth;
	    TimesDeriv(&d_qd,&d_VgsVth,-co4v15*WLCox);
            goto finished;
        }


    }

}   
finished:       /* returning Values to Calling Routine */
       
      
        

	    /*
	     * the above has set up (DrCur) and (the node q's)
	     * and (their derivatives upto third order wrt vgs, vbs, and
	     * vds)
	     */
            /*
             *  
             */
                /*
                 *  charge storage elements
                 *
                 *   bulk-drain and bulk-source depletion capacitances
                 *  czbd : zero bias drain junction capacitance
                 *  czbs : zero bias source junction capacitance
                 * czbdsw:zero bias drain junction sidewall capacitance
                 * czbssw:zero bias source junction sidewall capacitance
                 */

                czbd  = model->B1unitAreaJctCap * DrainArea;
                czbs  = model->B1unitAreaJctCap * SourceArea;
                czbdsw= model->B1unitLengthSidewallJctCap * DrainPerimeter;
                czbssw= model->B1unitLengthSidewallJctCap * SourcePerimeter;
                PhiB = model->B1bulkJctPotential;
                PhiBSW = model->B1sidewallJctPotential;
                MJ = model->B1bulkJctBotGradingCoeff;
                MJSW = model->B1bulkJctSideGradingCoeff;

                /* Source Bulk Junction */
                if( vbs < 0 ) {  
                    arg = 1 - vbs / PhiB;
                    argsw = 1 - vbs / PhiBSW;
                    sarg = exp(-MJ*log(arg));
                    sargsw = exp(-MJSW*log(argsw));
                    /* *(ckt->CKTstate0 + here->B1qbs) =
                        PhiB * czbs * (1-arg*sarg)/(1-MJ) + PhiBSW * 
                    czbssw * (1-argsw*sargsw)/(1-MJSW); */
                    capbs1 = czbs * sarg + czbssw * sargsw ;
		    capbs2 = (czbs * MJ * sarg / (PhiB*arg)
				+ czbssw * MJSW * sargsw /(PhiBSW*argsw))/2.0;
		    capbs3 = (czbs * (MJ) * (MJ + 1.) * sarg /((PhiB*arg)*(PhiB*arg))
				+ czbssw * MJSW * (MJSW + 1.) * sargsw / (PhiBSW*argsw*PhiBSW*argsw))/6.0;
                } else {  
                    /* *(ckt->CKTstate0+here->B1qbs) =
                        vbs*(czbs+czbssw)+ vbs*vbs*(czbs*MJ*0.5/PhiB 
                        + czbssw * MJSW * 0.5/PhiBSW); */
                    capbs1 = czbs + czbssw + vbs *(czbs*MJ/PhiB+
                        czbssw * MJSW / PhiBSW );
		    capbs2 = (czbs*MJ/PhiB+czbssw * MJSW / PhiBSW )*0.5;
		    capbs3 = 0.0;
                }

                /* Drain Bulk Junction */
                if( vbd < 0 ) {  
                    arg = 1 - vbd / PhiB;
                    argsw = 1 - vbd / PhiBSW;
                    sarg = exp(-MJ*log(arg));
                    sargsw = exp(-MJSW*log(argsw));
                    /* *(ckt->CKTstate0 + here->B1qbd) =
                        PhiB * czbd * (1-arg*sarg)/(1-MJ) + PhiBSW * 
                    czbdsw * (1-argsw*sargsw)/(1-MJSW); */
                    capbd1 = czbd * sarg + czbdsw * sargsw ;
		    capbd2 = (czbd * MJ * sarg / (PhiB*arg)
				+ czbdsw * MJSW * sargsw /(PhiBSW*argsw))*0.5;
		    capbd3 = (czbd * (MJ) * (MJ + 1.) * sarg /((PhiB*arg)*(PhiB*arg))
				+ czbdsw * MJSW * (MJSW + 1.) * sargsw / (PhiBSW*argsw*PhiBSW*argsw))/6.0;
                } else {  
                    /* *(ckt->CKTstate0+here->B1qbd) =
                        vbd*(czbd+czbdsw)+ vbd*vbd*(czbd*MJ*0.5/PhiB 
                        + czbdsw * MJSW * 0.5/PhiBSW); */
                    capbd1 = czbd + czbdsw + vbd *(czbd*MJ/PhiB+
                        czbdsw * MJSW / PhiBSW );
		    capbd2 = 0.5*(czbs*MJ/PhiB+czbssw * MJSW / PhiBSW );
		    capbd3 = 0.0;
                }





line755:
             
    /*
    qgd = GateDrainOverlapCap * (vgs - vds);
    qgs = GateSourceOverlapCap * vgs;
    qgb = GateBulkOverlapCap * (vgs -vbs);
    *qGatePointer = *qGatePointer + qgd + qgs + qgb;
    *qBulkPointer = *qBulkPointer - qgb;
    *qDrainPointer = *qDrainPointer - qgd;
    *qSourcePointer = -(*qGatePointer + *qBulkPointer + *qDrainPointer);
    */
    d_qg.d1_p += GateDrainOverlapCap + GateSourceOverlapCap + GateBulkOverlapCap;
    d_qg.d1_q += -GateBulkOverlapCap;
    d_qg.d1_r += -GateDrainOverlapCap;
    d_qb.d1_p += -GateBulkOverlapCap;
    d_qb.d1_q += GateBulkOverlapCap + capbs1 + capbd1;
    d_qb.d1_r += -capbd1;
    d_qd.d1_p += - GateDrainOverlapCap;
    d_qd.d1_q += -capbd1;
    d_qd.d1_r += GateDrainOverlapCap + capbd1;
    /*
    d[23]_qg_d[vlgbds23] += 0.0;
    d2_qb_dvgs2 += 0.0;
    d3_qb_dvgs3 += 0.0
    d[23]_qb_dvgs[dvbds23] += 0.0
    */
    d_qb.d2_q2 += 2* ( capbd2 + capbs2);
    d_qb.d3_q3 += 6*(capbd3 + capbs3);
    d_qb.d2_qr += -2*capbd2;
    d_qb.d3_q2r += -capbd3*6;
    d_qb.d3_qr2 += capbd3*6;
    d_qb.d2_r2 += 2*capbd2;
    d_qb.d3_r3 += -6*capbd3;
    /*
    d[23]_qd_dp[dvbds23] += 0.0
    */
    d_qd.d2_q2 -= 2*capbd2;
    d_qd.d3_q3 -= 6*capbd3;
    d_qd.d2_r2 -= 2*capbd2;
    d_qd.d3_r3 -= -6*capbd3;
    d_qd.d2_qr -= -2*capbd2;
    d_qd.d3_q2r -= -6*capbd3;
    d_qd.d3_qr2 -= 6*capbd3;


	/* get all the coefficients and adjust for mode and type */
	if (here->B1mode == 1)
		{
		/* normal mode - no source-drain interchange */
here->qg_x = d_qg.d1_p;
here->qg_y = d_qg.d1_q;
here->qg_z = d_qg.d1_r;
here->qg_x2 = d_qg.d2_p2;
here->qg_y2 = d_qg.d2_q2;
here->qg_z2 = d_qg.d2_r2;
here->qg_xy = d_qg.d2_pq;
here->qg_yz = d_qg.d2_qr;
here->qg_xz = d_qg.d2_pr;
here->qg_x3 = d_qg.d3_p3;
here->qg_y3 = d_qg.d3_q3;
here->qg_z3 = d_qg.d3_r3;
here->qg_x2z = d_qg.d3_p2r;
here->qg_x2y = d_qg.d3_p2q;
here->qg_y2z = d_qg.d3_q2r;
here->qg_xy2 = d_qg.d3_pq2;
here->qg_xz2 = d_qg.d3_pr2;
here->qg_yz2 = d_qg.d3_qr2;
here->qg_xyz = d_qg.d3_pqr;

here->qb_x = d_qb.d1_p;
here->qb_y = d_qb.d1_q;
here->qb_z = d_qb.d1_r;
here->qb_x2 = d_qb.d2_p2;
here->qb_y2 = d_qb.d2_q2;
here->qb_z2 = d_qb.d2_r2;
here->qb_xy = d_qb.d2_pq;
here->qb_yz = d_qb.d2_qr;
here->qb_xz = d_qb.d2_pr;
here->qb_x3 = d_qb.d3_p3;
here->qb_y3 = d_qb.d3_q3;
here->qb_z3 = d_qb.d3_r3;
here->qb_x2z = d_qb.d3_p2r;
here->qb_x2y = d_qb.d3_p2q;
here->qb_y2z = d_qb.d3_q2r;
here->qb_xy2 = d_qb.d3_pq2;
here->qb_xz2 = d_qb.d3_pr2;
here->qb_yz2 = d_qb.d3_qr2;
here->qb_xyz = d_qb.d3_pqr;

here->qd_x = d_qd.d1_p;
here->qd_y = d_qd.d1_q;
here->qd_z = d_qd.d1_r;
here->qd_x2 = d_qd.d2_p2;
here->qd_y2 = d_qd.d2_q2;
here->qd_z2 = d_qd.d2_r2;
here->qd_xy = d_qd.d2_pq;
here->qd_yz = d_qd.d2_qr;
here->qd_xz = d_qd.d2_pr;
here->qd_x3 = d_qd.d3_p3;
here->qd_y3 = d_qd.d3_q3;
here->qd_z3 = d_qd.d3_r3;
here->qd_x2z = d_qd.d3_p2r;
here->qd_x2y = d_qd.d3_p2q;
here->qd_y2z = d_qd.d3_q2r;
here->qd_xy2 = d_qd.d3_pq2;
here->qd_xz2 = d_qd.d3_pr2;
here->qd_yz2 = d_qd.d3_qr2;
here->qd_xyz = d_qd.d3_pqr;

here->DrC_x = d_DrCur.d1_p;
here->DrC_y = d_DrCur.d1_q;
here->DrC_z = d_DrCur.d1_r;
here->DrC_x2 = d_DrCur.d2_p2;
here->DrC_y2 = d_DrCur.d2_q2;
here->DrC_z2 = d_DrCur.d2_r2;
here->DrC_xy = d_DrCur.d2_pq;
here->DrC_yz = d_DrCur.d2_qr;
here->DrC_xz = d_DrCur.d2_pr;
here->DrC_x3 = d_DrCur.d3_p3;
here->DrC_y3 = d_DrCur.d3_q3;
here->DrC_z3 = d_DrCur.d3_r3;
here->DrC_x2z = d_DrCur.d3_p2r;
here->DrC_x2y = d_DrCur.d3_p2q;
here->DrC_y2z = d_DrCur.d3_q2r;
here->DrC_xy2 = d_DrCur.d3_pq2;
here->DrC_xz2 = d_DrCur.d3_pr2;
here->DrC_yz2 = d_DrCur.d3_qr2;
here->DrC_xyz = d_DrCur.d3_pqr;

here->gbs1 = lgbs1;
here->gbs2 = lgbs2;
here->gbs3 = lgbs3;

here->gbd1 = lgbd1;
here->gbd2 = lgbd2;
here->gbd3 = lgbd3;

} else {
		/*
		 * inverse mode - source and drain interchanged
		 * inversion equations for realqg and realqb are the 
		 * same; a minus is added for the
		 * realDrCur equation; 
		 * realqd = -(realqb + realqg + fakeqd)
		 */

here->qg_x = -(-d_qg.d1_p);
here->qg_y = -(-d_qg.d1_q);
here->qg_z = -(d_qg.d1_p + d_qg.d1_q + d_qg.d1_r);
here->qg_x2 = -(-d_qg.d2_p2);
here->qg_y2 = -(-d_qg.d2_q2);
here->qg_z2 = -(-(d_qg.d2_p2 + d_qg.d2_q2 + d_qg.d2_r2 + 2*(d_qg.d2_pq + d_qg.d2_pr + d_qg.d2_qr)));
here->qg_xy = -(-d_qg.d2_pq);
here->qg_yz = -(d_qg.d2_pq + d_qg.d2_q2 + d_qg.d2_qr);
here->qg_xz = -(d_qg.d2_p2 + d_qg.d2_pq + d_qg.d2_pr);
here->qg_x3 = -(-d_qg.d3_p3);
here->qg_y3 = -(-d_qg.d3_q3);
here->qg_z3 = -(d_qg.d3_p3 + d_qg.d3_q3 + d_qg.d3_r3 + 3*(d_qg.d3_p2q + d_qg.d3_p2r + d_qg.d3_pq2 + d_qg.d3_q2r + d_qg.d3_pr2 + d_qg.d3_qr2) + 6*d_qg.d3_pqr );
here->qg_x2z = -(d_qg.d3_p3 + d_qg.d3_p2q + d_qg.d3_p2r);
here->qg_x2y = -(-d_qg.d3_p2q);
here->qg_y2z = -(d_qg.d3_pq2 + d_qg.d3_q3 + d_qg.d3_q2r);
here->qg_xy2 = -(-d_qg.d3_pq2);
here->qg_xz2 = -(-(d_qg.d3_p3 + 2*(d_qg.d3_p2q + d_qg.d3_p2r + d_qg.d3_pqr) + d_qg.d3_pq2 + d_qg.d3_pr2));
here->qg_yz2 = -(-(d_qg.d3_q3 + 2*(d_qg.d3_pq2 + d_qg.d3_q2r + d_qg.d3_pqr) + d_qg.d3_p2q + d_qg.d3_qr2));
here->qg_xyz = -(d_qg.d3_p2q + d_qg.d3_pq2 + d_qg.d3_pqr);

here->qb_x = -(-d_qb.d1_p);
here->qb_y = -(-d_qb.d1_q);
here->qb_z = -(d_qb.d1_p + d_qb.d1_q + d_qb.d1_r);
here->qb_x2 = -(-d_qb.d2_p2);
here->qb_y2 = -(-d_qb.d2_q2);
here->qb_z2 = -(-(d_qb.d2_p2 + d_qb.d2_q2 + d_qb.d2_r2 + 2*(d_qb.d2_pq + d_qb.d2_pr + d_qb.d2_qr)));
here->qb_xy = -(-d_qb.d2_pq);
here->qb_yz = -(d_qb.d2_pq + d_qb.d2_q2 + d_qb.d2_qr);
here->qb_xz = -(d_qb.d2_p2 + d_qb.d2_pq + d_qb.d2_pr);
here->qb_x3 = -(-d_qb.d3_p3);
here->qb_y3 = -(-d_qb.d3_q3);
here->qb_z3 = -(d_qb.d3_p3 + d_qb.d3_q3 + d_qb.d3_r3 + 3*(d_qb.d3_p2q + d_qb.d3_p2r + d_qb.d3_pq2 + d_qb.d3_q2r + d_qb.d3_pr2 + d_qb.d3_qr2) + 6*d_qb.d3_pqr );
here->qb_x2z = -(d_qb.d3_p3 + d_qb.d3_p2q + d_qb.d3_p2r);
here->qb_x2y = -(-d_qb.d3_p2q);
here->qb_y2z = -(d_qb.d3_pq2 + d_qb.d3_q3 + d_qb.d3_q2r);
here->qb_xy2 = -(-d_qb.d3_pq2);
here->qb_xz2 = -(-(d_qb.d3_p3 + 2*(d_qb.d3_p2q + d_qb.d3_p2r + d_qb.d3_pqr) + d_qb.d3_pq2 + d_qb.d3_pr2));
here->qb_yz2 = -(-(d_qb.d3_q3 + 2*(d_qb.d3_pq2 + d_qb.d3_q2r + d_qb.d3_pqr) + d_qb.d3_p2q + d_qb.d3_qr2));
here->qb_xyz = -(d_qb.d3_p2q + d_qb.d3_pq2 + d_qb.d3_pqr);

here->qd_x= -here->qg_x - here->qb_x +(-d_qd.d1_p);
here->qd_y= -here->qg_y - here->qb_y +(-d_qd.d1_q);
here->qd_z= -here->qg_z - here->qb_z +(d_qd.d1_p + d_qd.d1_q + d_qd.d1_r);
here->qd_x2 = -here->qg_x2 - here->qb_x2 +(-d_qd.d2_p2);
here->qd_y2 = -here->qg_y2 - here->qb_y2 +(-d_qd.d2_q2);
here->qd_z2 = -here->qg_z2 - here->qb_z2 +(-(d_qd.d2_p2 + d_qd.d2_q2 + d_qd.d2_r2 + 2*(d_qd.d2_pq + d_qd.d2_pr + d_qd.d2_qr)));
here->qd_xy = -here->qg_xy - here->qb_xy +(-d_qd.d2_pq);
here->qd_yz = -here->qg_yz - here->qb_yz +(d_qd.d2_pq + d_qd.d2_q2 + d_qd.d2_qr);
here->qd_xz = -here->qg_xz - here->qb_xz +(d_qd.d2_p2 + d_qd.d2_pq + d_qd.d2_pr);
here->qd_x3 = -here->qg_x3 - here->qb_x3 +(-d_qd.d3_p3);
here->qd_y3 = -here->qg_y3 - here->qb_y3 +(-d_qd.d3_q3);
here->qd_z3 = -here->qg_z3 - here->qb_z3 +(d_qd.d3_p3 + d_qd.d3_q3 + d_qd.d3_r3 + 3*(d_qd.d3_p2q + d_qd.d3_p2r + d_qd.d3_pq2 + d_qd.d3_q2r + d_qd.d3_pr2 + d_qd.d3_qr2) + 6*d_qd.d3_pqr );
here->qd_x2z = -here->qg_x2z - here->qb_x2z +(d_qd.d3_p3 + d_qd.d3_p2q + d_qd.d3_p2r);
here->qd_x2y = -here->qg_x2y - here->qb_x2y +(-d_qd.d3_p2q);
here->qd_y2z = -here->qg_y2z - here->qb_y2z +(d_qd.d3_pq2 + d_qd.d3_q3 + d_qd.d3_q2r);
here->qd_xy2 = -here->qg_xy2 - here->qb_xy2 +(-d_qd.d3_pq2);
here->qd_xz2 = -here->qg_xz2 - here->qb_xz2 +(-(d_qd.d3_p3 + 2*(d_qd.d3_p2q + d_qd.d3_p2r + d_qd.d3_pqr) + d_qd.d3_pq2 + d_qd.d3_pr2));
here->qd_yz2 = -here->qg_yz2 - here->qb_yz2 +(-(d_qd.d3_q3 + 2*(d_qd.d3_pq2 + d_qd.d3_q2r + d_qd.d3_pqr) + d_qd.d3_p2q + d_qd.d3_qr2));
here->qd_xyz = -here->qg_xyz - here->qb_xyz +(d_qd.d3_p2q + d_qd.d3_pq2 + d_qd.d3_pqr);

here->DrC_x = -d_DrCur.d1_p;
here->DrC_y = -d_DrCur.d1_q;
here->DrC_z = d_DrCur.d1_p + d_DrCur.d1_q + d_DrCur.d1_r;
here->DrC_x2 = -d_DrCur.d2_p2;
here->DrC_y2 = -d_DrCur.d2_q2;
here->DrC_z2 = -(d_DrCur.d2_p2 + d_DrCur.d2_q2 + d_DrCur.d2_r2 + 2*(d_DrCur.d2_pq + d_DrCur.d2_pr + d_DrCur.d2_qr));
here->DrC_xy = -d_DrCur.d2_pq;
here->DrC_yz = d_DrCur.d2_pq + d_DrCur.d2_q2 + d_DrCur.d2_qr;
here->DrC_xz = d_DrCur.d2_p2 + d_DrCur.d2_pq + d_DrCur.d2_pr;
here->DrC_x3 = -d_DrCur.d3_p3;
here->DrC_y3 = -d_DrCur.d3_q3;
here->DrC_z3 = d_DrCur.d3_p3 + d_DrCur.d3_q3 + d_DrCur.d3_r3 + 3*(d_DrCur.d3_p2q + d_DrCur.d3_p2r + d_DrCur.d3_pq2 + d_DrCur.d3_q2r + d_DrCur.d3_pr2 + d_DrCur.d3_qr2) + 6*d_DrCur.d3_pqr ;
here->DrC_x2z = d_DrCur.d3_p3 + d_DrCur.d3_p2q + d_DrCur.d3_p2r;
here->DrC_x2y = -d_DrCur.d3_p2q;
here->DrC_y2z = d_DrCur.d3_pq2 + d_DrCur.d3_q3 + d_DrCur.d3_q2r;
here->DrC_xy2 = -d_DrCur.d3_pq2;
here->DrC_xz2 = -(d_DrCur.d3_p3 + 2*(d_DrCur.d3_p2q + d_DrCur.d3_p2r + d_DrCur.d3_pqr) + d_DrCur.d3_pq2 + d_DrCur.d3_pr2);
here->DrC_yz2 = -(d_DrCur.d3_q3 + 2*(d_DrCur.d3_pq2 + d_DrCur.d3_q2r + d_DrCur.d3_pqr) + d_DrCur.d3_p2q + d_DrCur.d3_qr2);
here->DrC_xyz = d_DrCur.d3_p2q + d_DrCur.d3_pq2 + d_DrCur.d3_pqr;

here->gbs1 = lgbd1;
here->gbs2 = lgbd2;
here->gbs3 = lgbd3;

here->gbd1 = lgbs1;
here->gbd2 = lgbs2;
here->gbd3 = lgbs3;
}

/* now to adjust for type and multiply by factors to convert to Taylor coeffs. */

here->qg_x2 = 0.5*model->B1type*here->qg_x2;
here->qg_y2 = 0.5*model->B1type*here->qg_y2;
here->qg_z2 = 0.5*model->B1type*here->qg_z2;
here->qg_xy = model->B1type*here->qg_xy;
here->qg_yz = model->B1type*here->qg_yz;
here->qg_xz = model->B1type*here->qg_xz;
here->qg_x3 = here->qg_x3/6.;
here->qg_y3 = here->qg_y3/6.;
here->qg_z3 = here->qg_z3/6.;
here->qg_x2z = 0.5*here->qg_x2z;
here->qg_x2y = 0.5*here->qg_x2y;
here->qg_y2z = 0.5*here->qg_y2z;
here->qg_xy2 = 0.5*here->qg_xy2;
here->qg_xz2 = 0.5*here->qg_xz2;
here->qg_yz2 = 0.5*here->qg_yz2;

here->qb_x2 = 0.5*model->B1type*here->qb_x2;
here->qb_y2 = 0.5*model->B1type*here->qb_y2;
here->qb_z2 = 0.5*model->B1type*here->qb_z2;
here->qb_xy = model->B1type*here->qb_xy;
here->qb_yz = model->B1type*here->qb_yz;
here->qb_xz = model->B1type*here->qb_xz;
here->qb_x3 = here->qb_x3/6.;
here->qb_y3 = here->qb_y3/6.;
here->qb_z3 = here->qb_z3/6.;
here->qb_x2z = 0.5*here->qb_x2z;
here->qb_x2y = 0.5*here->qb_x2y;
here->qb_y2z = 0.5*here->qb_y2z;
here->qb_xy2 = 0.5*here->qb_xy2;
here->qb_xz2 = 0.5*here->qb_xz2;
here->qb_yz2 = 0.5*here->qb_yz2;

here->qd_x2 = 0.5*model->B1type*here->qd_x2;
here->qd_y2 = 0.5*model->B1type*here->qd_y2;
here->qd_z2 = 0.5*model->B1type*here->qd_z2;
here->qd_xy = model->B1type*here->qd_xy;
here->qd_yz = model->B1type*here->qd_yz;
here->qd_xz = model->B1type*here->qd_xz;
here->qd_x3 = here->qd_x3/6.;
here->qd_y3 = here->qd_y3/6.;
here->qd_z3 = here->qd_z3/6.;
here->qd_x2z = 0.5*here->qd_x2z;
here->qd_x2y = 0.5*here->qd_x2y;
here->qd_y2z = 0.5*here->qd_y2z;
here->qd_xy2 = 0.5*here->qd_xy2;
here->qd_xz2 = 0.5*here->qd_xz2;
here->qd_yz2 = 0.5*here->qd_yz2;

here->DrC_x2 = 0.5*model->B1type*here->DrC_x2;
here->DrC_y2 = 0.5*model->B1type*here->DrC_y2;
here->DrC_z2 = 0.5*model->B1type*here->DrC_z2;
here->DrC_xy = model->B1type*here->DrC_xy;
here->DrC_yz = model->B1type*here->DrC_yz;
here->DrC_xz = model->B1type*here->DrC_xz;
here->DrC_x3 = here->DrC_x3/6.;
here->DrC_y3 = here->DrC_y3/6.;
here->DrC_z3 = here->DrC_z3/6.;
here->DrC_x2z = 0.5*here->DrC_x2z;
here->DrC_x2y = 0.5*here->DrC_x2y;
here->DrC_y2z = 0.5*here->DrC_y2z;
here->DrC_xy2 = 0.5*here->DrC_xy2;
here->DrC_xz2 = 0.5*here->DrC_xz2;
here->DrC_yz2 = 0.5*here->DrC_yz2;

here->gbs2 = model->B1type*here->gbs2;
here->gbd2 = model->B1type*here->gbd2;
        }   /* End of Mosfet Instance */

    }       /* End of Model Instance */
    return(OK);
}