/**********
Copyright 1990 Regents of the University of California.  All rights reserved.
Author: 1985 Hong J. Park, Thomas L. Quarles
**********/

#include "ngspice/ngspice.h"
#include "ngspice/smpdefs.h"
#include "ngspice/cktdefs.h"
#include "bsim2def.h"
#include "ngspice/const.h"
#include "ngspice/sperror.h"
#include "ngspice/suffix.h"

/* ARGSUSED */
int
B2temp(GENmodel *inModel, CKTcircuit *ckt)
        /* load the B2 device structure with those pointers needed later 
         * for fast matrix loading 
         */

{
    B2model *model = (B2model*) inModel;
    B2instance *here;
    struct bsim2SizeDependParam *pSizeDependParamKnot, *pLastKnot;
    double  EffectiveLength;
    double EffectiveWidth;
    double CoxWoverL, Inv_L, Inv_W, tmp;
    int Size_Not_Found;

    NG_IGNORE(ckt);

    /*  loop through all the B2 device models */
    for( ; model != NULL; model = model->B2nextModel ) {
    
/* Default value Processing for B2 MOSFET Models */
        /* Some Limiting for Model Parameters */
        if( model->B2bulkJctPotential < 0.1)  {
            model->B2bulkJctPotential = 0.1;
        }
        if( model->B2sidewallJctPotential < 0.1)  {
            model->B2sidewallJctPotential = 0.1;
        }

        model->B2Cox = 3.453e-13/(model->B2tox * 1.0e-4);/*in F/cm**2 */
	model->B2vdd2 = 2.0 * model->B2vdd;
	model->B2vgg2 = 2.0 * model->B2vgg;
	model->B2vbb2 = 2.0 * model->B2vbb;
	model->B2Vtm = 8.625e-5 * (model->B2temp + 273.0);
	model->pSizeDependParamKnot = NULL;
	pLastKnot = NULL;

        /* loop through all the instances of the model */
        for (here = model->B2instances; here != NULL ;
                here=here->B2nextInstance) {

	    pSizeDependParamKnot = model->pSizeDependParamKnot;
	    Size_Not_Found = 1;

	    while ((pSizeDependParamKnot != NULL) && Size_Not_Found)
	    {   if ((here->B2l == pSizeDependParamKnot->Length)
		    && (here->B2w == pSizeDependParamKnot->Width))
                {   Size_Not_Found = 0;
		    here->pParam = pSizeDependParamKnot;
		}
		else
		{   pLastKnot = pSizeDependParamKnot;
		    pSizeDependParamKnot = pSizeDependParamKnot->pNext;
		}
            }

	    if (Size_Not_Found)
	    {   here->pParam = TMALLOC(struct bsim2SizeDependParam, 1);
                if (pLastKnot == NULL)
		    model->pSizeDependParamKnot = here->pParam;
                else
		    pLastKnot->pNext = here->pParam;
                here->pParam->pNext = NULL;

		EffectiveLength = here->B2l - model->B2deltaL * 1.0e-6;
		EffectiveWidth = here->B2w - model->B2deltaW * 1.0e-6;

                if(EffectiveLength<=0)
		{  IFuid namarray[2];
                   namarray[0] = model->B2modName;
                   namarray[1] = here->B2name;
                   SPfrontEnd->IFerror (ERR_FATAL,
                    "B2: mosfet %s, model %s: Effective channel length <=0",
                    namarray);
                   return(E_BADPARM);
                }

                if(EffectiveWidth <= 0)
		{  IFuid namarray[2];
                   namarray[0] = model->B2modName;
                   namarray[1] = here->B2name;
                   SPfrontEnd->IFerror (ERR_FATAL,
                    "B2: mosfet %s, model %s: Effective channel width <=0",
                    namarray);
                   return(E_BADPARM);
                }

		Inv_L = 1.0e-6 / EffectiveLength;
		Inv_W = 1.0e-6 / EffectiveWidth;
		here->pParam->Width = here->B2w;
		here->pParam->Length = here->B2l;
                here->pParam->B2vfb = model->B2vfb0 + model->B2vfbW * Inv_W
			       + model->B2vfbL * Inv_L;
                here->pParam->B2phi = model->B2phi0 + model->B2phiW * Inv_W
			       + model->B2phiL * Inv_L;
                here->pParam->B2k1 = model->B2k10 + model->B2k1W * Inv_W
			       + model->B2k1L * Inv_L;
                here->pParam->B2k2 = model->B2k20 + model->B2k2W * Inv_W
			       + model->B2k2L * Inv_L;
                here->pParam->B2eta0 = model->B2eta00 
			       + model->B2eta0W * Inv_W
			       + model->B2eta0L * Inv_L;
                here->pParam->B2etaB = model->B2etaB0 + model->B2etaBW 
			       * Inv_W + model->B2etaBL * Inv_L;
                here->pParam->B2beta0 = model->B2mob00;
                here->pParam->B2beta0B = model->B2mob0B0 
			       + model->B2mob0BW * Inv_W
			       + model->B2mob0BL * Inv_L;
                here->pParam->B2betas0 = model->B2mobs00 
			       + model->B2mobs0W * Inv_W
			       + model->B2mobs0L * Inv_L;
                if (here->pParam->B2betas0 < 1.01 * here->pParam->B2beta0)
	            here->pParam->B2betas0 = 1.01 * here->pParam->B2beta0;
                here->pParam->B2betasB = model->B2mobsB0 
			       + model->B2mobsBW * Inv_W
			       + model->B2mobsBL * Inv_L;
                tmp = (here->pParam->B2betas0 - here->pParam->B2beta0
		               - here->pParam->B2beta0B * model->B2vbb);
                if ((-here->pParam->B2betasB * model->B2vbb) > tmp)
	            here->pParam->B2betasB = -tmp / model->B2vbb;
                here->pParam->B2beta20 = model->B2mob200 
			      + model->B2mob20W * Inv_W
			      + model->B2mob20L * Inv_L;
                here->pParam->B2beta2B = model->B2mob2B0 
			      + model->B2mob2BW * Inv_W
			      + model->B2mob2BL * Inv_L;
                here->pParam->B2beta2G = model->B2mob2G0 
			      + model->B2mob2GW * Inv_W
			      + model->B2mob2GL * Inv_L;
                here->pParam->B2beta30 = model->B2mob300 
			      + model->B2mob30W * Inv_W
			      + model->B2mob30L * Inv_L;
                here->pParam->B2beta3B = model->B2mob3B0 
			      + model->B2mob3BW * Inv_W
			      + model->B2mob3BL * Inv_L;
                here->pParam->B2beta3G = model->B2mob3G0 
			      + model->B2mob3GW * Inv_W
			      + model->B2mob3GL * Inv_L;
                here->pParam->B2beta40 = model->B2mob400 
			      + model->B2mob40W * Inv_W
			      + model->B2mob40L * Inv_L;
                here->pParam->B2beta4B = model->B2mob4B0 
			      + model->B2mob4BW * Inv_W
			      + model->B2mob4BL * Inv_L;
                here->pParam->B2beta4G = model->B2mob4G0 
			      + model->B2mob4GW * Inv_W
			      + model->B2mob4GL * Inv_L;

		CoxWoverL = model->B2Cox * EffectiveWidth / EffectiveLength;

	        here->pParam->B2beta0 *= CoxWoverL;
	        here->pParam->B2beta0B *= CoxWoverL;
	        here->pParam->B2betas0 *= CoxWoverL;
	        here->pParam->B2betasB *= CoxWoverL;
	        here->pParam->B2beta30 *= CoxWoverL;
	        here->pParam->B2beta3B *= CoxWoverL;
	        here->pParam->B2beta3G *= CoxWoverL;
	        here->pParam->B2beta40 *= CoxWoverL;
	        here->pParam->B2beta4B *= CoxWoverL;
	        here->pParam->B2beta4G *= CoxWoverL;

                here->pParam->B2ua0 = model->B2ua00 + model->B2ua0W * Inv_W
			       + model->B2ua0L * Inv_L;
                here->pParam->B2uaB = model->B2uaB0 + model->B2uaBW * Inv_W
			       + model->B2uaBL * Inv_L;
                here->pParam->B2ub0 = model->B2ub00 + model->B2ub0W * Inv_W
			       + model->B2ub0L * Inv_L;
                here->pParam->B2ubB = model->B2ubB0 + model->B2ubBW * Inv_W
			       + model->B2ubBL * Inv_L;
                here->pParam->B2u10 = model->B2u100 + model->B2u10W * Inv_W
			       + model->B2u10L * Inv_L;
                here->pParam->B2u1B = model->B2u1B0 + model->B2u1BW * Inv_W
			       + model->B2u1BL * Inv_L;
                here->pParam->B2u1D = model->B2u1D0 + model->B2u1DW * Inv_W
			       + model->B2u1DL * Inv_L;
                here->pParam->B2n0 = model->B2n00 + model->B2n0W * Inv_W
			       + model->B2n0L * Inv_L;
                here->pParam->B2nB = model->B2nB0 + model->B2nBW * Inv_W
			       + model->B2nBL * Inv_L;
                here->pParam->B2nD = model->B2nD0 + model->B2nDW * Inv_W
			       + model->B2nDL * Inv_L;
                if (here->pParam->B2n0 < 0.0)
	            here->pParam->B2n0 = 0.0;

                here->pParam->B2vof0 = model->B2vof00 
			       + model->B2vof0W * Inv_W
			       + model->B2vof0L * Inv_L;
                here->pParam->B2vofB = model->B2vofB0 
			       + model->B2vofBW * Inv_W
			       + model->B2vofBL * Inv_L;
                here->pParam->B2vofD = model->B2vofD0 
			       + model->B2vofDW * Inv_W
			       + model->B2vofDL * Inv_L;
                here->pParam->B2ai0 = model->B2ai00 + model->B2ai0W * Inv_W
			       + model->B2ai0L * Inv_L;
                here->pParam->B2aiB = model->B2aiB0 + model->B2aiBW * Inv_W
			       + model->B2aiBL * Inv_L;
                here->pParam->B2bi0 = model->B2bi00 + model->B2bi0W * Inv_W
			       + model->B2bi0L * Inv_L;
                here->pParam->B2biB = model->B2biB0 + model->B2biBW * Inv_W
			       + model->B2biBL * Inv_L;
                here->pParam->B2vghigh = model->B2vghigh0 
			       + model->B2vghighW * Inv_W
			       + model->B2vghighL * Inv_L;
                here->pParam->B2vglow = model->B2vglow0 
			       + model->B2vglowW * Inv_W
			       + model->B2vglowL * Inv_L;

		here->pParam->CoxWL = model->B2Cox * EffectiveLength
			       * EffectiveWidth * 1.0e4;
		here->pParam->One_Third_CoxWL = here->pParam->CoxWL / 3.0;
		here->pParam->Two_Third_CoxWL = 2.0 
			       * here->pParam->One_Third_CoxWL;
                here->pParam->B2GSoverlapCap = model->B2gateSourceOverlapCap 
			       * EffectiveWidth;
                here->pParam->B2GDoverlapCap = model->B2gateDrainOverlapCap 
			       * EffectiveWidth;
                here->pParam->B2GBoverlapCap = model->B2gateBulkOverlapCap 
			       * EffectiveLength;
                here->pParam->SqrtPhi = sqrt(here->pParam->B2phi);
		here->pParam->Phis3 = here->pParam->SqrtPhi
		               * here->pParam->B2phi;
                here->pParam->Arg = here->pParam->B2betasB
			       - here->pParam->B2beta0B - model->B2vdd
	                       * (here->pParam->B2beta3B - model->B2vdd
	                       * here->pParam->B2beta4B);


             }


            /* process drain series resistance */
            if( (here->B2drainConductance=model->B2sheetResistance *
                    here->B2drainSquares) != 0.0 ) {   
                here->B2drainConductance = 1. / here->B2drainConductance ;
            } 
                   
            /* process source series resistance */
            if( (here->B2sourceConductance=model->B2sheetResistance *
                    here->B2sourceSquares) != 0.0 ) { 
                here->B2sourceConductance = 1. / here->B2sourceConductance ;
            }


            here->pParam->B2vt0 = here->pParam->B2vfb 
			  + here->pParam->B2phi 
			  + here->pParam->B2k1 * here->pParam->SqrtPhi 
			  - here->pParam->B2k2 * here->pParam->B2phi;
            here->B2von = here->pParam->B2vt0; /* added for initialization*/
        }
    }
    return(OK);
}  



