Main.f90 - Main modul
From CFD-Wiki
(Difference between revisions)
(fix the bug) |
|||
Line 3: | Line 3: | ||
!Sample program for solving Smith-Hutton Test using different schemes | !Sample program for solving Smith-Hutton Test using different schemes | ||
!of covective terms approximation - Main modul | !of covective terms approximation - Main modul | ||
- | !Copyright (C) 2005 Michail | + | !Copyright (C) 2005 Michail Kiričkov |
!This program is free software; you can redistribute it and/or | !This program is free software; you can redistribute it and/or | ||
Line 30: | Line 30: | ||
Call Ini | Call Ini | ||
Call Grid_rectangular | Call Grid_rectangular | ||
- | |||
Call Geom | Call Geom | ||
Call Init_all | Call Init_all |
Latest revision as of 17:54, 20 April 2012
!Sample program for solving Smith-Hutton Test using different schemes !of covective terms approximation - Main modul !Copyright (C) 2005 Michail Kiričkov !This program is free software; you can redistribute it and/or !modify it under the terms of the GNU General Public License !as published by the Free Software Foundation; either version 2 !of the License, or (at your option) any later version. !This program is distributed in the hope that it will be useful, !but WITHOUT ANY WARRANTY; without even the implied warranty of !MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the !GNU General Public License for more details. !You should have received a copy of the GNU General Public License !along with this program; if not, write to the Free Software !Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. !*************************************************************************** Program Main include 'icomm_1.f90' Dimension F_out(nx,ny) Character Filename*10 Call Ini Call Grid_rectangular Call Geom Call Init_all Call Control_gem Do 100 Niter=1,100 write(*,*) '------------------------',Niter,'------------------------',Niter Call Coef_1(5) Call Check_koeff ! Alfa = 0.95 ! Call Under_Relaxation(Alfa,5) Call TDMA_1(5) Call Bound_condition Call Convergence_Criteria(5) 100 continue !---------------------------------------------------------------- NImax = NXmaxp NJmax = NYmaxp F_out = F(:,:,1) Filename ='1_U_s.txt' ! Call Out_array(F_out,NImax,NJmax,Filename) !------------------------------------------------------------------- !---------------------------------------------------------------- NImax = NXmaxp NJmax = NYmaxp F_out = F(:,:,2) Filename ='1_V_s.txt' ! Call Out_array(F_out,NImax,NJmax,Filename) !------------------------------------------------------------------- !---------------------------------------------------------------- NImax = NXmaxp NJmax = NYmaxp F_out = F(:,:,5) Filename ='1_T_s.txt' ! Call Out_array(F_out,NImax,NJmax,Filename) !------------------------------------------------------------------- Call output_all WRITE(*,*) 'PRIVET' STOP END