<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body 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>
</body>
</html>