#include "global_var.h" 
void default_setting()
{ 
  int i;
  nfmax=10;              
  np=11;
  nrho=12;
  ngam=13;
    /****************
	     nf[i][j][0]~到f[i][j][nfmax]是這樣分配的:
		 nf[i][j][1]: x方向的速度u(sorry,現在我的認識水平只停留在直角坐標系上） 
         nf[i][j][2]: y方向的速度v
         nf[i][j][3]: 壓力修正值p'
         nf[i][j][4]: 溫度t
		 nf[i][j][5]~nf[i][j][10]: 其它φ變量
         nf[i][j][11]: 壓力p
		 nf[i][j][12]: 密度ρ
         nf[i][j][13]: 擴散系數Γ
   *********************/
  lstop=0; //bool .false.
  for (i=1;i<=nfmax;i++) 
  {
	  lsolve[i]=0; //不求解nf[i][j][0]~到f[i][j][nfmax]            
	  lblk[i]=1;   //進行塊修正                                    
	  ntimes[i]=1; //在solve子程序中求解變量f[i][j][nf]的掃描一次  
  }
  for (i=1;i<=13;i++)
  {
	  lprint[i]=0;  //不打印中間結果
	  relax[i]=1;   //relax=1其實就是不松弛了                      
  }
  mode=1;           //直角坐標系
  last=5;           //迭代5次                                      
  time=0;           //開始計算的時間點		 
  dt=1.0e10;        //時間步長                                     	 
  ipref=1;
  jpref=1;          //壓力參考點的(i,j)為(1,1)
  rhocon=1;         //不可壓縮流的密度為1
}