Rocstar  1.0
Rocstar multiphysics simulation application
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
RFLO_PrintUserInput.F90
Go to the documentation of this file.
1 ! *********************************************************************
2 ! * Rocstar Simulation Suite *
3 ! * Copyright@2015, Illinois Rocstar LLC. All rights reserved. *
4 ! * *
5 ! * Illinois Rocstar LLC *
6 ! * Champaign, IL *
7 ! * www.illinoisrocstar.com *
8 ! * sales@illinoisrocstar.com *
9 ! * *
10 ! * License: See LICENSE file in top level of distribution package or *
11 ! * http://opensource.org/licenses/NCSA *
12 ! *********************************************************************
13 ! *********************************************************************
14 ! * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, *
15 ! * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES *
16 ! * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND *
17 ! * NONINFRINGEMENT. IN NO EVENT SHALL THE CONTRIBUTORS OR *
18 ! * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER *
19 ! * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, *
20 ! * Arising FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE *
21 ! * USE OR OTHER DEALINGS WITH THE SOFTWARE. *
22 ! *********************************************************************
23 !******************************************************************************
24 !
25 ! Purpose: write out user input for checking purposes.
26 !
27 ! Description: none.
28 !
29 ! Input: regions = user input.
30 !
31 ! Output: to standard output.
32 !
33 ! Notes: none.
34 !
35 !******************************************************************************
36 !
37 ! $Id: RFLO_PrintUserInput.F90,v 1.16 2009/08/12 04:15:58 mtcampbe Exp $
38 !
39 ! Copyright: (c) 2001 by the University of Illinois
40 !
41 !******************************************************************************
42 
43 SUBROUTINE rflo_printuserinput( regions )
44 
45  USE moddatatypes
46  USE moddatastruct, ONLY : t_region
47  USE modglobal, ONLY : t_global
48 #ifdef PLAG
50 #endif
51 #ifdef PEUL
53 #endif
54 #ifdef PERI
56 #endif
57 #ifdef RADI
59 #endif
60 #ifdef SPEC
62 #endif
63 #ifdef TURB
65 #endif
66 #ifdef INRT
68 #endif
69  USE moderror
70  USE modparameters
71  IMPLICIT NONE
72 
73 ! ... parameters
74  TYPE(t_region), POINTER :: regions(:)
75 
76 ! ... loop variables
77  INTEGER :: ireg, iproc
78 
79 ! ... local variables
80  TYPE(t_global), POINTER :: global
81 
82 !******************************************************************************
83 
84  global => regions(1)%global
85 
86  CALL registerfunction( global,'RFLO_PrintUserInput',&
87  'RFLO_PrintUserInput.F90' )
88 
89 ! start -----------------------------------------------------------------------
90 
91  WRITE(stdout,1000) solver_name
92 
93 ! region mapping
94 
95  WRITE(stdout,1005) solver_name//' Region mapping:'
96  IF(global%nRegionsProc > 0) THEN
97  WRITE(stdout,1010) solver_name//' nRegions/Processor:',global%nRegionsProc
98  ENDIF
99  IF (global%verbLevel.gt.verbose_med) then
100  DO iproc=0,global%nProcAlloc-1
101  WRITE(stdout,'(/,A,I5,A)',advance='no') solver_name//' proc. ',iproc,': '
102  DO ireg=1,global%nRegions
103  IF (regions(ireg)%procid == iproc) THEN
104  WRITE(stdout,'(1X,I5)',advance='no') ireg
105  ENDIF
106  ENDDO
107  ENDDO
108  WRITE(stdout,'(/)')
109  ENDIF
110 ! reference values
111 
112  WRITE(stdout,1005) solver_name//' Reference values:'
113  WRITE(stdout,1020) solver_name//' absvel ',global%refVelocity
114  WRITE(stdout,1020) solver_name//' press ',global%refPressure
115  WRITE(stdout,1020) solver_name//' dens ',global%refDensity
116  WRITE(stdout,1020) solver_name//' cp ',global%refCp
117  WRITE(stdout,1020) solver_name//' gamma ',global%refGamma
118  WRITE(stdout,1020) solver_name//' length ',global%refLength
119  WRITE(stdout,1020) solver_name//' renum ',global%refREnum
120  WRITE(stdout,1020) solver_name//' visc ',global%refVisc
121  WRITE(stdout,1020) solver_name//' prlam ',global%prLam
122  WRITE(stdout,1020) solver_name//' prturb ',global%prTurb
123  WRITE(stdout,1020) solver_name//' scnlam ',global%scnLam
124  WRITE(stdout,1020) solver_name//' scnturb',global%scnTurb
125 
126 ! acceleration
127 
128  WRITE(stdout,1005) solver_name//' Acceleration terms:'
129  IF (global%accelOn) THEN
130  WRITE(stdout,1030) solver_name//' terms on'
131  WRITE(stdout,1020) solver_name//' accel-x ',global%accelX
132  WRITE(stdout,1020) solver_name//' accel-y ',global%accelY
133  WRITE(stdout,1020) solver_name//' accel-z ',global%accelZ
134  ELSE
135  WRITE(stdout,1030) solver_name//' terms off'
136  ENDIF
137 
138 ! materials
139 
140 #ifdef INRT
141  CALL inrt_printmaterialinput(global)
142 #endif
143 
144 ! time stepping
145 
146  WRITE(stdout,1005) solver_name//' Time stepping:'
147  IF (global%solverType == solv_explicit) THEN
148  WRITE(stdout,1030) solver_name//' explicit scheme'
149  ELSE
150  WRITE(stdout,1030) solver_name//' implicit scheme'
151  ENDIF
152  IF (global%flowType == flow_unsteady) THEN
153  WRITE(stdout,1020) solver_name//' timestep',global%dtImposed
154  WRITE(stdout,1020) solver_name//' time ',global%timeStamp
155  WRITE(stdout,1020) solver_name//' maxtime ',global%maxTime
156  WRITE(stdout,1020) solver_name//' writime ',global%writeTime
157  WRITE(stdout,1020) solver_name//' prntime ',global%printTime
158  IF (global%solverType == solv_implicit) THEN
159  WRITE(stdout,1015) solver_name//' order ',global%tstepOrder
160  WRITE(stdout,1015) solver_name//' subiter ',global%maxSubIter
161  WRITE(stdout,1020) solver_name//' tolsub ',global%tolSubIter
162  IF (global%predictSol) THEN
163  WRITE(stdout,1030) solver_name//' predict = yes'
164  ELSE
165  WRITE(stdout,1030) solver_name//' predict = no'
166  ENDIF
167  IF (global%dtFixed) THEN
168  WRITE(stdout,1030) solver_name//' fixed dt = yes'
169  ELSE
170  WRITE(stdout,1030) solver_name//' fixed dt = no'
171  ENDIF
172  ENDIF
173  WRITE(stdout,1015) solver_name//' rkScheme',global%rkScheme
174  ELSE
175  WRITE(stdout,1015) solver_name//' iter ',global%currentIter
176  WRITE(stdout,1015) solver_name//' maxiter',global%maxIter
177  WRITE(stdout,1020) solver_name//' restol ',global%resTol
178  WRITE(stdout,1015) solver_name//' wriiter',global%writeIter
179  WRITE(stdout,1015) solver_name//' prniter',global%printIter
180  ENDIF
181 
182 ! grid motion
183 
184  WRITE(stdout,1005) solver_name//' Grid motion:'
185  IF (global%moveGridScheme == movegrid_blocks) THEN
186  WRITE(stdout,1030) solver_name//' block-TFI scheme'
187  ELSEIF (global%moveGridScheme == movegrid_global) THEN
188  WRITE(stdout,1030) solver_name//' block-wght-Laplacian scheme'
189  WRITE(stdout,1015) solver_name//' niter ',global%moveGridNiter
190  WRITE(stdout,1020) solver_name//' power ',global%moveGridPower
191  ELSEIF (global%moveGridScheme == movegrid_frame) THEN
192  WRITE(stdout,1030) solver_name//' global-wght-Laplacian scheme'
193  WRITE(stdout,1015) solver_name//' niter ',global%moveGridNiter
194  WRITE(stdout,1020) solver_name//' amplifx ',global%moveGridAmplifX
195  WRITE(stdout,1020) solver_name//' amplify ',global%moveGridAmplifY
196  WRITE(stdout,1020) solver_name//' amplifz ',global%moveGridAmplifZ
197  WRITE(stdout,1020) solver_name//' power ',global%moveGridPower
198  WRITE(stdout,1015) solver_name//' neighbors ',global%moveGridNbour
199  WRITE(stdout,1015) solver_name//' orthodir ',global%moveGridOrthDir
200  WRITE(stdout,1020) solver_name//' orthowghtx',global%moveGridOrthWghtX
201  WRITE(stdout,1020) solver_name//' orthowghtY',global%moveGridOrthWghtY
202  WRITE(stdout,1020) solver_name//' orthowghtZ',global%moveGridOrthWghtZ
203  WRITE(stdout,1015) solver_name//' nsurfmatch',global%moveGridNsmatch
204  ELSEIF (global%moveGridScheme == movegrid_foms) THEN
205  WRITE(stdout,1030) solver_name//' global-orthogonal-Laplacian scheme'
206  WRITE(stdout,1015) solver_name//' niter ',global%moveGridNiter
207  WRITE(stdout,1020) solver_name//' weight ',global%moveGridWeight
208  WRITE(stdout,1020) solver_name//' amplifx ',global%moveGridAmplifX
209  WRITE(stdout,1020) solver_name//' amplify ',global%moveGridAmplifY
210  WRITE(stdout,1020) solver_name//' amplifz ',global%moveGridAmplifZ
211  WRITE(stdout,1020) solver_name//' power ',global%moveGridPower
212  WRITE(stdout,1015) solver_name//' neighbors ',global%moveGridNbour
213  WRITE(stdout,1015) solver_name//' orthodir ',global%moveGridOrthDir
214  WRITE(stdout,1020) solver_name//' orthowghtx',global%moveGridOrthWghtX
215  WRITE(stdout,1020) solver_name//' orthowghtY',global%moveGridOrthWghtY
216  WRITE(stdout,1020) solver_name//' orthowghtZ',global%moveGridOrthWghtZ
217  WRITE(stdout,1020) solver_name//' weight ',global%moveGridWeight
218  WRITE(stdout,1020) solver_name//' orthocell ',global%moveGridOrthCell
219  WRITE(stdout,1015) solver_name//' nsurfmatch',global%moveGridNsmatch
220  ELSEIF (global%moveGridScheme == movegrid_elglobal) THEN
221  WRITE(stdout,1030) solver_name//' block-elliptic-PDE scheme'
222  WRITE(stdout,1015) solver_name//' niter ',global%moveGridNiter
223  WRITE(stdout,1015) solver_name//' viter ',global%moveGridViter
224  WRITE(stdout,1015) solver_name//' siter ',global%moveGridSiter
225  WRITE(stdout,1020) solver_name//' power ',global%moveGridPower
226  ELSEIF (global%moveGridScheme == movegrid_elframe) THEN
227  WRITE(stdout,1030) solver_name//' global-elliptic-PDE scheme'
228  WRITE(stdout,1015) solver_name//' niter ',global%moveGridNiter
229  WRITE(stdout,1015) solver_name//' viter ',global%moveGridViter
230  WRITE(stdout,1015) solver_name//' siter ',global%moveGridSiter
231  WRITE(stdout,1020) solver_name//' amplifx ',global%moveGridAmplifX
232  WRITE(stdout,1020) solver_name//' amplify ',global%moveGridAmplifY
233  WRITE(stdout,1020) solver_name//' amplifz ',global%moveGridAmplifZ
234  WRITE(stdout,1020) solver_name//' power ',global%moveGridPower
235  WRITE(stdout,1015) solver_name//' neighbors ',global%moveGridNbour
236  WRITE(stdout,1015) solver_name//' orthodir ',global%moveGridOrthDir
237  WRITE(stdout,1020) solver_name//' orthowghtx',global%moveGridOrthWghtX
238  WRITE(stdout,1020) solver_name//' orthowghtY',global%moveGridOrthWghtY
239  WRITE(stdout,1020) solver_name//' orthowghtZ',global%moveGridOrthWghtZ
240  WRITE(stdout,1015) solver_name//' nsurfmatch',global%moveGridNsmatch
241  ELSEIF (global%moveGridScheme == movegrid_vms) THEN
242  WRITE(stdout,1030) solver_name//' block-VMS scheme'
243  WRITE(stdout,1015) solver_name//' niter ',global%moveGridNiter
244  WRITE(stdout,1020) solver_name//' weight',global%moveGridWeight
245  ENDIF
246 
247 #ifdef STATS
248 ! statistics
249 
250  IF (global%doStat == active) THEN
251  IF (global%reStat == active) THEN
252  WRITE(stdout,1005) solver_name//' Continued statistics:'
253  ELSE
254  WRITE(stdout,1005) solver_name//' New statistics:'
255  ENDIF
256  IF (global%mixtNStat > 0) THEN
257  WRITE(stdout,*) solver_name//' mixture NStat =',global%mixtNStat
258  WRITE(stdout,*) solver_name//' mixture StatID =', &
259  global%mixtStatID(1,:)*10+global%mixtStatID(2,:)
260  ENDIF
261  IF (global%turbNStat > 0) THEN
262  WRITE(stdout,*) solver_name//' TURB NStat =',global%turbNStat
263  WRITE(stdout,*) solver_name//' TURB StatID =', &
264  global%turbStatID(1,:)*10+global%turbStatID(2,:)
265  ENDIF
266  ENDIF
267 #endif
268 
269 #ifdef PERI
270  IF (regions(1)%periInput%flowKind /= off) THEN
271  CALL peri_printuserinput( regions(1) )
272  ENDIF
273 #endif
274 
275 ! repeat for all regions ------------------------------------------------------
276 
277 
278  IF (global%verbLevel.gt.verbose_med) then
279  DO ireg=1,global%nRegions
280  WRITE(stdout,1025) solver_name,ireg
281 
282  ! - dimensions
283 
284  WRITE(stdout,1005) solver_name//' Dimensions:'
285  WRITE(stdout,1015) solver_name//' icells ',regions(ireg)%levels(1)%grid%ipc
286  WRITE(stdout,1015) solver_name//' jcells ',regions(ireg)%levels(1)%grid%jpc
287  WRITE(stdout,1015) solver_name//' kcells ',regions(ireg)%levels(1)%grid%kpc
288  WRITE(stdout,1010) solver_name//' dumcells',regions(ireg)%nDumCells
289 
290  ! - flow model
291 
292  WRITE(stdout,1005) solver_name//' Flow type, turbulence:'
293  IF (regions(ireg)%mixtInput%flowModel == flow_euler) THEN
294  WRITE(stdout,1030) solver_name//' flow model = Euler'
295  ELSE
296  WRITE(stdout,1030) solver_name//' flow model = Navier-Stokes'
297  ENDIF
298  IF (regions(ireg)%mixtInput%moveGrid) THEN
299  WRITE(stdout,1030) solver_name//' moving grid = yes'
300  ELSE
301  WRITE(stdout,1030) solver_name//' moving grid = no'
302  ENDIF
303 
304  IF (regions(ireg)%mixtInput%flowModel == flow_navst) THEN
305  IF (regions(ireg)%mixtInput%turbModel == turb_model_none) THEN
306  WRITE(stdout,1030) solver_name//' turb. model = none'
307  ELSE
308 #ifdef TURB
309  CALL turb_printuserinput( regions(ireg) )
310 #endif
311  ENDIF
312  ENDIF
313 
314  ! - viscosity model
315 
316  IF (regions(ireg)%mixtInput%computeTv) THEN
317  WRITE(stdout,1005) solver_name//' Viscosity model:'
318 
319  IF (regions(ireg)%mixtInput%viscModel == visc_suthr) THEN
320  WRITE(stdout,1030) solver_name//' viscosity model = Sutherland'
321  WRITE(stdout,1020) solver_name//' reference viscosity ',&
322  regions(ireg)%mixtInput%refVisc
323  WRITE(stdout,1020) solver_name//' reference temperature ',&
324  regions(ireg)%mixtInput%refTemp
325  WRITE(stdout,1020) solver_name//' Sutherland coefficient',&
326  regions(ireg)%mixtInput%suthCoef
327 
328  ELSEIF (regions(ireg)%mixtInput%viscModel == visc_fixed) THEN
329  WRITE(stdout,1030) solver_name//' viscosity model = Fixed'
330  WRITE(stdout,1020) solver_name//' viscosity value ',&
331  regions(ireg)%mixtInput%refVisc
332 
333  ELSEIF (regions(ireg)%mixtInput%viscModel == visc_antib) THEN
334  WRITE(stdout,1030) solver_name//' viscosity model = Antibes'
335  WRITE(stdout,1020) solver_name//' reference viscosity',&
336  regions(ireg)%mixtInput%refVisc
337 
338  ENDIF ! viscModel
339  ENDIF ! computeTv
340 
341  ! - multi-physics modules
342 
343  WRITE(stdout,1005) solver_name//' Multi-physics modules:'
344 
345  IF (regions(ireg)%mixtInput%gasModel == gas_model_tcperf) THEN
346  WRITE(stdout,1030) solver_name//' species model = none'
347  ELSE
348 #ifdef SPEC
349  CALL spec_printuserinput( regions(ireg) )
350 #endif
351  ENDIF
352 
353  IF (global%peulUsed) THEN
354 #ifdef PEUL
355  WRITE(stdout,1030) solver_name//' con.part. used = yes'
356  CALL peul_printuserinput( regions(ireg) )
357 #endif
358  ELSE
359  WRITE(stdout,1030) solver_name//' con.part. used = no'
360  ENDIF
361 
362  IF (global%plagUsed) THEN
363 #ifdef PLAG
364  WRITE(stdout,1030) solver_name//' dis.part. used = yes'
365  CALL plag_printuserinput( regions(ireg) )
366 #endif
367  ELSE
368  WRITE(stdout,1030) solver_name//' dis.part. used = no'
369  ENDIF
370 
371  IF (regions(ireg)%mixtInput%radiUsed) THEN
372 #ifdef RADI
373  CALL radi_printuserinput( regions(ireg) )
374 #endif
375  ELSE
376  WRITE(stdout,1030) solver_name//' radiation used = no'
377  ENDIF
378 
379 #ifdef INRT
380  IF (global%inrtUsed) THEN
381  CALL inrt_printuserinput( regions(ireg) )
382  ENDIF
383 #endif
384 
385  ! - numerics
386 
387  WRITE(stdout,1005) solver_name//' Numerics:'
388  WRITE(stdout,1010) solver_name//' levels ',regions(ireg)%nGridLevels
389  IF (regions(ireg)%mixtInput%timeScheme == tst_hyb5rk) THEN
390  WRITE(stdout,1030) solver_name//' timedis = explicit multistage'
391  ELSE IF (regions(ireg)%mixtInput%timeScheme == tst_std4rk) THEN
392  WRITE(stdout,1030) solver_name//' timedis = explicit classical Runge-Kutta'
393  ENDIF
394  WRITE(stdout,1020) solver_name//' CFL ',regions(ireg)%mixtInput%cfl
395  WRITE(stdout,1020) solver_name//' smoocf ',regions(ireg)%mixtInput%smoocf
396  WRITE(stdout,1010) solver_name//' discr ',regions(ireg)%mixtInput%spaceDiscr
397  IF (regions(ireg)%mixtInput%spaceDiscr == discr_cen_scal) THEN
398  WRITE(stdout,1020) solver_name//' k2 ',regions(ireg)%mixtInput%vis2
399  WRITE(stdout,1020) solver_name//' 1/k4 ',1./regions(ireg)%mixtInput%vis4
400  IF (regions(ireg)%mixtInput%pSwitchType == pswitch_std) THEN
401  WRITE(stdout,1030) solver_name//' pswitch = standard'
402  ELSE
403  WRITE(stdout,1030) solver_name//' pswitch = mixed TVD and standard'
404  WRITE(stdout,1020) solver_name//' pswOmega',regions(ireg)%mixtInput%pSwitchOmega
405  ENDIF
406  ELSE
407  WRITE(stdout,1010) solver_name//' order ',regions(ireg)%mixtInput%spaceOrder
408  WRITE(stdout,1020) solver_name//' limfac ',regions(ireg)%mixtInput%limfac
409  WRITE(stdout,1020) solver_name//' entropy ',regions(ireg)%mixtInput%epsentr
410  ENDIF
411  IF (regions(ireg)%mixtInput%faceEdgeAvg == fe_avg_uniform) THEN
412  WRITE(stdout,1030) solver_name//' f/e avg = uniform'
413  ELSE
414  WRITE(stdout,1030) solver_name//' f/e avg = grid-dependent linear'
415  ENDIF
416 
417  ENDDO ! iReg
418  ENDIF
419 ! finish ----------------------------------------------------------------------
420 
421  WRITE(stdout,1035) solver_name
422 
423  CALL deregisterfunction( global )
424 
425 1000 FORMAT(/,a,1x,80('-'))
426 1005 FORMAT(/,a)
427 1010 FORMAT(a,' = ',i2)
428 1015 FORMAT(a,' = ',i8)
429 1020 FORMAT(a,' = ',e12.5)
430 1025 FORMAT(/,a,' Region ',i6,':')
431 1030 FORMAT(a)
432 1035 FORMAT(/,a,1x,80('-'),/)
433 
434 END SUBROUTINE rflo_printuserinput
435 
436 !******************************************************************************
437 !
438 ! RCS Revision history:
439 !
440 ! $Log: RFLO_PrintUserInput.F90,v $
441 ! Revision 1.16 2009/08/12 04:15:58 mtcampbe
442 ! Major update, bugfix from Abe development, more propagation compatibility,
443 ! some Rocstar IO changes, Ju's temporary clipping fix for turbulence. A bug
444 ! fix for initialization IO.
445 !
446 ! Revision 1.15 2008/12/06 08:44:27 mtcampbe
447 ! Updated license.
448 !
449 ! Revision 1.14 2008/11/19 22:17:38 mtcampbe
450 ! Added Illinois Open Source License/Copyright
451 !
452 ! Revision 1.13 2006/05/09 23:33:44 wasistho
453 ! added DTFIXED in timestep section
454 !
455 ! Revision 1.12 2006/04/27 01:55:14 wasistho
456 ! added screenprint neigbors
457 !
458 ! Revision 1.11 2006/03/18 13:26:45 wasistho
459 ! added orthDir and orthWghtX,Y,Z
460 !
461 ! Revision 1.10 2006/03/08 06:34:38 wasistho
462 ! added movegrid_elglobal and elframe
463 !
464 ! Revision 1.9 2005/11/18 07:21:01 wasistho
465 ! rearranged gridmotion screen print
466 !
467 ! Revision 1.8 2005/10/31 21:09:36 haselbac
468 ! Changed specModel and SPEC_MODEL_NONE
469 !
470 ! Revision 1.7 2005/10/28 22:47:24 wasistho
471 ! print FOMS orthocell
472 !
473 ! Revision 1.6 2005/10/28 05:42:44 wasistho
474 ! printout FOMS gridmotion
475 !
476 ! Revision 1.5 2005/08/28 23:49:05 wasistho
477 ! added orthoWght for block orthogonality of RFLO global-gridmotion
478 !
479 ! Revision 1.4 2005/08/18 19:49:19 wasistho
480 ! added print NSURFMATCH
481 !
482 ! Revision 1.3 2005/06/04 01:02:20 wasistho
483 ! distinguished to AMPLIFX,Y,Z
484 !
485 ! Revision 1.2 2005/06/02 22:59:54 wasistho
486 ! added moveGridAmplif and moveGridPower
487 !
488 ! Revision 1.1 2004/11/29 20:51:39 wasistho
489 ! lower to upper case
490 !
491 ! Revision 1.34 2004/11/17 16:30:44 haselbac
492 ! Added printing of rkScheme
493 !
494 ! Revision 1.33 2004/09/02 02:59:42 wasistho
495 ! screen output face-edge averaging option
496 !
497 ! Revision 1.32 2004/03/05 22:09:02 jferry
498 ! created global variables for peul, plag, and inrt use
499 !
500 ! Revision 1.31 2004/03/02 21:49:22 jferry
501 ! Added inrtUsed flag to mixture data structure
502 !
503 ! Revision 1.30 2003/11/20 16:40:40 mdbrandy
504 ! Backing out RocfluidMP changes from 11-17-03
505 !
506 ! Revision 1.26 2003/10/01 23:52:10 jblazek
507 ! Corrected bug in moving noslip wall BC and grid speeds.
508 !
509 ! Revision 1.25 2003/09/26 21:44:28 fnajjar
510 ! Modified ModInterfaces calls to proper physical modules
511 !
512 ! Revision 1.24 2003/08/28 20:37:43 wasistho
513 ! excluced ModInterfacesTurbulence,Radiation,Periodic from ModInterfaces
514 !
515 ! Revision 1.23 2003/08/28 20:05:39 jblazek
516 ! Added acceleration terms.
517 !
518 ! Revision 1.22 2003/08/11 21:51:18 jblazek
519 ! Added basic global grid smoothing scheme.
520 !
521 ! Revision 1.21 2003/07/03 21:48:45 jblazek
522 ! Implemented dual-time stepping.
523 !
524 ! Revision 1.20 2003/05/15 02:57:04 jblazek
525 ! Inlined index function.
526 !
527 ! Revision 1.19 2003/04/10 23:31:20 fnajjar
528 ! Added printouts for viscosity models
529 !
530 ! Revision 1.18 2003/03/29 03:30:18 wasistho
531 ! install ROCPERI
532 !
533 ! Revision 1.17 2003/03/24 23:27:01 jferry
534 ! converted PrintMaterialInput to INRT_PrintMaterialInput
535 !
536 ! Revision 1.16 2003/03/11 16:04:19 jferry
537 ! Enclosed USE statements for multi-physics routines within ifdefs
538 !
539 ! Revision 1.15 2003/03/04 22:12:35 jferry
540 ! Initial import of Rocinteract
541 !
542 ! Revision 1.14 2002/11/04 22:30:37 wasistho
543 ! Put stats print userinput within ifdef STATS
544 !
545 ! Revision 1.13 2002/11/02 01:57:21 wasistho
546 ! Added TURB statistics
547 !
548 ! Revision 1.12 2002/09/17 22:51:23 jferry
549 ! Removed Fast Eulerian particle type
550 !
551 ! Revision 1.11 2002/09/05 17:40:21 jblazek
552 ! Variable global moved into regions().
553 !
554 ! Revision 1.10 2002/07/25 00:40:25 jblazek
555 ! Option for TVD type pressure switch.
556 !
557 ! Revision 1.9 2002/07/16 21:34:37 jblazek
558 ! Prefixed screen output with SOLVER_NAME.
559 !
560 ! Revision 1.8 2002/06/22 01:13:38 jblazek
561 ! Modified interfaces to BC routines.
562 !
563 ! Revision 1.7 2002/02/27 18:38:20 jblazek
564 ! Changed extrapol. to dummy cells at injection boundaries and slip walls.
565 !
566 ! Revision 1.6 2002/02/25 22:36:53 jblazek
567 ! Simplified solver initialization routine.
568 !
569 ! Revision 1.5 2002/02/21 23:25:06 jblazek
570 ! Blocks renamed as regions.
571 !
572 ! Revision 1.4 2002/01/11 17:20:19 jblazek
573 ! Added time stamp or iteration number to file names.
574 !
575 ! Revision 1.3 2002/01/08 22:09:17 jblazek
576 ! Added calculation of face vectors and volumes.
577 !
578 ! Revision 1.2 2001/12/19 23:09:22 jblazek
579 ! Added routines to read grid and solution.
580 !
581 ! Revision 1.1 2001/12/11 21:59:29 jblazek
582 ! memory allocation added.
583 !
584 !******************************************************************************
585 
586 
587 
588 
589 
590 
591 
subroutine plag_printuserinput(region)
subroutine registerfunction(global, funName, fileName)
Definition: ModError.F90:449
subroutine rflo_printuserinput(regions)
subroutine turb_printuserinput(region)
subroutine inrt_printmaterialinput(global)
void int int REAL * x
Definition: read.cpp:74
subroutine spec_printuserinput(region)
subroutine peul_printuserinput(region)
subroutine inrt_printuserinput(region)
subroutine deregisterfunction(global)
Definition: ModError.F90:469
subroutine peri_printuserinput(region)
RT a() const
Definition: Line_2.h:140