_RL ssq0, ssq1, ssq2 ! constant used in surface specific humidity
DATA ssq0, ssq1, ssq2
ssq = ssq0*exp( lath*(ssq1-ssq2/Tsf) ) / p0
_RL ssq0, ssq1, ssq2
& ssq0, ssq1, ssq2,
ssqt=ssq0*EXP( lath*(ssq1-ssq2/tta) ) / p0
ssqa = ssq0*EXP( lath*(ssq1-ssq2/local)) / p0
ssqa = ssq0*EXP( lath*(ssq1-ssq2/local) ) / p0
ssq= ssq0*exp( lath*(ssq1-ssq2/toa) ) / p0
ssqt = ssq0*exp( lath*(ssq1-ssq2/ttt) ) / pt
ssq1 = 7.93252 _d -6