# -----------------------------------------------------------
# MODIS/Aqua sensor-specific atmospheric correction data
# -----------------------------------------------------------

#
# Number of bands
#
Nbands = 9

#
# Wavelengths (um)
#
Lambda(1) = 412
Lambda(2) = 443
Lambda(3) = 488
Lambda(4) = 531
Lambda(5) = 551
Lambda(6) = 667
Lambda(7) = 678
Lambda(8) = 748
Lambda(9) = 869

#
# Extraterrestrial Solar Irradiance (mW/cm^2/um/sr)
# Thuillier, G. Meister, 3-Nov-2003
#
F0(1) = 172.718
F0(2) = 187.688
F0(3) = 194.973
F0(4) = 185.820
F0(5) = 186.594
F0(6) = 152.163
F0(7) = 148.097
F0(8) = 128.123
F0(9) = 95.849

#
# Rayleigh optical thickness (M. Wang, Oct 2004)
#
Tau_r(1) = 0.3109
Tau_r(2) = 0.2375
Tau_r(3) = 0.1596
Tau_r(4) = 0.1131
Tau_r(5) = 0.0994
Tau_r(6) = 0.0446
Tau_r(7) = 0.0416
Tau_r(8) = 0.0286
Tau_r(9) = 0.0155

#
# Ozone absorption (/cm) (HRG, 20-May-02) [ave]
#
k_oz(1) = 0.00244
k_oz(2) = 0.00383
k_oz(3) = 0.0225
k_oz(4) = 0.0663
k_oz(5) = 0.0836
k_oz(6) = 0.0485
k_oz(7) = 0.0395
k_oz(8) = 0.0119
k_oz(9) = 0.00263

#
# co2 transmittance
# Z. Amhad, May 2006
#
t_co2(1) = 1.0
t_co2(2) = 1.0
t_co2(3) = 1.0
t_co2(4) = 1.0
t_co2(5) = 1.0
t_co2(6) = 1.0
t_co2(7) = 1.0
t_co2(8) = 1.0
t_co2(9) = 1.0

#
# no2 absorption coefficient (per molecule)
# Z. Amhad, May 2006
#
k_no2(1) = 5.7317e-19
k_no2(2) = 4.9163e-19
k_no2(3) = 2.8023e-19
k_no2(4) = 1.4372e-19
k_no2(5) = 1.1241e-19
k_no2(6) = 0.08e-19
k_no2(7) = 0.085e-19
k_no2(8) = 0.015e-19
k_no2(9) = 0.0

#
# h2o absorption transittance function
# Z. Amhad, May 2006
#
a_h2o(1) = 1.0
a_h2o(2) = 1.0
a_h2o(3) = 1.0
a_h2o(4) = 1.0
a_h2o(5) = 1.0
a_h2o(6) = 1.0
a_h2o(7) = 1.0
a_h2o(8) = 1.0
a_h2o(9) = 1.0

b_h2o(1) = 0.0
b_h2o(2) = 0.0
b_h2o(3) = 0.0
b_h2o(4) = 0.0
b_h2o(5) = 0.0
b_h2o(6) = 0.0
b_h2o(7) = 0.0
b_h2o(8) = 0.0
b_h2o(9) = 0.0

c_h2o(1) = 0.0
c_h2o(2) = 0.0
c_h2o(3) = 0.0
c_h2o(4) = 0.0
c_h2o(5) = 0.0
c_h2o(6) = 0.0
c_h2o(7) = 0.0
c_h2o(8) = 0.0
c_h2o(9) = 0.0

d_h2o(1) = 0.0
d_h2o(2) = 0.0
d_h2o(3) = 0.0
d_h2o(4) = 0.0
d_h2o(5) = 0.0
d_h2o(6) = 0.0
d_h2o(7) = 0.0
d_h2o(8) = 0.0
d_h2o(9) = 0.0

e_h2o(1) = 0.0
e_h2o(2) = 0.0
e_h2o(3) = 0.0
e_h2o(4) = 0.0
e_h2o(5) = 0.0
e_h2o(6) = 0.0
e_h2o(7) = 0.0
e_h2o(8) = 0.0
e_h2o(9) = 0.0

