Product Suite: Surface Reflectance --Polarization correction, de-trending and gain adjustments of Terra and Aqua L1B --Improvement to the aerosol retrieval and correction algorithm, and use of new aerosol retrieval LUTs--Refinements to the internal snow, cloud, and cloud shadow detection algorithms. Uses BRDF database to better constraint the different thresholds used.--Process ocean bands and provide QA datasets for these bands--Improved discrimination of salt pans from cloud and snow and flag salt pan in band QA.

Product Suite: Thermal Anomalies and Active Fire--Refinements to internal cloud mask, which sometimes flags heavy smoke as clouds--Fix for frequent false alarms in the Amazon that are caused by small (~1 km2) clearings within forests--Fix to correct a bug that causes incorrect assessment of cloud and water pixels adjacent to fire pixels near the scan edge--Detect small fires using dynamic thresholding--Process ocean and coastline pixels to detect fire from oil rigs.

Product Suite: Sea-ice --No change to the science algorithm.

Product Suite: Land Surface Temperature and Emissivity--Remove cloud-contaminated LSTs not only from L3 LST products but also from L2 LST products.--Update the coefficient LUT for the split-window algorithm with comprehensive regression analysis of MODIS simulation data in bands 31 and 32 over wide range of surface and atmospheric conditions, especially extending the upper boundary for (LST -- Ts-air) in arid and semi-arid regions and increasing the overlapping between various sub-ranges in order to reduce the sensitivity of the algorithm to the uncertainties in the input data (i.e., column water vapor and air surface temperature from MxD07).--Make minor adjustments in the classification-based surface emissivity values, especially for land-cover type of bare soil and rocks.--Tune the day/night algorithm by adjusting weights to improve its performance in desert regions where the incorporated split-window algorithm may not work well.--Generate new LST products for 8-day and monthly at 6 km grids (in response to user community requests).

Product Suite: Vegetation Indices--Generates 16-day composite VI using two 8-day composite reflectance product (MxD09A1) in the 16-day period.--This Surface Reflectance Input is based on the Minimum Blue compositing approach used to generate the 8 day LSR product--Product format is consistent with the C5 version of the product generated using the L2G daily surface reflectance product.--Update to the Long term Global CMG Average Vegetation Index Product database used to fill the gaps in the CMG product suite

Product Suite: BRDF/Albedo--Generated daily, still using the 16-day composite algorithm.--Observations are weighted to estimate the BRDF/Albedo at the first day of the 2nd 8-day period of the 16-day period.--Uses improved back up database -- pixel based updated from latest full inversion as opposed to the land cover based DB used in C5.--Improved quality and more retrieval at high latitudes from use of all available observations. C5 used only four observations per day.--Uses current day snow status instead of the majority snow/no-snow status from the 16-day period.--CMG Albedo is estimated through retrieval of albedo using all the clear sky observations within 1 km grid for MCD43D as opposed to aggregating from 500 m albedo--Uses L2G-lite surface reflectance as input

Product Suite: LAI/FPAR--Uses daily L2G-lite surface reflectance as input as opposed to MODAGAGG used in C5--Products are generated at native resolution of 500 m--Uses improved multi-year land cover product

Product Suite: GPP/NPP--Uses updated BPLUT and updated version of the daily GMAO--Generate product at native resolution of 500 m using the 8-day composite LAI/FPAR at 500 m resolution

Product Suite: Gridding and L2G products--Better geolocation accuracy

Product Suite: Snow Cover--Generates C5 snow as input to the downstream processes: L2 snow for LST and L3 daily snow for BRDF/Albedo

Product Suite: Burned Area Intermediate--Improvement to the quality filtering algorithm and additional quality flags

Product Suite: VCF Intermediate--No changes

Product Suite: L1 Subsets--AERONET and FLUXNET subset updated to include additional sites--New set of subsets generated for over 400 sites called the "BELMANIP" sites.

