/*
* exam6.c : example 6 for SIMPLE
*/
/*
* 注意C與FORTRAN的差別：
* 1. C中的數組下標是從0開始的,而FORTRAN則是從1開始的
*   因此本程序中的一些量(如L1,L2,L3,M1,M2,M3,NFMAX,nGam,nP,nRho,IPref,JPref)比FORTRAN程序中的小1
*   所有數組的下標從0開始,程序中數組下標中出現的數字比FORTRAN中的小1
* 2. C的數組是按行優先讀的,而FORTRAN的則是按列讀的
*   因此本程序中的F的定義不同,是F[NFMAX][I][J]
* 3. C用指針來實現FORTRAN的EQUIVALENCE功能。
*   如定義double (*U)[NSIZE],U=F[0],則U[i][j]等同于F[0][i][j];
*/

#include "staincl.h"
#include "extern.h"

#define PI 3.14159
/*global variable for user */
	double (*T)[NSIZE];
	double Tin;
	double Tw;
	double TRef;
	double Vin;
	double Vout;
	double Pr;
	double Amu;
	double Amup;
	double RhoRef;
	double RhoT;
	double Flowin;
	double Fl;

/*******************************************************************************
* functon name:	fnGrid
* description:	in this file you should do:
*			--initialize variable you defined
*			--set the filename of the output file, the step can be ignored 
*			--set the title of the variable that will be calculated
*			--set which variable will be calc and print
*			--set axis mode-nMode: 1-cartesian, 2-axisymmetric, 3-polar
*			--set the number of steps
*			--set the number of grids
*			--devide grids
********************************************************************************/
void fnGrid()
{
	int i;

	T=F[3];
	szTitle[0]="U";
	szTitle[1]="V";
	szTitle[2]="Str Fun";
	szTitle[3]="temprature";
	szTitle[10]="pressure";
	szTitle[11]="Density";
	Relax[0]=0.5;
	Relax[1]=0.5;
	Relax[10]=0.8;
	for( i=0;i<=3;i++)
	{
		bSolve[i]=TRUE;
		bPrint[i]=TRUE;
	}
	bPrint[10]=TRUE;
	bPrint[11]=TRUE;
	iLast=30;
	XL=0.5;
	YL=2.0;
	L1=6;//
	M1=11;//

	fnUGrid();
}

/*******************************************************************************
* functon name:	fnStart
* description:	in this file you should do:
*			--initialize the 
********************************************************************************/
void fnStart()
{
	int i,j;
	Tin=500;
	Tw=300;
	Vin=100;
	Vout=Vin*XCV[L2]/X[L1]*Tw/Tin;
	for( i=0;i<=L1;i++)
	{
		for( j=0;j<=M1;j++)
		{
			U[i][j]=0;
			V[i][j]=Vout;
			V[i][1]=0;
			V[0][j]=0;
			T[i][j]=Tw;
		}
	}
	V[L2][1]=Vin;
	T[L2][0]=Tin;
	Pr=0.7;
	Amu=1;
	Amup=Amu/Pr;
	TRef=300;
	RhoRef=1;
	RhoT=RhoRef*TRef;
}

/*******************************************************************************
* functon name:	fnDense
* description:	in this file you should do:
*			--calc the density if needed
********************************************************************************/
void fnDense()
{
	int i,j;
	for( i=0;i<=L1;i++)
		for( j=0;j<=M1;j++)
			Rho[i][j]=RhoT/T[i][j];

}

/*******************************************************************************
* functon name:	fnbound
* description:	in this file you should do:
*			--calc the boundary
********************************************************************************/
void fnBound()
{
	int i;
	double fac;
	if(iIter==0)
		Flowin=Rho[L2][0]*V[L2][1]*XCV[L2];
	Fl=0;
	for( i=1;i<=L2;i++)
		Fl+=Rho[i][M1]*V[i][M2]*XCV[i];
	fac=Flowin/Fl;
	for(i=1;i<=L2;i++)
	{
		V[i][M1]=V[i][M2]*fac;
		T[i][M1]=T[i][M2];
	}
}

/*******************************************************************************
* functon name:	fnOutput
* description:	in this file you should do:
*			--print calc info to terminal or write to file
********************************************************************************/
void fnOutPut()
{
	if(iIter==0)
		fprintf(pFileOut,"\n    Iter    SMax    SSum      V[4][7]      T[4][7]");
	fprintf(pFileOut,"\n   %4d  %10.3e  %10.3e  %10.3e  %10.3e ",iIter,SMax,SSum,V[3][6],T[3][6]);

	if(iIter==iLast)
		fnPrint();
}

/*******************************************************************************
* functon name:	fnGamsor
* description:	in this file you should do:
*			--set the gam 
*			--set the sc,sp
********************************************************************************/
void fnGamSor()
{
	int i,j;
	for( i=0;i<=L1;i++)
	{
		for( j=0;j<=M1;j++)
		{
			Gam[i][j]=Amu;
			if(nF==3)
				Gam[i][j]=Amup;
			if(nF!=0)
				Gam[L1][j]=0;
			Gam[i][M1]=0;
		}
	}
}