diff options
Diffstat (limited to 'Dragon/src/LIBE3R.f')
| -rw-r--r-- | Dragon/src/LIBE3R.f | 22 |
1 files changed, 4 insertions, 18 deletions
diff --git a/Dragon/src/LIBE3R.f b/Dragon/src/LIBE3R.f index c644c52..045db16 100644 --- a/Dragon/src/LIBE3R.f +++ b/Dragon/src/LIBE3R.f @@ -51,7 +51,7 @@ DOUBLE PRECISION SUM PARAMETER (IOUT=6,MAXR2=12) PARAMETER (KDECAY=1,KFISSP=2,KCAPTU=3,KN2N=4,KN3N=5,KN4N=6) - CHARACTER RECNAM*80,RECNAM2*80,HSMG*131,NMDEPA(MAXR2)*6 + CHARACTER RECNAM*80,HSMG*131,NMDEPA(MAXR2)*6 *---- * ALLOCATABLE ARRAYS *---- @@ -141,6 +141,7 @@ * MAIN LOOP OVER ISOTOPES *---- NDFP2=0 + BPAX(:NBESP,:NEL+MAXR,:NEL)=0.0 DO ISO=1,NISOT ITZEA(ISO)=IZ(ISO)*10000+IA(ISO)*10 II=LEN(TRIM(NAMES(ISO))) @@ -171,22 +172,8 @@ ENDDO DEALLOCATE(LIST) ENDIF - WRITE(RECNAM,'(10HIsotopeXS/,A,8H/Energy/)') TRIM(NAMES(ISO)) - IF(hdf5_group_exists(IPAP1,TRIM(RECNAM))) THEN - IF(KPAX(NEL+KFISSP,ISO).EQ.1) THEN - WRITE(RECNAM2,'(A,16HFISS/EnergyValue)') TRIM(RECNAM) - CALL hdf5_read_data(IPAP1,TRIM(RECNAM2),VALUE) - BPAX(:,NEL+KFISSP,ISO)=VALUE - ENDIF - IF(KPAX(NEL+KCAPTU,ISO).EQ.1) THEN - WRITE(RECNAM2,'(A,18HMT-102/EnergyValue)') TRIM(RECNAM) - CALL hdf5_read_data(IPAP1,TRIM(RECNAM2),VALUE) - BPAX(:,NEL+KCAPTU,ISO)=VALUE - ENDIF - ENDIF IF(IMPX.GT.2) THEN - WRITE(IOUT,100) NAMES(ISO),BPAX(1,NEL+KDECAY,ISO), - 1 BPAX(1,NEL+KFISSP,ISO),BPAX(1,NEL+KCAPTU,ISO) + WRITE(IOUT,100) NAMES(ISO),BPAX(1,NEL+KDECAY,ISO) WRITE(IOUT,110) (NMDEPA(I),KPAX(NEL+I,ISO),I=1,MAXR) ENDIF *---- @@ -291,8 +278,7 @@ RETURN * 100 FORMAT(/44H LIBE3R: DECAY AND BURNOUT DATA FOR ISOTOPE=,A/ - 1 5X,6HDECAY=,1P,E12.5,7H E-8 /S,16H FISSION ENERGY=,E12.5,4H MEV, - 1 16H CAPTURE ENERGY=,E12.5,4H MEV) + 1 5X,6HDECAY=,1P,E12.5,7H E-8 /S) 110 FORMAT(5X,12(A6,2H= ,I1,2X)) 120 FORMAT(5X,14HPARENT NAMES: ,12A8/(19X,12A8)) END |