From 77ef0c30c3b436b1918a21bfe8eb1b760f22b8e8 Mon Sep 17 00:00:00 2001 From: Nick Wogan Date: Tue, 24 Oct 2023 18:09:07 -0700 Subject: [PATCH] prevent double counting of water continuum --- src/radtran/clima_radtran_types_create.f90 | 16 ++++++++++++++++ 1 file changed, 16 insertions(+) diff --git a/src/radtran/clima_radtran_types_create.f90 b/src/radtran/clima_radtran_types_create.f90 index e2e8633..b6b8b73 100644 --- a/src/radtran/clima_radtran_types_create.f90 +++ b/src/radtran/clima_radtran_types_create.f90 @@ -467,6 +467,22 @@ module function create_OpticalProperties(datadir, optype, species_names, & !!! Continuum !!! !!!!!!!!!!!!!!!!! if (allocated(sop%water_continuum)) then + + ! Make sure water continuum is not already accounted for with CIA + if (allocated(sop%cia)) then + do i = 1,op%ncia + if (trim(sop%cia(i)) == 'H2O-H2O') then + err = 'Optical property "water-continuum" is on, but water '// & + 'continuum is already accounted for with H2O-H2O CIA.' + return + elseif (trim(sop%cia(i)) == 'H2O-N2') then + err = 'Optical property "water-continuum" is on, but water '// & + 'continuum is already accounted for with H2O-N2 CIA.' + return + endif + enddo + endif + allocate(op%cont) filename = datadir//"/water_continuum/"//trim(sop%water_continuum)//".h5" op%cont = create_WaterContinuum(sop%water_continuum, filename, species_names, op%wavl, err)