Grid sub.f90
From CFD-Wiki
| Line 1: | Line 1: | ||
| + | Sample program for solving Smith-Hutton Test using different schemes | ||
| + | of covective terms approximation - Grid's computing 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 Grid_rectangular | Subroutine Grid_rectangular | ||
| Line 16: | Line 35: | ||
! SLx = 1. | ! SLx = 1. | ||
| - | |||
! SLy = 1. | ! SLy = 1. | ||
! Xbeg = 0. | ! Xbeg = 0. | ||
| - | |||
! Ybeg = 0. | ! Ybeg = 0. | ||
NXmax = 121 | NXmax = 121 | ||
| - | |||
NYmax = 61 | NYmax = 61 | ||
| - | |||
NXmaxP = NXmax + 1 | NXmaxP = NXmax + 1 | ||
| - | |||
NYmaxP = NYmax + 1 | NYmaxP = NYmax + 1 | ||
| - | + | ||
!------------------------------------------------------------------------------------------- | !------------------------------------------------------------------------------------------- | ||
DO 2 I=1, NXmax | DO 2 I=1, NXmax | ||
| - | |||
DO 2 J=1, NYmax | DO 2 J=1, NYmax | ||
| Line 43: | Line 56: | ||
DO 3 J=1, NYmax | DO 3 J=1, NYmax | ||
| - | + | DO 3 I=1, NXmax | |
| - | + | ||
Y(I,J) = Ybeg + (SLy /(NYmax-1))* (j-1) | Y(I,J) = Ybeg + (SLy /(NYmax-1))* (j-1) | ||
3 continue | 3 continue | ||
| + | |||
| + | !------------------------------------------------------------------------------------------- | ||
open (22,file='GRID.dat') | open (22,file='GRID.dat') | ||
| Line 57: | Line 71: | ||
DO 1 J=1, NYmax | DO 1 J=1, NYmax | ||
| - | |||
DO 1 I=1, NXmax | DO 1 I=1, NXmax | ||
| Line 70: | Line 83: | ||
DO 4 J=1, NYmax | DO 4 J=1, NYmax | ||
| - | |||
DO 4 I=1, NXmax | DO 4 I=1, NXmax | ||
| Line 80: | Line 92: | ||
Return | Return | ||
| - | |||
End | End | ||
Revision as of 20:57, 20 September 2005
Sample program for solving Smith-Hutton Test using different schemes of covective terms approximation - Grid's computing 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 Grid_rectangular
include 'icomm_1.f90'
x=0.
y=0.
SLx = 2.
SLy = 1.
Xbeg = -1.
Ybeg = 0.
! SLx = 1. ! SLy = 1.
! Xbeg = 0. ! Ybeg = 0.
NXmax = 121 NYmax = 61
NXmaxP = NXmax + 1 NYmaxP = NYmax + 1
!-------------------------------------------------------------------------------------------
DO 2 I=1, NXmax DO 2 J=1, NYmax
X(I,J) = Xbeg + (SLx /(NXmax-1))* (i-1)
2 continue
DO 3 J=1, NYmax
DO 3 I=1, NXmax
Y(I,J) = Ybeg + (SLy /(NYmax-1))* (j-1)
3 continue
!-------------------------------------------------------------------------------------------
open (22,file='GRID.dat')
WRITE(22,*)'VARIABLES = "X", "Y" '
WRITE (22,*)' ZONE I=' ,NXmax, ', J=', NYmax, ', F=POINT'
DO 1 J=1, NYmax DO 1 I=1, NXmax
WRITE (22,*) X(I,J), Y(I,J)
1 continue
open (21,file='GRID.txt')
write(21,*) NYmax,NXmax
DO 4 J=1, NYmax DO 4 I=1, NXmax
WRITE (21,*) X(I,J), Y(I,J)
4 continue
close(21)
Return End
