Fluid_preHdfOutput.F90

Go to the documentation of this file.

! *********************************************************************
! * Rocstar Simulation Suite                                          *
! * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *
! *                                                                   *
! * Illinois Rocstar LLC                                              *
! * Champaign, IL                                                     *
! * www.illinoisrocstar.com                                           *
! * sales@illinoisrocstar.com                                         *
! *                                                                   *
! * License: See LICENSE file in top level of distribution package or *
! * http://opensource.org/licenses/NCSA                               *
! *********************************************************************
! *********************************************************************
! * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
! * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
! * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
! * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
! * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *
! * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
! * Arising FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
! * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
! *********************************************************************
!******************************************************************************
!
! Purpose: treatments before solution dump to Hdf files.
!
! Description: none.
!
! Input: globalGenx = global data structure (contains name of the window)
!
! Output: none.
!
! Notes: none.
!
!******************************************************************************
!
! $Id: Fluid_preHdfOutput.F90,v 1.7 2008/12/06 08:44:00 mtcampbe Exp $
!
! Copyright: (c) 2002 by the University of Illinois
!
!******************************************************************************

SUBROUTINE fluid_prehdfoutput( globalGenx )

  USE moddatatypes
  USE modrocstar, ONLY       : t_globalgenx
  USE modglobal, ONLY     : t_global
  USE moddatastruct, ONLY : t_region
  USE moderror
  USE modparameters

  IMPLICIT NONE
  include 'roccomf90.h'

! ... parameters
  TYPE(t_globalgenx), POINTER  :: globalgenx
  TYPE(t_global), POINTER      :: global
  TYPE(t_region), POINTER      :: regions(:)

! ... loop variables
  INTEGER :: ireg, istat

! ... local variables
  INTEGER :: ilev
  REAL(RFREAL) :: eps

!******************************************************************************

  global  => globalgenx%global

#ifdef RFLO
  regions => globalgenx%regions
#endif
#ifdef RFLU
  regions => globalgenx%levels(1)%regions
#endif

  CALL registerfunction( global,'Fluid_preHdfOutput',&
  'Fluid_preHdfOutput.F90' )

! set tav from accumulated values ---------------------------------------------
! to be moved to a wrapper

  eps = 100._rfreal*epsilon( 1._rfreal )

#ifdef STATS

#ifdef RFLO
  DO ireg=1,global%nRegions
    IF (regions(ireg)%procid==global%myProcid .AND. &   ! region active and
        regions(ireg)%active==active) THEN              ! on my processor
      ilev = regions(ireg)%currLevel
#endif
#ifdef RFLU
  DO ireg=1,global%nRegionsLocal
#endif

      IF ((global%flowType==flow_unsteady) .AND. (global%doStat==active)) THEN
        IF (global%integrTime > eps) THEN
          IF (global%mixtNStat > 0) THEN
            DO istat=1,global%mixtNStat

#ifdef RFLO
              regions(ireg)%levels(ilev)%mixt%tav(istat,:) = &
              regions(ireg)%levels(ilev)%mixt%tav(istat,:)/global%integrTime
#endif
#ifdef RFLU
              regions(ireg)%mixt%tav(istat,:) = &
              regions(ireg)%mixt%tav(istat,:)/global%integrTime
#endif

            ENDDO
          ENDIF  ! mixtNstat

#ifdef TURB
          IF ((global%turbActive .EQV. .true.) .AND. (global%turbNStat > 0)) THEN
            DO istat=1,global%turbNStat

#ifdef RFLO
              regions(ireg)%levels(ilev)%turb%tav(istat,:) = &
              regions(ireg)%levels(ilev)%turb%tav(istat,:)/global%integrTime
#endif
#ifdef RFLU
              regions(ireg)%turb%tav(istat,:) = &
              regions(ireg)%turb%tav(istat,:)/global%integrTime
#endif

            ENDDO
          ENDIF  ! turbNstat
#endif

        ENDIF    ! integrTime
      ENDIF      ! unsteady and dostat

#ifdef RFLO
    ENDIF        ! region on this processor and active
  ENDDO          ! iReg
#endif
#ifdef RFLU
  ENDDO          ! iReg
#endif

#endif

  CALL deregisterfunction( global )

END SUBROUTINE fluid_prehdfoutput

!******************************************************************************
!
! RCS Revision history:
!
! $Log: Fluid_preHdfOutput.F90,v $
! Revision 1.7  2008/12/06 08:44:00  mtcampbe
! Updated license.
!
! Revision 1.6  2008/11/19 22:17:14  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.5  2006/01/07 10:19:46  wasistho
! added Rocflu treatment
!
! Revision 1.4  2005/12/08 19:56:44  wasistho
! modified purpose (cosmetics)
!
! Revision 1.3  2005/12/08 07:03:43  wasistho
! added ifdef TURB
!
! Revision 1.2  2005/12/08 02:04:57  wasistho
! added register deregister function
!
! Revision 1.1  2005/12/08 00:22:04  wasistho
! initial import
!
!
!******************************************************************************