Problems with "collision-partners-broadening": true #7
Description
I tried to run an optical depth only calculation, in the range 1230-1344cm-1, using TAPE3_aer_v_3.8.1_ex_ir as per the release_example case (~all molecules included).
This worked and reproduced an old LBLRTM calculation (with some small differences in optical depth values).
However, when I then tried to do the same calculation but with {"collision-partners-broadening": true}, I had the following problems:
- MXBRDMOL needed to be increased from value of 7 set in clblm_LineData.f90, to avoid an array bounds error (I set it to 15)
- Further array bounds errors then occur reference AVRAT in clblm_ODLAY_LineF4.f90 and clblm_ODLAY.f90 , when one of other of the widths used to define ZETA is negative, resulting in ZETA>1 and then IZ>101.
I have attempted to fix this by adding lines as follows in clblm_ODLAY.f90 (line 1886)
if (alfl .lt. 0) alfl=0.
if (alfad .lt. 0) alfad=0.
(and similar in clblm_ODLAY.f90)
The code then runs to completion without crashing/reporting errors (giving results very similar to the run without collision-partners-broadening) however probably it should be fixed in a more correct way...