LunarLambertEmperical not working properly in Photomet

[volcanopele]

ISIS version(s) affected: 3.9

Description
Using Photomet on a batch of CaSSIS Images. Photometric model used is LunarLambertEmpirical. Parameters in PVL files generated by phoempglobal . I get several PROGRAMMER ERRORS:

**PROGRAMMER ERROR** Evaluate() - Unable to evaluate the function at the point a = 0.0.
**PROGRAMMER ERROR** InsideDomain() - Unable to compute domain boundaries.
**PROGRAMMER ERROR** DomainMinimum() - Unable to calculate the domain minimum for the data set.
**PROGRAMMER ERROR** ValidateDataSet() - cspline-clamped interpolation requires a minimum of 3 data points - currently have 1.

How to reproduce

Code used, Can provide images for testing:

phoempglobal to=NIR_lunarlambert.pvl phtname=HAPKEHEN wh=0.72 hh=0.005 b0=0.755 theta=30 hg1=0.213 hg2=1 nph=40 emamax_pcoeff=1 model=LUNARLAMBERT
phoempglobal to=RED_lunarlambert.pvl phtname=HAPKEHEN wh=0.723 hh=0.011 b0=0.70 theta=30 hg1=0.213 hg2=1 nph=40 emamax_pcoeff=1 model=LUNARLAMBERT 
phoempglobal to=PAN_lunarlambert.pvl phtname=HAPKEHEN wh=0.69 hh=0.017 b0=0.55 theta=30 hg1=0.213 hg2=1 nph=40 emamax_pcoeff=1 model=LUNARLAMBERT
phoempglobal to=BLU_lunarlambert.pvl phtname=HAPKEHEN wh=0.185 hh=0.062 b0=1.00 theta=30 hg1=0.145 hg2=1 nph=40 emamax_pcoeff=1 model=LUNARLAMBERT

photomet from=\$1D1.cub to=\$1P1.cub frompvl=/data/cassis/users/perry/data/photometrytest/NIR_lunarlambert.pvl anglesource=backplane normname=albedoatm phtname=LUNARLAMBERTEMPIRICAL atmname=hapkeatm1 nulneg=no wha=0.925 hga=0.679 tau=$TAU tauref=0 hnorm=0.003 incref=0 -batchlist=nir_image_name.lis
photomet from=\$1D1.cub to=\$1P1.cub frompvl=/data/cassis/users/perry/data/photometrytest/RED_lunarlambert.pvl anglesource=backplane normname=albedoatm phtname=LUNARLAMBERTEMPIRICAL atmname=hapkeatm1 nulneg=no wha=0.922 hga=0.675 tau=$TAU tauref=0 hnorm=0.003 incref=0 -batchlist=red_image_name.lis
photomet from=\$1D1.cub to=\$1P1.cub frompvl=/data/cassis/users/perry/data/photometrytest/PAN_lunarlambert.pvl anglesource=backplane normname=albedoatm phtname=LUNARLAMBERTEMPIRICAL atmname=hapkeatm1 nulneg=no wha=0.940 hga=0.653 tau=$TAU tauref=0 hnorm=0.003 incref=0 -batchlist=pan_image_name.lis
photomet from=\$1D1.cub to=\$1P1.cub frompvl=/data/cassis/users/perry/data/photometrytest/BLU_lunarlambert.pvl anglesource=backplane normname=albedoatm phtname=LUNARLAMBERTEMPIRICAL atmname=hapkeatm1 nulneg=no wha=0.724 hga=0.746 tau=$BLUETAU tauref=0 hnorm=0.003 incref=0 -batchlist=blu_image_name.lis

Possible Solution

Additional context

[krlberry]

@volcanopele Could you send images for testing to me at [email protected]? I can't reproduce this with some arbitrary CaSSIS images I tried to test with. Also, what values are $TAU and $BLUETAU in your example above?

[krlberry]

Data received, thanks!

[krlberry]

#3614 should fix the problem described in this ticket, but I just noticed that MinnaertEmpirical has the same failure mode, and could be fixed the same way. @jlaura does it make sense to add this fix as part of this ticket, or to create a new ticket?

[jlaura]

@kberryUSGS We should open a new ticket and label as missions. Can we close this issue with #3614 being merged?