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 | !       Modified         : D. Poitou (sept2011)                         !
  13 | !                                                                       !
  14 | !       ================================================================!
  15 | 
  16 |         USE mod_inout
  17 |         USE mod_prissma
  18 |         USE mod_slave
  19 | 
  20 |         IMPLICIT NONE
  21 | 
  22 |         include 'dom_constants.h'
  23 | 
  24 | !       IN
  25 | 
  26 | !       LOCAL
  27 |         DOM_INT       :: i, j, k, ibnd
  28 | 
  29 |         CHARACTER*80  :: GFile, SrFile, HFile, QwFile, QrFile, Qpfile
  30 |         CHARACTER*80  :: FMT, file
  31 | 
  32 |         GFile  = trim(outpath)//'/G.out'
  33 |         SrFile = trim(outpath)//'/Sr.out'
  34 |         QwFile = trim(outpath)//'/Qw.out'
  35 |         QrFile = trim(outpath)//'/Qr.out'
  36 |         HFile  = trim(outpath)//'/H.out'
  37 |         Qpfile = trim(outpath)//'/Qr_probe.out'
  38 | 
  39 |         OPEN(FILE_G ,FILE=Gfile ,FORM='UNFORMATTED')
  40 |         OPEN(FILE_Sr,FILE=SrFile,FORM='UNFORMATTED')
  41 |         OPEN(FILE_H ,FILE=HFile ,FORM='UNFORMATTED')
  42 |         OPEN(FILE_Qw,FILE=QwFile,FORM='UNFORMATTED')
  43 |         OPEN(FILE_Qr,FILE=QrFile,FORM='UNFORMATTED')
  44 | 
  45 |         DO j=1,i_dom_nnodes
  46 | 
  47 |           WRITE(FILE_G ) j,Gtot(j)
  48 |           WRITE(FILE_Sr) j,Srtot(j)
  49 |           WRITE(FILE_Qr) j,(Q_rtot(i,j),i=1,3)
  50 | 
  51 |         ENDDO
  52 | 
  53 |         file = path(1:len_trim(path))//'/CLdata.in'
  54 |         OPEN(FILE_CLDAT, FILE=file , FORM='UNFORMATTED')
  55 | 
  56 |         DO ibnd = 1,is_ntot_bfaces
  57 | 
  58 |           READ(FILE_CLDAT) j, i
  59 |           WRITE(FILE_Qw) i,j,Qw(ibnd)
  60 |           WRITE(FILE_H ) i,j,Htot(ibnd)
  61 | 
  62 |         ENDDO
  63 |         CLOSE(FILE_CLDAT)
  64 | 
  65 |         IF(i_dom_nprobes.gt.0) THEN
  66 |           OPEN(FILE_QRP,FILE=Qpfile,FORM='FORMATTED')
  67 |           WRITE(FILE_QRP,*) i_dom_nprobes,'pts Inc_X Inc_Y Inc_Z'
  68 |           FMT ='(I6,2X,I9,2X,3(1X,E15.6))'
  69 | 
  70 |           DO i = 1,i_dom_nprobes
  71 |             WRITE(FILE_QRP,FMT) i,is_norm_probe(i),(Q_ptot(k,i),k=1,3)
  72 |           ENDDO
  73 | 
  74 |         ENDIF
  75 | 
  76 |         CLOSE(FILE_G )
  77 |         CLOSE(FILE_Sr)
  78 |         CLOSE(FILE_H )
  79 |         CLOSE(FILE_Qw)
  80 |         CLOSE(FILE_Qr)
  81 | 
  82 |         PRINT*
  83 |         PRINT*, "    =============================="
  84 |         IF(Srvtot.gt.1.E3) THEN
  85 |           PRINT*, "    -> Total Radiative Heat Loss =", (Srvtot/1000), "kW"
  86 |         ELSE
  87 |           PRINT*, "    -> Total Radiative Heat Loss =", Srvtot, "W"
  88 |         ENDIF
  89 |         PRINT*, "    =============================="
  90 |         PRINT*
  91 | 
  92 |         WRITE(*,*) "    -------------"
  93 |         WRITE(*,*) "    *-> H max  = ", MAXVAL(Htot)
  94 |         WRITE(*,*) "    *-> H min  = ", MINVAL(Htot)
  95 |         WRITE(*,*) "    -------------"
  96 |         WRITE(*,*) "    *-> Qw max = ", MAXVAL(Qw)
  97 |         WRITE(*,*) "    *-> Qw min = ", MINVAL(Qw)
  98 |         WRITE(*,*) "    -------------"
  99 |         WRITE(*,*) "    *-> Sr max = ", MAXVAL(Srtot)
 100 |         WRITE(*,*) "    *-> Sr min = ", MINVAL(Srtot)
 101 |         WRITE(*,*) "    -------------"
 102 |         WRITE(*,*)
 103 | 
 104 |       END SUBROUTINE outprocessing