<div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Hi Vanucia,</div><div><br></div><div>You are doing nothing "wrong."</div><div><br></div><div>The '<b>Golden Rule</b>' data processing is '<b>Know Your Data</b>'</div><div><br></div><div>A very useful function to explore a variable is <b> <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml">stat_dispersion</a></b><br></div></div><div dir="ltr"><br></div><div>See attached: <br></div><div><br></div><div>%> ncl <b>MOD11A1.ncl_explore</b><br></div><div><br></div><div>Not only will you get a feel for the data, the output can guide your contour settings.</div><div><br></div><div>%> ncl <b>MOD11A1.ncl</b></div><div><b><br></b></div><div>Good luck<b><br></b></div><div><b>D<br></b></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Mar 2, 2019 at 2:22 PM Vanúcia Schumacher <<a href="mailto:vanucia-schumacher@hotmail.com">vanucia-schumacher@hotmail.com</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,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
I'm trying to read MODIS data for daily mean temperature, but I think I'm doing something wrong because Maximo value is going up to 336 K. Any tips if I am reading correctly. Thanks<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
SCRIPT:</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span>;---HDF file name, also with HDFEOS suffix<br>
</span>
<div> fname_hdf = systemfunc("ls MOD11A1.A2*.hdf")<br>
</div>
<div> fname_he1 = "MOD11A1.A2015359.mosaic.006.2019043165512.psmcrpgscs_000501303977.LST_Day_1km.hdf"<br>
</div>
<div> fname_he2 = fname_he1 + ".he2"<br>
</div>
<div><br>
</div>
<div>;---Read data to be plotting as HDF lst_dayiable (not HE2).<br>
</div>
<div> a = addfiles(fname_hdf,"r")<br>
</div>
<div> ListSetType(a,"join") <br>
</div>
<div> lst_day = short2flt(a[:]->LST_Day_1km)<br>
</div>
<div><br>
</div>
<div>;---Change the _FillValue, which is equal to 0.0 (NCL doesn't like 0.0 as a missing value)<br>
</div>
<div> lst_day@_FillValue = default_fillvalue(typeof(lst_day))<br>
</div>
<div><br>
</div>
<div>;---Read lat/lon grid as HDFEOS lst_day<br>
</div>
<div> agrid = addfile(fname_he2,"r")<br>
</div>
<div> lat = agrid->GridLat_MODIS_Grid_Daily_1km_LST<br>
</div>
<div> lon = agrid->GridLon_MODIS_Grid_Daily_1km_LST<br>
</div>
<div><br>
</div>
<div>lst_day!0 = "time"<br>
</div>
<div>lst_day!1 = "lat"<br>
</div>
<div>lst_day!2 = "lon"</div>
<div>lst_day@lat = lat<br>
</div>
<div>lst_day@lon = lon<br>
</div>
<span></span><br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span>Variable: lst_day<br>
</span>
<div>Type: float<br>
</div>
<div>Total Size: 6896628 bytes<br>
</div>
<div> 1724157 values<br>
</div>
<div>Number of Dimensions: 3<br>
</div>
<div>Dimensions and sizes: [ncl_join | 3] x [YDim_MODIS_Grid_Daily_1km_LST | 1003] x [XDim_MODIS_Grid_Daily_1km_LST | 573]<br>
</div>
<div>Coordinates: <br>
</div>
<div>Number Of Attributes: 12<br>
</div>
<div> NCL_converted_from_type : ushort<br>
</div>
<div> hdf_name : LST_Day_1km<br>
</div>
<div> calibrated_nt : 5<br>
</div>
<div> add_offset_err : 0<br>
</div>
<div> scale_factor_err : 0<br>
</div>
<div> LST : LST data * scale_factor<br>
</div>
<div> Number_Type : uint16<br>
</div>
<div> units : K<br>
</div>
<div> long_name : Daily daytime 1km grid Land-surface Temperature<br>
</div>
<div> _FillValue_original : 0<br>
</div>
<div> _FillValue : 9.96921e+36<br>
</div>
<div> valid_range : ( 150, -0.02 )<br>
</div>
<div>(0) Daily daytime 1km grid Land-surface Temperature (K) : <span style="color:rgb(200,38,19)">
min=270.34 max=336.6</span><br>
</div>
<div>(0) latitude (degrees_north) : min=-40.03057991847813 max=-30.96683235691647<br>
</div>
<div>(0) longitude (degrees_east) : min=-72.35461066157526 max=-65.65836999384004<br>
</div>
<span></span><br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<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>