hyben.cpp

海洋油污扩散模型的计算机实现。简单了点

// HyBen.cpp: implementation of the CHyBen class.
//
//////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "Pollution.h"
#include "HyBen.h"
#include "Condation.h"
#include "math.h"
#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////
//CCondation condation;

extern CCondation condation;

CHyBen::CHyBen()
{
	//SYF20050305 m_fBenXDiffse=0.001;
	//SYF20050305 m_fBenYDiffse=0.002;
	m_fBenLeak=20;
	m_fLogitude=30;
	m_fDistanceMax=50;
	m_fDistanceMin=5;
	m_fSpreyDistance=5;
	m_fTime=0;
	//SYF20050305 m_fSlope=0.01;

	m_fBenXDiffse=0.05*condation.m_fWaterSpeed*condation.m_fWaterDepth/condation.m_fWaterSlope;
    m_fBenYDiffse=0.068871*pow(0.75,(condation.m_fWaterWidth*condation.m_fWaterDepth))*sqrt(condation.m_fWaterSlope);





}


CHyBen::~CHyBen()
{

}


float CHyBen::GetBenDensity(int x, int y, float m_fBenLeak, 
							float m_fLogitude, float m_fDistanceMax, float m_fDistanceMin)
{
    float _fDensity=0.25*m_fBenLeak*m_fLogitude*(m_fDistanceMax-m_fDistanceMin)/
		             (condation.m_fWaterDepth*m_fSpreyDistance*sqrt(condation.m_fEarthRotSpeed*m_fBenXDiffse*m_fBenYDiffse)*m_fTime);
	_fDensity*=exp(-(x-condation.m_fWaterSpeed)*(x-condation.m_fWaterSpeed)/(4*m_fBenXDiffse*m_fTime));
	return _fDensity;
}

void CHyBen::PrintBen(int nBeginX, int nBeginY, int m_nTime, CClientDC *pClientDC)

{

    float nColor[900][200],Max,Min,fTemp;
          Max=Min=fTemp=0;
	for (int i=0;i<900;i++)
	   for(int j=0;j<200;j++)
	   {
		   
		   nColor[i][j]=GetBenDensity(i+m_nTime,j+m_nTime,200,30.0,50,5);
		   if (nColor[i][j]>Max)
			   Max=nColor[i][j];
		   if(nColor[i][j]<Min)
			   Min=nColor[i][j];
		   
		   //pDC->SetPixel(nBeginY+i,nBeginZ+j,(255,0,0));
	   }
	   for (int p=0;p<900;p++)
	   for(int q=0;q<200;q++)
	   {
		   
		  fTemp=250-nColor[p][q]*255/(Max-Min);
		  if(fTemp<120&&fTemp>0) 
		  pClientDC->SetPixel(nBeginX+p-20, nBeginY+q-20, ((int)fTemp,(int)fTemp,(int)fTemp));
	   }


}