<div dir="ltr"><div>[1] If you do:<br><br>%> ncl_filedump xMOD06_L2.A2011349.0550.006.2015056134430.hdf | less<br><br>You will see:<br><br>[snip]<br>Variable: f<br>Type: file<br>filename: xMOD06_L2.A2011349.0550.006.2015056134430<br>path: xMOD06_L2.A2011349.0550.006.2015056134430.hdf<br> file global attributes:<br> HDFEOSVersion : HDFEOS_V2.17<br> StructMetadata_0 : GROUP=SwathStructure<br> GROUP=SWATH_1<br> SwathName="mod06"<br>[snip]<br><br>So, these files contain satellite SWATH data<br><br>[2] For convenience, when HDFEOS_V2 data are encountered, it is convenient to append the file suffix .he2 or, equivalently .hdfeos in the addfile(s) function. This makes for a 'cleaner' look at the file. Further, in some cases, value added information is added. [Side point: for hdf5-eos use .he5)<br><br>%> ncl_filedump xMOD06_L2.A2011349.0550.006.2015056134430.hdf.he2 | less<br><br>You will see a different look. More like a classic netCDF file. Value added features include adding a 'coordinates' attribute to each variable.<br><br>[3] While the dimension sizes and dimension names are the same for each swath file:<br><br> dimensions:<br> Cell_Across_Swath_5km_mod06 = 270<br> Cell_Along_Swath_5km_mod06 = 406<br><br>The locations [ Latitude (Latitude_mod06) and Longitude (Longitude_mod06)] are different for each file.<br><br>[4] Your use of<br><br>f=addfiles(fils,"r")<br>ListSetType(f,"cat") ; this is the default<br>cl = f[:]->Cloud_Top_Pressure<br><br>is not appropriate. 'cat' mode is more appropriate when time(s) are changing but the physical locations for each file are the same. In your case, each swath has a different location imprint.<br><br>In your case "join" would be more appropriate.<br><br>f=addfiles(fils,"r")<br>ListSetType(f,"join") ; this is the default<br>cl = f[:]->Cloud_Top_Pressure<br><br>You will see a different look. More like a classic netCDF file. Value added features include adding a 'coordinates' attribute to each variable.<br><br>[3] While the dimension sizes and dimension names are the same for each swath file:<br><br> dimensions:<br> Cell_Across_Swath_5km_mod06 = 270<br> Cell_Along_Swath_5km_mod06 = 406<br><br>The locations [ Latitude (Latitude_mod06) and Longitude (Longitude_mod06)] are different for each file.<br><br>[4] Your use of<br><br>f=addfiles(fils,"r")<br>ListSetType(f,"cat") ; this is the default<br>cl = f[:]->Cloud_Top_Pressure<br><br>is not appropriate. 'cat' mode is more appropriate when time(s) are changing but the physical locations for each file are the same. In your case, each swath has a different location imprint.<br><br>In your case "join" would be more appropriate.<br><br>f=addfiles(fils,"r")<br>ListSetType(f,"join") ; this is the default<br>cl = f[:]->Cloud_Top_Pressure<br><br>printVarSummary(cl) ; ALWAYS look at your data<br><br>Variable: cl<br>Type: short<br>Total Size: 657720 bytes<br> 328860 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [ncl_join | 3] x [Cell_Along_Swath_5km_mod06 | 406] x [Cell_Across_Swath_5km_mod06 | 270]<br> scale_factor : 0.10000000<br> add_offset : 0<br><br>[snip]<br><br>[5] The user must unpack the variable of type short. You could use<br><br><a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml">https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml</a><br><br>======<br>Really, I would not use addfiles for your application. I would loop over each file and overlay each swath file/swath individually. Then, advance the frame.<br><br>See attached<br>======<br><br></div>Good Luck<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Dec 30, 2017 at 2:54 AM, Krishnakumar AP via ncl-talk <span dir="ltr"><<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Dear NCL Community,<br><br></div>This is a re-post of my previous question, regarding distortions appear when merging multiple modis hdf files into a single plot. <br clear="all"><div><div><div>I am using NCL 6.4.0 and data used from Modis Cloud Products(06_L2)</div><div><br></div><div>I have uploaded few of data to the ftp (<a href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a>), file names are <br></div><div>xMOD06_L2.A2011349.0550.006.<wbr>2015056134430.hdf</div><div>xMOD06_L2.A2011349.0730.006.<wbr>2015056133545.hdf</div><div>xMOD06_L2.A2011349.0735.006.<wbr>2015056132747.hdf</div><div><br></div><div>NCL code, <br></div><div><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>begin<br><br>;:::::::::::::::::::::::::::::<wbr>:::::<br>fils=systemfunc("ls ./xMOD06_L2*.hdf")<br>f=addfiles(fils,"r")<br>;:::::::::::::::::::::::::::::<wbr>:::::<br><br>ListSetType(f,"cat")<br>cl = f[:]->Cloud_Top_Pressure<br><br>;:::::::::::::::::::::::::::::<wbr>:::::<br><br>printVarSummary(cl)<br><br>cl@lat2d = f[:]->Latitude<br>cl@lon2d = f[:]->Longitude <br><br>wks=gsn_open_wks("png","1_<wbr>modis_ctp_349") <br>;:::::::::::::::::::::::::::::<wbr>::::::::::::::<br><br>res=True ; plot mods desired<br>res@gsnMaximize =True ; make plot large<br>res@cnLinesOn =False ; turn off contour lines<br>;res@cnLineLabelsOn =True ; turn off contour line labels<br>res@gsnSpreadColors =True ; use the entire color spectrum<br>res@cnFillOn =True ; enable contour fill<br>res@gsnMaximize =True ; make plot large<br>res@<wbr>mpGeophysicalLineThicknessF= 2. ; double the thickness of geophysical boundaries<br>res@cnMissingValFillPattern = 0<br>res@cnFillMode = "RasterFill"<br><br>;::::::ZOOM:::::::::::;<br>res@mpMinLatF = 20<br>res@mpMaxLatF = 32<br>res@mpMinLonF = 68<br>res@mpMaxLonF = 95<br><br>;:::::::::::::Title&Plotting::<wbr>::::::::::::::<br> res@tiMainString="Cloud Top Pressure" ; create title<br> plot=gsn_csm_contour_map_ce(<wbr>wks,cl,res) ; plot on world map first<br> delete(res@tiMainString) ; delete title for subsequent plots<br><br>;:::::::::::::Clean up resources used:::::::::::::::::::::<br> delete(plot) <br> delete(wks)<br> delete(res)<br>end<br></div></blockquote><div> </div><div>Individually each file plotted properly, but when it get merged, distortions or clippings occurred. <br></div><div><br></div><div>As per Rick Brownrigg's suggestion, res@cnFillMode = "RasterFill" added, but could not get any better. Thanks for your kind response. <br></div><div>I have attached the ncl_filedump of three example files, PrintVarSummary and the output.<br></div><div><br></div><div>Please help me on this. Thanks in advance. <br></div><span class="HOEnZb"><font color="#888888"><div></div><div><br></div><div><br></div><div><br>-- <br><div class="m_209394707849348050gmail-m_3908748141857064765gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><br></div><div>Best regards<br>Krishnakumar.AP<br></div><div>Research Fellow <br></div><div>Banaras Hindu University, India.<br></div></div></div></div></div>
</div></font></span></div></div></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>