1 | include(dom.inc)
   2 | 
   3 |        SUBROUTINE OUTPROCESSING(ncells,nbounds,xc,yc,zc,G,Sr,H,Lo,Q_r,  &
   4 |      &                          Lb,Kp,epsil,nfacemax,path,outpath)
   5 | 
   6 |        IMPLICIT NONE
   7 | 
   8 |        include 'dom_constants.h'
   9 | 
  10 |        DOM_INT :: ncells, nbounds,nfacemax
  11 | 
  12 |        CHARACTER*80 :: path, outpath
  13 |        CHARACTER*80 :: gfile, srfile, hfile, kpfile
  14 |        CHARACTER*80 :: qwfile, qrfile, cfacesfile
  15 | 
  16 |        DOM_REAL :: Qw
  17 |        DOM_REAL , DIMENSION (nfacemax) :: xfc,yfc,zfc
  18 |        DOM_REAL , DIMENSION (ncells) :: G,xc,yc,zc,Lb,Sr,Kp
  19 |        DOM_REAL , DIMENSION (ncells,nfacemax) :: H,Lo,epsil
  20 |        DOM_REAL , DIMENSION (ncells,3) :: Q_r
  21 |   
  22 |        DOM_INT :: i, j, icell, l, infcelt
  23 | 
  24 |        gfile      = outpath(1:len_trim(outpath))//'/G.out'
  25 |        srfile     = outpath(1:len_trim(outpath))//'/Sr.out'
  26 |        hfile      = outpath(1:len_trim(outpath))//'/H.out'
  27 |        qwfile     = outpath(1:len_trim(outpath))//'/Qw.out'
  28 |        qrfile     = outpath(1:len_trim(outpath))//'/Qr.out'
  29 |        kpfile     = outpath(1:len_trim(outpath))//'/Kp.out'
  30 |        cfacesfile = path(1:len_trim(path))//'/Centerfaces.in'
  31 | 
  32 |        OPEN(UNIT=10,FILE=gfile,FORM='unformatted')
  33 |        OPEN(UNIT=20,FILE=srfile,FORM='unformatted')
  34 |        OPEN(UNIT=30,FILE=hfile,FORM='unformatted')
  35 |        OPEN(UNIT=31,FILE=qwfile,FORM='unformatted')
  36 |        OPEN(UNIT=32,FILE=qrfile,FORM='unformatted')
  37 |        OPEN(UNIT=33,FILE=kpfile,FORM='unformatted')
  38 |        OPEN(UNIT=40,FILE=cfacesfile,FORM='unformatted')
  39 | 
  40 |        DO i=1,ncells
  41 | 
  42 | !         print*, " Sr(",i,")=", Sr(i)
  43 | !         print*, "           Lbtot:", Lb(i)
  44 | !         print*, "           Gtot :", G(i)
  45 | 
  46 |           READ(40) icell,infcelt,(xfc(l),yfc(l),zfc(l),l=1,infcelt)
  47 |           WRITE(10) xc(icell),yc(icell),zc(icell),G(icell)
  48 |           WRITE(20) xc(icell),yc(icell),zc(icell),Sr(icell)
  49 |           WRITE(33) xc(icell),yc(icell),zc(icell),Kp(icell)
  50 |           WRITE(32) xc(icell),yc(icell),zc(icell),(Q_r(icell,l),l=1,3)
  51 | 
  52 |           DO j=1,infcelt
  53 |  
  54 |              IF (epsil(icell,j)/=-1.) THEN
  55 |  
  56 |                 WRITE(30) xfc(j),yfc(j),zfc(j),H(icell,j)
  57 | 
  58 |                 Qw=epsil(icell,j)*H(icell,j)-pi*Lo(icell,j)
  59 | 
  60 |                 WRITE(31) xfc(j),yfc(j),zfc(j),Qw                 
  61 | 
  62 |              ENDIF
  63 |      
  64 |           ENDDO
  65 | 
  66 |        ENDDO
  67 | 
  68 |        CLOSE(10)
  69 |        CLOSE(20)
  70 |        CLOSE(30)
  71 |        CLOSE(31)
  72 |        CLOSE(32)
  73 |        CLOSE(33)
  74 |        CLOSE(40)
  75 | 
  76 |        END SUBROUTINE OUTPROCESSING