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         : T. PEDOT (November) Add Probes               !
  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, l, k, iface, ibnd
  28 | 
  29 |         CHARACTER*80  :: GFile, SrFile, HFile, QwFile, QrFile, Qpfile
  30 |         CHARACTER*80  :: FMT
  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 |         DO ibnd = 1,is_nbfaces
  54 |  
  55 |           j = is_bcell(ibnd)
  56 |           i = is_bface(ibnd)
  57 |           WRITE(FILE_Qw) i,j,Qw(ibnd)
  58 |           WRITE(FILE_H ) i,j,Htot(ibnd)
  59 | 
  60 |         ENDDO
  61 | 
  62 |         IF(i_dom_nprobes.gt.0) THEN
  63 |           OPEN(FILE_QRP,FILE=Qpfile,FORM='FORMATTED')
  64 |           WRITE(FILE_QRP,*) i_dom_nprobes,'pts Inc_X Inc_Y Inc_Z'
  65 |           FMT ='(I6,2X,I9,2X,3(1X,E15.6))'
  66 | 
  67 |           DO i = 1,i_dom_nprobes
  68 |             WRITE(FILE_QRP,FMT) i,is_norm_probe(i),(Q_ptot(k,i),k=1,3)
  69 |           ENDDO
  70 | 
  71 |         ENDIF
  72 | 
  73 |         CLOSE(FILE_G )
  74 |         CLOSE(FILE_Sr)
  75 |         CLOSE(FILE_H )
  76 |         CLOSE(FILE_Qw)
  77 |         CLOSE(FILE_Qr)
  78 | 
  79 | 
  80 | 
  81 |         WRITE(*,*) "    -------------"
  82 |         WRITE(*,*) "    *-> H max  = ", MAXVAL(Htot)
  83 |         WRITE(*,*) "    *-> H min  = ", MINVAL(Htot)
  84 |         WRITE(*,*) "    -------------"
  85 |         WRITE(*,*) "    *-> Qw max = ", MAXVAL(Qw)
  86 |         WRITE(*,*) "    *-> Qw min = ", MINVAL(Qw)
  87 |         WRITE(*,*) "    -------------"
  88 |         WRITE(*,*) "    *-> Sr max = ", MAXVAL(Srtot)
  89 |         WRITE(*,*) "    *-> Sr min = ", MINVAL(Srtot)
  90 |         WRITE(*,*) "    -------------"
  91 |         WRITE(*,*)
  92 | 
  93 |       END SUBROUTINE outprocessing