volpfv.F File Reference

#include "implicit_f.inc"
#include "com06_c.inc"
#include "com08_c.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	volpfv (ivolu, rvolu, vol, fsav, npc, tf, pmain, wfext, python, nfunct)

Function/Subroutine Documentation

volpfv()

subroutine volpfv	(	integer, dimension(*)	ivolu,
			rvolu,
			vol,
			fsav,
		integer, dimension(*)	npc,
			tf,
		integer	pmain,
		double precision, intent(inout)	wfext,
		type(python_), intent(inout)	python,
		integer, intent(in)	nfunct )

Definition at line 31 of file volpfv.F.

         USE python_funct_mod
         USE finter_mixed_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com06_c.inc"
#include      "com08_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NPC(*),IVOLU(*), PMAIN
      my_real tf(*),rvolu(*),fsav(*)
      DOUBLE PRECISION,INTENT(INOUT) :: WFEXT
      TYPE(PYTHON_), intent(inout) :: PYTHON
      INTEGER, intent(in) :: NFUNCT
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER ITFUN,ITYPFUN
      my_real vol,v0,area,veps,vinc,pold,pext,vold,pres,q,qold,dv,scale,xfun
C-----------------------------------------------
C
C CALCUL DE LA PRESSION
C
      itfun  =ivolu(7)
      itypfun=ivolu(19)
      pext   =rvolu(3)
      v0     =rvolu(4)
      vinc   =rvolu(5)
      pold   =rvolu(12)
      vold   =rvolu(16)
      veps   =rvolu(17)
      area   =rvolu(18)
      qold   =rvolu(23)
      scale  =rvolu(15)
      vol = vol + veps
      v0  =  v0
      SELECT CASE (itypfun)
      CASE(0)
         xfun = (v0-vinc)/(vol-vinc)
      CASE(1, 3)
         xfun = tt
      CASE(2)
         xfun = (vol-vinc) / (v0-vinc)
      CASE DEFAULT
         
      END SELECT
      pres = scale * finter_mixed(python,nfunct,itfun,xfun,npc,tf)
      IF (itypfun == 3) THEN
         pres = pres * (v0-vinc) / (vol-vinc)
      ENDIF
      q = zero
      IF (ispmd+1==pmain) THEN
        dv=vol-vold
        wfext=wfext+(half*(q+qold+pres+pold))*dv
        fsav(1)=zero
        fsav(2)=vol
        fsav(3)=pres
        fsav(4)=area
        fsav(5)=zero
        fsav(6)=zero
        fsav(7)=zero
        fsav(8)=zero
        fsav(9)=zero
        fsav(10)=zero
        fsav(11)=zero
        fsav(12)=zero
      ENDIF
C
      rvolu(16)=vol
      rvolu(12)=pres
      rvolu(23)=q
C
      RETURN