RFLU_PostProcessRegions.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: Wrapper for postprocessing results without merging the regions.
!
! Description: None.
!
! Input: 
!   levels      Level data
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************
!
! $Id: RFLU_PostProcessRegions.F90,v 1.26 2008/12/06 08:45:06 mtcampbe Exp $
!
! Copyright: (c) 2004-2005 by the University of Illinois
!
! ******************************************************************************

SUBROUTINE rflu_postprocessregions(levels)

  USE moddatatypes
  USE modglobal, ONLY: t_global 
  USE moderror
  USE moddatastruct, ONLY: t_level,t_region
  USE modbndpatch, ONLY: t_patch
  USE modparameters

  USE rflu_modplottingvars

#ifndef NO_TECPLOT  
  USE rflu_modtecplot
#endif

#ifdef PLAG
  USE plag_modsurfstats
#endif

  USE modinterfaces, ONLY: rflu_closepostinfo, &
                           rflu_destroygrid, &
                           rflu_openpostinfo, &
                           rflu_postprocessregionscommon1, & 
                           rflu_postprocessregionscommon2, &                           
                           rflu_readpostinfo, & 
                           rflu_setpatchplotflags
    
  IMPLICIT NONE

! ******************************************************************************
! Declarations and definitions
! ******************************************************************************

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_level), POINTER :: levels(:)
  
! ==============================================================================
! Local variables
! ==============================================================================

  LOGICAL :: postinfofileexists
  CHARACTER(CHRLEN) :: casename,rcsidentstring
  INTEGER :: errorflag,ipatch,ireg
  TYPE(t_region), POINTER :: pregion
  TYPE(t_global), POINTER :: global
  TYPE(t_patch), POINTER :: ppatch

! ******************************************************************************
! Start
! ******************************************************************************

  rcsidentstring = '$RCSfile: RFLU_PostProcessRegions.F90,v $ $Revision: 1.26 $'

  global => levels(1)%regions(1)%global

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

! ******************************************************************************  
! Open files
! ******************************************************************************

  CALL rflu_openpostinfo(global,file_status_old,postinfofileexists) 

! ******************************************************************************
! Open files for plotting. NOTE for TECPLOT, initialize interface and need to 
! count number of plotting variables first.
! ******************************************************************************

  SELECT CASE ( global%postOutputFormat ) 
#ifndef NO_TECPLOT
    CASE ( post_output_format_tecplot )  
      CALL rflu_tec_init(global)

      IF ( global%nRegions == 1 ) THEN ! single region
        pregion => levels(1)%regions(0)
      ELSE 
        pregion => levels(1)%regions(1)
      END IF ! global%nRegions

      CALL rflu_countplottingvars(pregion)
      CALL rflu_createplottingvarmaps(pregion)
      CALL rflu_buildplottingvarmaps(pregion)
      CALL rflu_printplottingvarsinfo(pregion)

      CALL rflu_tec_openfilefield(pregion)

      CALL rflu_destroyplottingvarmaps(pregion)

      IF ( global%patchCoeffFlag .EQV. .true. ) THEN
        CALL rflu_tec_openfilepatch(pregion)
      END IF ! global%patchCoeffFlag  
    
#ifdef PLAG
      IF ( global%plagUsed .EQV. .true. ) THEN
        CALL rflu_tec_openfilepnt(pregion)
        CALL rflu_tec_openfilepatchstats(pregion)
      END IF ! global%plagUsed
#endif
#endif
    CASE default
      CALL errorstop(global,err_reached_default,__line__)
  END SELECT ! global%postOutputFormat  

! ******************************************************************************  
! Loop over regions
! ******************************************************************************

  DO ireg = 1,global%nRegionsLocal
    IF ( global%nRegions == 1 ) THEN ! single region
      pregion => levels(1)%regions(0)
    ELSE 
      pregion => levels(1)%regions(ireg)
    END IF ! global%nRegions

    IF ( postinfofileexists .EQV. .true. ) THEN 
      CALL rflu_readpostinfo(pregion,infofile_readmode_flag)    
    END IF ! postInfoFileExists  

! ==============================================================================  
!   Check for activation status and write files
! ==============================================================================  

    IF ( pregion%postActiveFlag .EQV. .true. ) THEN 
      CALL rflu_postprocessregionscommon1(pregion,postinfofileexists)
   
      SELECT CASE ( global%postOutputFormat )   
#ifndef NO_TECPLOT   
        CASE ( post_output_format_tecplot )  
          IF ( global%postPlotVolFlag .EQV. .true. ) THEN      
            CALL rflu_tec_builddatafieldvol(pregion)
            CALL rflu_tec_writefilefieldvol(pregion)
            CALL rflu_tec_destroydatafieldvol(pregion)                       
          END IF ! global%postPlotVolFlag

          CALL rflu_setpatchplotflags(pregion)

          DO ipatch = 1,pregion%grid%nPatches
            ppatch => pregion%patches(ipatch)

            IF ( ppatch%plotFlag .EQV. .true. ) THEN 
              CALL rflu_tec_builddatafieldsurf(pregion,ppatch)               
              CALL rflu_tec_writefilefieldsurf(pregion,ppatch)
              CALL rflu_tec_destroydatafieldsurf(pregion,ppatch)
            END IF ! pPatch%plotFlag
          END DO ! iPatch

          IF ( global%patchCoeffFlag .EQV. .true. ) THEN 
            DO ipatch = 1,pregion%grid%nPatches
              ppatch => pregion%patches(ipatch)

              CALL rflu_tec_builddatapatch(pregion,ppatch)  
              CALL rflu_tec_writefilepatch(pregion,ppatch)
              CALL rflu_tec_destroydatapatch(pregion,ppatch) 
            END DO ! iPatch 
          END IF ! global%patchCoeffFlag                                    

