b1dset.c

ngspice又一个电子CAD仿真软件代码.功能更全

	    TimesDeriv(&d_qb,&d_Argl7,-1.0);
	    d_qb.value += 0.5;
	    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;
      }

#if 0
      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);
#endif

      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);
}