Check koeff.f90 - Performing checking of koefficients
From CFD-Wiki
(Difference between revisions)
(2 intermediate revisions not shown) | |||
Line 1: | Line 1: | ||
- | + | <pre> | |
+ | ! 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' | include 'icomm_1.f90' | ||
Dimension F_out(nx,ny) | Dimension F_out(nx,ny) | ||
- | |||
Character Filename*10 | Character Filename*10 | ||
- | |||
!---------------------------------------------------------------- | !---------------------------------------------------------------- | ||
- | |||
NImax = NXmaxp | NImax = NXmaxp | ||
- | |||
NJmax = NYmaxp | NJmax = NYmaxp | ||
Line 20: | Line 35: | ||
Call Out_array(F_out,NImax,NJmax,Filename) | Call Out_array(F_out,NImax,NJmax,Filename) | ||
!------------------------------------------------------------------- | !------------------------------------------------------------------- | ||
- | |||
!---------------------------------------------------------------- | !---------------------------------------------------------------- | ||
- | |||
NImax = NXmaxp | NImax = NXmaxp | ||
- | |||
NJmax = NYmaxp | NJmax = NYmaxp | ||
Line 32: | Line 44: | ||
Call Out_array(F_out,NImax,NJmax,Filename) | Call Out_array(F_out,NImax,NJmax,Filename) | ||
- | |||
!------------------------------------------------------------------- | !------------------------------------------------------------------- | ||
- | |||
!---------------------------------------------------------------- | !---------------------------------------------------------------- | ||
- | |||
NImax = NXmaxp | NImax = NXmaxp | ||
- | |||
NJmax = NYmaxp | NJmax = NYmaxp | ||
F_out = Aw | F_out = Aw | ||
- | |||
Filename ='A_3ww.txt' | Filename ='A_3ww.txt' | ||
Call Out_array(F_out,NImax,NJmax,Filename) | Call Out_array(F_out,NImax,NJmax,Filename) | ||
- | |||
!------------------------------------------------------------------- | !------------------------------------------------------------------- | ||
- | |||
!---------------------------------------------------------------- | !---------------------------------------------------------------- | ||
- | |||
NImax = NXmaxp | NImax = NXmaxp | ||
- | |||
NJmax = NYmaxp | NJmax = NYmaxp | ||
Line 62: | Line 65: | ||
Call Out_array(F_out,NImax,NJmax,Filename) | Call Out_array(F_out,NImax,NJmax,Filename) | ||
- | |||
!------------------------------------------------------------------- | !------------------------------------------------------------------- | ||
- | |||
!---------------------------------------------------------------- | !---------------------------------------------------------------- | ||
- | |||
NImax = NXmaxp | NImax = NXmaxp | ||
- | |||
NJmax = NYmaxp | NJmax = NYmaxp | ||
Line 76: | Line 75: | ||
Call Out_array(F_out,NImax,NJmax,Filename) | Call Out_array(F_out,NImax,NJmax,Filename) | ||
- | |||
!------------------------------------------------------------------- | !------------------------------------------------------------------- | ||
- | |||
- | |||
Return | Return | ||
+ | End | ||
- | + | </pre> |
Latest revision as of 14:42, 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