function ChargeDensity_deVries()
{
	this.cf = new CommonFun();
	
	this.Debugging = false;
	this.N = 0;
	this.Nref = 0;
	this.ShortReferences = [];
	this.FullReferences =[];
	this.Z = -1;
	this.A = -1;
	this.indexes = [];
	this.density = [];
//----------------------------------------------------------------------	
// все индексы в имеющемся списке, соответствующие данному ядру
this.getIndex = function(z,a)
{
    var n = 0,i;
	this.Z = z;
	this.A = a;	
    this.indexes = [];
    for(i=0;i<this.N;i++)
    {
        if(z>0 && z==this.density[i].Z && a>0 && a==this.density[i].A) 
		{
			if(this.Debugging) console.log(i+" "+this.density[i].ShortReference);
			this.indexes[n]=i; n++;
		}
        if(z<0 && a>0 && a==this.density[i].A) {this.indexes[n]=i; n++;}
        if(a<0 && z>0 && z==this.density[i].Z) {this.indexes[n]=i; n++;}
    }
    if(n==0) this.indexes = null;
}
//-----------------------------------------------------------------------------
// содержимое файла density-deVries.txt во входном параметре
this.readDensity = function(deVries)
{
    var z,a,i,ii,
        r2,c_a,z_alpha,w,qmin,qmax,rp,
        Ri = [],
        Qi = [],
        Ai = [],
    	ref,
		type,
		str1,str2,str3;
	if(this.Debugging) console.log("Reading the density-deVries.txt");
	//-----------------------------------------------------------------------
	str1 = deVries.split("\n");// разбили построчно
	str2 = str1[0].trim();     // анализируем 1ую строку чтобы узнать число записей
	str3 = str2.split(" ");
	this.N = str3[0]; if(this.Debugging) console.log("Число записей в deVries: " + this.N);
	if(this.N <= 0) return false;
	
	this.density = [];	
	for(i=0;i<this.N;i++)
	{
		this.density[i] = new NuclearChargeDensity();
		
		str2 = this.cf.fullTrim(str1[i+2]);
		str3 = str2.split(" ");
		
		z = str3[0];
		a = str3[1];
		type = str3[2];;

		this.density[i].setNucleus(z,a);
		this.density[i].setType(type);
		switch(this.density[i].Type)
		{
			default: break;
			case 0: // Harmonic oscilator
			case 1: // Modified Harmonic oscilator
			case 5: // 2 parameters Fermi model
			case 6: // 3 parameters Fermi model
			case 7: // 3 parameters Gauss model
				r2 = str3[3];
				c_a = str3[4];
				z_alpha = str3[5];
				w = str3[6];
				qmin = str3[7];
				qmax = str3[8];
				ref = str3[9];
					
				this.density[i].setUniformParameters(r2,c_a,z_alpha,w,qmin,qmax);
				this.density[i].setReference(ref);
               	break;
            case 2: // Model independent evaluation
            case 3: // Fourier-Bessel
				r2 = str3[3];
				rp = str3[4];
				for(ii=0;ii<17;ii++) Ai[ii] = str3[5+ii];
				ref = str3[22];				
               	this.density[i].setFourierBesselParameters(r2,rp,Ai);
               	this.density[i].setReference(ref);
				break;
            case 4: // Gauss sum
				r2 = str3[3];
				rp = str3[4];
				for(ii=0;ii<12;ii++)
				{
					Ri[ii] = str3[5+2*ii];
					Qi[ii] = str3[6+2*ii];
				}
				ref = str3[29];				
				this.density[i].setGaussParameters(r2,rp,Ri,Qi);
				this.density[i].setShortReference(ref);
				break;
           	case 8: // Uniform Gauss model
				break;
        }
	}
	if(this.Debugging) console.log("Densities were read!");
	return true;
}
this.readReferences = function(refFile)
{
	var i,ii,j,str1,str2,str3;

	str1 = refFile.split("\n");// разбили построчно
	str2 = str1[0].trim();     // анализируем 1ую строку чтобы узнать число записей
	str3 = str2.split(" ");
	this.Nref = str3[0]; if(this.Debugging) console.log("Number of references: " + this.Nref);
	if(this.Nref <= 0) return false;
	
	this.ShortReferences = [];
	this.FullReferences = [];

	if(this.Debugging) console.log("\nStrat to read references ...");
	for(i=1;i<=this.Nref;i++)
	{			
		str2 = this.cf.fullTrim(str1[i]);
		str3 = str2.split(" ");

		this.ShortReferences[i] = str3[0]; 
		this.FullReferences[i] = "";
		for(ii=1;ii<str3.length;ii++) this.FullReferences[i] += str3[ii]+" ";
	}
	if(this.Debugging) console.log("... references were read");
	//-------------------------------------------------------------------
	if(this.density != null)
	{
		for(i=0;i<this.N;i++)
		{
			this.density[i].IndexRef = -1;
			for(j=0;j<this.Nref;j++) 
				if(this.ShortReferences[j] == this.density[i].ShortReference) 
				{
					this.density[i].Reference = this.FullReferences[j];
					this.density[i].IndexRef = j;
					break;
				}
			if(this.density[i].IndexRef < 0)
				console.log("Unknown reference: "+this.density[i].ShortReference);
		}
	}
	//-------------------------------------------------------------------
	return true;
}
}