-------OpenCFD-EC  (c) by Li Xinliang, lixl@imech.ac.cn-------------------

编译方法： 
         主程序编译，运行make即可 （如有问题，请修改makefile).
            注： 需要MPI库及Fortran 90 编译器支持；
         
         CGNS网格转换器编译:  ifort -O3 -o convert-cgns-inp-1.1a.out convert-cgns-inp-1.1a.f90  -lcgns
            注： 需要CGNS库， 下载地址 http://cgns.sourceforge.net/
                安装方法（以cgnslib_3.1.4为例）： 
                 1) 下载原始文件 （例如cgnslib_3.1.4.tar.gz), 解压
                 2） 进入 ./cgnslib_3.1.4/src
                 3) ./configure --with-fortran
                 4) make  (或make all, 编译全部工具)
                 5) sudo make install (需要超级用户口令）
                   （该命令将头文件(*.h)拷贝到./usr/local/include 文件夹下面； 将库文件 （libcgns.a）拷贝到/usr/local/lib文件夹下）
                   如无超级用户口令，可手工将 *.h及 libcgns.a 拷贝到用户指定的目录。

运行方法：
          1) 准备好网格文件Mesh3d.dat (PLOT3D格式， format与unformat 均可； 即Gridgen的 Generic格式）
          2）准备好边界条件文件 bc3d.inp  (Gridgen 的Generic 格式）
          3) 如使用CGNS格式的网格文件 (grid.cgns)，可用convert-cgns-1.1a.out 将其转换为Mesh3d.dat及bc3d.inp 
          4）编辑修改控制文件control.ec 设置 流动参数及计算方法等；
             其中Ma,Re 分别为Mach数及Reynolds数；  
             AoA 为攻角 ， AoS为侧滑角 ; (单位为度)

          5）运行 mpirun -np xx ./opencfd-1.04.out （或用bsub 等方式提交）
            xx 可为1,2,3,4, ... 但不能超过网格的总块数。
          
          注1：
           如无wall_dist.dat (记录每个网格点到最近壁面的距离)， 则程序开始时先生成该文件，需要消耗一些时间；
           运行时输出残差到Residual.dat （各方程的最大残差即均方根残差）
           运行时输出整个物体承受的气动力6分量到force.log 
              各列含义： 计算步， CL, CD, CS, CMx, CMy, CMz 
			  (升力系数、阻力系数、侧向力系数，三个方向的力矩系数， Ver1.04起，均按照通常定义无量纲化）。

           对于定常问题，请根据残差（或气动力）的收敛情况，判断是否可结束计算。

         注2：
            如使用OpenMP, 则需要：
                1）在makefile中打开openmp开关（设定 opt=-O3 -openmp），然后重新编译。         
                2）在控制文件中设定NUM_THREADS的数目 （例如NUM_THREADS=4， 设定为1为单线程串行）。


后处理方法： 运行 readflow3d-ver1.7.f90； 该程序将以tecplot格式输出壁面及对称面上的二维场，以及全流场的三维场。 

CGNS网格处理方法：
      本版通过转换器本支持CGNS网格，具体方法：
       1） 准备好CGNS格式的网格文件 grid.cgns
       2） 运行 ./convert-cgns-1.1a.out (读入grid.cgns, 生成Mesh3d.dat 及bc3d.inp)


---------------------相关输入数据-----------------------------------------------------
 以下文件夹中包含了一些示例计算，为了减少数据量，文件夹中并不包含网格文件，网格文件请到指定位置下载。
 
 1） ./case/DLR-F4:   DLR-F4 翼身组合体 （粗网格）
     具体使用方法见该目录下的readme.txt
   
 2） ./case/M6-wing:   ONEAR M6-wing 三维翼 （粗网格）
      control.ec 主控文件
      bc3d.inp 边界条件文件  
      网格文件(Mesh3d.dat) 请到作者的网盘下载 （https://skydrive.live.com/?cid=1cc0dcbff560c149, 公开->OpenCFD-EC->算例->M6-Wing）
      该网格为4块，因而MPI并行运行，最多不超过4进程 （推荐2进程）。



!-----------------------------
 版本更新记录：
  0.98c: 改写了SA模型部分，采用类似无量纲计算（含Re）， 参考《CFL3D手册》 （与CFL3D无量纲化有所区别，因而不包含Mach数）；
  0.99: 改写了滤波部分
  1.00: 支持自定义边界条件
  1.01: 增加了边界格式功能（允许边界使用低阶精度格式）
  1.02: 增加了WENO7格式
  1.03: 对控制文件的格式进行了修改支持 NAMELIST格式
  1.04: 对输出进行了修改（报错信息，气动力、矩信息）
  1.05: 对差分格式进行了修改； 去掉了WENO3格式，补充了 MUSCL2 with k=-1 (迎风型); 原MUSCL2 是k=1的中心型MUSSCL
  1.05a: 对压力、密度的增量进行了限制(类似CFL3D的做法)，以增强稳定性；
  1.05b: 修改了角点处理方式 (角点薄层近似)
  1.05d: 在压力、密度限制方面进行了修改，（局部平均）
  1.06:  采用限制措施（限制最低、最高密度、压力及速度）
  1.06a: 对SA等扩展方程也进行了限制；
  1.07:  对OMUSCL2进行了修改， 添加了New-SA模型；
  1.07a: new SA模型的系数可在控制文件中修改；
  1.1:  增加了叶轮机械（压气机、涡轮等）计算功能，使用旋转坐标系 （含源项）
  1.1a: 修正了背压定义的Bug  （压力应当用动压无量纲化，而不是用来流总压）
  1.11：修正了与X,Y,Z坐标定义及计算CL,CD的Bug, 使得对于Z-坐标轴向上的网格也能输出正确的气动力系数 （需设置Cood_Y_UP=0)
  1.12: 恢复对差分-有限体积混合方法的支持；
  1.13: 增加了内流模式；与叶轮机械模式类似(但不含旋转);
  1.13b: 修改了Roe 格式中的Bug;
  1.13c：添加了时间相关入口边界条件 （读取时间序列）， 修正了Runge-Kutta 3格式中OpenMP bug， 支持时间平均
  1.14:  增加了数据文件output方式（可写入不同的文件名）；
  1.14a: 添加了含吹吸气扰动的壁面边界条件 （作为用户自定义边界条件）
  1.15： A Bug in comput_force( ) is removed;  修正了后处理粘性力计算的Bug （i-, j-, k- 面 cf应为负！）
  1.15a: 修正了call output_vt() 中的Bug;  修正了readflow3d-ver2.4.f90中的Bug (read wall_dist.dat)；
         对边界条件进行了修改，  Bc%bc .eq. BC_Inflow .and. IF_InnerFlow .eq. 0 时， 采用完全给定来流的边界条件 （不区分超、亚声速）；
  1.16： 提升程序鲁棒性。  如发现某块出现物理量超限（如负温度），则强制该块采用1阶迎风格式计算，该块时间步降为原先的1/10， 该块不计算粘性项；
   1.16a:  sub_turbulence_SST 中的Bug 被修正；