;Function for including geometry effects in a phase function. This ;particular one puts in the secant of the view angle. function geometry_factor,ssa,cva,sza,_extra=extra return,1.0/cos(cva*!const.dtor) end