PGE26 v6.0.2. Correct the angle SDSs attribute. Add MOD_PRL3CRS process for coarse product M*D13A3C.

PGE27 v6.0.5. Change perl script to handle the new long term 5year CMG DB ancillary of gap filling.

PGE28 v6.0.0. Updated PGE28 to port to C6. Added the new Metadata files.

1. Pass QA to function calc_albedo() (line 2344 in funcs.c), in which the band QA and mandatory QA can be altered to fill value if the calculated albedo is invalid.

2. Add function init_qa() to util.c (line 153).

3. Pass starting row and ending row number to the function main() to be able to split a tile into multiple simultaneous processes.

PGE12 v6.0.13, PGE13 v6.0.32, PGE14 v6.0.6: L2G process updated to work with changes in L2. PGE13 updated to generate the L2G-lite product MxD10GA work with a new snow specific sorting scheme.

PGE43 v6.0.11, 45 v6.0.2: L3 daily and 8-day PGEs updated to work with the L2G-lite input and new science changes made in L2.

Ported the science version of the code for operational production of the product at MODAPS with following changes

Generate 500m resolution product using 500m resolution input.

Use GMAO for daily meteorological data.

Change all climate variable inputs from binary to HDF format.

Use the C6 DEM replacing the tiled SRTM data.

Use C51 multi-year land cover product MCDLCHKM in place of the C4 1km resolution land cover.

PGE21 v6.0.13 Update to MOD_PRSS.

PGE23 v6.0.38 Update to MOD_PRSS.

PGE29 v6.0.7 Update to MOD_PRSS.

PGE31 v6.3.3 Update to MOD_PRSS.

PGE34 6.0.8 New MOD_PRSS products.

PGE37 v6.0.7 Update to MOD_PRSS.

PGE130 v6.0.2. Updated the operational version of the generic L2G code to process the L2 swath product and generate the L2G product (PGE 130).

PGE120 v6.0.29 - updated projection code to coincide with SIPS versions

PGE121 v6.0.16 - modified output ESDT to have an _t.

PGE125 v6.0.28 - fix for partial polar tiles; fix for compiler warnings.

PGE126 v6.0.18 - fix typo from copy+paste. A variable had the wrong name.

Updated MOD_PRLCAT - Bug fix was made to MOD_PRAlg17/alg17.cpp in the marine low cloud section of the algorithm. There was a "one-off" indexing problem based on the month of input data ('m').

- To handle end of the year and days from next cycle

- Evaluate alternative pixels as opposed to the MVC (dictated by the 500m selection)

- Apply a thresholds test based on NDVI and RED/NIR ratios

- Change the logic for composite to bin optimize NDVI under 30VZA for good quality data

- Apply new narrow-to-broadband coefficients without band 6 if only band 6 is fill value due to non-functional or noisy detectors and the other bands are valid full or magnitude retrievals.

- Polarization correction

- Residual trending and MODIS Terra-to-Aqua cross calibration

- Retrieve surface BRDF from accumulated multi-angle set of observations

- Fix the code reading the LUNs for processing parameters

- Change code to measure height in km

- Initialize variable “count” to fix error estimate

- Change cloud processing to choose c95 instead of c66

- Set QC bits 1 and 0 to 10 when LST > 400K

- Use different PWV threshold and smoothing scale when Tgmean < 295K

- Clamp gi to (-2.0 3.0)

- Change ASTER minimum from 0.7 to 0.8

PGE132 v6.0.15 -- updated PGE*.pl to correct extra single quotation mark; updated LoaderModule to check for MaxProfile and run for PM1M files.

PGE130 v6.0.3 - Use of PGE12/MOD_PRMGR/HISTORY.txt (v6.0.21); Fixed the data type, units and valid range for MOD21 dataset "PWV" in lh2g_lst21.h and lhdf_lst21.h files to reflect the new PWV sds in MOD21 L2.

