ADAS Subroutine d6mfll
SUBROUTINE D6MFLL ( NTDIM , IZDIMD , IPDIMD , IMDIMD , & NMDIM , & NSTAGE , NPART , & ACDA , SCDA , CCDA , QCDA , XCDA , & DENS , DENSH , & ITEM , & A & ) C----------------------------------------------------------------------- C C ****************** FORTRAN 77 SUBROUTINE: D6MFLL ******************** C C PURPOSE: FILLS MATRIX WITH RECOMBINATION, IONISATION AND METASTABLE C CROSS-COUPLING COEFFICIENTS READY FOR EIGENVECTOR SOLUTION C C CALLING PROGRAM: D6MPOP C C SUBROUTINE: C C C INPUT : (I*4) NTDIM = MAXIMUM NUMBER OF DTEV/DDENS PAIRS C INPUT : (I*4) IZDIMD = MAXIMUM NUMBER OF CHARGE STATES C IN ISONUCLEAR MASTER FILES C INPUT : (I*4) IPDIMD = MAXIMUM NUMBER OF METASTABLES FOR EACH C IONISATION STAGE C INPUT : (I*4) IMDIMD = MAXIMUM NUMBER OF POPULATIONS C C INPUT : (I*4) NMDIM = MAX. NUMBER OF POPULATIONS (FROM D6MPOP) C C INPUT : (I*4) NSTAGE = NUMBER OF IONISATION STATES (EXCL. C EXTRA ONE (BARE NUCLEUS) ADDED AT END C INPUT : (I*4) NPART() = METASTABLE PARTITION. I.E. NUMBER OF C METASTABLES FROM CHARGE STATE IZ1MIN-1 TO C IZ1MAX ON INPUT C C INPUT : (R*8) ACDA(,,,) = INTERPOLATION OF ACD COEFFICIENT (CM3 S-1) C 1ST DIM: TEMPERATURE INDEX C 2ND DIM: CHARGE STATE INDEX C 3RD DIM: RECOMBINING METASTABLE INDEX C 4TH DIM: RECOMBINED METASTABLE INDEX C INPUT : (R*8) SCDA(,,,) = INTERPOLATION OF SCD COEFFICIENT (CM3 S-1) C 1ST DIM: TEMPERATURE INDEX C 2ND DIM: CHARGE STATE INDEX C 3RD DIM: RECOMBINING METASTABLE INDEX C 4TH DIM: RECOMBINED METASTABLE INDEX C INPUT : (R*8) CCDA(,,,) = INTERPOLATION OF CCD COEFFICIENT (CM3 S-1) C 1ST DIM: TEMPERATURE INDEX C 2ND DIM: CHARGE STATE INDEX C 3RD DIM: RECOMBINING METASTABLE INDEX C 4TH DIM: RECOMBINED METASTABLE INDEX C INPUT : (R*8) QCDA(,,,) = INTERPOLATION OF QCD COEFFICIENT (CM3 S-1) C 1ST DIM: TEMPERATURE INDEX C 2ND DIM: CHARGE STATE INDEX C 3RD DIM: FIRST METASTABLE INDEX C 4TH DIM: SECOND METASTABLE INDEX C INPUT : (R*8) XCDA(,,,) = INTERPOLATION OF XCD COEFFICIENT (CM3 S-1) C 1ST DIM: TEMPERATURE INDEX C 2ND DIM: CHARGE STATE INDEX C 3RD DIM: FIRST PARENT METASTABLE INDEX C 4TH DIM: SECOND PARENT METASTABLE INDEX C C INPUT : (R*8) DENS() = ELECTRON DENSITIES FOR MODEL C INPUT : (R*8) DENSH() = NEUTRAL HYDROGEN DENSITIES FOR MODEL C C INPUT : (I*4) ITEM = CURRENT TEMP/DENSITY INDEX C C OUTPUT : (R*8) A(,) = RECOMB/IONIS COLL. RAD. MATRIX C C (I*4) I = GENERAL INDEX C (I*4) IGRD = GENERAL INDEX C (I*4) IND = GENERAL INDEX C (I*4) IND1 = GENERAL INDEX C (I*4) IPRT = GENERAL INDEX C (I*4) IZ = IONISATION STAGE COUNTER C (I*4) I = GENERAL INDEX C (I*4) JGRD = GENERAL INDEX C C C ROUTINES: C ROUTINE SOURCE BRIEF DESCRIPTION C ------------------------------------------------------------ C C C AUTHOR : H. P. SUMMERS, JET C K1/1/57 C JET EXT. 4941 C C DATE : 27/07/94 C UNIX-IDL PORT: C C AUTHOR: WILLIAM OSBORN (TESSELLA SUPPORT SERVICES PLC) C C DATE: 07/06/96 C C VERSION: 1.1 DATE:07/06/96 C MODIFIED: WILLIAM OSBORN C - FIRST VERSION C C VERSION: 1.2 DATE:28/03/07 C MODIFIED: ALLAN WHITEFORD C - REMOVED VERSIONING INFORMATION FROM COLUMNS 73-80 C C----------------------------------------------------------------------- INTEGER IMDIMD, IPDIMD, ITEM, IZDIMD INTEGER NMDIM, NPART(IZDIMD), NSTAGE INTEGER NTDIM REAL*8 A(NMDIM,NMDIM) REAL*8 ACDA(NTDIM,IZDIMD,IPDIMD,IPDIMD) REAL*8 CCDA(NTDIM,IZDIMD,IPDIMD,IPDIMD) REAL*8 DENS(NTDIM), DENSH(NTDIM) REAL*8 QCDA(NTDIM,IZDIMD,IPDIMD,IPDIMD) REAL*8 SCDA(NTDIM,IZDIMD,IPDIMD,IPDIMD) REAL*8 XCDA(NTDIM,IZDIMD,IPDIMD,IPDIMD)