diff options
| author | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
|---|---|---|
| committer | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
| commit | 7dfcc480ba1e19bd3232349fc733caef94034292 (patch) | |
| tree | 03ee104eb8846d5cc1a981d267687a729185d3f3 /Dragon/src/LIBD10.F | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Dragon/src/LIBD10.F')
| -rw-r--r-- | Dragon/src/LIBD10.F | 62 |
1 files changed, 62 insertions, 0 deletions
diff --git a/Dragon/src/LIBD10.F b/Dragon/src/LIBD10.F new file mode 100644 index 0000000..c98de35 --- /dev/null +++ b/Dragon/src/LIBD10.F @@ -0,0 +1,62 @@ +*DECK LIBD10 +#if defined(HDF5_LIB) + SUBROUTINE LIBD10 (MAXDIL,NAMFIL,HNISOR,NDIL,DILUT) +* +*----------------------------------------------------------------------- +* +*Purpose: +* Find the dilutions corresponding to a resonant isotope within a +* library in Apolib-3 format. +* +*Copyright: +* Copyright (C) 2022 Ecole Polytechnique de Montreal +* This library is free software; you can redistribute it and/or +* modify it under the terms of the GNU Lesser General Public +* License as published by the Free Software Foundation; either +* version 2.1 of the License, or (at your option) any later version +* +*Author(s): A. Hebert +* +*Parameters: input +* MAXDIL maximum number of dilutions. +* NAMFIL name of the APOLIB-3 file. +* HNISOR library name of the resonant isotope. +* +*Parameters: output +* NDIL number of finite dilutions. +* DILUT dilutions. +* +*----------------------------------------------------------------------- +* + USE GANLIB + USE hdf5_wrap +*---- +* SUBROUTINE ARGUMENTS +*---- + INTEGER MAXDIL,NDIL + CHARACTER HNISOR*12,NAMFIL*(*) + REAL DILUT(MAXDIL) +*---- +* LOCAL VARIABLES +*---- + TYPE(C_PTR) IPAP + CHARACTER RECNAM*80 + REAL, ALLOCATABLE, DIMENSION(:) :: BGXS +*---- +* OPEN AND PROBE THE APOLIB-3 FILE. +*---- + CALL hdf5_open_file(NAMFIL, IPAP, .TRUE.) + WRITE(RECNAM,'(9HIsotopes/,A,15H/HomoRates/BgXS)') TRIM(HNISOR) + CALL hdf5_read_data(IPAP, RECNAM, BGXS) + NDIL=SIZE(BGXS) + IF(NDIL.GT.MAXDIL) CALL XABORT('LIBD10: MAXDIL IS TOO SMALL.') + DILUT(:NDIL)=BGXS(:) + NDIL=NDIL-1 + DEALLOCATE(BGXS) + CALL hdf5_close_file(IPAP) +#else + SUBROUTINE LIBD10 + CALL XABORT('LIBD10: THE HDF5 API IS NOT AVAILABLE.') +#endif /* defined(HDF5_LIB) */ + RETURN + END |