| |
| | [No title] |
 | | subroutine deriv(t,npart,x,y,z,vx,vy,vz,u,rho,h, 1 dxdt,dydt,dzdt, 2 fx,fy,fz,du,drho,dh,tdel) c****************************************************************** c * c This subroutine is the driver for all right-hand-side terms * c for all differential equations. |
 | | All processors read message file one at a time do j = 0,shmem_n_pes()-1 if (shmem_my_pe().eq.j) then c c--open message file first c open(idisk3,file='message', & form='formatted') c c--read message file (one mes. |
 | | call deriv (tfull,npart,dumx,dumy,dumz,dumvx, 1 dumvy,dumvz,dumu,dumrho,dumh, 2 f2x,f2y,f2z,f2vx, 3 f2vy,f2vz,f2u,f2rho,f2h,dt) call barrier() c--check if damage derivative not too steep, if so set damage c and cancel derivative c nset=0 do i=1,npart dtfrac1=1e20 dtfrac2=1e20 dtfrac=min(dtfrac1,dtfrac2) c if(vsound(i).lt.1e-13) write(*,*) 'vsound', vsound(i) dthydro=0.01*h(i)/vsound(i) flag(i)=.true. |
| www.ucolick.org /~randi/esnew.f (3327 words) |
|