Commit d095a38e authored by Paul Stryck's avatar Paul Stryck
Browse files

two profiles

parent 5360c6c7
<log>
<WARNING outeriter="0000000000" >
Preparing grid without edgeparents, this might cause problems with cellsearch. Otherwise use -we option to generate edges.
</WARNING>
found NEW ICEM file ../../grids/project1_1
number of cells and faces 30940 123002
<WARNING outeriter="0000000000" >
Number of families smaller than boundary conditions,
probably BC for points and other objects in cgns file.
Skipping exeeding entries.
Proceed with fingers crossed!
</WARNING>
load distribution 30940 30940
load distribution 0.100000E+01 0.100000E+01
</log>
<?xml version="1.0"?>
<controls>
<general>
<name>Bent Pipe Flow</name>
<description>FDS-Team</description>
<caseid>CFD2 Tutorial</caseid>
</general>
<grids>
<grid name="grid">
<fromFile>
<gridFilePath>../../grids</gridFilePath>
<gridFileName>project1_1</gridFileName>
<restartFileName/>
<restartFilePath>.</restartFilePath>
</fromFile>
</grid>
</grids>
<timeLoop>
<unsteady>false</unsteady>
</timeLoop>
<outerLoop>
<maxIteration>10000</maxIteration>
<convergenceTolerance>1e-8</convergenceTolerance>
</outerLoop>
<equations>
<equation name="mome">
<EQMomentum>
<solve_z>false</solve_z>
<convectiveFluxDiscretisation>
<UD1_CD_BLEND>
<blendingFactor>0.7</blendingFactor>
</UD1_CD_BLEND>
</convectiveFluxDiscretisation>
<relaxationFactor>0.7</relaxationFactor>
</EQMomentum>
</equation>
<equation name="pres">
<EQPressure>
<pressureReferencePosition>0.5 0.5 0.5</pressureReferencePosition>
<pressureReference>1.e5</pressureReference>
<maxIteration>1000</maxIteration>
</EQPressure>
</equation>
<equation name="turb">
<EQTurbulence>
<turbulenceModel>KOMEGA_STANDARD</turbulenceModel>
</EQTurbulence>
</equation>
</equations>
<initial>
<velocity userCode="true">0. 0. 0.</velocity>
</initial>
<materials>
<backgroundMaterial>Air</backgroundMaterial>
<material>
<materialEntry name="Air">
<fluid>
<viscosityMolecular>1.8e-5</viscosityMolecular>
<density>1.2</density>
</fluid>
</materialEntry>
</material>
</materials>
<boundaries>
<family name="INLET">
<BCVelocity>
<velocity userCode="false">1.2 0.0 0.0</velocity>
<!-- <velocity userCode="false">6.0 0.0 0.0</velocity> -->
</BCVelocity>
</family>
<family name="OUTLET">
<BCOutflow/>
</family>
<family name="YMIN">
<BCWall>
<useWallFunction>LOW_RN</useWallFunction>
<!--useWallFunction>HIGH_RN</useWallFunction-->
</BCWall>
</family>
<family name="YMAX">
<BCWall>
<useWallFunction>LOW_RN</useWallFunction>
<!--useWallFunction>HIGH_RN</useWallFunction-->
</BCWall>
</family>
<family name="ZMIN">
<BCSymmetryPlane/>
</family>
<family name="ZMAX">
<BCSymmetryPlane/>
</family>
</boundaries>
<output>
<event name="default">
<saveTimeStep>1</saveTimeStep>
</event>
<event name="init">
<saveTimeStep>1</saveTimeStep>
<filePrefix>init</filePrefix>
<eventType>BEGIN_COMPUTATION</eventType>
<forceGridWrite>true</forceGridWrite>
</event>
</output>
</controls>
<readerControls>
<!--File generated by FreSCo+ after saving files for an event
File generation can be turned off by setting <generateMetaFile>false<generateMetaFile>
in controls section of the controls.xml file-->
<!-- event name init -->
<gridFileName>.//init</gridFileName>
<dataFileName>.//init</dataFileName>
<numberOfProcesses> 1</numberOfProcesses>
<firstTimeStep> 1</firstTimeStep>
<lastTimeStep> 1</lastTimeStep>
<timeStepIntervall> 1</timeStepIntervall>
<unsteady>false</unsteady>
</readerControls>
!***********************************************************************
!
! Hamburg Ship Model Basin (HSVA)
! TU Hamburg-Harburg (TUHH)
! Maritime Research Institute of Netherlands (MARIN)
!
! FreSCo
!
! $Id: post.F90 21 2005-06-08 12:50:20Z dieke $
!
! (c) 2004-2006 HSVA,TUHH,MARIN
! Proprietary data. Unauthorized use, distribution,
! or duplication is prohibited. All rights reserved.
!
!***********************************************************************
!
!! module post
! this module can be used to set up special post-processing routines. Use
! -post_outit to modyfy values in each outer iteration
! -post_timestep to modyfy values in each timesteo
! -post_final to set values after the solution
! algorithm is finished or stopped
!
!***********************************************************************
MODULE post
USE common_user_code
USE controls
USE equations
USE fielddata
USE fielddata_tools
USE fielddata_pointers
USE forces
USE geometry
USE logging
USE math
USE movegrid
USE parameters
USE topology
USE tracing
USE grids
USE parallel
USE counters
USE files
IMPLICIT NONE
PRIVATE
INTEGER :: unit_pressure
PUBLIC post_init_set_initial_fields_6,&
post_init_setup_lists_3,&
post_init_register_fields_4,&
post_init_after_grid_resize_5,&
post_timestep,&
post_outit,&
post_adoutit,&
post_final,&
post_exit_default_1
CONTAINS
!======================================================================
SUBROUTINE post_init_set_initial_fields_6
!=====================================================================
! use this to set initital values. To resolve dependencies between
! several *_inital modules this is called several times. Thus
! do not put ALLOCATE statements here.
!----------------------------------------------------------------------
CALL tracing_trace_begin(trace_post_init_set_initial_fields_6)
CALL tracing_trace_end(trace_post_init_set_initial_fields_6)
END SUBROUTINE post_init_set_initial_fields_6
!======================================================================
SUBROUTINE post_init_setup_lists_3
!======================================================================
! use this routine to set up linked lists, to get sizes or to allocate fields
!----------------------------------------------------------------------
CALL tracing_trace_begin(trace_post_init_setup_lists_3)
CALL files_register_file( fileBaseName = 'pressure',&
fileExtension = 'dat',&
timeStepFiles = .FALSE.,&
subDomainFiles = .FALSE.,&
appendOnRestart = .FALSE.,&
fileunit = unit_pressure)
CALL files_open_file(fileBaseName = 'pressure')
CALL tracing_trace_end(trace_post_init_setup_lists_3)
END SUBROUTINE post_init_setup_lists_3
!======================================================================
SUBROUTINE post_init_register_fields_4
!=====================================================================
! use *_init_register_fields_4 module for fielddata_init_register_fields_4,
! equation_system_register_equation_system ... statements
!----------------------------------------------------------------------
CALL tracing_trace_begin(trace_post_init_register_fields_4)
CALL tracing_trace_end(trace_post_init_register_fields_4)
END SUBROUTINE post_init_register_fields_4
!======================================================================
SUBROUTINE post_init_after_grid_resize_5
!=====================================================================
! use *_init_after_grid_resize_5 module for all *_get_pointer statements
!----------------------------------------------------------------------
CALL tracing_trace_begin(trace_post_init_after_grid_resize_5)
CALL tracing_trace_end(trace_post_init_after_grid_resize_5)
END SUBROUTINE post_init_after_grid_resize_5
!======================================================================
SUBROUTINE post_timestep(timestep,simultime)
!======================================================================
! this routine is executed after each timestep
!----------------------------------------------------------------------
IMPLICIT NONE
INTEGER :: timestep
REAL(dp) :: simultime
CALL tracing_trace_begin(trace_post_timestep)
timestep=timestep
simultime = simultime
CALL tracing_trace_end(trace_post_timestep)
END SUBROUTINE post_timestep
!======================================================================
SUBROUTINE post_outit(out_iter)
!======================================================================
! this routine is executed after each outer iteration
!----------------------------------------------------------------------
IMPLICIT NONE
INTEGER :: out_iter
INTEGER :: icell1, icell2, iface
REAL(dp) :: pressure1
REAL(dp) :: area1
REAL(dp) :: pressure2
REAL(dp) :: area2
REAL(dp) :: x1
REAL(dp) :: eps1
REAL(dp) :: y2
REAL(DP) :: eps2
CALL tracing_trace_begin(trace_post_outit)
out_iter = out_iter
x1 = 0
eps1 = 0.01
y2 = 0
eps2 = 0.01
pressure1 = 0
area1 = 0
pressure2 = 0
area2 = 0
! Interpolate pressure to the faces
p_f = fielddata_interpolate_f(p_c)
! Loop all boundary faces
DO iface = 1, nIntFaces
! Select Faces at the start of the curve
IF (( ABS(cent_f_3(iface,1) - x1) .LE. eps1 ) .AND. &
(cent_f_3(iface, 2) - y2) .GT. eps1 ) THEN
pressure1 = pressure1 + p_f(iface) * surf_f(iface)
area1 = area1 + surf_f(iface)
! Select Faces at the end of the curve
ELSE IF (( ABS(cent_f_3(iface,1) - x1) .LE. eps1 ) .AND. &
(cent_f_3(iface, 2) - y2) .LT. eps1 ) THEN
pressure2 = pressure2 + p_f(iface) * surf_f(iface)
area2 = area2 + surf_f(iface)
END IF
END DO
! Write pressures to file
WRITE(unit_pressure,*) out_iter, pressure1 / area1, pressure2 / area2
CALL tracing_trace_end(trace_post_outit)
END SUBROUTINE post_outit
!======================================================================
SUBROUTINE post_adoutit(adout_iter)
!======================================================================
! this routine is executed after each adjoint outer iteration
!----------------------------------------------------------------------
IMPLICIT NONE
INTEGER :: adout_iter
CALL tracing_trace_begin(trace_post_adoutit)
adout_iter = adout_iter
CALL tracing_trace_end(trace_post_adoutit)
END SUBROUTINE post_adoutit
!======================================================================
SUBROUTINE post_final
!======================================================================
! this routine is executed only once after the computation has stopped
! due to max no of iteration, truncation error or stopfile. It is not
! executed when killfile is used
!----------------------------------------------------------------------
CALL tracing_trace_begin(trace_post_final)
CALL tracing_trace_end(trace_post_final)
END SUBROUTINE post_final
!======================================================================
SUBROUTINE post_exit_default_1
!=====================================================================
! use *_exit_default_1 module for DEALLOCATE statements etc.
!----------------------------------------------------------------------
CALL tracing_trace_begin(trace_post_exit_default_1)
close(unit=400)
CALL tracing_trace_end(trace_post_exit_default_1)
END SUBROUTINE post_exit_default_1
END MODULE post
GFORTRAN module version '10' created from ./post.F90
MD5:fee54ef7b71ce46c817226aa4ef167e5 -- If you edit this, you'll get what you deserve.
(() () () () () () () () () () () () () () () () () () () () () () ()
() () () ())
()
()
(('mpi_fortran_argv_null' 2 0 0 '') ('h5i_flags' 3 0 0 '') ('h5f_flags'
4 0 0 '') ('h5fd_flags' 5 0 0 '') ('h5fd_hid_flags' 6 0 0 '') (
'h5o_flags' 7 0 0 '') ('h5g_flags' 8 0 0 '') ('h5generic_flags' 9 0 0 '')
('h5r_flags' 10 0 0 '') ('mpi_fortran_bottom' 11 0 0 '') ('h5l_flags' 12
0 0 '') ('h5lib_flags' 13 0 0 '') ('mpi_fortran_errcodes_ignore' 14 0 0 '')
('mpi_fortran_in_place' 15 0 0 '') ('mpi_fortran_status_ignore' 16 0 0 '')
('mpi_fortran_statuses_ignore' 17 0 0 '') ('mpi_fortran_unweighted' 18 0
0 '') ('h5s_flags' 19 0 0 '') ('mpi_fortran_weights_empty' 20 0 0 '') (
'h5s_hsize_flags' 21 0 0 '') ('floating_types' 22 0 0 '') ('h5p_flags'
23 0 0 '') ('h5p_flags_int' 24 0 0 '') ('h5t_flags' 25 0 0 '') (
'h5d_flags' 26 0 0 '') ('h5z_flags' 27 0 0 '') ('integer_types' 28 0 0 '')
('mpi_fortran_argvs_null' 29 0 0 '') ('h5e_hid_flags' 30 0 0 '') (
'predefined_types' 31 0 0 '') ('h5d_size_flags' 32 0 0 '') ('h5e_flags'
33 0 0 ''))
()
()
(34 'post_adoutit' 'post' '' 1 ((PROCEDURE UNKNOWN-INTENT MODULE-PROC
DECL UNKNOWN 0 0 SUBROUTINE) (UNKNOWN 0 0 0 0 UNKNOWN ()) 35 0 (36) () 0
() () () 0 0)
37 'post_exit_default_1' 'post' '' 1 ((PROCEDURE UNKNOWN-INTENT
MODULE-PROC DECL UNKNOWN 0 0 SUBROUTINE) (UNKNOWN 0 0 0 0 UNKNOWN ()) 0
0 () () 0 () () () 0 0)
38 'post_final' 'post' '' 1 ((PROCEDURE UNKNOWN-INTENT MODULE-PROC DECL
UNKNOWN 0 0 SUBROUTINE) (UNKNOWN 0 0 0 0 UNKNOWN ()) 0 0 () () 0 () () ()
0 0)
39 'post_init_after_grid_resize_5' 'post' '' 1 ((PROCEDURE
UNKNOWN-INTENT MODULE-PROC DECL UNKNOWN 0 0 SUBROUTINE) (UNKNOWN 0 0 0 0
UNKNOWN ()) 0 0 () () 0 () () () 0 0)
40 'post_init_register_fields_4' 'post' '' 1 ((PROCEDURE UNKNOWN-INTENT
MODULE-PROC DECL UNKNOWN 0 0 SUBROUTINE) (UNKNOWN 0 0 0 0 UNKNOWN ()) 0
0 () () 0 () () () 0 0)
41 'post_init_set_initial_fields_6' 'post' '' 1 ((PROCEDURE
UNKNOWN-INTENT MODULE-PROC DECL UNKNOWN 0 0 SUBROUTINE) (UNKNOWN 0 0 0 0
UNKNOWN ()) 0 0 () () 0 () () () 0 0)
42 'post_init_setup_lists_3' 'post' '' 1 ((PROCEDURE UNKNOWN-INTENT
MODULE-PROC DECL UNKNOWN 0 0 SUBROUTINE) (UNKNOWN 0 0 0 0 UNKNOWN ()) 0
0 () () 0 () () () 0 0)
43 'post_outit' 'post' '' 1 ((PROCEDURE UNKNOWN-INTENT MODULE-PROC DECL
UNKNOWN 0 0 SUBROUTINE) (UNKNOWN 0 0 0 0 UNKNOWN ()) 44 0 (45) () 0 () ()
() 0 0)
46 'post_timestep' 'post' '' 1 ((PROCEDURE UNKNOWN-INTENT MODULE-PROC
DECL UNKNOWN 0 0 SUBROUTINE) (UNKNOWN 0 0 0 0 UNKNOWN ()) 47 0 (48 49) ()
0 () () () 0 0)
2 'mpi_argv_null' 'mpi' '' 1 ((VARIABLE UNKNOWN-INTENT UNKNOWN-PROC
UNKNOWN UNKNOWN 0 0 DIMENSION IN_COMMON) (CHARACTER 1 0 0 0 CHARACTER (
(CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '1'))) 0 0 () (1 0 EXPLICIT (
CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '1') (CONSTANT (INTEGER 4 0 0 0
INTEGER ()) 0 '1')) 0 () () () 0 0)
3 'h5i_flags' 'h5global' '' 1 ((VARIABLE UNKNOWN-INTENT UNKNOWN-PROC
UNKNOWN IMPLICIT-SAVE 0 0 DIMENSION IN_COMMON) (INTEGER 4 0 0 0 INTEGER
()) 0 0 () (1 0 EXPLICIT (CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '1') (
CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '7')) 0 () () () 0 0)
4 'h5f_flags' 'h5global' '' 1 ((VARIABLE UNKNOWN-INTENT UNKNOWN-PROC
UNKNOWN IMPLICIT-SAVE 0 0 DIMENSION IN_COMMON) (INTEGER 4 0 0 0 INTEGER
()) 0 0 () (1 0 EXPLICIT (CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '1') (
CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '19')) 0 () () () 0 0)
5 'h5fd_flags' 'h5global' '' 1 ((VARIABLE UNKNOWN-INTENT UNKNOWN-PROC
UNKNOWN IMPLICIT-SAVE 0 0 DIMENSION IN_COMMON) (INTEGER 4 0 0 0 INTEGER
()) 0 0 () (1 0 EXPLICIT (CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '1') (
CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '11')) 0 () () () 0 0)
6 'h5fd_hid_flags' 'h5global' '' 1 ((VARIABLE UNKNOWN-INTENT
UNKNOWN-PROC UNKNOWN IMPLICIT-SAVE 0 0 DIMENSION IN_COMMON) (INTEGER 4 0
0 0 INTEGER ()) 0 0 () (1 0 EXPLICIT (CONSTANT (INTEGER 4 0 0 0 INTEGER
()) 0 '1') (CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '7')) 0 () () () 0 0)
7 'h5o_flags' 'h5global' '' 1 ((VARIABLE UNKNOWN-INTENT UNKNOWN-PROC
UNKNOWN IMPLICIT-SAVE 0 0 DIMENSION IN_COMMON) (INTEGER 4 0 0 0 INTEGER
()) 0 0 () (1 0 EXPLICIT (CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '1') (
CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '27')) 0 () () () 0 0)
8 'h5g_flags' 'h5global' '' 1 ((VARIABLE UNKNOWN-INTENT UNKNOWN-PROC
UNKNOWN IMPLICIT-SAVE 0 0 DIMENSION IN_COMMON) (INTEGER 4 0 0 0 INTEGER
()) 0 0 () (1 0 EXPLICIT (CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '1') (
CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '12')) 0 () () () 0 0)
9 'h5generic_flags' 'h5global' '' 1 ((VARIABLE UNKNOWN-INTENT
UNKNOWN-PROC UNKNOWN IMPLICIT-SAVE 0 0 DIMENSION IN_COMMON) (INTEGER 4 0
0 0 INTEGER ()) 0 0 () (1 0 EXPLICIT (CONSTANT (INTEGER 4 0 0 0 INTEGER
()) 0 '1') (CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '9')) 0 () () () 0 0)
10 'h5r_flags' 'h5global' '' 1 ((VARIABLE UNKNOWN-INTENT UNKNOWN-PROC
UNKNOWN IMPLICIT-SAVE 0 0 DIMENSION IN_COMMON) (INTEGER 4 0 0 0 INTEGER
()) 0 0 () (1 0 EXPLICIT (CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '1') (
CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '2')) 0 () () () 0 0)
11 'mpi_bottom' 'mpi' '' 1 ((VARIABLE UNKNOWN-INTENT UNKNOWN-PROC
UNKNOWN UNKNOWN 0 0 IN_COMMON) (INTEGER 4 0 0 0 INTEGER ()) 0 0 () () 0
() () () 0 0)
12 'h5l_flags' 'h5global' '' 1 ((VARIABLE UNKNOWN-INTENT UNKNOWN-PROC
UNKNOWN IMPLICIT-SAVE 0 0 DIMENSION IN_COMMON) (INTEGER 4 0 0 0 INTEGER
()) 0 0 () (1 0 EXPLICIT (CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '1') (
CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '6')) 0 () () () 0 0)
13 'h5lib_flags' 'h5global' '' 1 ((VARIABLE UNKNOWN-INTENT UNKNOWN-PROC
UNKNOWN IMPLICIT-SAVE 0 0 DIMENSION IN_COMMON) (INTEGER 4 0 0 0 INTEGER
()) 0 0 () (1 0 EXPLICIT (CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '1') (
CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '2')) 0 () () () 0 0)
14 'mpi_errcodes_ignore' 'mpi' '' 1 ((VARIABLE UNKNOWN-INTENT
UNKNOWN-PROC UNKNOWN UNKNOWN 0 0 DIMENSION IN_COMMON) (INTEGER 4 0 0 0
INTEGER ()) 0 0 () (1 0 EXPLICIT (CONSTANT (INTEGER 4 0 0 0 INTEGER ())
0 '1') (CONSTANT (INTEGER 4 0 0 0 INTEGER ()) 0 '1')) 0 () () () 0 0)
15 'mpi_in_place' 'mpi' '' 1 ((VARIABLE UNKNOWN-INTENT UNKNOWN-PROC
UNKNOWN UNKNOWN 0 0 IN_COMMON) (INTEGER 4 0 0 0 INTEGER ()) 0 0 () () 0
() () () 0 0)