PGE131 v6.0.16 - Slight change to the QA legend to match the UG.

Update pre_process_*VI.c , timesat_*VI.f90, timesate_*VI_2nd.f90, gapfill_*VI.c. do_gapfill.csh

1) Fixed LAI wrong parameters MxD15LAIA2.txt and MxD15LAIA3.txt

2) Updated LUTs: MxD12IGBP_lut.xml and MxD44WM_lut.xml

3) Filling empty tiles and missing values with water color for MxD12IGBP and MxD44WM (MODIS_Sinusoidal_L3_2_LL.sh, sds2png_xml_c, added sds2png_xml_emptytile.c)

4) Changed background color of MxD43NBR from blue to gray (MODIS_Sinusoidal_L3_2_LL.sh)

5) Added GDAL option to fix gaps in polar projections (MODIS_CM_EASE_L3_2_polar.sh)

1. Create a new PGE (PGE157) for temporal gap-filling of all bad quality or out of range LAI-FPAR (MxD15A2H) input values.

2. Assign new ESDTs (MOD15A2HGF/MYD15A2HGF) for these new year-end gap-filled LAI-FPAR products.

3. Set up new production rules and recipes to employ PGE157 to generate the MxD15A2HGF products at the end of each year.

1. All land PGEs updated to reflect collection specific metadata changes and correct SDS attributes scale, offset, fill value, unit and valid range

2. L2 PGEs updated to read and process C61 version of the L1B and geolocation product.

3. Snow PGEs updated to reduce ice/now commission errors and developed a new Cloud Gap filled daily L3 product (MxD10A1F)

4. Fix known issues in DSR/PAR (MCD18) process, generate the product at 1km resolution and generate new CMG products

5. Update the JPL-LST process (MxD21) to use GEOS data in place of MERRA2 as ancillary, and develop new process to generated daily, 8-day and monthly CMG products.

6. Update the VCF intermediate 16-day composite process to fix a known bug in the LWM which is also put in C6 forward processing starting Jan 2019

7. Update VI code to fix known product issues and improve product quality. The change was put in C6 forward processing starting day 2017254

8. Update BRDF PGE to take into account the dead detector in band 5, while estimating the narrow to broad band coefficients

9. L2 Surface Reflectance updated to use a new approach to flag the aerosol quantiy in the L2 swath product.

1. Introduced correct initialization of input metadata structure array in the PGE13 main function. This metadata structure array is used to keep track of each input L2G file and read in all file level metadata from each input file. Incorrect initialization of this array was causing errors, when less than expected number of input files were present, giving rise to missing RANGEBEGININGDATE values.

2. The incorrect SPSOPARAMETER value, seen in some of the L2G files output from PGE12, was fixed by fixing the buffer setting of the "l2g_spso" variable in "L2G_main" function. This array buffer was not set correctly in C6 to reflect several new additional L2G files that were introduced in C6.

1. Define new GIBS MAIAC ESDTs and update the GIBS ICD for MAIAC STD and NRT images

2. Include new color LUTs for MAIAC MCD19A1/A2 products.

3. Updated code to distinguish between the MAIAC STD and NRT GIBS ESDTs.

4. Added loader module specific changes to load MAIAC inputs

1. At the PGE level, a condition was enforced in PGE12, whereby the PGE would throw a fatal error and fail if it failed to open an input file that has been specified in the PCF.

2. An additional check was implemented in the PGE12/PGE13/PGE14/PGE15/PGE130 and PGE140 loader module that would force the inputs to be copied to the working directory instead of being linked from /tis/lads/data/

1. Changed "if(all->wData->HD_snow_qa[row][col] != 255)" to "if(all->wData->HD_snow_qa[row][col] < 2)"

1. Split MxD09CMG1 into two files, MxD09CMG1 and MxD09CMB1

2. Updated M{OY}D09CM{AG}1.fs;

1. Fix fatal exit on bad pixel

