// de_drag.cpp: implementation of the de_drag class.
//
//////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "2level.h"
#include "de_dragmodel.h"

#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////

de_dragmodel::de_dragmodel()
{


}

de_dragmodel::~de_dragmodel()
{

}
double de_dragmodel::de_drag(int m,double rH)
{
	double z0;
	double z1;
	double z2;
	double z3;
	double dd;
	double de_drag;
	//double z[9][4];
	
   //z[m][j]=(65.150000,00.935500,00.009200,00.230000,
   //        82.260000,02.019000,00.008200,00.230000,
	//		 65.150000,00.935500,00.009200,00.230000,
	//		 63.713500,01.112670,00.711273,00.444762,
	//		 76.398600,01.000000,00.267355,00.316020,
	//		 71.237100,01.000000,00.620330,00.425969,
	//		 33.413800,01.000000,00.431080,00.562023,
	//		 68.767700,01.000000,01.039400,00.443062,
	//		 84.151700,00.991309,00.922673,00.395288);
	if (m=0)
	{
		z0=65.150000;
		z1=00.935500;
		z2=00.009200;
		z3=00.230000;
	}
	if (m=1)
	{
		z0=82.260000;
		z1=02.019000;
		z2=00.008200;
		z3=00.230000;
	}
	if (m=2)
	{
		z0=65.150000;
		z1=00.935500;
		z2=00.009200;
		z3=00.230000;
	}
	if (m=3)
	{
		z0=63.713500;
		z1=01.112670;
		z2=00.711273;
		z3=00.44476;
	}
	
	
	
	
	
	//		 84.151700,00.991309,00.922673,00.395288);
	if (m=4)
	{
		z0=76.398600;
		z1=01.000000;
		z2=00.267355;
		z3=00.316020;
	}
	if (m=5)
	{
		z0=71.237100;
		z1=01.000000;
		z2=00.431080;
		z3=00.425969;
	}
	if (m=6)
	{
		z0=33.413800;
		z1=01.000000;
		z2=00.431080;
		z3=00.562023;
	}
	if (m=7)
	{
		z0=68.767700;
		z1=01.000000;
		z2=01.039400;
		z3=00.443062;
	}
	if (m=8)
	{
		z0=84.151700;
		z1=00.991309;
		z2=00.922673;
		z3=00.395288;
	}
	
	
	dd=z3*log(z0*(log10((z1+rH)/(z1-rH))+z2));
	de_drag=exp(dd)*100;	
	return de_drag;



}