RFLU_WriteIntegrals1245OLES.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: Write optimal LES integrals in ASCII ROCFLU format.
!
! Description: None.
!
! Input:
!   pRegion     Pointer to region
!
! Output: None.
!
! Notes: None.
!
!******************************************************************************
!
! $Id: RFLU_WriteIntegrals1245OLES.F90,v 1.7 2008/12/06 08:44:30 mtcampbe Exp $
!
! Copyright: (c) 2001 by the University of Illinois
!
!******************************************************************************

SUBROUTINE rflu_writeintegrals1245oles(pRegion)

  USE modglobal, ONLY: t_global
  USE moddatatypes
  USE modparameters
  USE moderror
  USE modbndpatch, ONLY: t_patch
  USE modgrid, ONLY: t_grid
  USE moddatastruct, ONLY: t_region    
  USE modmpi
  
  USE rflu_modoles

  IMPLICIT NONE
  
! ******************************************************************************
! Declarations and definitions
! ******************************************************************************  
   
! ==============================================================================
! Local variables
! ==============================================================================

  INTEGER :: errorflag,hloc,i,icmp,ic1l,ic2l,ic3l,ic4l,ifcp,ifile,ifun,j,k, & 
             l,m,ncells,vloc
  CHARACTER(CHRLEN) :: ifilename,sectionstring,rcsidentstring
  TYPE(t_grid), POINTER :: pgrid
  TYPE(t_global), POINTER :: global
  
! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region), POINTER :: pregion  
  
! ******************************************************************************
! Start, set pointer
! ******************************************************************************

  rcsidentstring = '$RCSfile: RFLU_WriteIntegrals1245OLES.F90,v $ $Revision: 1.7 $'

  pgrid => pregion%grid

! ******************************************************************************
! Open file
! ******************************************************************************

  global => pregion%global

  CALL registerfunction(global,'RFLU_WriteIntegrals1245OLES',&
  'RFLU_WriteIntegrals1245OLES.F90')
  
  ncells = SIZE(pgrid%fsOLES,1)

  IF ( global%myProcid == masterproc .AND. &
       global%verbLevel > verbose_none ) THEN 
    WRITE(stdout,'(A)') solver_name  
    WRITE(stdout,'(A,1X,A)') solver_name,'Writing optimal LES integrals...'
  END IF ! global%verbLevel

  WRITE(ifilename,'(A,I1)') trim(global%outDir)// & 
                            trim(global%casename)//'.int_',ncells

  ifile = if_integ_oles
  OPEN(ifile,file=ifilename,form="FORMATTED",status="UNKNOWN", & 
       iostat=errorflag)
  global%error = errorflag        
  IF ( global%error /= err_none ) THEN 
    CALL errorstop(global,err_file_open,__line__,ifilename)
  END IF ! global%error  
      
! ==============================================================================
! Integral 1
! ==============================================================================

  IF ( global%myProcid == masterproc .AND. &
       global%verbLevel > verbose_low ) THEN 
    WRITE(stdout,'(A,3X,A)') solver_name,'Integral 1...'
  END IF ! global%verbLevel

  sectionstring = '# Integral 1'
  WRITE(ifile,'(A)') trim(sectionstring)
  
!  DO icmp = 1,3
!    DO ifcp = 1,3
!      DO i = 1,3*nCells
!        WRITE(iFile,'(E15.6)') pGrid%int1OLES(icmp,ifcp,i)
!      END DO ! i
!    END DO ! ifcp
!  END DO ! icmp
  
  DO ifcp = 1,3
    DO ic1l = 1,ncells
      DO ifun = 1,9
        CALL rflu_mapk2ij(ifun,l,i)
        vloc = rflu_geti1posoles(l,ic1l)
        
        WRITE(ifile,'(E15.6,3X,I2,3X,I2,2X,2(1X,I2),3X,I3)') & 
              pgrid%int1OLES(i,ifcp,vloc),ifcp,ic1l,l,i,vloc
      END DO ! iFun
    END DO ! ic1l
  END DO ! ifcp
  
! ==============================================================================
! Integral 2
! ==============================================================================

  IF ( global%myProcid == masterproc .AND. &
       global%verbLevel > verbose_low ) THEN 
    WRITE(stdout,'(A,3X,A)') solver_name,'Integral 2...'
  END IF ! global%verbLevel

  sectionstring = '# Integral 2'
  WRITE(ifile,'(A)') trim(sectionstring)
  
!  DO ifcp = 1,3
!    DO i = 1,3*nCells
!      DO j = 1,3*nCells
!        WRITE(iFile,'(2(E15.6,1X))') pGrid%int20OLES(ifcp,i,j), & 
!                                     pGrid%int21OLES(ifcp,i,j)
!      END DO ! j
!    END DO ! i
!  END DO ! ifcp
  
  DO ifcp = 1,3
    DO ic1l = 1,ncells
      DO ic2l = 1,ncells
        DO ifun = 1,9
          CALL rflu_mapk2ij(ifun,l,j)   
          vloc = rflu_geti1posoles(l,ic1l)
          hloc = rflu_getlposoles(j,ic2l)   
          
          WRITE(ifile,1) & 
                pgrid%int20OLES(ifcp,vloc,hloc), & 
                pgrid%int21OLES(ifcp,vloc,hloc), & 
                ifcp,ic1l,ic2l,l,j,vloc,hloc                                  
        END DO ! iFun
      END DO ! ic2l    
    END DO ! ic1l    
  END DO ! ifcp
  
