flowoutput.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	flowoutput (itg, ieg, ntg, nteq, bc, lakon, ntmat_, v, shcon, nshcon, ipkon, kon, co, nflow, dtime, ttime, time, ielmat, prop, ielprop, nactdog, nacteq, iin, physcon, camt, camf, camp, rhcon, nrhcon, vold, jobnamef, set, istartset, iendset, ialset, nset, mi, iaxial, istep, iit)

Function/Subroutine Documentation

flowoutput()

subroutine flowoutput	(	integer, dimension(*)	itg,
		integer, dimension(*)	ieg,
		integer	ntg,
		integer	nteq,
		real*8, dimension(*)	bc,
		character*8, dimension(*)	lakon,
		integer	ntmat_,
		real*8, dimension(0:mi(2),*)	v,
		real*8, dimension(0:3,ntmat_,*)	shcon,
		integer, dimension(*)	nshcon,
		integer, dimension(*)	ipkon,
		integer, dimension(*)	kon,
		real*8, dimension(3,*)	co,
		integer	nflow,
		real*8	dtime,
		real*8	ttime,
		real*8	time,
		integer, dimension(mi(3),*)	ielmat,
		real*8, dimension(*)	prop,
		integer, dimension(*)	ielprop,
		integer, dimension(0:3,*)	nactdog,
		integer, dimension(0:3,*)	nacteq,
		integer	iin,
		real*8, dimension(*)	physcon,
		real*8, dimension(*)	camt,
		real*8, dimension(*)	camf,
		real*8, dimension(*)	camp,
		real*8, dimension(0:1,ntmat_,*)	rhcon,
		integer, dimension(*)	nrhcon,
		real*8, dimension(0:mi(2),*)	vold,
		character*132, dimension(*)	jobnamef,
		character*81, dimension(*)	set,
		integer, dimension(*)	istartset,
		integer, dimension(*)	iendset,
		integer, dimension(*)	ialset,
		integer	nset,
		integer, dimension(*)	mi,
		integer	iaxial,
		integer	istep,
		integer	iit
	)

!     
      implicit none
!     
      logical identity
!
      character*8 lakon(*)
      character*81 set(*)
      character*132 jobnamef(*),fnnet
!     
      integer mi(*),itg(*),ieg(*),ntg,nflow,ielmat(mi(3),*),i,
     &     nrhcon(*),node,iaxial,ider,idirf(8),ieq,imat,kflag,
     &     ntmat_,nteq,nshcon(*),nelem,index,ipkon(*),kon(*),iin,
     &     nactdog(0:3,*),nacteq(0:3,*),ielprop(*),node1,nodem,node2,
     &     istartset(*),iendset(*),ialset(*),nset,nodef(8),numf,
     &     istep,iit
!     
      real*8 physcon(*),v(0:mi(2),*),shcon(0:3,ntmat_,*),co(3,*),
     &     prop(*),dtime,ttime,time,xflow,camt(*),camf(*),camp(*),
     &     rhcon(0:1,ntmat_,*),vold(0:mi(2),*),uamt,uamf,uamp,eta,
     &     bc(*),cp,dvi,df(8),gastemp,f,g(3),r,rho,ts1,ts2,tg1,tg2
!     
      do i=1,132
         if(jobnamef(1)(i:i).eq.' ') exit
      enddo
      i=i-1
      fnnet=jobnamef(1)(1:i)//'.net'
      open(1,file=fnnet,status='unknown')
!
      kflag=3
!
      do i=1,nflow
         nelem=ieg(i)
!     
!        output for gas networks
!
         if((lakon(nelem)(2:5).ne.'LIPI').and.
     &      (lakon(nelem)(2:5).ne.'LICH')) then
!     
            index=ipkon(nelem)
            node1=kon(index+1)
            nodem=kon(index+2)
            node2=kon(index+3)
!
            xflow=v(1,nodem)
!
            if(lakon(nelem)(2:3).ne.'LP') then
!
!              incompressible
!
               if(node1.eq.0) then
c                  tg1=v(0,node2)
c                  tg2=tg1
                  ts1=v(3,node2)
                  ts2=ts1
               elseif(node2.eq.0) then
c                  tg1=v(0,node1)
c                  tg2=tg1
                  ts1=v(3,node1)
                  ts2=ts1
               else
c                  tg1=v(0,node1)
c                  tg2=v(0,node2)
                  ts1=v(3,node1)
                  ts2=v(3,node2)
               endif
               gastemp=(ts1+ts2)/2.d0
            else
!
!              compressible
!
               if(xflow.gt.0) then
c                  tg1=v(0,node1)
c                  ts1=v(3,node1)
c                  tg2=v(0,node2)
c                  ts2=v(3,node2)
c                  gastemp=ts1
                  gastemp=v(3,node1)
               else
c                  tg2=v(0,node1)
c                  ts2=v(3,node1)
c                  tg1=v(0,node2)
c                  ts1=v(3,node2)
c                  gastemp=ts2
                  gastemp=v(3,node2)
               endif
            endif
!
            imat=ielmat(1,nelem)
!
            call materialdata_tg(imat,ntmat_,gastemp,shcon,nshcon,cp,r,
     &         dvi,rhcon,nrhcon,rho)
!
            if(nacteq(2,nodem).ne.0) then
               ieq=nacteq(2,nodem)
c               xflow=v(1,nodem)
!
!              dummy set number
!
               numf=1
!
               call flux(node1,node2,nodem,nelem,lakon,kon,ipkon,
     &              nactdog,identity,
     &              ielprop,prop,kflag,v,xflow,f,nodef,idirf,df,
     &              cp,r,rho,physcon,g,co,dvi,numf,vold,set,shcon,
     &              nshcon,rhcon,nrhcon,ntmat_,mi,ider,ttime,time,
     &              iaxial)
            endif
         endif
!            
         if(lakon(ieg(i))(2:5).eq.'LICH') then
            if((lakon(ieg(i))(6:7).eq.'SG').or.
     &           (lakon(ieg(i))(6:7).eq.'WE').or.
     &           (lakon(ieg(i))(6:7).eq.'DS')) then
               index=ipkon(ieg(i))
               node=kon(index+2)
               if(nactdog(3,node).eq.0) cycle
               index=ielprop(ieg(i))
               if(lakon(ieg(i))(6:7).eq.'SG') then
                  eta=prop(index+4)
                  nelem=int(prop(index+7))      
               elseif(lakon(ieg(i))(6:7).eq.'WE') then
                  eta=prop(index+4)
                  nelem=int(prop(index+7))      
               elseif(lakon(ieg(i))(6:7).eq.'DS') then
                  eta=prop(index+7)
                  nelem=int(prop(index+9))      
               endif
               if(nelem.ne.0) then
                  write(*,*) '     *INFO in flowoutput: hydraulic jump'
                  write(*,*) '           in element ',nelem,'.'
                  write(*,*) '           relative location:',eta
                  write(*,*)
               endif
            endif
         endif
      enddo
!
      close(1)
!
      return