SUBROUTINE QPOTCG ( I )
Argument Definitions (+ indicates altered content)
INTEGER            I
Description
  Copyright (C) 1978 Charlotte Froese Fisher 
  Copyright (C) 2012 Jacek Kobus 
 
 
The qpotwb subroutine fills for orbital i the vrel array with the 
values of the quasirelativistic correction due to R.D.Cowan and 
D.C.Griffin (J.Opt.Soc.AM. 66, 1010 (1976))
Source file:qpotcg.f
I/O Operations:
Unit ID  Unit No       Access  Form   Operation
      *                   SEQ  FMTD           W  
    OUC                   SEQ  FMTD           W  

Operation codes A=rewind,B=backspace,C=close,E=endfile
                I=inquire,O=open,R=read,W=write
Intrinsic Functions Called
GENERIC*8          ABS,        SQRT
External Functions and Subroutines Called
REAL*8             BB
SUBROUTINE         DIFFS
Parameter Variables Used
INTEGER            MAXNO              (MAXNO = 29999)
INTEGER            MAXORB             (MAXORB = 20)
REAL*8             DERCUTTOFF         (DERCUTTOFF = 1.D12)
Local Variables (+ indicates altered content)
INTEGER           +J
REAL*8            +A(MAXNO),  +A1,        +A2,        +AR(MAXNO)
REAL*8            +ARR(MAXNO),            +DLI,       +EE
REAL*8            +P1,        +P2,        +PP(MAXNO), +RT
REAL*8            +RT2,        V,         +VV,        +XK
REAL*8             XY1
Referenced Equivalenced Variables (+ indicates altered content)

EQUIV              REAL*8            +A()        Local Var
                   REAL*8            +ARR()      Local Var
Referenced Common Block Variables (+ indicates altered content)
CONST              REAL*8             CVEL,      PRECIS
DEBUG              INTEGER            IDBG(100)
INOUT              INTEGER            OUC
MCOMMBLOCK         INTEGER            L(MAXORB)
MCOMMBLOCK         REAL*8             P(MAXORB,MAXNO)
MCOMMBLOCK         REAL*8             R(MAXNO),  R2(MAXNO)
PARAM              INTEGER            NO
PARAM              REAL*8             Z
PARAM0             INTEGER            IONE
PARAM1             REAL*8             D1,        D2
PARAM1             REAL*8             D4,        D5
REL                INTEGER            IQPOT,     IUPD
REL                REAL*8             CVELDAMP, +DLR(MAXORB)
REL                REAL*8             FSC2,     +Q0(MAXORB)
REL                REAL*8            +QM1(MAXORB)
REL                REAL*8            +QM2(MAXORB)
REL                REAL*8            +VREL(MAXNO)
REL                REAL*8             VXCH(MAXNO)
REL                REAL*8             Z1,       +ZRL(MAXORB)
STATE              REAL*8             QC(MAXORB,40)
WAVE               REAL*8             E(MAXORB,MAXORB)