This line https://github.com/sandialabs/MATLAB_PV_LIB/blob/684553ac1d21c59752240295ef225af7ccc1ca82/pvl_louche.m#L45 should restrict zenith to [0, 90]