f_h2o(1) = 0.0
f_h2o(2) = 0.0
f_h2o(3) = 0.0
f_h2o(4) = 0.0
f_h2o(5) = 0.0
f_h2o(6) = 0.0
f_h2o(7) = 0.0
f_h2o(8) = 0.0
f_h2o(9) = 0.0

g_h2o(1) = 0.0
g_h2o(2) = 0.0
g_h2o(3) = 0.0
g_h2o(4) = 0.0
g_h2o(5) = 0.0
g_h2o(6) = 0.0
g_h2o(7) = 0.0
g_h2o(8) = 0.0
g_h2o(9) = 0.0

#
# White-cap albedo (Frouin, May 1999)
#
awhite(1) = 1.0
awhite(2) = 1.0
awhite(3) = 1.0
awhite(4) = 1.0
awhite(5) = 1.0
awhite(6) = 0.889225
awhite(7) = 0.889225
awhite(8) = 0.760046
awhite(9) = 0.644950

#
# Band-averaged water absorption (Pope & Fry 1997 + Kou et. al. 1993)
#
aw(1) = 0.02027
aw(2) = 0.01110
aw(3) = 0.01677
aw(4) = 0.04496
aw(5) = 0.05492
aw(6) = 0.4920
aw(7) = 0.5118
aw(8) = 2.801
aw(9) = 4.727

#
# Band-averaged backscatter (Smith and Baker 1981)
#
bbw(1) = 0.003242
bbw(2) = 0.002454
bbw(3) = 0.001622
bbw(4) = 0.001132
bbw(5) = 0.0009890
bbw(6) = 0.0004262
bbw(7) = 0.0003968
bbw(8) = 0.0002671
bbw(9) = 0.0001404 

#
# Out-of-band water-vapor functions (Yang algorithm)
# from HRG modisoob2.cal.coeff
#
oobwv01(1) =  9.999e-01
oobwv02(1) = -1.359e-05
oobwv03(1) = -1.298e+00
oobwv04(1) =  3.715e-03
oobwv05(1) = -4.213e-05
oobwv06(1) = -1.325e-05
oobwv07(1) =  1.093e-02
oobwv08(1) =  3.093e-03
oobwv09(1) =  5.127e-06
oobwv10(1) =  1.572e-06
oobwv11(1) = -1.116e-03
oobwv12(1) = -4.346e-04

oobwv01(2) =  1.000e+00
oobwv02(2) = -5.484e-06
oobwv03(2) =  5.386e-01
oobwv04(2) =  2.282e-03
oobwv05(2) = -2.380e-05
oobwv06(2) = -4.496e-06
oobwv07(2) =  8.961e-03
oobwv08(2) =  5.385e-04
oobwv09(2) =  3.269e-06
oobwv10(2) =  5.677e-07
oobwv11(2) = -1.594e-03
oobwv12(2) =  1.034e-04

oobwv01(3) =  1.000e+00
oobwv02(3) = -1.121e-05
oobwv03(3) =  1.415e+00
oobwv04(3) =  2.382e-03
oobwv05(3) = -4.739e-05
oobwv06(3) = -2.105e-05
oobwv07(3) =  9.637e-03
oobwv08(3) =  2.045e-03
oobwv09(3) =  4.952e-06
oobwv10(3) =  2.428e-06
oobwv11(3) = -1.341e-03
oobwv12(3) = -2.221e-04

oobwv01(4) =  1.000e+00
oobwv02(4) = -7.868e-06
oobwv03(4) =  1.250e+00
oobwv04(4) =  1.015e-03
oobwv05(4) = -3.302e-05
oobwv06(4) = -1.263e-05
oobwv07(4) =  5.191e-03
oobwv08(4) =  9.931e-04
oobwv09(4) =  3.590e-06
oobwv10(4) =  1.464e-06
oobwv11(4) = -6.739e-04
oobwv12(4) = -8.700e-05

oobwv01(5) =  1.000e+00
oobwv02(5) = -4.530e-06
oobwv03(5) =  4.389e+00
oobwv04(5) = -1.887e-04
oobwv05(5) = -2.210e-05
oobwv06(5) = -9.093e-06
oobwv07(5) =  1.497e-03
oobwv08(5) =  1.225e-03
oobwv09(5) =  2.605e-06
oobwv10(5) =  1.002e-06
oobwv11(5) =  1.395e-04
oobwv12(5) = -2.259e-04

