<div dir="ltr"><div>[1]</div><div>If you use:</div><div><br></div><div>%> <a href="https://www.ncl.ucar.edu/Document/Tools/ncl_filedump.shtml"><b>ncl_filedump </b></a>MCD19A2.A2019071.h10v11.006.2019073033901.<b>hdf|</b> less</div><div><br></div><div>Variable: f</div>Type: file<br>filename: <b>MCD19A2.A2019071.</b>h10v11.006.2019073033901<br>path: MCD19A2.A2019071.h10v11.006.2019073033901.hdf<br> file global attributes:<br> <b>HDFEOSVersion : HDFEOS_V2.19<br> StructMetadata_0 : GROUP=SwathStructure<br></b><div><b>END_GROUP=SwathStructure</b></div><div><b><br></b></div><div>[SNIP]<b><br></b></div><div><br></div><div>This indicates the file use the<b> HDF-EOS </b>'convention' <br></div><div>HDF-EOS supports three geospatial data types: grid, point, and swath, <br></div><div>These data types are <b>not</b> accessible via standard HDF library calls.</div><div>NCL has a very nice simple feature that will return additional [value-added] information.</div><div>For example: There are no latitude/longitude variables on the original .<b>hdf</b> file. <br></div><div><br></div><div>NCL recognizes addtional 'virtual' file extensions: .<b>he2</b> [<b>.he5</b>] or <b>hdfeos</b>.</div><div><br></div><div>%> <a href="https://www.ncl.ucar.edu/Document/Tools/ncl_filedump.shtml"><b>ncl_filedump </b></a>MCD19A2.A2019071.h10v11.006.2019073033901.<b>hdf.<span style="color:rgb(0,0,255)">he2</span>|</b> less</div><div><br></div><div>The file will be 'cleaned-up' . Further, it will create classic latitude/longitude arrays. EG<br></div><div><br></div><div> double <b>GridLat_grid1km ( YDim_grid1km, XDim_grid1km )</b><br> projection : Sinusoidal<br> corners : ( -20.00416666486723, -20.00416666486723, -29.99583333063962, -29.99583333063962 )<br> long_name : latitude<br> units : degrees_north<br><br> double <b>GridLon_grid1km ( YDim_grid1km, XDim_grid1km )</b><br> projection : Sinusoidal<br> corners : ( -85.13204127176456, -74.49885022173193, -80.83045538778872, -92.36735390691214 )<br> long_name : longitude<br> units : degrees_east</div><div><br></div><div>These can be read and used by an NCL script.</div><div><br></div><div>----</div><div>The attached <b>MCD19A2.JasonHerd.ncl </b>should do what you want. <b>It uses the virtual .he2 file extension.</b> It is a slight modification of:<br></div><div><br></div><div><div><a href="https://www.ncl.ucar.edu/Applications/binning.shtml" target="_blank"><b>https://www.ncl.ucar.edu/Applications/binning.shtml</b></a></div><div>Example 1</div><div><br></div><div>There is much about Example 1 that is a bit out of date ...but it works!</div></div><div>Note that this is 'read-and-plot' script. No filtering by data quality flags has been performed. That is your job.</div><div><br></div><div><b>I do not have a complete set of MCD19A2 files for March 12, 2019. </b>However, the attached plot gives the flavor of what you will get.<br></div><div><br></div><div>===</div><div>The golden rule of data processing is look at tour data. The are multiple orbits per file.</div><div><br></div><div>Good luck<br></div><div><b></b></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jun 14, 2019 at 12:48 PM Herb, Jason <<a href="mailto:jherb@albany.edu">jherb@albany.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Hello,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I am working with the MCD19A2 aerosol data. To get an idea as to the overall coverage extent and quantity of measurements for each 1x1km grid box, I am attempting to make all measured data to be " 1 " and then add with each orbit pass in the file (Ex. the location
being studied has a max of 6 passes for a day.) Assuming the full number of passes over the studied area and say fully clear skies to give no missing data when I map the replaced 1 into each grid point and sum the data points the map should show a uniform
map of 6.</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
That being said how can I replace any of the AOD data that falls in to the range of 0.00-1.00 to become 1 and sum up for from the array data?</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Thanks,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Jason</div>
<br>
</div>
</div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>