Most of the processes for generation of these gap-filled products for C61 are similar to C6 (ref: PCR-19-018, 19-019). The following are the additional improvements for C61:

1. qa_map.h: Changed instances of "aerosol quantity" to "aerosol quantity:level of uncertainty in aerosol correction".

2. versions.h: Updated the PROCESSVERSIONNUMBER to 6.4.9

Most of the processes for generation of these gap-filled products for C61 are similar to C6 (ref: PCR-19-018, 19-019). The following are the additional improvements for C61:

1. Added support of reading in and processing of the gap-filled daily BRDF/Albedo input (MCD43GFbsasw) product.

2. Enable the 8-day process to exit with FATAL ERROR when the number of daily data days processed, exceeds the standard number of days of 8 (5 or 6 for the last 8-day of the year).

3. Added daily GMAO and annual average air temperature file names to the INPUTPOINTER of dailyMOD16A1 or MOD16A1GF HDF file.

4. Change the scale_factor for "LE_500m" and "PLE_500", from 10000 to 0.01 and increase their units from J/m^2/day to MJ/m^2/day for more reasonable magnitude of the final physical values.

5. Enable the archiving of the daily rolling intermediate files at every "n-day" interval.

1. Updated LM (RealTimeDelay of 'DAYS=2' for M1_GAC AVHRR inputs)(Delay running PGE199 -- run PGE199 on data-day X on actual-day X+2)

1. Updated a solar zenith threshold from 105 to 110 degrees in no_EPH_prescan() (in file PGE199/GAPS/src/navigate.c). Several other non-science changes (metadata changes) for outputs were also made (in postprocessing programs), to reflect the contents of the SDSs more accurately.

- Added lat/lon prescreening to the "Projection into grid" portion of main() (gaps.c) -- if geolocation has problems (zeros in it (lt/gt +/-0.005)), replace such values with -999.0.

1) Updated LoaderModule to include $startdate/$enddate to handle end of year for 8-day products.

1. Correction for PGE23N_LoaderModule to properly skip the tiles that have previously been flagged as unnecessary because they are not produced for MOD10A1.

1) Changed background color from grey to black for MAIAC RGB products.

2) Updated color LUT MCDWDL3FLD_lut.xml (in the Paramfiles folder) to set 'NODATA' color from transparent to grey;

3) Special LoaderModule for AMPM GL3 recipe to skip Profiles 3&4 for forward and NRT processing until PGE113 is added.

1) Added sanity checks for geolocation data when determining lat/lon extent of data in a granule (routine NO_EPH_prescan() in navigate.c) and when interpolating the lat/lon data read from the AVHRR file itself (routine interpol_noaa_geom() in file calib_data.c).

CountPGE('${PreviousPGE[$i]}','$startdate','$enddate',{Tile=>undef,Profile=>1}) >= $TileTotal[$i]

To

CountPGE('${PreviousPGE[$i]}','$startdate','$enddate',{Tile=>undef,Profile=>1}) >= $TileTotal[$i]

To

1. New version of code in python2/3.

2. Adds header to output file: "latitude,longitude,Mid-IR,scan,track,acq_date,acq_time,satellite,confidence,version,Long-IR,frp,daynight".

3. Header only files are deleted in the COMB/PGE603.pl file by checking m5sum and file size.

Updates to PGE145:
- Added new PGE145L for leading edge process
- PGE145 will use all data from previous year + current year + next year
- PGE145L will use previous year 12 months + current year 12 months
- Updates to PGE145 Perl script to add the met file inside the tar ball

Updates to PGE146:
- Updates to Perl script to add the met file inside the tar ball
- Use only current year PGE145 outputs as PGE145 contains data for all three years in a single tar file.
- Update the valid range values
- Update the unit values
- Update to add the scale_factor for EVI_Area, EVI_Minimum and EVI_Amplitude
- Update to change the string "fill_value" with "_FillValue
- Updated MOD_PR12Q2 to V6.1.1
- The function MCD12Q2C6_AnnualPhenologyFunctions.R is updated for the following:

