crival = 2
omitot = []
gometot = []
omistat = []
gomestat = []
omitemp = np.zeros([3])
gometemp = np.zeros([3])
res = 0.5
for iyear in years:
for imon in months:
for i, iday in enumerate(days):
# latmod, lonmod, aodmod0 = gridMODIS(iyear, imon, iday, ROI, res, data)
latgome, longome, aaigome0 = gridGOME2(iyear, imon, iday, ROI, res)
latomi, lonomi, aaiomi0 = gridOMAERO(iyear, imon, iday, ROI,
jetlag, 'AAI', res)
latomi, lonomi, aodomi0 = gridOMAERO(iyear, imon, iday, ROI,
jetlag, 'AOD', res)
latcal, loncal, bs532, z = trackCALIOP(iyear, imon, iday, ROI,
'L2')
# aerosite = AERONET_day(iyear, imon, iday, ROIaer)
## plumemsk = np.zeros(aaiomi0.shape)
plumemsk = np.logical_or(np.isnan(aaiomi0), np.isnan(aaigome0))
plumemsk = np.logical_or(plumemsk, aaigome0 < crival)
plumemsk = np.logical_or(plumemsk, aaiomi0 < crival)
#
## plumemsk = np.logical_or(aaigome0 < crival, aodmod0 < 0.)
## plumemsk = np.logical_or(plumemsk, aodomi0 < 0.)
lat, lon, aai, aod, aodstd, sza, saa, vza, vaa, ssa, Ps, plumemskomi, alb = OMAERO4cfg(
iyear, imon, iday, ROI, jetlag, plumemsk, crival, res)
alhval = np.array([
alhest[dd] - 0.5, alhest[dd] - 0.25, alhest[dd],
alhest[dd] + 0.25, alhest[dd] + 0.5
])
altval = np.array(
[altest[dd] - 0.25, altest[dd], altest[dd] + 0.25])
factor[dd, :] = factorval
alh[dd, :] = alhval
alt[dd, :] = altval
latmod, lonmod, aodmod0 = gridMODIS(iyear, imon, iday, ROIval[dd],
res, '3K')
latgome, longome, aaigome0 = gridGOME2(iyear, imon, iday,
ROIval[dd], res)
latomi, lonomi, aaiomi0 = gridOMAERO(iyear, imon, iday, ROIval[dd],
jetlag, 'AAI', res)
plumemsk = np.logical_or(aodmod0 < 0., aaigome0 < crival)
lat, lon, aai, aod, aodstd, sza, saa, vza, vaa, ssa, Ps, plumemskomi, alb = OMAERO4cfg(
iyear, imon, iday, ROIval[dd], jetlag, plumemsk, crival, res)
aaiomi = np.ma.masked_array(aaiomi0, plumemskomi)
aaigome = np.ma.masked_array(aaigome0, plumemskomi)
aodmod = np.ma.masked_array(aodmod0, plumemskomi)
for ff, ifactor in enumerate(factorval):
for hh, ialh in enumerate(alhval):
for tt, ialt in enumerate(altval):
if ialh > ialt:
testname = '/fine_LUT_alh-%1.2f_alt-%1.2f_ni-%1.2f' % (