! ==============================================================================
! Integral 4
! ==============================================================================

  IF ( global%myProcid == masterproc .AND. &
       global%verbLevel > verbose_low ) THEN 
    WRITE(stdout,'(A,3X,A)') solver_name,'Integral 4...'
  END IF ! global%verbLevel

  sectionstring = '# Integral 4'
  WRITE(ifile,'(A)') trim(sectionstring)
  
!  DO icmp = 1,3
!    DO ifcp = 1,3
!      DO i = 1,9*nCells*nCells
!        WRITE(iFile,'(3(E15.6,1X))') pGrid%int40OLES(icmp,ifcp,i), & 
!                                     pGrid%int41OLES(icmp,ifcp,i), &
!                                     pGrid%int42OLES(icmp,ifcp,i)
!      END DO ! i
!    END DO ! ifcp
!  END DO ! icmp  
  
  DO ifcp = 1,3
    DO ic1l = 1,ncells
      DO ic2l = 1,ncells
        DO ifun = 1,27
          CALL rflu_mapl2ijk(ifun,l,m,i)    
          vloc = rflu_geti4posoles(l,m,ic1l,ic2l,ncells)   
          
          WRITE(ifile,2) & 
                pgrid%int40OLES(i,ifcp,vloc), & 
                pgrid%int41OLES(i,ifcp,vloc), & 
                pgrid%int42OLES(i,ifcp,vloc), & 
                ifcp,ic1l,ic2l,l,m,i,vloc                             
        END DO ! iFun
      END DO ! ic2l    
    END DO ! ic1l    
  END DO ! ifcp  
  
! ==============================================================================
! Integral 5
! ==============================================================================

  IF ( global%myProcid == masterproc .AND. &
       global%verbLevel > verbose_low ) THEN 
    WRITE(stdout,'(A,3X,A)') solver_name,'Integral 5...'
  END IF ! global%verbLevel

  sectionstring = '# Integral 5'
  WRITE(ifile,'(A)') trim(sectionstring)
  
!  DO ifcp = 1,3
!    DO i = 1,9*nCells*nCells
!      DO j = 1,9*nCells*nCells
!       WRITE(iFile,'(3(E15.6,1X))') pGrid%int50OLES(ifcp,i,j), & 
!                                    pGrid%int51OLES(ifcp,i,j), &
!                                    pGrid%int52OLES(ifcp,i,j)
!      END DO ! j
!    END DO ! i
!  END DO ! ifcp    

  DO ifcp = 1,3
    DO ic1l = 1,ncells
      DO ic2l = 1,ncells
        DO ic3l = 1,ncells
          DO ic4l = 1,ncells
            DO ifun = 1,81
              CALL rflu_mapm2ijkl(ifun,l,m,j,k)
              vloc = rflu_getqposoles(j,k,ic3l,ic4l,ncells)
              hloc = rflu_geti4posoles(l,m,ic1l,ic2l,ncells) 

              WRITE(ifile,3) & 
                    pgrid%int50OLES(ifcp,vloc,hloc), & 
                    pgrid%int51OLES(ifcp,vloc,hloc), & 
                    pgrid%int52OLES(ifcp,vloc,hloc), & 
                    ifcp,ic1l,ic2l,ic3l,ic4l,l,m,j,k,vloc,hloc                                       
            END DO ! iFun
          END DO ! ic4l
        END DO ! ic3l
      END DO ! ic2l    
    END DO ! ic1l    
  END DO ! ifcp

! ==============================================================================
! Write footer
! ==============================================================================

  sectionstring = '# End'
  WRITE(ifile,'(A)') trim(sectionstring)  
  
! ==============================================================================
! Close file
! ==============================================================================

  CLOSE(ifile,iostat=errorflag)
  global%error = errorflag   
  IF ( global%error /= err_none ) THEN 
    CALL errorstop(global,err_file_close,__line__,ifilename)
  END IF ! global%error
   
  IF ( global%myProcid == masterproc .AND. &
       global%verbLevel > verbose_none ) THEN   
    WRITE(stdout,'(A,1X,A)') solver_name,'Writing optimal LES integrals done.'
  END IF ! global%verbLevel

  CALL deregisterfunction(global)
    
1 FORMAT(2(e15.6,1x),2x,i2,2x,2(1x,i2),2x,2(1x,i2),2x,2(1x,i3))
2 FORMAT(3(e15.6,1x),2x,i2,2x,2(1x,i2),2x,3(1x,i2),3x,i3)
3 FORMAT(3(e15.6,1x),2x,i2,2x,4(1x,i2),2x,4(1x,i2),2x,2(1x,i3))        
       
! ******************************************************************************
! End
! ******************************************************************************
 
END SUBROUTINE rflu_writeintegrals1245oles


! ******************************************************************************
!
! RCS Revision history:
!
!   $Log: RFLU_WriteIntegrals1245OLES.F90,v $
!   Revision 1.7  2008/12/06 08:44:30  mtcampbe
!   Updated license.
!
!   Revision 1.6  2008/11/19 22:17:43  mtcampbe
!   Added Illinois Open Source License/Copyright
!
!   Revision 1.5  2006/12/15 13:27:02  haselbac
!   Fixed bug in format statement, found by ifort
!
!   Revision 1.4  2004/01/22 16:04:34  haselbac
!   Changed declaration to eliminate warning on ALC
!
!   Revision 1.3  2002/10/08 15:49:29  haselbac
!   {IO}STAT=global%error replaced by {IO}STAT=errorFlag - SGI problem
!
!   Revision 1.2  2002/10/05 19:35:41  haselbac
!   GENX integration: added outDir to file name
!
!   Revision 1.1  2002/09/09 16:28:02  haselbac
!   Initial revision
!
! ******************************************************************************