-- Change QA_Detailed to not redundantly report QA_Overall.
-- Then 32767 can be maintained as the fill value for all SDS and QA_Detailed needs only 0-16383 range.
-- The function "DefaultPhenoParameters" is updated to include new default values
-- The function "GetSegs" is modified to allow for the specification of both absolute and relative amplitude thresholds.

Updates to PGE150:
- Update class_name_list.txt to change types for UMD, LCCS2 and LCCS3
- See MOD_PRBIN2HDFLC HISTORY.txt 6.1.2
- Patch PGE150.pl to used correct class_name_list.txt file (2022-05-25)

Updates to PGE151:
- Updates perl script PGE151.pl by changing $MOD_PR12CMGA_cmd through adding if($LC_list eq "LC2") so that "-u 1" argument works for UMD classification case to drop two types 11 and 14 off the output binary type fraction layers in BIP band interleave.

Updates to PGE41:
- In quarterly.c, comment out the lines between #5314 and 5341 under this statement: /* Override LCCS3 legend with the correct legend defined here */
- Comment out the lines between #5499 and 5516 under this statement:/* Set the correct legend of LC_Prop3 here */
- In quarterly.c, comment out #5478 to 5516 not to Set the correct legend of LC_Prop2and LC_Prop3
- Set the correct legend of LC_Prop2 and LC_Prop3
- Patched PGE41.pl to export coarse products (2022-06-03)

C61 Specific Algorithm Change/Improvements:

- Update MOD_PR44B3/MOD44B3.c MOD_PR44B3/make_grid.c to update valid_range from 253 to 100 for sds Percent_Tree_Cover, Percent_NonTree_Vegetation, Percent_NonVegetated.

- Update MOD_PR44B3/MOD44B3.c small change that populates the outputs with fill values when there are fill values in the input files.

- Changed DFACC_RDWR to DFACC_READ in MOD_PR44B3/MOD44B3.c to read MOD44W input.

- Added the water value as 200 and the scale factor as 0.01 to the metadata of the appropriate layers.

For all science changes, refer to section 8 above. In addition, the following are to be noted:

1. New version of code in python2/3.
2. Adds header to output file: "latitude,longitude,Mid-IR,scan,track,acq_date,acq_time,satellite,confidence,version,Long-IR,frp,daynight".
3. Header only files are deleted in the COMB/PGE603.pl file by checking m5sum and file size.
4. Added code to fireloc2.py to deal with inconsistent handling of single element arrays. i.e., the error that we see in scan/track when there is only one fire pixel entry in the input files.

The full set of C61 MAIAC changes have already been summarized in the PCR# 22-004. This PCR only highlights the key changes made for NRT processing.

1. Updated PGE113N loader module Perl script to do the following:

1. Added MCD19A3DN.C61.MCF file for MCD19A3DN NRT files.
2. Process polar tiles following the same logic used for PGE113F loader module.
3. Added M[OY]DMAIACOqkm inputs.

• New PGE (PGE63) that can be run in production.
• Update post processing ancillary data.
• New seven class layer.
• Production of Antarctica tiles.

1. Integrated science code delivery. The C6-C7 blending is now carried out from 2024-03-31 to 2024-04-01.
2. Fixed the problem related to incorrect InputPointer attribute in MOD02PC*KM.
3. Integrated auto-formatting tools.
4. Enabled dynamic staging of coeff files in LM script. If PGE is run for Terra, the LM script stages two sets of xcal coeff files.

1. set 1: v1.13
2. the most recent version available in the database

Only major changes have been shown here for brevity:

1. Fixed a bug in MODIS MAIAC C61 code, which will affect snow fraction output.
2. Changed MAIAC_File.cpp

