(OAD) Document for Ozone Mapping and Profiler Suite ... - NASA

OAD – OMPS TC EDR D37031
Rev C5
Page 19
intnor: performs table look-ups for the Northern Hemisphere
intsou: performs table look-ups for the Southern Hemisphere
plint: performs polynomial interpolations.
2.1.2.9 Subroutine residue.f
Compute two sets of residues and sensitivities for wavelengths comprising from two sets of
profiles (low and mid for |lat| < 45�, mid and high otherwise). Also, calculate residues for the
336 nm and 377 nm wavelengths.
1) For each of the two latitude bands and wavelengths, use �0 to determine I0, I1 (from
Z1), I2 (from Z2), T, and Sb for the 1.0, 0.7, 0.4, and 0.1 atm tables. Using f*, Rt*, Rc* and
reversing the calculations in subroutine oznot, determine Ncalc.
2) Calculate residues:
meas calc N N r � �
3) Calculate the OMPS “aerosol index” (R336-377).
Subroutines called by residue:
prfind: determines profile index (the pointer that indexes the correct profile in the table);
qvalue: determines table values It calls the following subroutines:
intnor: performs table look-ups for the Northern Hemisphere
intsou: performs table look-ups for the Southern Hemisphere
plint: performs polynomial interpolations.
2.1.2.10 Subroutine ozone.f
Compute a new ozone estimate for each of the 3 wavelength triplets.
1) Determine which three triplets to use by calculating the optical path length:
2) For each of the 3 triplets:
�
� est
somega
�
(sec 0
� � sec�
)
1000
a. Compute two separate values of � at two latitudes using the two sets of
residues and sensitivities obtained in subroutine residue and the following formula:
r1��2
� r2��1
� profile � �0
�
; � �i � �i
� �R
�1��2
��2��1
where subscripts refer to two selected wavelengths, � is the sensitivity, and r is the
residue.
b. To account for the change in calculated ozone value, adjust residues using the
following formula:
r
adj
�r � 1�
f ) � r � f �� ( � � � ��
� )
prf 1 ( m prf m profile 0