Check koeff.f90 - Performing checking of koefficients
From CFD-Wiki
(Difference between revisions)
Line 20: | Line 20: | ||
Subroutine Check_koeff | 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 40: | 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 52: | 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 82: | 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 96: | Line 75: | ||
Call Out_array(F_out,NImax,NJmax,Filename) | Call Out_array(F_out,NImax,NJmax,Filename) | ||
- | |||
!------------------------------------------------------------------- | !------------------------------------------------------------------- | ||
- | |||
- | |||
Return | Return | ||
- | |||
End | End | ||
- | + | </pre> |
Revision as of 20:39, 20 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