From 1b127b999bb6af837584075f8b3a97cf730cc583 Mon Sep 17 00:00:00 2001 From: Emily Sandford Date: Mon, 11 May 2026 15:08:37 +0200 Subject: [PATCH] fixed incorrect metal mass fraction array for Freedman lowT opacity tables --- kap/preprocessor/src/freedman.f90 | 14 +++++++------- kap/public/kap_def.f90 | 4 ++-- 2 files changed, 9 insertions(+), 9 deletions(-) diff --git a/kap/preprocessor/src/freedman.f90 b/kap/preprocessor/src/freedman.f90 index 338e67a73..82b2b9d79 100644 --- a/kap/preprocessor/src/freedman.f90 +++ b/kap/preprocessor/src/freedman.f90 @@ -43,19 +43,19 @@ subroutine get_Freedman_fname(data_dir, Z, fname) integer :: iz iz = floor(Z*1d5 + 0.1d0) select case (iz) - case (1000) + case (700) fname = trim(data_dir)//'/m0.3.txt' - case (2000) + case (1500) fname = trim(data_dir)//'/p0.0.txt' - case (4000) + case (2900) fname = trim(data_dir)//'/p0.3.txt' - case (10000) + case (7000) fname = trim(data_dir)//'/p0.7.txt' - case (20000) + case (13000) fname = trim(data_dir)//'/p1.0.txt' - case (63000) + case (32100) fname = trim(data_dir)//'/p1.5.txt' - case (100000) + case (42900) fname = trim(data_dir)//'/p1.7.txt' case default write (*, *) 'get_Freedman_fname: unexpected Z value for Freedman data', Z diff --git a/kap/public/kap_def.f90 b/kap/public/kap_def.f90 index 3eaeb69a9..66ae19f19 100644 --- a/kap/public/kap_def.f90 +++ b/kap/public/kap_def.f90 @@ -618,7 +618,7 @@ subroutine kap_def_init(kap_cache_dir_in) kap_lowT_Xs(1:num_kap_lowT_Xs(i), i) = [ 0.00d0 ] num_kap_lowT_Zs(i) = 7 kap_lowT_Zs(1:num_kap_lowT_Zs(i), i) = & - [ 0.01d0, 0.02d0, 0.04d0, 0.100d0, 0.200d0, 0.63d0, 1.00d0 ] + [ 0.007d0, 0.015d0, 0.029d0, 0.07d0, 0.13d0, 0.321d0, 0.429d0 ] num_kap_lowT_Xs_for_this_Z(1:num_kap_lowT_Zs(i), i) = num_kap_lowT_Xs(i) case DEFAULT num_kap_lowT_Xs(i) = 10 @@ -831,4 +831,4 @@ end subroutine free_co_table end subroutine do_Free_Kap_Tables -end module kap_def +end module kap_def \ No newline at end of file