Main.f90 - Main modul

From CFD-Wiki

Revision as of 03:42, 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