Esope source files

   1 : C P501      SOURCE    CHAT      05/01/13    02:08:37     5004           
   2 :       SUBROUTINE P501(D,U,G,R,N,DT,G1)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 :       REAL*8 D(N),U(N),G(N),R(N),G1(N)
   6 :       DO 1 I=1,N
   7 :       R(I)=(U(I)*D(I)-DT*G(I))/(D(I)+DT*G1(I))
   8 :  1    CONTINUE
   9 :       RETURN
  10 :       END
  11 :