Use of Hold in DoNewGlobalFit

function ngf(base0,base1,templist,cfstr0,cfvar0,cfstr1,cfvar1,cfstr2,cfvar2,cfstr3,cfvar3,cfstr4,cfvar4)

// base0 and base1 are the basenames of the two data sets (without temperature)

string base0, base1, templist
string cfstr0
variable cfvar0
string cfstr1
variable cfvar1
string cfstr2
variable cfvar2
string cfstr3
variable cfvar3
string cfstr4
variable cfvar4

variable i, numdatapts

// name of data
	make/o/n=(2,2)/t datasets
	DataSets[0][0]=base0+"_"+stringfromlist(i,templist)
	DataSets[1][0]=base1+"_"+stringfromlist(i,templist)
	DataSets[0][1]="_calculated_"
	DataSets[1][1]="_calculated_"
	
// extract number of points in data wave	
wave sampledataset = $datasets[0][0]
numdatapts = numpnts(sampledataset)

// format: [row][column]
// a row for each data set
// # columns: (max # of coefficients used by any fit functions) + 4
	make/o/n=(2,8) CoefDataSetLinkage
	
	
// fit functions
	make/o/t/n=2 fitfuncnames
	fitfuncnames[0]="sigma1ff"
	fitfuncnames[1]="eps1ff"
	
// column 0: contains index in the fitfuncnames where the name of a fitting function will be found
	CoefDataSetLinkage[0][0]=0
// column 1: point # within concatenated data set where this data set begins 	
	CoefDataSetLinkage[0][1]=0
// column 2: point # where a given data set ends (this number is not used...)	
	CoefDataSetLinkage[0][2]=numdatapts-1
// column 3: how many fit coefficents are used by the fit function for the corresponding data set
	CoefDataSetLinkage[0][3]=4
// column 4: index to the master coefficient wave 	
	CoefDataSetLinkage[0][4]=0
	CoefDataSetLinkage[0][5]=1
	CoefDataSetLinkage[0][6]=2
	CoefDataSetLinkage[0][7]=3
	
// column 0: contains index in the fitfuncnames where the name of a fitting function will be found
	CoefDataSetLinkage[1][0]=0
// column 1: point # within concatenated data set where this data set begins 	
	CoefDataSetLinkage[1][1]=numdatapts
// column 2: point # where a given data set ends (this number is not used...)	
	CoefDataSetLinkage[1][2]=2*numdatapts-1
// column 3: how many fit coefficents are used by the fit function for the corresponding data set
	CoefDataSetLinkage[1][3]=4
// column 4: index to the master coefficient wave 	
	CoefDataSetLinkage[1][4]=4
	CoefDataSetLinkage[1][5]=1
	CoefDataSetLinkage[1][6]=2
	CoefDataSetLinkage[1][7]=3


// coefwave contains initial guesses
// inidices established in coefdatasetlinkage are repeated here
// i.e. the length of coefwave will be the total number of variables, not just the number
// used by any particular fit function
	make/o/n=(5,2) coefwave
	setdimlabel 1, -1, Hold, coefwave
	coefwave[0][0]=cfvar0
	coefwave[1][0]=cfvar1
	coefwave[2][0]=cfvar2
	coefwave[3][0]=cfvar3
	coefwave[4][0]=cfvar4
// set coefwave[n][1]=0 to not hold n, 1 to hold n	
	coefwave[][1]=1
	

	
	
// coefnames provides strings for the variable names
// same index convention as coefwave and coefdatasetlinkage	
	make/o/t/n=5 coefnames
	coefnames[0]=cfstr0
	coefnames[1]=cfstr1
	coefnames[2]=cfstr2
	coefnames[3]=cfstr3
	coefnames[4]=cfstr4

donewglobalfit(FitFuncNames, DataSets, CoefDataSetLinkage, CoefWave, CoefNames, $"",NewGFOptionAPPEND_RESULTS, 200, 1)
	
end

function eps1ff(w,x) : FitFunc
	wave w
	variable x

	//CurveFitDialog/ These comments were created by the Curve Fitting dialog. Altering them will
	//CurveFitDialog/ make the function less convenient to work with in the Curve Fitting dialog.
	//CurveFitDialog/ Independent Variables 1
	//CurveFitDialog/ x
	//CurveFitDialog/ Coefficients 4
	//CurveFitDialog/ w[0] = eps_inf 
	//CurveFitDialog/ w[1] = sj1
	//CurveFitDialog/ w[2] = fj1
	//CurveFitDialog/ w[3] = gj1

//	return real(8.854*4*pi*2*pi*x*(w[0]-1+w[1]*(2*pi*w[2])^2 / ((2*pi*w[2])^2-(2*pi*x)^2-sqrt(-1)*w[3]*x*(2*pi)^2)+ w[4]*(2*pi*w[5])^2 / ((2*pi*w[5])^2-(2*pi*x)^2-sqrt(-1)*w[6]*x*(2*pi)^2)  )  / (100*4* pi * sqrt(-1)))
       return real(w[0]-1+w[1]*(2*pi*w[2])^2 / ((2*pi*w[2])^2-(2*pi*x)^2-sqrt(-1)*w[3]*x*(2*pi)^2))
end


function sigma1ff(w,x) : FitFunc
	wave w
	variable x

	//CurveFitDialog/ These comments were created by the Curve Fitting dialog. Altering them will
	//CurveFitDialog/ make the function less convenient to work with in the Curve Fitting dialog.
	//CurveFitDialog/ Independent Variables 1
	//CurveFitDialog/ x
	//CurveFitDialog/ Coefficients 4
	//CurveFitDialog/ w[0] = sigma0
	//CurveFitDialog/ w[1] = sj1
	//CurveFitDialog/ w[2] = fj1
	//CurveFitDialog/ w[3] = gj1	

	return real( 8.854*4*pi*2*pi*x *(w[0]+w[1]*(2*pi*w[2])^2 / ((2*pi*w[2])^2-(2*pi*x)^2-sqrt(-1)*w[3]*x*(2*pi)^2)))
end