subroutine vueos(
C Read only -
     *     nblock, 
     *     jElem, kIntPt, kLayer, kSecPt, 
     *     stepTime, totalTime, dt, cmname,
     *     nstatev, nfieldv,
     *     nprops, props, 
     *     tempOld, tempNew, fieldOld, fieldNew,
     *     stateOld,
     *     charLength, coordMp,
     *     densityMean, 
     *     refDensity, 
     *     density, dkk,
     *     Em,
C Write only -
     *     press, dpdrho,dpdem,
     *     stateNew )
C
      include 'vaba_param.inc'
C
      dimension props(nprops), 
     *  tempOld(nblock),
     *  fieldOld(nblock,nfieldv), 
     *  stateOld(nblock,nstatev), 
     *  tempNew(nblock),
     *  fieldNew(nblock,nfieldv),
     *  press(nblock), density(nblock),
     *  refDensity(nblock),charLength(nblock),
     *  coordMp(nblock,*),dkk(nblock), 
     *  dpdrho(nblock), densityMean(nblock),
     *  dpdem(nblock),
     *  stateNew(nblock,nstatev)
C
      character*80 cmname
C
	    parameter ( zero=0.d0, one = 1.d0, half=0.5d0,three=3.0d0 )
C
      rK    = props(1) ! K

      do k = 1, nblock
         evol = log (refDensity(k)/density(k))   
         densityInv = one / density(k)

*
         press(k)= -rK*evol
         dpdem(k) = zero

         dpdrho(k)=  rK*densityInv


      end do
C
      return
      end