Line data    Source code

       1             :       MODULE m_stpot1
       2             :       CONTAINS
       3          84 :       SUBROUTINE stpot1(
       4          84 :      >                  msh,n,z,rad,
       5          84 :      <                  vr1)
       6             : c     **********************************************************
       7             : c     create a starting potential
       8             : c     **********************************************************
       9             : 
      10             :       IMPLICIT NONE
      11             : C     ..
      12             :       INTEGER,INTENT (IN) :: n,msh
      13             :       REAL,INTENT (IN)    :: z,rad(msh)
      14             :       REAL,INTENT (OUT)   :: vr1(msh)
      15             : C     ..
      16             : C     .. Local Scalars ..
      17             :       INTEGER i,i1
      18             :       REAL    d1,h1
      19             : C     ..
      20          84 :       IF (z.LT.1.0) THEN
      21           0 :         DO i = 1,n
      22           0 :            vr1(i) = 1.e-3
      23             :         ENDDO
      24             :         RETURN
      25             :       ENDIF
      26             : C
      27          84 :       d1 = 0.9e0
      28          84 :       h1 = d1* (z-1)**0.4e0
      29          84 :       i1 = 0.75e0*n
      30       72758 :       DO i = 1,n
      31       72758 :          vr1(i) = -1.0
      32             :       ENDDO
      33       54564 :       DO i = 1,i1
      34       54564 :          vr1(i) = vr1(i)* ((z-1)/ (h1* (exp(rad(i)/d1)-1)+1)+1)
      35             :       ENDDO
      36             :       RETURN
      37             :       END SUBROUTINE
      38             :       END