PGE03 v6.0.47 -- Use destripe_config_terra.dat.v5nrt in NRT.

PGE06 v51.0.18 -- Added destripe_config_terra.dat.v5nrt for nrt. Matched with the logic used by C6 for chosing destripe_config file.

PGE04 6.1.3 - Updated MOD_PR04DB to version 6.1.3

PGE06 6.1.2 - Updated MOD_PR04DB to version 6.1.3; Use of new destripe_config_terra files (revision 2)

PGE85 6.1.0 - Release of Collection 6.1 version

PGE93 6.1.0 - Produces products with filenames and metadata appropriate for Collection 6.1

PGE90 6.1.0_5 - Updated metadata for Collection 6.1; Corrected to resolve Bug 1674.

PGE108 6.1.0 - Release of Collection 6.1 version

PGE83 6.1.0 - Release of Collection 6.1 version

PGE69 6.2.0 - Release of Collection 6.1 version

PGE56 6.2.0 - Release of Collection 6.1 version

PGE57 6.2.1 - Updated template file to version 620

PGE70 6.2.0 - Release of Collection 6.1 version

PGE55 6.1.0 - Release of Collection 6.1 version

PGE81 6.1.0 - Release of Collection 6.1 version

PGE89 6.1.0 - Release of Collection 6.1 version

Updated MOD_PR06OD to V6.1.1

1. A new release of C61 version of PGE98, v6.1.0.
2. The new PGE98 contains the science code revisions (fixed the ice concentration over inland), new coefficient file, and updated input requirements.
3. Resampled the new input daily OISSTD21_AVHRR files and produced REYNSST equivalent REYNSSD21 data.
4. The new REYNSSTD21 files is used for downstream PGE03, PGE06, and PGE85 products.
5. This new daily oisst file will be used for reprocessing since January 2023 and forward processing.

For all science changes, refer to section 8 above. In addition, the following are to be noted:

1. New version of code in python2/3.
2. Adds header to output file: "latitude,longitude,Mid-IR,scan,track,acq_date,acq_time,satellite,confidence,version,Long-IR,frp,daynight".
3. Header only files are deleted in the COMB/PGE603.pl file by checking m5sum and file size.
4. Added code to fireloc2.py to deal with inconsistent handling of single element arrays. i.e., the error that we see in scan/track when there is only one fire pixel entry in the input files

Initial version of the science processing algorithm to make the MEaSUREs GEO-LEO DT Aerosol product was developed by DT Aerosol science team. Loader module (LM) and production rule updates and additions were made by MODAPS STIG to complete the transition of this process to a MODAPS PGE.

Comparing with C61 MODIS DT product, Specific Algorithm Change/Improvements for this DT Package:

1. A new release of Collection 1 version of PGE104, v1.0.0.
2. Uses all L1b data, geolocation, angles, terrain height, Land Sea mask interpolated or congridded to 1km resolution.
3. Final output retrieval is for 10 KM resolution, so 10*10 pixels are used for averaging overboard.
4. Ocean Cloud Masking uses 3x3 spatial variability test at 0.66 �m (500m).
5. Land Cloud Masking and NDVI Water Masking: Cirrus cloud testing is at 1.37 (1km) �m band, and NDVI water masking is done at 1km.
6. All masks are done for whole granule so there is continuity.
7. Change the data type of the bounding coordinates from string to float.

1. The crosstalk correction algorithm described in #8 has been implemented in Terra MODIS PGE02 V6.2.0.

2. A penalty to the uncertainty is added to bands 27-30 based on the magnitude of the corrected signal.

3. The detector_quality_flag_value table has been rederived from the start of the mission to better reflect the QA status for detectors in bands 27-30.

1. The processing version number has been updated to 6.2.1.

I was instructed to increase the maximum sensor zenith angle to 70°, just to be safe.

Major improvements implemented in the C7 L1 products include:
Geolocation:

