ADAS Subroutine rqvnew
FUNCTION RQVNEW(Z,N11,N,PHI,ZP, AMSIMP,TP,VDISP) C IMPLICIT REAL*8(A-H,O-Z) C C----------------------------------------------------------------------- C C ****************** FORTRAN77 FUNCTION: RQVNEW *********************** C C----------------------------------------------------------------------- C PURPOSE: EVALUATES ION IMPACT RATE COEFFICIENTS OF C VAINSHTEIN ET AL 1981 C C (ALTERNATIVE TO RQVAIN WITH BETTER MAXWELL AVERAGING BUT SLOWER) C C ROUTINE MUST RETURN EXCITATION RATE COEFFICIENT IF N11.LT.N AND C DEXCITATION RATE COEFFICIENT IF N11.GT.N C C NOTE THAT THE RELATION BETWEEN INVERSE PROCESSES IS DETERMINED BY THE C TEMPERATURE TP AND THE SPEED VDISP C THE TREATMENT GIVEN IS APPROXIMATE EXCEPT IN THE LIMITS VDISP=0 OR C VDISP >> DSQRT(2*TP/AMSIMP) C C INPUT C Z = TARGET ION CHARGE+1 C N11 = PRINCIPAL QUANTUM NUMBER OF INITIAL LEVEL C N = PRINCIPAL QUANTUM NUMBER OF FINAL LEVEL C PHI = (IH/EIJ)F(N ---> N'') C ZP = PROJECTILE CHARGE C AMSIMP = PROJECTILE MASS (PROTON UNITS) C TP = PROJECTILE ION TEMPERATURE (K) C VDISP = CONSTANT MEAN SPEED SHIFT FOR THE COLLISION (CM/SEC) C (DESCRIBES BEAM PLASMA SITUATIONS) C C OUTPUT C RQVNEW = RATE COEFFICIENT (CM**3 SEC-1) C C C ********* H.P. SUMMERS, JET 2 JUL 1991 ************ C C NOTES: THIS ROUTINE IS NOT YET PROPERLY ANNOTATED C C UNIX-IDL PORT: C C VERSION: 1.1 DATE: 16-1-96 C MODIFIED: TIM HAMMOND (TESSELLA SUPPORT SERVICES PLC) C - FIRST VERSION C C VERSION: 1.2 DATE: 08-02-96 C MODIFIED: TIM HAMMOND (TESSELLA SUPPORT SERVICES PLC) C - REMOVED SUPERFLUOUS VARIABLES C C VERSION: 1.3 DATE: 16-05-07 C MODIFIED: Allan Whiteford C - Updated comments as part of subroutine documentation C procedure. C----------------------------------------------------------------------- C C----------------------------------------------------------------------- C INTEGER N, N11 REAL*8 AMSIMP, PHI, TP, VDISP REAL*8 Z, ZP