#ifdef PLAG
          IF ( global%plagUsed .EQV. .true. ) THEN
            CALL rflu_tec_writefilepnt(pregion)

            DO ipatch = 1,pregion%grid%nPatches
              ppatch => pregion%patches(ipatch)

              IF ( ppatch%plotStatsFlag .EQV. .true. ) THEN
                CALL rflu_tec_builddatapatchstats(pregion,ppatch)  
                CALL rflu_tec_writefilepatch(pregion,ppatch)
                CALL rflu_tec_destroydatapatch(pregion,ppatch) 
              END IF ! pPatch%plotStatsFlag
            END DO ! iPatch 
          END IF ! global%plagUsed
#endif
#endif                     
        CASE default
          CALL errorstop(global,err_reached_default,__line__)
      END SELECT ! global%postOutputFormat

      CALL rflu_postprocessregionscommon2(pregion)
    END IF ! pRegion%postActiveFlag                       
  END DO ! iReg
  
! ******************************************************************************
! Close files
! ******************************************************************************  
 
  SELECT CASE ( global%postOutputFormat )   
#ifndef NO_TECPLOT   
    CASE ( post_output_format_tecplot )  
      CALL rflu_tec_closefilefield(global)

      IF ( global%patchCoeffFlag .EQV. .true. ) THEN
        CALL rflu_tec_closefilepatch(global)
      END IF ! global%patchCoeffFlag 

#ifdef PLAG
      IF ( global%plagUsed .EQV. .true. ) THEN
        CALL rflu_tec_closefilepnt(global)  
        CALL rflu_tec_closefilepatchstats(global)
      END IF ! global%plagUsed
#endif
#endif      
    CASE default
      CALL errorstop(global,err_reached_default,__line__)
  END SELECT ! global%postOutputFormat
  
  IF ( postinfofileexists .EQV. .true. ) THEN 
    CALL rflu_closepostinfo(global) 
  END IF ! postInfoFileExists  
  
! ******************************************************************************
! End
! ******************************************************************************

  CALL deregisterfunction(global)

END SUBROUTINE rflu_postprocessregions

! ******************************************************************************
!
! RCS Revision history:
!
! $Log: RFLU_PostProcessRegions.F90,v $
! Revision 1.26  2008/12/06 08:45:06  mtcampbe
! Updated license.
!
! Revision 1.25  2008/11/19 22:18:17  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.24  2007/03/19 21:44:41  haselbac
! Adapted to changes related to plotting variables
!
! Revision 1.23  2006/04/07 15:19:26  haselbac
! Removed tabs
!
! Revision 1.22  2005/12/01 23:38:55  haselbac
! Bug fix: Added IFs on plagUsed
!
! Revision 1.21  2005/10/05 20:53:00  haselbac
! Fixed missing module bugs
!
! Revision 1.20  2005/10/05 20:22:40  haselbac
! Partial rewrite as part of adding ENSIGHT
!
! Revision 1.19  2005/10/05 16:20:11  haselbac
! Bug fix: Added missing USE RFLU_ModStencilCells
!
! Revision 1.18  2005/10/05 14:32:18  haselbac
! Adapted to changes in stencil modules, added use of vertex list module
!
! Revision 1.17  2005/09/23 19:02:51  haselbac
! Added writing of patch stats files
!
! Revision 1.16  2005/08/18 18:48:16  haselbac
! Added postVortCoreFlag to IF statements
!
! Revision 1.15  2005/08/10 00:36:13  haselbac
! Adapted to changes in RFLU_ModPlottingVars
!
! Revision 1.14  2005/08/09 01:10:49  haselbac
! Added IFs for patch coeffs, writing serial files, plotting patches
!
! Revision 1.13  2005/07/25 12:23:30  haselbac
! Added vorticity plotting variables
!
! Revision 1.12  2005/07/01 21:29:01  haselbac
! Bug fix: Added counting of plotting vars before creating them
!
! Revision 1.11  2005/05/31 16:44:58  haselbac
! Enclosed point TEC routines within ifdef PLAG
!
! Revision 1.10  2005/05/18 22:26:41  fnajjar
! ACH: Adapted point files to multiple regions
!
! Revision 1.9  2005/05/09 20:35:22  haselbac
! Bug fix: Added calls for PLAG patch data back in
!
! Revision 1.8  2005/05/03 20:39:34  haselbac
! Only compute plotting vars if postDiscFlag is TRUE
!
! Revision 1.7  2005/05/01 14:22:50  haselbac
! Added postprocessing of plotting vars
!
! Revision 1.6  2005/04/29 12:53:15  haselbac
! Added calls to compute errors at probe locations
!
! Revision 1.5  2005/04/15 15:08:35  haselbac
! Changed call to RFLU_SetVars
!
! Revision 1.4  2005/01/30 22:04:55  haselbac
! Bug fix for no interpolation of data
!
! Revision 1.3  2005/01/20 14:53:27  haselbac
! Cosmetics only
!
! Revision 1.2  2005/01/03 16:05:49  haselbac
! Adapted to changes in RFLU_ModStencils
!
! Revision 1.1  2004/12/29 20:58:54  haselbac
! Initial revision
!
! ******************************************************************************