CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Wiki > Main.f90 - Main modul

Main.f90 - Main modul

From CFD-Wiki

(Difference between revisions)
Jump to: navigation, search
Line 1: Line 1:
-
Subroutine Input_all
+
Program Main
include 'icomm_1.f90'
include 'icomm_1.f90'
-
open(21,file='GRID.txt')  
+
Dimension F_out(nx,ny)
-
+
-
read (21,*) NYmax,NXmax
+
-
DO 4 J=1, NYmax
+
Character  Filename*10
-
DO 4 I=1, NXmax
+
Call Ini
-
read (21,*) X(I,J), Y(I,J) 
+
Call Grid_rectangular
-
  4    continue
+
Call input_all
-
close(21)
+
Call Geom
 +
Call Init_all
-
!--------------------------------------------------------------------
+
    Call Control_gem
-
+
-
open(23,file='input_control.txt')
+
-
+
-
write (23,*) NYmax,NXmax
+
-
DO 5 J=1, NYmax
+
Do 100 Niter=1,100
-
DO 5 I=1, NXmax
+
write(*,*) '------------------------',Niter,'------------------------',Niter
-
write (23,*) X(I,J), Y(I,J) 
+
Call Coef_1(5)
-
  5    continue
+
Call Check_koeff
-
close(22)
+
!  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'
-
Return
+
STOP
-
End
+
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

My wiki