ADAS Subroutine a5data
SUBROUTINE A5DATA( DSFULL , INDXREF , TITLE , CAMETH , Z0 , Z , & Z1 , NIGRP , EMIN , CIA , NSHELA, & EIONA , IZETAA , NRGRP , CRA , NRESOA, & ENERA , WGHTA , ICT , XA , YA ,APA, & ITOUT , TOA , YOA , YOAP , & ISTDIM , IREAD , NA , LA & ) C----------------------------------------------------------------------- C C ********************* FORTRAN77 SUBROUTINE A5DATA ******************* C C PURPOSE: TO REFRESH A DATA INDEX FROM AN ADAS105 ARCHIVE. READS C IN THE INDEX CODE A-ADAS, B-BURGESS AND THE THE REST OF C THE DATA AS APPROPRIATE. C C CALLING PROGRAM: C ADAS105.FOR C C INPUT: C (C*80) DSFULL - THE USERS' CHOSEN ARCHIVE FILE NAME. C (I*4) INDXREF - THE INDEX NUMBER TO REFRESH FROM. C (C*4) CAMETH - THE TAG TO DISTINGUISH BETWEEN THE C TWO TYPES OF ANALYSIS. C A - ADAS, B- BURGESS C (I*4) ISTDIM = ARRAY DIMENSIONS : MAX. NO OF VALUES THAT C CAN BE READ IN C (I*4) IREAD = THE INPUT UNIT C C OUTPUTS: C (C*40) TITLE - THE INFORMATION LINE IN THE ARCHIVE C FILE. C (R*8) Z0 = NUCLEAR CHARGE OF ION C (R*8) Z = INITIAL ION CHARGE C (R*8) Z1 = FINAL ION CHARGE C (I*4) NIGRP = NO. OF SHELL GROUPS C (R*8) EMIN = MINIMUM ENERGY (?) C (R*8) CIA() = SCALING PARAMETERS FOR SHELL GROUPS C (I*4) NSHELA()=NO. OF ENTRIES FOR EACH SHELL GROUP C (I*4) NA(,) = SHELL GROUP DATA : N C (I*4) LA(,) = SHELL GROUP DATA : L C (R*8) EIONA(,)=SHELL GROUP DATA : EION(RYD) C (I*4) IZETAA(,)=SHELL GROUP DATA : IZETA C (I*4) NRGRP = NO. OF RESONANCE GROUPS C (R*8) CRA() = SCALING PARAMETERS FOR RESONANCE GROUPS C (I*4) NRESOA()=NO. OF ENTRIES FOR EACH RESONANCE GROUP C (R*8) ENERA(,)=RESONANCE GROUP DATA : ENERGY(RYD) C (R*8) WGHTA(,)=RESONANCE GROUP DATA : WEIGHT C (R*8) ICT = NO. OF ENERGY / X-SECTION PAIRS C (R*8) XA() = X, THRESHOLD PARAMETER RELATIVE TO FIRST IONIS. POT. C (R*8) YA() = Q/(1-1/X), Q=CROSS-SECTION ? C (R*8) APA() = QEM/(1-1/X), QEM=APPROX. X-SECTION ) C (R*4) YPA() = Q/QEM C (I*4) ITOUT = NO. OF TEMPS. C (R*8) TOA() = TEMP (KELVIN) C (R*8) YOA() = S, MAXWELL AVERAGED IONISATION RATE COEFF.(CM^3 S^-1) C (R*8) YOAP() = SEM, APPROXIMATE RATE COEFF. C C ROUTINES: NONE C C AUTHOR: WILLIAM OSBORN, TESSELLA SUPPORT SERVICES PLC., 28TH AUG 1996 C C VERSION 1.1 DATE: 28-08-96 C MODIFIED: WILLIAM OSBORN C - FIRST RELEASE C C VERSION 1.2 DATE: 08-10-96 C MODIFIED: WILLIAM OSBORN C - REMOVED FILE READ FOR NRGRP=0 CASE C - ADDED READ OF NA AND LA FROM ARCHIVE C----------------------------------------------------------------------- CHARACTER*4 CAMETH CHARACTER*80 DSFULL CHARACTER*40 TITLE INTEGER ICT, INDXREF, IREAD, ISTDIM INTEGER ITOUT, IZETAA(6,2), LA(6,2), NA(6,2) INTEGER NIGRP, NRESOA(2), NRGRP INTEGER NSHELA(2) REAL*8 APA(ISTDIM), CIA(2), CRA(2) REAL*8 EIONA(6,2), EMIN, ENERA(6,2) REAL*8 TOA(ISTDIM), WGHTA(6,2), XA(ISTDIM) REAL*8 YA(ISTDIM), YOA(ISTDIM), YOAP(ISTDIM) REAL*8 Z, Z0, Z1