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

Contr gem1.f90

From CFD-Wiki

(Difference between revisions)
Jump to: navigation, search
 
(2 intermediate revisions not shown)
Line 1: Line 1:
-
''This subroutine make only output of all geometrical date. (it was be done in order to control all geometry)''
+
''This subroutine make only output of all geometrical date. (it was done in order to control all geometry)''
-
Subroutine Control_gem
+
<pre>
 +
!Sample program for solving Smith-Hutton Test using different schemes
 +
!of covective terms approximation - Geometry control 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 Control_gem
include 'icomm_1.f90'
include 'icomm_1.f90'
Dimension F_out(nx,ny)
Dimension F_out(nx,ny)
-
 
Character  Filename*10
Character  Filename*10
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmax
NImax = NXmax
-
 
NJmax = NYmax
NJmax = NYmax
Line 20: Line 37:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!-------------------------------------------------------------------
!-------------------------------------------------------------------
-
 
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmax
NImax = NXmax
-
 
NJmax = NYmax
NJmax = NYmax
Line 34: Line 47:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!-------------------------------------------------------------------
!-------------------------------------------------------------------
-
 
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmaxP
NImax = NXmaxP
-
 
NJmax = NYmaxP
NJmax = NYmaxP
Line 48: Line 57:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!-------------------------------------------------------------------
!-------------------------------------------------------------------
-
 
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmaxP
NImax = NXmaxP
-
 
NJmax = NYmaxP
NJmax = NYmaxP
Line 62: Line 67:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!-------------------------------------------------------------------
!-------------------------------------------------------------------
-
 
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmaxP
NImax = NXmaxP
-
 
NJmax = NYmaxP
NJmax = NYmaxP
Line 76: Line 77:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!-------------------------------------------------------------------
!-------------------------------------------------------------------
-
 
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmaxP
NImax = NXmaxP
-
 
NJmax = NYmaxP
NJmax = NYmaxP
Line 90: Line 87:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!-------------------------------------------------------------------
!-------------------------------------------------------------------
-
 
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmaxP
NImax = NXmaxP
-
 
NJmax = NYmaxP
NJmax = NYmaxP
Line 104: Line 97:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!-------------------------------------------------------------------
!-------------------------------------------------------------------
-
 
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmaxP
NImax = NXmaxP
-
 
NJmax = NYmaxP
NJmax = NYmaxP
Line 118: Line 107:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmaxP
NImax = NXmaxP
-
 
NJmax = NYmaxP
NJmax = NYmaxP
Line 130: Line 116:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!-------------------------------------------------------------------
!-------------------------------------------------------------------
-
 
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmaxP
NImax = NXmaxP
-
 
NJmax = NYmaxP
NJmax = NYmaxP
Line 144: Line 126:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!-------------------------------------------------------------------
!-------------------------------------------------------------------
-
 
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmaxP
NImax = NXmaxP
-
 
NJmax = NYmaxP
NJmax = NYmaxP
Line 158: Line 136:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!-------------------------------------------------------------------
!-------------------------------------------------------------------
-
 
!----------------------------------------------------------------
!----------------------------------------------------------------
-
 
NImax = NXmaxP
NImax = NXmaxP
-
 
NJmax = NYmaxP
NJmax = NYmaxP
Line 172: Line 146:
     Call  Out_array(F_out,NImax,NJmax,Filename)
     Call  Out_array(F_out,NImax,NJmax,Filename)
-
 
!-------------------------------------------------------------------
!-------------------------------------------------------------------
Return
Return
 +
End
-
End
+
</pre>

Latest revision as of 14:46, 21 September 2005

This subroutine make only output of all geometrical date. (it was done in order to control all geometry)


!Sample program for solving Smith-Hutton Test using different schemes 
!of covective terms approximation - Geometry control 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 Control_gem
include 'icomm_1.f90'

Dimension F_out(nx,ny)
Character  Filename*10

!----------------------------------------------------------------
	NImax = NXmax
	NJmax = NYmax

	F_out     = X

 	Filename  ='X_ooo.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!-------------------------------------------------------------------
!----------------------------------------------------------------
	NImax = NXmax
	NJmax = NYmax

	F_out     = Y

 	Filename  ='Y_ooo.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!-------------------------------------------------------------------
!----------------------------------------------------------------
	NImax = NXmaxP
	NJmax = NYmaxP

	F_out     = Xp

 	Filename  ='X_poo.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!-------------------------------------------------------------------
!----------------------------------------------------------------
	NImax = NXmaxP
	NJmax = NYmaxP

	F_out     = Yp

 	Filename  ='Y_poo.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!-------------------------------------------------------------------
!----------------------------------------------------------------
	NImax = NXmaxP
	NJmax = NYmaxP

	F_out     = Del_X_xi

 	Filename  ='DX_xi.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!-------------------------------------------------------------------
!----------------------------------------------------------------
	NImax = NXmaxP
	NJmax = NYmaxP

	F_out     = Del_Y_xi

 	Filename  ='DY_xi.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!-------------------------------------------------------------------
!----------------------------------------------------------------
	NImax = NXmaxP
	NJmax = NYmaxP

	F_out     = Del_X_et

 	Filename  ='DX_et.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!-------------------------------------------------------------------
!----------------------------------------------------------------
	NImax = NXmaxP
	NJmax = NYmaxP

	F_out     = Del_Y_et

 	Filename  ='DY_et.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!----------------------------------------------------------------
	NImax = NXmaxP
	NJmax = NYmaxP

	F_out     = X_et

 	Filename  ='7X_et.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!-------------------------------------------------------------------
!----------------------------------------------------------------
	NImax = NXmaxP
	NJmax = NYmaxP

	F_out     = Y_et

 	Filename  ='7Y_et.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!-------------------------------------------------------------------
!----------------------------------------------------------------
	NImax = NXmaxP
	NJmax = NYmaxP

	F_out     = X_xi

 	Filename  ='7X_xi.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!-------------------------------------------------------------------
!----------------------------------------------------------------
	NImax = NXmaxP
	NJmax = NYmaxP

	F_out     = Y_xi

 	Filename  ='7Y_xi.txt' 

    Call  Out_array(F_out,NImax,NJmax,Filename)
!-------------------------------------------------------------------

Return
End

My wiki