1 | include(dom.inc)
   2 | 
   3 |       SUBROUTINE outprocessing
   4 | 
   5 | !       ================================================================!
   6 | !                                                                       !
   7 | !       outprocessing.F  : Creates output files.                        !
   8 | !                                                                       !
   9 | !       out              :                                              !
  10 | !                                                                       !
  11 | !       author           : J. AMAYA (october 2007)                      !
  12 | !                                                                       !
  13 | !       ================================================================!
  14 | 
  15 |         USE mod_inout
  16 |         USE mod_prissma
  17 |         USE mod_slave
  18 | 
  19 |         IMPLICIT NONE
  20 | 
  21 |         include 'dom_constants.h'
  22 | 
  23 | !       IN
  24 | 
  25 | !       LOCAL
  26 |         DOM_INT       :: i, j, l, k, iface, ibnd
  27 | 
  28 |         CHARACTER*80  :: GFile, SrFile, HFile, QwFile, QrFile
  29 | 
  30 | !       ---------------------------------------------------------!
  31 | !       Version 2.04 output files with data in the cell's center !
  32 | !       ---------------------------------------------------------!
  33 | 
  34 |         GFile  = trim(outpath)//'/G.out'
  35 |         SrFile = trim(outpath)//'/Sr.out'
  36 |         QwFile = trim(outpath)//'/Qw.out'
  37 |         QrFile = trim(outpath)//'/Qr.out'
  38 |         HFile  = trim(outpath)//'/H.out'
  39 | 
  40 |         OPEN(FILE_G ,FILE=Gfile ,FORM='UNFORMATTED')
  41 |         OPEN(FILE_Sr,FILE=SrFile,FORM='UNFORMATTED')
  42 |         OPEN(FILE_H ,FILE=HFile ,FORM='UNFORMATTED')
  43 |         OPEN(FILE_Qw,FILE=QwFile,FORM='UNFORMATTED')
  44 |         OPEN(FILE_Qr,FILE=QrFile,FORM='UNFORMATTED')
  45 | 
  46 |         DO j=1,i_dom_nnodes
  47 | 
  48 |           WRITE(FILE_G ) j,Gtot(j)
  49 |           WRITE(FILE_Sr) j,Srtot(j)
  50 |           WRITE(FILE_Qr) j,(Q_rtot(i,j),i=1,3)
  51 | 
  52 |         ENDDO
  53 | 
  54 |         DO ibnd = 1,is_nbfaces
  55 |  
  56 |           j = is_bcell(ibnd)
  57 |           i = is_bface(ibnd)
  58 |           WRITE(FILE_Qw) i,j,Qw(ibnd)
  59 |           WRITE(FILE_H ) i,j,Htot(ibnd)
  60 | 
  61 |         ENDDO
  62 | 
  63 |         CLOSE(FILE_G )
  64 |         CLOSE(FILE_Sr)
  65 |         CLOSE(FILE_H )
  66 |         CLOSE(FILE_Qw)
  67 |         CLOSE(FILE_Qr)
  68 | 
  69 |         WRITE(*,*) " -------------"
  70 |         WRITE(*,*) " *-> H max  = ", MAXVAL(Htot)
  71 |         WRITE(*,*) " *-> H min  = ", MINVAL(Htot)
  72 |         WRITE(*,*) " -------------"
  73 |         WRITE(*,*) " *-> Qw max = ", MAXVAL(Qw)
  74 |         WRITE(*,*) " *-> Qw min = ", MINVAL(Qw)
  75 |         WRITE(*,*) " -------------"
  76 |         WRITE(*,*) " *-> Sr max = ", MAXVAL(Srtot)
  77 |         WRITE(*,*) " *-> Sr min = ", MINVAL(Srtot)
  78 |         WRITE(*,*) " -------------"
  79 | 
  80 |       END SUBROUTINE outprocessing