PLAG_WriteTecplotAscii.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 solution data to plot file in ASCII format.
!
! Description: none.
!
! Input: iReg    = region number
!        iLev    = current level
!        region  = region data (dimensions, plag variables)
!
! Output: to plot file.
!
! Notes: none.
!
!******************************************************************************
!
! $Id: PLAG_WriteTecplotAscii.F90,v 1.3 2008/12/06 08:45:07 mtcampbe Exp $
!
! Copyright: (c) 2003 by the University of Illinois
!
!******************************************************************************

SUBROUTINE plag_writetecplotascii( iReg, iLev, region )

  USE moddatatypes
  USE moderror
  USE modpartlag, ONLY    : t_plag, t_plag_input
  USE moddatastruct, ONLY : t_region
  USE modglobal, ONLY     : t_global

  USE modmpi
  USE modparameters
  USE plag_modparameters
  IMPLICIT NONE

! ... parameters
  INTEGER :: ireg, ilev

  TYPE(t_region) :: region

! ... loop variables
  INTEGER :: i

! ... local variables
  CHARACTER(CHRLEN)   :: rcsidentstring
  CHARACTER(CHRLEN+4) :: fname

  INTEGER :: errorflag, ndimplag
  INTEGER :: pidini,regini,icell,indexi,indexj,indexk
  INTEGER, POINTER :: cvplagmass(:)
  INTEGER, POINTER :: aiv(:,:)

  REAL(RFREAL)          :: x,y,z,diam,evapor,temp,up,vp,wp
  REAL(RFREAL), POINTER :: cv(:,:), dv(:,:)
  REAL(RFREAL)          :: massp,muem,taup

  TYPE(t_global), POINTER :: global

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

  rcsidentstring = '$RCSfile: PLAG_WriteTecplotAscii.F90,v $ $Revision: 1.3 $'

  global => region%global

  CALL registerfunction( global,'PLAG_WriteTecplotAscii', 'PLAG_WriteTecplotAscii.F90' )

! set parameters --------------------------------------------------------------

  ndimplag = region%levels(ilev)%plag%nPcls
  
!  PRINT*,'PLAG_writeTecplotAscii: iReg, nDimPlag = ',iReg,nDimPlag

! open file and write the header ----------------------------------------------

  IF (ireg == 1) THEN
    WRITE(fname,'(A,1PE11.5,A)') &
       trim(global%casename)//'.plag_',global%currentTime,'.dat'
    OPEN(if_plot,file=fname,status='unknown',form='formatted',iostat=errorflag)
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global,err_file_open,__line__,fname )

    IF (global%currentTime <= 0._rfreal) THEN
      WRITE(if_plot,1005,err=10) trim(global%casename),global%timeStamp
    ELSE
      WRITE(if_plot,1005,err=10) trim(global%casename),global%currentTime
    ENDIF
    
    WRITE(if_plot,1010,err=10) 'x y z diam Evapor temp up vp wp tau pidini regini icellp ip jp kp'

  ENDIF   ! iReg=1

! for zero nDimPlag exit ------------------------------------------------------
 
  IF ( ndimplag <= 0 ) goto 1999

! write zone header -----------------------------------------------------------

  WRITE(if_plot,1015) ireg, ndimplag

! write data ------------------------------------------------------------------
! pointer to variables

  aiv => region%levels(ilev)%plag%aiv
  cv  => region%levels(ilev)%plag%cv
  dv  => region%levels(ilev)%plag%dv
  cvplagmass  => region%levels(ilev)%plag%cvPlagMass

  DO i = 1, ndimplag
    x      = cv(cv_plag_xpos,     i)
    y      = cv(cv_plag_ypos,     i)
    z      = cv(cv_plag_zpos,     i)
    diam   = dv(dv_plag_diam,     i)
    evapor = cv(cv_plag_enervapor,i)
    temp   = dv(dv_plag_temp,     i)
    up     = dv(dv_plag_uvel,     i)
    vp     = dv(dv_plag_vvel,     i)
    wp     = dv(dv_plag_wvel,     i)

    pidini = aiv(aiv_plag_pidini, i)
    regini = aiv(aiv_plag_regini, i)
    icell  = aiv(aiv_plag_icells, i)
    indexi = aiv(aiv_plag_indexi, i)
    indexj = aiv(aiv_plag_indexj, i)
    indexk = aiv(aiv_plag_indexk, i)

    massp = sum( cv(cvplagmass(:),i) )
    muem  = 3.6e-04_rfreal
    taup  = massp/(3.0_rfreal*global%pi*muem*diam)

    WRITE(if_plot,1020,err=10) x,y,z,diam,evapor,temp,up,vp,wp,taup, &
       pidini,regini,icell,indexi,indexj,indexk

  ENDDO ! i 

! close file, handle errors ---------------------------------------------------

  IF (ireg == global%nRegions) THEN
    CLOSE(if_plot,iostat=errorflag)
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global,err_file_close,__line__,fname )
  ENDIF

1999 CONTINUE
  CALL deregisterfunction( global )
  goto 999

10   CONTINUE
  CALL errorstop( global,err_file_write,__line__,fname )

! formats ---------------------------------------------------------------------

1005 FORMAT('TITLE="',a,'. Time: ',1pe11.5,'."')
1010 FORMAT('VARIABLES= ',a)
1015 FORMAT('ZONE T="',i5.5,'", I=',i10,', F=POINT')
1020 FORMAT(1p,10(1x,e13.6),1x,i10,5(1x,i8))

999  CONTINUE
END SUBROUTINE plag_writetecplotascii

!******************************************************************************
!
! RCS Revision history:
!
! $Log: PLAG_WriteTecplotAscii.F90,v $
! Revision 1.3  2008/12/06 08:45:07  mtcampbe
! Updated license.
!
! Revision 1.2  2008/11/19 22:18:18  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.1  2004/12/01 22:00:47  fnajjar
! Initial revision after changing case
!
! Revision 1.7  2004/11/18 22:06:30  fnajjar
! Modified formatting to properly write high particle count index
!
! Revision 1.6  2004/11/13 22:03:55  fnajjar
! Included particle timescale in Tecplot output file
!
! Revision 1.5  2004/06/29 14:07:03  fnajjar
! Removed mixture properties from Tecplot file since mixture data not allocated and read
!
! Revision 1.4  2004/05/24 14:23:24  fnajjar
! Added variables to Tecplot file
!
! Revision 1.3  2004/03/05 22:52:00  fnajjar
! Added particle temperature to Tecplot output file and changed extension to .dat
!
! Revision 1.2  2004/03/02 21:51:15  jferry
! Added output of vapor energy to rplagpost output file
!
! Revision 1.1.1.1  2003/05/06 16:14:38  fnajjar
! Import of postprocessing tool for Rocpart
!
!******************************************************************************