<div dir="ltr">Thanks a lot Dennis and NCL community. That was really helpful. Thanks again. <br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 3, 2018 at 7:58 AM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@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>[1]  If you do:<br><br>%> ncl_filedump xMOD06_L2.A2011349.0550.006.<wbr>2015056134430.hdf | less<br><br>You will see:<br><br>[snip]<br>Variable: f<br>Type: file<br>filename:       xMOD06_L2.A2011349.0550.006.<wbr>2015056134430<br>path:   xMOD06_L2.A2011349.0550.006.<wbr>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.<wbr>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")          <wbr>      ; this is the default<span class=""><br>cl = f[:]->Cloud_Top_Pressure<br><br></span>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")         <wbr>       ; this is the default<span class=""><br>cl = f[:]->Cloud_Top_Pressure<br><br></span>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")          <wbr>      ; this is the default<span class=""><br>cl = f[:]->Cloud_Top_Pressure<br><br></span>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")         <wbr>       ; this is the default<span class=""><br>cl = f[:]->Cloud_Top_Pressure<br><br></span>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" target="_blank">https://www.ncl.ucar.edu/<wbr>Document/Functions/<wbr>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"><span class="">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></span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><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.20<wbr>15056134430.hdf</div><div>xMOD06_L2.A2011349.0730.006.20<wbr>15056133545.hdf</div><div>xMOD06_L2.A2011349.0735.006.20<wbr>15056132747.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_modi<wbr>s_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@mpGeophysicalLineThickness<wbr>F= 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(wk<wbr>s,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="m_-535740472223946084HOEnZb"><font color="#888888"><div></div><div><br></div><div><br></div><div><br>-- <br><div class="m_-535740472223946084m_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></div></div><span class="">______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></span></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><br></div><div>Best regards<br>Krishnakumar.AP<br>   </div></div></div></div></div>
</div>