Search Site | Contact Details | FAQ

ADAS Subroutine b8wrps

      subroutine b8wrps( open15 , iunt15 , open13 , iunt13 ,
     &                   dsninc , dsfull , dsnexp , ibsela ,
     &                   titled , date   , user   ,
     &                   ndlev  , ndtem  , ndden  , ndmet  , ndtrn ,
     &                   lnorm  ,
     &                   iz     , iz0    , iz1    ,
     &                   il     , nmet   , nord   ,
     &                   maxt   , maxd   , icntr  , icnti  , icnth ,
     &                   isa    , ila    , xja    ,
     &                   cstrga , wa     ,
     &                   icnte  ,
     &                   ie1a   , ie2a   , aa     ,
     &                   imetr  , iordr  , teva   , densa  ,
     &                   npl    , nplr   , npli   , npl3   ,
     &                   lrsel  , lisel  , lhsel  , liosel ,
     &                   lpsel  , lzsel  , lnsel  ,
     &                   wvls   , wvll   , avlt   ,
     &                   stvr   , stvi   , stvh   ,
     &                   ratpia , ratmia , stack  ,
     &                   fvionr , sgrda  ,
     &                   lsseta , lss04a
     &                 )

C-----------------------------------------------------------------------
C
C  ******************** FORTRAN77 SUBROUTINE: B8WRPS *******************
C
C  PURPOSE:  To output data to pec and sxb passing files.
C
C  CALLING PROGRAM: ADAS208
C
C
C  SUBROUTINE:
C
C  INPUT : (I*4)  IUNT15    = UNIT NUMBER FOR PECS.
C  INPUT : (L*4)  OPEN15    =.TRUE. IF SEC OUTPUT IS REQUIRED
C  INPUT : (I*4)  IUNT13    = UNIT NUMBER FOR SXBS.
C  INPUT : (L*4)  OPEN13    =.TRUE. IF SXB OUTPUT IS REQUIRED
C  INPUT : (C*44) DSNINC    = INPUT ADF04 DATA SET NAME
C  INPUT : (C*80) DSFULL    = INPUT ADF07 DATA SET NAME
C  INPUT : (C*80) DSNEXP    = INPUT EXPANSION FILE
C  INPUT : (I*4)  IBSELA(,) = IONISATION DATA BLOCK SELECTION INDICES
C                              1ST DIMENSION - (Z) ION METASTABLE COUNTER
C                              2ND DIMENSION - (Z+1) ION METASTABLE COUNTER
C  INPUT : (C*3)  TITLED    = ELEMENT SYMBOL.
C  INPUT : (C*8)  DATE      = CURRENT DATE.
C  INPUT : (C*30) USER      = FULL NAME OF AUTHOR.
C
C  INPUT : (I*4)  NDLEV     = MAXIMUM NUMBER OF LEVELS ALLOWED
C  INPUT : (I*4)  NDTEM     = MAXIMUM NUMBER OF TEMPERATURES ALLOWED
C  INPUT : (I*4)  NDDEN     = MAXIMUM NUMBER OF DENSITIES ALLOWED
C  INPUT : (I*4)  NDMET     = MAXIMUM NUMBER OF METASTABLES ALLOWED
C  INPUT : (I*4)  NDTRN     = PARAMETER = MAX. NO. OF TRANSITIONS ALLOWED
C  INPUT : (L*4)  LNORM     =.TRUE.  => IF NMET=1 THEN VARIOUS
C                                       EMISSIVITY OUTPUT FILES
C                                       NORMALISED TO STAGE TOT.POPULATN.
C                                       (** NORM TYPE = T)
C                           =.FALSE. => OTHERWISE NORMALISE TO IDENTIFIED
C                                       METASTABLE POPULATIONS.
C                                        (** NORM TYPE = M)
C  INPUT : (I*4)  IZ        =  RECOMBINED ION CHARGE READ
C  INPUT : (I*4)  IZ0       =         NUCLEAR CHARGE READ
C  INPUT : (I*4)  IZ1       = RECOMBINING ION CHARGE READ
C                             (NOTE: IZ1 SHOULD EQUAL IZ+1)
C
C  INPUT : (I*4)  IL        = NUMBER OF ENERGY LEVELS
C  INPUT : (I*4)  NMET      = NUMBER OF METASTABLES LEVELS: 1<=NMET<=NDMET
C  INPUT : (I*4)  NORD      = NUMBER OF ORDINARY LEVELS ('IL' - 'NMET')
C  INPUT : (I*4)  MAXT      = NUMBER OF INPUT TEMPERATURES ( 1 -> 'NDTEM')
C  INPUT : (I*4)  MAXD      = NUMBER OF INPUT DENSITIES ( 1 -> 'NDDEN')
C  INPUT : (I*4)  ICNTR     = NUMBER OF FREE ELECTRON RECOMBINATIONS INPUT
C  INPUT : (I*4)  ICNTI     = NUMBER OF LOWER STAGE IONISATIONS      INPUT
C  INPUT : (I*4)  ICNTH     = NO. OF CHARGE EXCHANGE RECOMBINATIONS INPUT
C  INPUT : (I*4)  ISA()     = MULTIPLICITY FOR LEVEL 'IA()'
C                             NOTE: (ISA-1)/2 = QUANTUM NUMBER (S)
C  INPUT : (I*4)  ILA()     = QUANTUM NUMBER (L) FOR LEVEL 'IA()'
C  INPUT : (R*8)  XJA()     = QUANTUM NUMBER (J-VALUE) FOR LEVEL 'IA()'
C                             NOTE: (2*XJA)+1 = STATISTICAL WEIGHT
C  INPUT : (C*18) CSTRGA()  = NOMENCLATURE/CONFIGURATION FOR LEVEL 'IA()'
C  INPUT : (R*8)  WA()      = ENERGY RELATIVE TO LEVEL 1 (CM-1)
C                             DIMENSION: LEVEL INDEX
C  INPUT : (I*4)  ICNTE     = NUMBER OF ELECTRON IMPACT TRANSITIONS INPUT
C  INPUT : (I*4)  IE1A()    = ELECTRON IMPACT TRANSITION:
C                              LOWER ENERGY LEVEL INDEX
C  INPUT : (I*4)  IE2A()    = ELECTRON IMPACT TRANSITION:
C                              UPPER ENERGY LEVEL INDEX
C  INPUT : (R*8)  AA()      = ELECTRON IMPACT TRANSITION: A-VALUE (SEC-1)
C  INPUT : (I*4)  IMETR()   = INDEX OF METASTABLE IN COMPLETE LEVEL LIST
C  INPUT : (I*4)  IORDR()   = INDEX OF ORDINARY LEVELS IN COMPLETE LEVEL
C                             LIST.
C  INPUT : (R*8)  TEVA()    = ELECTRON TEMPERATURES (UNITS: KELVIN)
C  INPUT : (R*8)  DENSA()   = ELECTRON DENSITIES  (UNITS: CM-3)
C  INPUT : (I*4)  NPL       = NO. OF METASTABLES OF(Z+1) ION ACCESSED
C                                BY EXCITED STATE IONISATION IN COPASE
C                                FILE WITH IONISATION POTENTIALS GIVEN
C                                ON THE FIRST DATA LINE
C  INPUT : (I*4)  NPLR      = NO. OF ACTIVE METASTABLES OF (Z+1) ION
C  INPUT : (I*4)  NPLI      = NO. OF ACTIVE METASTABLES OF (Z+1) ION
C  INPUT : (I*4)  NPL3      = NO. OF ACTIVE METASTABLES OF (Z+1) ION
C  INPUT   (L*4)  LRSEL     = .TRUE.  => INCLUDE FREE ELECTRON
C                                        RECOMBINATION.
C                           = .FALSE. => DO NOT INCLUDE FREE ELECTRON
C                                        RECOMBINATION.
C  INPUT : (L*4)  LISEL     = .TRUE.  => INCLUDE ELECTRON IMPACT
C                                        IONISATION.
C                           = .FALSE. => DO NOT INCLUDE FREE ELECTRON
C                                        RECOMBINATION.
C  INPUT : (L*4)  LHSEL     = .TRUE.  => INCLUDE CHARGE TRANSFER FROM
C                                        NEUTRAL HYDROGREN.
C                           = .FALSE. => DO NOT INCLUDE CHARGE TRANSFER
C                                        FROM NEUTRAL HYDROGREN.
C  INPUT : (L*4)  LIOSEL    = .TRUE.  => INCLUDE IONISATION RATES
C                           = .FALSE. => DO NOT INCLUDE IONISATION RATES
C                              FOR RECOM AND  3-BODY
C  INPUT : (L*4)  LPSEL     = .TRUE.  => INCLUDE PROTON COLLISIONS
C                           = .FALSE. =>DO NOT INCLUDE PROTON COLLISIONS
C  INPUT : (L*4)  LZSEL     = .TRUE.  => SCALE PROTON COLLISIONS WITH
C                                        PLASMA Z EFFECTIVE'ZEFF'.
C                           = .FALSE. => DO NOT SCALE PROTON COLLISIONS
C                                         WITH PLASMA Z EFFECTIVE 'ZEFF'.
C                             (ONLY USED IF 'LPSEL=.TRUE.')
C  INPUT : (L*4)  LNSEL     = .TRUE.  => INCLUDE PROJECTED BUNDLE-N DATA
C                                         FROM DATAFILE IF AVAILABLE
C                           = .FALSE. => DO NOT INCLUDE PROJECTED
C                                         BUNDLE-N DATA
C  INPUT : (R*8)  WVLS      = SHORT WAVELENGTH LIMIT FOR PEC & SXB (A)
C  INPUT : (R*8)  WVLL      = LONG WAVELENGTH LIMIT FOR PEC & SXB  (A)
C  INPUT : (R*8)  AVLT      = LOWER LIMIT OF A-VALUES FOR PEC & SXB
C
C  INPUT : (R*4) STVR(,,,)  = FREE ELECTRON RECOMBINATION COEFFICIENTS
C                              1st DIMENSION: ORDINARY LEVEL INDEX
C                              2nd DIMENSION: TEMPERATURE INDEX
C                              3rd DIMENSION: DENSITY INDEX
C                              4TH DIMENSION: PARENT INDEX
C  INPUT : (R*4) STVI(,,,)  = ELECTRON IMPACT IONISATION  COEFFICIENTS
C                              1st DIMENSION: ORDINARY LEVEL INDEX
C                              2nd DIMENSION: TEMPERATURE INDEX
C                              3rd DIMENSION: DENSITY INDEX
C                              4TH DIMENSION: PARENT INDEX
C  INPUT : (R*4) STVH(,,,)  =  CHARGE EXCHANGE COEFFICIENTS
C                              1st DIMENSION: ORDINARY LEVEL INDEX
C                              2nd DIMENSION: TEMPERATURE INDEX
C                              3rd DIMENSION: DENSITY INDEX
C  INPUT : (R*8) RATPIA(,)  = RATIO ( N(Z+1)/N(Z)  STAGE ABUNDANCIES )
C                              1ST DIMENSION: TEMP/DENS INDEX
C                              2ND DIMENSION: PARENT INDEX
C  INPUT : (R*8) RATMIA(,)  = RATIO ( N(Z-1)/N(Z)  STAGE ABUNDANCIES )
C                              1ST DIMENSION: TEMP/DENS INDEX
C                              2ND DIMENSION: PARENT INDEX
C  INPUT : (R*4) STACK(,,,) = POPULATION DEPENDENCE
C                              1st DIMENSION: ORDINARY LEVEL INDEX
C                              2nd DIMENSION: METASTABLE INDEX
C                              3rd DIMENSION: TEMPERATURE INDEX
C                              4th DIMENSION: DENSITY INDEX
C  INPUT : (R*8) FVIONR(,,,)= GEN. COLL. RAD. IONIS. RATE COEFFTS.
C                             1ST DIMENSION: (Z) ION METASTABLE INDEX
C                             2ND DIMENSION: (Z+1) ION METASTABLE INDEX
C                             3rd DIMENSION: TEMPERATURE INDEX
C                             4th DIMENSION: DENSITY INDEX
C  INPUT : (R*8)  SGRDA(,,) = GROUND & METASTABLE IONISATION RATE
C                             COEFFICIENTS  FROM SZD FILES (CM3 SEC-1)
C                             1ST DIMENSION: TEMPERATURE INDEX
C                             2ND DIMENSION: (Z) ION METASTABLE INDEX
C                             3RD DIMENSION: (Z+1) ION METASTABLE INDEX
C  INPUT : (L*4) LSSETA(,)  = .TRUE.  - MET. IONIS RATE SET IN B8GETS
C                             .FALSE.- MET. IONIS RATE NOT SET IN B8GETS
C                              1ST DIMENSION: (Z) ION METASTABLE INDEX
C                              2ND DIMENSION: (Z+1) ION METASTABLE INDEX
C  INPUT : (L*4) LSS04A(,)  = .TRUE. => IONIS. RATE SET IN ADF04 FILE:
C                             .FALSE.=> NOT SET IN ADF04 FILE
C                               1ST DIM: LEVEL INDEX
C                               2ND DIM: PARENT METASTABLE INDEX
C
C          (I*4) NOTRN      = PARAMETER = MAXIMUM NUMBER OF TRANSITIONS
C          (I*4) NDPEC      = PARAMETER = MAXIMUM NUMBER OF PECS PER
C                                         METASTABLE FOR OUTPUT
C          (I*4) METCNT     = COUNTER OF PECS FOR EACH METASTABLE
C
C          (I*4) I4UNIT     = FUNCTION (SEE ROUTINE SELECTION BELOW)
C
C          (I*4) I          = GENERAL USE
C          (I*4) IP         = GENERAL USE
C          (I*4) J          = GENERAL USE
C          (I*4) K          = GENERAL USE
C          (I*4) L          = GENERAL USE
C
C          (R*8) DUM1       = GENERAL USE- DUMMY
C          (R*8) DUM2       = GENERAL USE- DUMMY
C          (R*8) DUM3       = GENERAL USE- DUMMY
C
C ROUTINES:
C          -------------------------------------------------------------
C          i4unit     ADAS      Fetch unit number for output of messages
C          b8norm     ADAS      Perform  stage population normalisation
C          b8corp     ADAS      'fixes' low te problem in rec. data of pecs
C          b8winf     ADAS      Determines ionis. source and writes
C                               comment block
C          xxeiam     ADAS      Get 2-character name of element.
C          xxslen     ADAS      Returns lenght of string.
C          xxordr     ADAS      Sorts a real*8 array.
C
C-----------------------------------------------------------------------
C
C NOTES: Based on b8wr11, b8wr12 and hapecf.
C
C
C VERSION  : 1.1
C DATE     : 15-02-2006
C MODIFIED : Martin O'Mullane
C              - First version.
C
C VERSION  : 1.2
C DATE     : 08-03-2006
C MODIFIED : Martin O'Mullane
C              - Make sure that CX pec data which is 0.0 is written 
C                as 1.0E-70 to avoid splining problems when reading
C                adf15 datasets.
C
C-----------------------------------------------------------------------
      CHARACTER*18        CSTRGA(NDLEV)
      CHARACTER*8         DATE
      CHARACTER*80        DSFULL,      DSNEXP,      DSNINC
      CHARACTER*3         TITLED
      CHARACTER*30        USER
      INTEGER             IBSELA(NDMET,NDMET),      ICNTE,       ICNTH
      INTEGER             ICNTI,       ICNTR,       IE1A(NDTRN)
      INTEGER             IE2A(NDTRN), IL,          ILA(NDLEV)
      INTEGER             IMETR(NDMET),             IORDR(NDLEV)
      INTEGER             ISA(NDLEV),  IUNT13,      IUNT15,      IZ
      INTEGER             IZ0,         IZ1,         MAXD,        MAXT
      INTEGER             NDDEN,       NDLEV,       NDMET,       NDTEM
      INTEGER             NDTRN,       NMET,        NORD,        NPL
      INTEGER             NPL3,        NPLI,        NPLR
      LOGICAL             LHSEL,       LIOSEL,      LISEL,       LNORM
      LOGICAL             LNSEL,       LPSEL,       LRSEL
      LOGICAL             LSS04A(NDLEV,NDMET),      LSSETA(NDMET,NDMET)
      LOGICAL             LZSEL,       OPEN13,      OPEN15
      REAL*8              AA(NDTRN),   AVLT,        DENSA(NDDEN)
      REAL*8              FVIONR(NDMET,NDMET,NDTEM,NDDEN)
      REAL*8              RATMIA(NDDEN,NDMET),      RATPIA(NDDEN,NDMET)
      REAL*8              SGRDA(NDTEM,NDMET,NDMET)
      REAL                STACK(NDLEV,NDMET,NDTEM,NDDEN)
      REAL                STVH(NDLEV,NDTEM,NDDEN,NDMET)
      REAL                STVI(NDLEV,NDTEM,NDDEN,NDMET)
      REAL                STVR(NDLEV,NDTEM,NDDEN,NDMET)
      REAL*8              TEVA(NDTEM), WA(NDLEV),   WVLL,        WVLS
      REAL*8              XJA(NDLEV)
© Copyright 1995-2018 The ADAS Project
Comments and questions to: adas-at-adas.ac.uk