oobwv01(6) =  9.993e-01
oobwv02(6) = -4.413e-04
oobwv03(6) = -1.662e-01
oobwv04(6) =  1.662e-02
oobwv05(6) = -1.815e-03
oobwv06(6) = -1.107e-03
oobwv07(6) =  7.503e-02
oobwv08(6) =  1.323e-02
oobwv09(6) =  1.808e-04
oobwv10(6) =  1.003e-04
oobwv11(6) = -8.489e-03
oobwv12(6) = -2.504e-03

oobwv01(7) =  9.989e-01
oobwv02(7) = -3.736e-04
oobwv03(7) = -4.620e-01
oobwv04(7) =  2.127e-02
oobwv05(7) = -8.808e-04
oobwv06(7) = -3.267e-04
oobwv07(7) =  8.673e-02
oobwv08(7) =  1.670e-02
oobwv09(7) =  9.649e-05
oobwv10(7) =  3.504e-05
oobwv11(7) = -9.963e-03
oobwv12(7) = -2.678e-03

oobwv01(8) =  9.983e-01
oobwv02(8) = -9.267e-04
oobwv03(8) =  3.304e+00
oobwv04(8) = -9.455e-03
oobwv05(8) = -3.692e-03
oobwv06(8) = -2.223e-03
oobwv07(8) = -2.516e-02
oobwv08(8) = -1.816e-02
oobwv09(8) =  3.295e-04
oobwv10(8) =  2.263e-04
oobwv11(8) =  2.104e-03
oobwv12(8) =  1.814e-03

oobwv01(9) =  9.995e-01
oobwv02(9) = -2.086e-04
oobwv03(9) =  3.197e+00
oobwv04(9) = -1.267e-03
oobwv05(9) = -8.618e-04
oobwv06(9) = -4.033e-04
oobwv07(9) =  5.332e-03
oobwv08(9) = -2.000e-03
oobwv09(9) =  9.147e-05
oobwv10(9) =  4.084e-05
oobwv11(9) = -8.560e-04
oobwv12(9) =  1.008e-04

#
# Full band to nominal band nLw functions (S. Bailey, September 2004)
#
ooblw01(1) =  9.691e-01
ooblw02(1) =  3.137e-02
ooblw03(1) = -3.279e-03

ooblw01(2) =  1.004e+00
ooblw02(2) =  2.283e-03
ooblw03(2) = -1.263e-03

ooblw01(3) =  1.002e+00
ooblw02(3) = -3.345e-03
ooblw03(3) = -3.033e-04

ooblw01(4) =  1.031e+00
ooblw02(4) = -2.499e-02
ooblw03(4) =  2.790e-03

ooblw01(5) =  1.041e+00
ooblw02(5) = -5.868e-02
ooblw03(5) =  6.910e-03

ooblw01(6) =  9.933e-01
ooblw02(6) =  8.146e-04
ooblw03(6) = -6.000e-04

ooblw01(7) =  9.995e-01
ooblw02(7) =  7.793e-04
ooblw03(7) = -4.220e-04

ooblw01(8) =  1.000e+00
ooblw02(8) =  0.000e+00
ooblw03(8) =  0.000e+00

ooblw01(9) =  1.000e+00
ooblw02(9) =  0.000e+00
ooblw03(9) =  0.000e+00

#
# Weights for Ed integration to ipar (Sean, 16 Aug 2004)
#
wed(1) = 26.51
wed(2) = 36.90
wed(3) = 45.85
wed(4) = 30.45
wed(5) = 103.65
wed(6) = 48.90
wed(7) = 0.0
wed(8) = 0.0
wed(9) = 0.0

#
# Weights for aph integration to arp (MODIS ATBD 20, Feb 2003)
#
waph(1) = 1.010440
waph(2) = 0.970991
waph(3) = 0.984819
waph(4) = 1.128290
waph(5) = 0.732397
waph(6) = 0.601099
waph(7) = 0.0
waph(8) = 0.0
waph(9) = 0.0