• Updated geolocation LUTs in pointing time series that correct for inaccurate predictions in the previous forward-processing Collection (and the slight differences in the errors measured by refresh chip library and old existing chip library). They are also extended to the end of the missions (2027).
• Updated Land Water Mask from a static database to year-by-year databases.
• Added a solar eclipse data field.

L1B Calibrated TOA radiance/reflectance:

• Mission-long crosstalk correction applied to the Aqua PV LWIR & MWIR bands and Terra MWIR bands.
• Improved characterization of the VIS bands gains and RVS, especially in Terra MODIS by correcting the polarization effects on the desert data.
• Enhanced offset and non-linear coefficient estimation for several MODIS TEB
• Time-dependent RVS corrections in the Terra MODIS SWIR bands.

500m Geolocation Layers:

• New 500m latitude/longitude layers were added to the L1/Geolocation product (MxD03). The PGE01 will retain the 1km and 500m geolocation and height layers.

Geolocation LUT Update:

• Updated C7/GEO LUT v7.0.11 to address the geolocation errors associated with C7/Terra L1 processing, for the latter half of 2000 and earlier portion of 2001.

Removal of Native Resolution Geolocation Layers:

• Updated PGE02 to remove native resolution geolocation layers from L1Bs that was previously added.The C7/L1B outputs will resemble the C61 version.

MODIS_Combine_V6.0.25 --

1) fix to calculation of $collection in MODIS_combine_global.sh

MODIS_Combine_V6.0.25 --

1) fix to calculation of $collection in MODIS_combine_global.sh

1) RR_combine_tile.sh - fix for partial polar tiles

2) composite_subsets.exe and sds2hdf.exe - changed "tweaks to remove compiler warnings" to more appropriate tweaks (and sync with PGE121 versions)

The changes for MODIS_Combine V6.0.31:

1) MODIS_combine_tile.sh - fix for partial polar tiles

Update the IDPS L2 SR-IP process for operational processing at Land SIPS with MODIS LSR C6 algorithm changes

1. Enable it to work with NASA L1B inputs in HDF5/netCDF4 format.

2. Include the land water mask (LWM) from the geolocation (VNP03/VJ103) products instead of the cloud mask, as was being done in AS 5000.

3. Put in a fix to ignore RSB bands during the night time processing

4. Put in a fix to include the inland water bodies as part of LWM that was ignored in #2.

1. Use of NASA Land SIPS variation of L1B/Geo code delivered by L1WG (V3.1.0.rc).

2. Changes to the GEO LUT approach to include:

a. VIGMU (VIIRS instrument geometric model update), LUT and code changes

b. Corrections for temporal pointing variations, unified LUT and code changes

c. Kalman filter “module” in L1A to improve attitude accuracy, thus geolocation accuracy.

d. Update EFL (effective focal length) in LUTs.

3. Changes to the L1B LUT, including:

a. Accounting for the whole mission period for S-NPP VIIRS

b. Improvements to RSB calibration factor estimations:

4. The major improvements on RSB are:

5. Apply 7 years of lunar F to calculate H_rta, from comparing with SD F(H_sdsm).

6. Use DCC and Libya4 data to determine H_rta SD positional dependence to remove the striping seen in the current version.

7. Solar azimuth dependence of H_rta for all bands.

8. Apply new DeltaC table for entire mission.

a. Improvements to DNB:

9. Update DN0 table to remove striping seen in the current version.

10. Update LGS gain and gain ratios with updated H-factor.

1. Updated num_day_375m_rows to prevent running on Night granules

2. Updated code for new filenames

3. Changed LoaderModule and production perl scripts to allow for different time intervals to run the code ex. 1 hour

PGE603 (FIRMS Fire table generation):

1. Updated to run at 1 hour intervals

2. Added code and logic for J1 FIRMS products .

PGE620 (Swath to tile mapping):

1. Updated to run at 1 hour intervals

2. Removed some of the output products for PGE620L1 (assorted VNR02 products)

