Main.f90 - Main modul

From CFD-Wiki

Revision as of 02:23, 19 September 2005 (view source) Michail (Talk | contribs) ← Older edit		Revision as of 02:25, 19 September 2005 (view source) Michail (Talk | contribs) Newer edit →
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