Calculates the optical depth of the Martian atmosphere as a function of latitude and season.
The optical depth varies spatially and temporally, with peaks during global dust storms.
Implements Equations 1 and 2 from Appelbaum, Landis & Sherman (1991).