From fa208dbd7397ee106566ac908a3d2adc8f4a10a0 Mon Sep 17 00:00:00 2001 From: Nick Wogan Date: Thu, 29 Feb 2024 11:16:32 -0800 Subject: [PATCH] bug with radtran initializaiton --- src/radtran/clima_radtran.f90 | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/src/radtran/clima_radtran.f90 b/src/radtran/clima_radtran.f90 index 9eed586..b790784 100644 --- a/src/radtran/clima_radtran.f90 +++ b/src/radtran/clima_radtran.f90 @@ -125,6 +125,8 @@ function create_Radtran_2(species_names, particle_names, s, star_f, & character(:), allocatable, intent(out) :: err type(Radtran) :: rad + + real(dp) :: photon_scale_factor if (nz < 1) then err = '"nz" can not be less than 1.' @@ -164,9 +166,14 @@ function create_Radtran_2(species_names, particle_names, s, star_f, & allocate(rad%surface_emissivity(rad%ir%nw)) rad%surface_emissivity(:) = 1.0_dp + if (s%planet_is_present) then + photon_scale_factor = s%photon_scale_factor + else + photon_scale_factor = 1.0_dp + endif ! photons hitting the planet allocate(rad%photons_sol(rad%sol%nw)) - call read_stellar_flux(star_f, rad%sol%nw, rad%sol%wavl, s%photon_scale_factor, rad%photons_sol, err) + call read_stellar_flux(star_f, rad%sol%nw, rad%sol%wavl, photon_scale_factor, rad%photons_sol, err) if (allocated(err)) return ! IR work arrays