Re: Vuhard

Hi Wendell, 

 Are you brazilian?
 

 please, try this and let me know if it'll be okay:
 

 

 subroutine vuhard(
C Read only -
* nblock,
* jElem, kIntPt, kLayer, kSecPt,
* lAnneal, stepTime, totalTime, dt, cmname,
* nstatev, nfieldv, nprops,
* props, tempOld, tempNew, fieldOld, fieldNew,
* stateOld,
* eqps, eqpsRate,
C Write only -
* yield, dyieldDtemp, dyieldDeqps,
* stateNew )

include 'vaba_param.inc'

dimension props(nprops), tempOld(nblock), tempNew(nblock),
1 fieldOld(nblock,nfieldv), fieldNew(nblock,nfieldv),
2 stateOld(nblock,nstatev), eqps(nblock), eqpsRate(nblock),
3 yield(nblock), dyieldDtemp(nblock), dyieldDeqps(nblock,2),
4 stateNew(nblock,nstatev), jElem(nblock)
C

 C
 real*8 :: val_K
real*8 :: val_n
real*8 :: val_m
real*8 :: Syield
real*8 :: Eel

 C
 C
 val_K = props(1)
val_n = props(2)
val_m = props(3)
Syield = props(4)
Eel = props(5)
 
character*80 cmname
C
 

 parameter ( one = 1.0d0)


 do km = 1,nblock


 yield(km)= Syield + val_K*((Eel+eqps(km))**val_n)*(eqpsRate(km)**val_m)


 dyieldDeqps(km,1)=val_K*val_n*((Eel + eqps(km))**(val_n - one)) * (eqpsRate(km)**val_m)


 dyieldDeqps(km,2)=val_K*val_m*((Eel+eqps(km))**val_n)*(eqpsRate(km)**(val_m - one))


 end do
C
return
end