檢視 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