Main.f90 - Main modul
From CFD-Wiki
Line 1: | Line 1: | ||
- | + | Program Main | |
include 'icomm_1.f90' | include 'icomm_1.f90' | ||
- | + | Dimension F_out(nx,ny) | |
- | + | ||
- | + | ||
- | + | Character Filename*10 | |
- | + | Call Ini | |
- | + | Call Grid_rectangular | |
- | + | Call input_all | |
- | + | 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 |
Revision as of 02:25, 19 September 2005
Program Main
include 'icomm_1.f90'
Dimension F_out(nx,ny)
Character Filename*10
Call Ini
Call Grid_rectangular
Call input_all
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