Search Site | Contact Details | FAQ

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)
© Copyright 1995-2018 The ADAS Project
Comments and questions to: adas-at-adas.ac.uk