WriteThrust.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 thrust history into a file.
!
! Description: none.
!
! Input: global%thrustTotal = total thrust
!
! Output: to file.
!
! Notes: none.
!
!******************************************************************************
!
! $Id: WriteThrust.F90,v 1.3 2008/12/06 08:44:10 mtcampbe Exp $
!
! Copyright: (c) 2003 by the University of Illinois
!
!******************************************************************************

SUBROUTINE writethrust( global )

  USE moddatatypes
  USE moderror
  USE modglobal, ONLY : t_global
  USE modmpi
  USE modparameters
  IMPLICIT NONE

! ... parameters
  TYPE(t_global), POINTER :: global

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

#ifdef MPI
  REAL(RFREAL) :: localthrust(2), globalthrust(2)
#endif

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

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

! sum up data from other processors

#ifdef MPI
  localthrust(1) = global%thrustMom
  localthrust(2) = global%thrustPress

  CALL mpi_reduce( localthrust,globalthrust,2,mpi_rfreal,mpi_sum,masterproc, &
                   global%mpiComm,global%mpierr )
  IF (global%mpierr /= 0) CALL errorstop( global,err_mpi_trouble,__line__ )

  global%thrustMom   = globalthrust(1)
  global%thrustPress = globalthrust(2)
#endif
  global%thrustTotal = global%thrustMom + global%thrustPress

! steady flow

  IF (global%flowType==flow_steady .AND. global%myProcid==masterproc) THEN
    WRITE(if_thrust,1000,err=10) global%currentIter,global%thrustMom, &
                                 global%thrustPress,global%thrustTotal

! unsteady flow

  ELSE IF (global%flowType==flow_unsteady .AND. global%myProcid==masterproc) THEN
    WRITE(if_thrust,2000,err=10) global%currentTime,global%thrustMom, &
                                 global%thrustPress,global%thrustTotal
  ENDIF

! close and open file (instead of fflush)

  IF (global%thrustOpenClose .AND. global%myProcid==masterproc) THEN
    WRITE(fname,'(A)') trim(global%outDir)//trim(global%casename)//'.thr'
    CLOSE(if_thrust)
    OPEN(if_thrust,file=fname,form='FORMATTED',status='OLD',position='APPEND')
  ENDIF

! finalize

  CALL deregisterfunction( global )
  goto 999

10   CONTINUE
  CALL errorstop( global,err_file_write,__line__,'Thrust history file.' )

1000 FORMAT(i6,1pe13.4,2e13.4)
2000 FORMAT(1pe12.5,3e13.4)

999  CONTINUE

END SUBROUTINE writethrust

!******************************************************************************
!
! RCS Revision history:
!
! $Log: WriteThrust.F90,v $
! Revision 1.3  2008/12/06 08:44:10  mtcampbe
! Updated license.
!
! Revision 1.2  2008/11/19 22:17:24  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.1  2004/12/01 16:52:27  haselbac
! Initial revision after changing case
!
! Revision 1.1  2003/06/02 17:12:00  jblazek
! Added computation of thrust.
!
!******************************************************************************