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

SUBROUTINE rflu_mergepostprocessregions(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_modallocatememory
  USE rflu_modboundlists
  USE rflu_modcellmapping
  USE rflu_moddeallocatememory
  USE rflu_moddimensions
  USE rflu_modensight
  USE rflu_modextractflowdata
  USE rflu_modfacelist
  USE rflu_modgeometry
  USE rflu_modinterpolation
  USE rflu_modmergeregions
  USE rflu_modpatchcoeffs
  USE rflu_modplottingvars
  USE rflu_modprobes
  USE rflu_modreadbcinputfile
  USE rflu_modreadwriteflow
  USE rflu_modreadwritegrid
  USE rflu_modrenumberings
  USE rflu_modstencilscells
  USE rflu_modstencilsutils
  USE rflu_modstencilsvert
  USE rflu_modvertexlists   
  USE rflu_modweights  

#ifndef NO_TECPLOT  
  USE rflu_modtecplot
#endif

#ifdef PLAG
  USE plag_moddimensions, ONLY: plag_printdimensions, &
                                plag_setdimensions
#endif

  USE modinterfaces, ONLY: rflu_allocmemsolwrapper, &
                           rflu_allocmemvertwrapper, & 
                           rflu_computeexactflowerror, &
                           rflu_computeexactflowprobeerror, & 
                           rflu_creategrid, & 
                           rflu_deallocmemsolwrapper, &
                           rflu_deallocmemvertwrapper, & 
                           rflu_decidebuildgeometry, &
                           rflu_destroygrid, &  
                           rflu_interpolatewrapper, &
                           rflu_printflowinfowrapper, &
                           rflu_printgridinfo, &
                           rflu_printheader, &
                           rflu_setpatchplotflags, & 
                           rflu_setvarinfowrapper, &
                           rflu_setvars, & 
                           rflu_setvarswrapper
  
  IMPLICIT NONE

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

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

  TYPE(t_level), POINTER :: levels(:)

! ==============================================================================
! Local variables
! ==============================================================================

  CHARACTER(CHRLEN) :: casename,rcsidentstring
  INTEGER :: errorflag,ipatchserial,ireg
  TYPE(t_region), POINTER :: pregion,pregionserial
  TYPE(t_global), POINTER :: global
  TYPE(t_patch), POINTER :: ppatchserial  

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

  rcsidentstring = & 
    '$RCSfile: RFLU_MergePostProcessRegions.F90,v $ $Revision: 1.29 $'

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

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

#ifndef NO_TECPLOT
! ******************************************************************************
! Initialize TECPLOT interface
! ******************************************************************************

  IF ( global%postOutputFormat == post_output_format_tecplot ) THEN 
    CALL rflu_tec_init(global)
  END IF ! global%postOutputFormat 
#endif

! ******************************************************************************
! Allocate memory for serial region
! ******************************************************************************

  pregionserial => levels(1)%regions(0)
  
  CALL rflu_readdimensionswrapper(pregionserial)
  CALL rflu_creategrid(pregionserial) 

  IF ( pregionserial%grid%nPatches > 0 ) THEN      
    CALL rflu_readbcinputfilewrapper(pregionserial)    
  END IF ! pRegionSerial%grid%nPatches  

  CALL rflu_createcellmapping(pregionserial)
  CALL rflu_readloc2globcellmapping(pregionserial)
  CALL rflu_buildglob2loccellmapping(pregionserial) 
    
  CALL rflu_allocmemsolwrapper(pregionserial) 
  CALL rflu_allocatememorygspeeds(pregionserial)     
  CALL rflu_setvarinfowrapper(pregionserial) 
  
  IF ( global%patchCoeffFlag .EQV. .true. ) THEN 
    CALL rflu_createpatchcoeffs(pregionserial)   
  END IF ! global%patchCoeffFlag

#ifdef PLAG
  IF ( global%plagUsed .EQV. .true. ) THEN 
    CALL plag_setdimensions(pregionserial,0)
  END IF ! global%plagUsed
#endif

! ******************************************************************************  
! Loop over regions and merge region data into serial region
! ******************************************************************************

  DO ireg = 1,global%nRegionsLocal
    pregion => levels(1)%regions(ireg)

    CALL rflu_readdimensionswrapper(pregion)
    CALL rflu_creategrid(pregion) 

    IF ( pregion%grid%nPatches > 0 ) THEN      
      CALL rflu_readbcinputfilewrapper(pregion)    
    END IF ! pRegion%grid%nPatches

    CALL rflu_readgridwrapper(pregion)  

    IF ( global%verbLevel > verbose_none ) THEN 
      CALL rflu_printgridinfo(pregion)
    END IF ! global%verbLevel

    CALL rflu_createcellmapping(pregion)
    CALL rflu_readloc2globcellmapping(pregion)
    CALL rflu_buildglob2loccellmapping(pregion) 

    CALL rflu_rnmb_createpc2scmap(pregion)
    CALL rflu_rnmb_createpv2svmap(pregion)
    CALL rflu_rnmb_createpbf2sbfmap(pregion)    
    CALL rflu_rnmb_readpxx2sxxmaps(pregion)

    CALL rflu_allocmemsolwrapper(pregion) 
    CALL rflu_setvarinfowrapper(pregion)
  
    IF ( global%patchCoeffFlag .EQV. .true. ) THEN   
      CALL rflu_createpatchcoeffs(pregion)
    END IF ! global%patchCoeffFlag 

    IF ( global%postPlotType == plot_grid_flow ) THEN
      CALL rflu_readflowwrapper(pregion)
      
      IF ( global%patchCoeffFlag .EQV. .true. ) THEN 
        CALL rflu_readpatchcoeffswrapper(pregion)
      END IF ! global%patchCoeffFlag
      
      IF ( global%verbLevel > verbose_none ) THEN 
        CALL rflu_printflowinfowrapper(pregion)
      END IF ! global%verbLevel       
    END IF ! global%postPlotType
            
    CALL rflu_merg_mergegrid(pregion,pregionserial)
    
    IF ( global%postPlotType == plot_grid_flow ) THEN
      CALL rflu_merg_mergesolwrapper(pregion,pregionserial)
      
      IF ( global%patchCoeffFlag .EQV. .true. ) THEN 
        CALL rflu_merg_mergepatchcoeffs(pregion,pregionserial)
      END IF ! global%patchCoeffFlag    
    END IF ! global%postPlotType

    CALL rflu_rnmb_destroypc2scmap(pregion)
    CALL rflu_rnmb_destroypv2svmap(pregion)
    CALL rflu_rnmb_destroypbf2sbfmap(pregion)

    IF ( global%patchCoeffFlag .EQV. .true. ) THEN 
      CALL rflu_destroypatchcoeffs(pregion)
    END IF ! global%patchCoeffFlag  
      
    CALL rflu_deallocmemsolwrapper(pregion)
    CALL rflu_destroycellmapping(pregion) 
    CALL rflu_destroygrid(pregion)
  END DO ! iReg    

! ******************************************************************************  
! Serial region is now available
! ******************************************************************************

  global%nRegions      = 1
  global%nRegionsLocal = 1

  IF ( global%verbLevel > verbose_none ) THEN 
    CALL rflu_printgridinfo(pregionserial)
  END IF ! global%verbLevel
  
! ==============================================================================
! Define dependent variables of merged region
! ==============================================================================

  IF ( global%postPlotType == plot_grid_flow ) THEN
    CALL rflu_setvarswrapper(pregionserial,1,pregionserial%grid%nCellsTot)  

    IF ( global%verbLevel > verbose_none ) THEN 
      CALL rflu_printflowinfowrapper(pregionserial)
    END IF ! global%verbLevel       

#ifdef PLAG
    IF ( global%plagUsed .EQV. .true. ) THEN 
      IF ( global%verbLevel > verbose_none ) THEN 
        CALL plag_printdimensions(pregionserial)
      END IF ! global%verbLevel
    END IF ! global%plagUsed
#endif
  END IF ! global%postPlotType

! ==============================================================================
! Write files for merged region
! ==============================================================================

  IF ( global%postWriteMergeFlag .EQV. .true. ) THEN
    CALL rflu_writedimensionswrapper(pregionserial,write_dimens_mode_force)
    CALL rflu_writeloc2globcellmapping(pregionserial)
    CALL rflu_writegridwrapper(pregionserial)

    IF ( global%postPlotType == plot_grid_flow ) THEN
      CALL rflu_writeflowwrapper(pregionserial)
    END IF ! global%postPlotType
  END IF ! global%postWriteMergeFlag

! ==============================================================================
! Build data structure and compute geometry of merged region
! ============================================================================== 

  CALL rflu_createbvertexlists(pregionserial)
  CALL rflu_buildbvertexlists(pregionserial)              

  CALL rflu_createfacelist(pregionserial)
  CALL rflu_buildfacelist(pregionserial)  
  CALL rflu_renumberbfacelists(pregionserial)      

  IF ( rflu_decidebuildgeometry(pregionserial%global) .EQV. .true. ) THEN        
    CALL rflu_creategeometry(pregionserial)
    CALL rflu_buildgeometry(pregionserial)
  END IF ! RFLU_DecideBuildGeometry

! ******************************************************************************
! Carry out specific operations on solution. NOTE this can only be done if have
! read solution.
! ******************************************************************************

  IF ( global%postPlotType == plot_grid_flow ) THEN

! ============================================================================== 
!   Create and compute plotting variables. NOTE always use order 1 for cell
!   gradients.
! ============================================================================== 
        
    IF ( rflu_decidecomputeplottingvars(pregionserial) .EQV. .true. ) THEN
      CALL rflu_countplottingvars(pregionserial)
      CALL rflu_createplottingvarmaps(pregionserial)
      CALL rflu_buildplottingvarmaps(pregionserial)
      CALL rflu_printplottingvarsinfo(pregionserial)
      CALL rflu_createplottingvars(pregionserial)          

      IF ( (global%postDiscFlag .EQV. .true.) .OR. & 
           (global%postVortFlag .EQV. .true.) .OR. & 
           (global%postVortCoreFlag .EQV. .true.) ) THEN
        CALL rflu_createvert2celllist(pregionserial)
        CALL rflu_buildvert2celllist(pregionserial)

        CALL rflu_setinfoc2cstencilwrapper(pregionserial,1)
        CALL rflu_createc2cstencilwrapper(pregionserial)
        CALL rflu_buildc2cstencilwrapper(pregionserial)

        CALL rflu_destroyvert2celllist(pregionserial)

        CALL rflu_createwtsc2cwrapper(pregionserial,1)
        CALL rflu_computewtsc2cwrapper(pregionserial,1)
      END IF ! RFLU_DecideComputePlottingVars

      CALL rflu_computeplottingvarswrapper(pregionserial)       

      IF ( (global%postDiscFlag .EQV. .true.) .OR. &
           (global%postVortFlag .EQV. .true.) .OR. &
           (global%postVortCoreFlag .EQV. .true.) ) THEN
        CALL rflu_destroywtsc2cwrapper(pregionserial)   
        CALL rflu_destroyc2cstencilwrapper(pregionserial)        
      END IF ! global%postDiscFlag
    END IF ! pRegionSerial%mixtInput%fluidModel

! ============================================================================== 
!   Compute errors
! ============================================================================== 
        
    IF ( global%initFlowFlag == initflow_fromhardcode ) THEN          
      CALL rflu_computeexactflowerror(pregionserial)  
    END IF ! global%initFlowFlag

! ==============================================================================
!   Compute errors at probe locations
! ==============================================================================

    IF ( (global%nProbes > 0) .AND. & 
         (global%postCompErrFlag .EQV. .true.) .AND. & 
         (global%flowType == flow_unsteady) .AND. &        
         (global%initFlowFlag == initflow_fromhardcode) ) THEN
      CALL rflu_createcell2facelist(pregionserial)
      CALL rflu_buildcell2facelist(pregionserial)
      CALL rflu_findprobecells(pregionserial) 
      CALL rflu_destroycell2facelist(pregionserial)
      CALL rflu_openprobefiles(pregionserial)
      CALL rflu_computeexactflowprobeerror(pregionserial)
      CALL rflu_closeprobefiles(pregionserial)
    END IF ! global%nProbes

! ==============================================================================
!   Extract data from solution
! ==============================================================================

    IF ( global%postExtractFlag .EQV. .true. ) THEN 
      CALL rflu_extractflowdata(pregionserial)
    END IF ! global%postExtractFlag
  END IF ! global%postPlotType

! ******************************************************************************
! Interpolate data from cell centers to vertices
! ******************************************************************************

  IF ( (global%postPlotType == plot_grid_flow) .AND. & 
       (global%postInterpType /= interp_type_none) ) THEN 
    CALL rflu_createvert2celllist(pregionserial)
    CALL rflu_buildvert2celllist(pregionserial)

    IF ( global%postInterpType == interp_type_proper ) THEN 
      CALL rflu_setinfostencilvert2cell(pregionserial,global%postInterpOrder)    
      CALL rflu_createstencilvert2cell(pregionserial)
      CALL rflu_buildstencilvert2cell(pregionserial)
    END IF ! global%postInterpType
      
    CALL rflu_allocmemvertwrapper(pregionserial)
    CALL rflu_interpolatewrapper(pregionserial)
    
    IF ( global%postInterpType == interp_type_proper ) THEN
      CALL rflu_destroystencilvert2cell(pregionserial) 
    END IF ! global%postInterpType
         
    CALL rflu_destroyvert2celllist(pregionserial)
  END IF ! global%postPlotType

! ******************************************************************************
! Write out files for postprocessing
! ******************************************************************************

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

      IF ( global%patchCoeffFlag .EQV. .true. ) THEN    
        CALL rflu_tec_openfilepatch(pregionserial)  
      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

      IF ( global%postPlotVolFlag .EQV. .true. ) THEN
        CALL rflu_tec_builddatafieldvol(pregionserial)      
        CALL rflu_tec_writefilefieldvol(pregionserial)         
        CALL rflu_tec_destroydatafieldvol(pregionserial)
      END IF ! global%postPlotVolFlag

      CALL rflu_setpatchplotflags(pregionserial)

      DO ipatchserial = 1,pregionserial%grid%nPatches
        ppatchserial => pregionserial%patches(ipatchserial)

        IF ( ppatchserial%plotFlag .EQV. .true. ) THEN 
          CALL rflu_tec_builddatafieldsurf(pregionserial,ppatchserial)             
          CALL rflu_tec_writefilefieldsurf(pregionserial,ppatchserial)
          CALL rflu_tec_destroydatafieldsurf(pregionserial,ppatchserial) 
        END IF ! pPatchSerial%plotFlag
      END DO ! iPatchSerial

      IF ( global%patchCoeffFlag .EQV. .true. ) THEN   
        DO ipatchserial = 1,pregionserial%grid%nPatches
          ppatchserial => pregionserial%patches(ipatchserial)

          CALL rflu_tec_builddatapatch(pregionserial,ppatchserial)  
          CALL rflu_tec_writefilepatch(pregionserial,ppatchserial)
          CALL rflu_tec_destroydatapatch(pregionserial,ppatchserial)
        END DO ! iPatchSerial
      END IF ! global%patchCoeffFlag

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

        DO ipatchserial = 1,pregionserial%grid%nPatches
          ppatchserial => pregionserial%patches(ipatchserial)

          IF ( ppatchserial%plotStatsFlag .EQV. .true. ) THEN
            CALL rflu_tec_builddatapatchstats(pregionserial,ppatchserial)
            CALL rflu_tec_writefilepatch(pregionserial,ppatchserial)
            CALL rflu_tec_destroydatapatch(pregionserial,ppatchserial)
          END IF ! pPatchSerial%plotStatsFlag
        END DO ! iPatchSerial
      END IF ! global%plagUsed
#endif

      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 ( post_output_format_ensight ) 
      pregion => levels(1)%regions(0)
      CALL rflu_ens_builddatainfo(pregion,1) 
      CALL rflu_ens_writefilecase(global,1) 
      
      CALL rflu_ens_openfilegeometry(global,1)
      CALL rflu_ens_openfileflowwrapper(global,1)        
              
      CALL rflu_ens_storepartnumber(global)   
      CALL rflu_ens_writegridwrapper(pregionserial)    
      CALL rflu_ens_writeflowwrapper(pregionserial)        
      
      CALL rflu_ens_closefilegeometry(global)
      CALL rflu_ens_closefileflowwrapper(global)      
      CALL rflu_ens_destroydatainfo(global)      
    CASE default
      CALL errorstop(global,err_reached_default,__line__)
  END SELECT ! global%postOutputFormat

! ******************************************************************************
! Deallocate memory
! ******************************************************************************

  IF ( rflu_decidecomputeplottingvars(pregionserial) .EQV. .true. ) THEN
    CALL rflu_destroyplottingvarmaps(pregionserial)
    CALL rflu_destroyplottingvars(pregionserial)
  END IF ! RFLU_DecideComputePlottingVars
  
  IF ( rflu_decidebuildgeometry(pregionserial%global) .EQV. .true. ) THEN   
    CALL rflu_destroygeometry(pregionserial)
  END IF ! RFLU_DecideBuildGeometry  

  CALL rflu_destroyfacelist(pregionserial)                        
  CALL rflu_destroybvertexlists(pregionserial)
  CALL rflu_destroycellmapping(pregionserial)            
  CALL rflu_deallocmemsolwrapper(pregionserial)
  CALL rflu_deallocatememorygspeeds(pregionserial)
  
  IF ( global%patchCoeffFlag .EQV. .true. ) THEN  
    CALL rflu_destroypatchcoeffs(pregionserial)
  END IF ! global%patchCoeffFlag     
    
  CALL rflu_destroygrid(pregionserial)   

  IF ( (global%postPlotType == plot_grid_flow) .AND. & 
       (global%postInterpType /= interp_type_none) ) THEN 
    CALL rflu_deallocmemvertwrapper(pregionserial)
  END IF ! global%postPlotType
  
! ******************************************************************************
! End
! ******************************************************************************

  CALL deregisterfunction(global)

END SUBROUTINE rflu_mergepostprocessregions

! ******************************************************************************
!
! RCS Revision history:
!
! $Log: RFLU_MergePostProcessRegions.F90,v $
! Revision 1.29  2008/12/06 08:45:06  mtcampbe
! Updated license.
!
! Revision 1.28  2008/11/19 22:18:16  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.27  2007/03/27 00:46:48  haselbac
! Completed merging of particle solution
!
! Revision 1.26  2007/03/19 21:42:09  haselbac
! Adapted to changes related to plotting variables
!
! Revision 1.25  2006/04/07 15:19:26  haselbac
! Removed tabs
!
! Revision 1.24  2006/01/24 21:16:20  mparmar
! Now pass pRegion instead of global into RFLU_ENS_BuildDataInfo
!
! Revision 1.23  2006/01/06 22:21:40  haselbac
! Adapted to name changes
!
! Revision 1.22  2005/12/10 16:56:03  haselbac
! Added RFLU_DecideComputePlottingVars and made accomp changes
!
! Revision 1.21  2005/11/10 02:48:41  haselbac
! Removed superfluous region pointers, now print grid info after merging
!
! Revision 1.20  2005/11/04 15:13:48  haselbac
! Added priniting of flow info after merging
!
! Revision 1.19  2005/10/27 19:22:30  haselbac
! Adapted to changes in stencil routine names
!
! Revision 1.18  2005/10/05 20:13:14  haselbac
! Added support for ENSIGHT
!
! Revision 1.17  2005/10/05 16:20:11  haselbac
! Bug fix: Added missing USE RFLU_ModStencilCells
!
! Revision 1.16  2005/10/05 14:27:44  haselbac
! Adapted to changes in stencil modules, added use of vertex list module
!
! Revision 1.15  2005/09/21 19:41:51  haselbac
! Bug fix: Incorrect region pointer for RFLU_SetPatchPlotFlags
!
! Revision 1.14  2005/08/18 18:48:16  haselbac
! Added postVortCoreFlag to IF statements
!
! Revision 1.13  2005/08/10 00:36:13  haselbac
! Adapted to changes in RFLU_ModPlottingVars
!
! Revision 1.12  2005/08/09 01:10:49  haselbac
! Added IFs for patch coeffs, writing serial files, plotting patches
!
! Revision 1.11  2005/07/25 12:23:30  haselbac
! Added vorticity plotting variables
!
! Revision 1.10  2005/06/11 20:35:03  haselbac
! Bug fix: Added IFs for postPlotType, now read dims for serial region at t = 0 always
!
! Revision 1.9  2005/05/03 20:39:34  haselbac
! Only compute plotting vars if postDiscFlag is TRUE
!
! Revision 1.8  2005/05/01 14:22:50  haselbac
! Added postprocessing of plotting vars
!
! Revision 1.7  2005/04/29 12:53:15  haselbac
! Added calls to compute errors at probe locations
!
! Revision 1.6  2005/04/15 15:08:32  haselbac
! Added merging of patch coeffs, modified call to RFLU_SetVarsWrapper
!
! Revision 1.5  2005/01/30 22:04:55  haselbac
! Bug fix for no interpolation of data
!
! Revision 1.4  2005/01/20 14:53:10  haselbac
! Adapted to changed RNMB routines
!
! Revision 1.3  2005/01/17 19:51:44  haselbac
! Bug fix: Now pass proper pointer to RFLU_AllocateMemoryGSpeeds
!
! 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
!
! ******************************************************************************