Subroutine updatetaunp

Subroutine updatetaunp is activated by the command $usedynataunrfunc. Below is an exmaple. User must modify this code.

subroutine updatetaunp
use nrtype
use global2d , only :me, totalnode,sprst               
implicit none
integer(I4B)  :: i,j,k,i1,i2 ,ksk
real(dp)   :: x,y,t
do i = 1 , totalnode
     t=sprst(i)
     x=me%x(i)
     y=me%y(i)
     me%taun(i)= 1.0e-30*me%n(i)
     me%taup(i)= 1.0e-30*me%p(i)
end do
print*,"you called the system's update-taunp function, but you should use 2D-ddcc-dyna.exe to run the program. Program stop!"
stop
end subroutine