3. Corrected VNR01 files to be VNR02 (in V2.0.2)

4. Added PGEVersion as a metadata element

5. Corrected TileDataDay to write proper grantime values

6. Made changes to find the NRT and J1 geolocation file names.

7. Removed _NRT from output ESDTs.

8. Added missing Geolocation Filename variations for STD runs (not relevant for NRT)

9. Fixed a missing entry looking for VNP geolocation Filenames in projection.c for STD run that was accidentally overwritten by an NRT filename search.

For VIGMU, the mag matrix as switched sign and the u_aft vector is flipped within the code. This also requires changing the alpha, beta, and gammaa angles within the LUT to accommodate the change.For the time-dependant attitude corrections, 3 new variables were added to the LUTs: a count, an array of times, and an array of roll/pitch/yaw correction angles. The arrays of times and corrections were assigned to be 1000 long, and the count variable tells the code which elements of the arrays are valid. This method allows the LUT to not change format/size while allowing essentially more elements to be added to the arrays. Currently only about 50 elements are used, so it is unlikely that more elements will be necessary when considering the lifetime of the mission.

1. Revision to daytime.c to correct flaw occassionally causing a '17' to appear in the fire mask over sun glint.

1. Update to num_day_375m_rows.c to use both sensor_mode flags and solar zenith angle to detect day granules.

PGE620 2.0.8 (Swath to tile mapping):

1. Updated day_night_rows.exe to remove day & night mixed rows by solar zenith angle.

2. Added check to PGE620*_LoaderModule.pl to see if granule to be processed is in the DataHiddenFromProduction table, and if it is, go to NoRun.

Initial version of the PGE641/642 was developed at SCF in Univ. Of Maryland, followed by suitable loader module (LM) and production rule updates and additions by LandSIPS STIG to complete the transition to a LandSIPS PGE. The upstream L1B gridding process (PGE635) was developed by STIG along the lines of its MODIS equivalent.

Specific Algorithm Change/Improvements, following the initial delivery:

• Integration of all necessary LM and Production rules.
• Integration of browse products for export to the DAACs.
• Several bug fixes by STIG and SCF to avoid certain memory errors on specific production hosts and to prevent failures when trying to open a read-only file.

Initial version of the Global Reservoir suite of PGEs, PGE687, PGE688, and PGE698, developed at SCF in Texas A&M University, followed by suitable loader module (LM) and production rule updates and additions by LandSIPS STIG to complete the transition to a LandSIPS PGE.

Specific Algorithm Change/Improvements, following the initial delivery:

• Integration of all necessary LM and Production rules.
• Several bug fixes by STIG and SCF to avoid certain memory errors on specific production hosts and to prevent failures when trying to open a read-only file.

The PGE510 V2.0.13 bit-mask definitions were corrected to set the the QA bits properly. In addition, the following maintenance revisions were made in this update:

1. corrected several malloc() calls for which the return value was cast to the incorrect data type;
2. replaced malloc() with calloc() in main processing module to ensure that all internal bit fields were initially clear;
3. increased netCDF compression level to reduce final product size;
4. corrected bug in input module that could potentially use an uninitalized array to check for invalid reflectance values;
5. minor changes to "fire mask" data-set legend for clarity. Changes to the science code were made by the current PI Louis Giglio after review of C2 datasets. Changes to PGE510.pl to correct missflagging of outputs implemented by SSTG.

Revision to the PGE554N_LoaderModule.pl include:

1. PGE554N runs hourly. If not processing the 2300 hour of the day, check if PGE554G has already finished processing the next hour, and if so, did it generate data for the tile to be processed. If it did, put current RI for that tile/hour to NoRun and go to the next hour for processing.
2. If RI start time is not 0000, then check if for the previous hour, of PGE554N is not completed yet (Status is not "RunOk", or "NoRun"). If still waiting to run, then immediately go to TryAgainLater.