Check koeff.f90 - Performing checking of koefficients (output for visual control)
From CFD-Wiki
(Difference between revisions)
Latest revision as of 17:20, 21 September 2005
! Sample program for solving Smith-Hutton Test using different schemes ! of covective terms approximation - Coefficients cheking modul ! Copyright (C) 2005 Michail Kirichkov ! 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. !*********************************************************************** Subroutine Check_koeff include 'icomm_1.f90' Dimension F_out(nx,ny) Character Filename*10 !---------------------------------------------------------------- NImax = NXmaxp NJmax = NYmaxp F_out = Ap Filename ='A_ppp.txt' Call Out_array(F_out,NImax,NJmax,Filename) !------------------------------------------------------------------- !---------------------------------------------------------------- NImax = NXmaxp NJmax = NYmaxp F_out = Ae Filename ='A_4ee.txt' Call Out_array(F_out,NImax,NJmax,Filename) !------------------------------------------------------------------- !---------------------------------------------------------------- NImax = NXmaxp NJmax = NYmaxp F_out = Aw Filename ='A_3ww.txt' Call Out_array(F_out,NImax,NJmax,Filename) !------------------------------------------------------------------- !---------------------------------------------------------------- NImax = NXmaxp NJmax = NYmaxp F_out = An Filename ='A_2nn.txt' Call Out_array(F_out,NImax,NJmax,Filename) !------------------------------------------------------------------- !---------------------------------------------------------------- NImax = NXmaxp NJmax = NYmaxp F_out = As Filename ='A_1ss.txt' Call Out_array(F_out,NImax,NJmax,Filename) !------------------------------------------------------------------- Return End