<div dir="ltr"><div class="gmail_default" style="font-size:small">There are two changes you need to make.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">[1] Your tmp1 longitude values go from -180 to 179.375, but you are trying to subscript it with 180 to 255. To do this, you first need to "flip" your longitudes to go from 0 to 360:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">tmp1 = lonFlip(tmp1)</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">[2] Your "tmp1" variable is on a rectilinear grid, which means you can use coordinate subscripting to subscript it. You need to use curly braces ({ and }) to subscript based on lat/lon:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px"> tmp1 = in1->SWNETTOA(:,{lat_min:lat_</span><wbr style="font-size:12.8px"><span style="font-size:12.8px">max},</span><span style="font-size:12.8px">{l</span><span style="font-size:12.8px">on_min:lon_max</span><span style="font-size:12.8px">}</span><span style="font-size:12.8px">)</span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px"> tmp = in->TMP_P0_L103_GGA0(:,</span><span style="font-size:12.8px">{l</span><span style="font-size:12.8px">at_</span><wbr style="font-size:12.8px"><span style="font-size:12.8px">min:lat_max</span><span style="font-size:12.8px">}</span><span style="font-size:12.8px">,</span><span style="font-size:12.8px">{l</span><span style="font-size:12.8px">on_min:lon_max</span><span style="font-size:12.8px">}</span><span style="font-size:12.8px">)</span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px"><br></span></div><div class="gmail_default"><span style="font-size:12.8px">The way you can tell that tmp1 is on a rectilinear grid is from your printVarSummary output. Under the "Coordinates" section, it listed three arrays, and each one of these is a coordinate array.</span></div><div class="gmail_default"><span style="font-size:12.8px"><br></span></div><div class="gmail_default"><span style="font-size:12.8px"><div class="gmail_default"><font face="monospace, monospace">Variable: tmp1</font></div><div class="gmail_default"><font face="monospace, monospace">Type: float</font></div><div class="gmail_default"><font face="monospace, monospace">Total Size: 19961856 bytes</font></div><div class="gmail_default"><font face="monospace, monospace"> 4990464 values</font></div><div class="gmail_default"><font face="monospace, monospace">Number of Dimensions: 3</font></div><div class="gmail_default"><font face="monospace, monospace">Dimensions and sizes: [time | 24] x [lat | 361] x [lon | 576]</font></div><div class="gmail_default"><font face="monospace, monospace">Coordinates: </font></div><div class="gmail_default"><font face="monospace, monospace"> time: [0..1380] <------ time is a coordinate array</font></div><div class="gmail_default"><font face="monospace, monospace"> lat: [ -90.. 90]</font><span style="font-family:monospace,monospace;font-size:12.8px"> <------ lat is a coordinate array</span></div></span><span style="font-size:12.8px"><div class="gmail_default"><font face="monospace, monospace"> lon: [-180..179.375]</font><span style="font-family:monospace,monospace;font-size:12.8px"> <------ lon is a coordinate array</span></div></span><div class="gmail_default" style="font-size:12.8px"><br></div><div class="gmail_default" style="font-size:12.8px">If you do not have any arrays listed under "Coordinates:", then you will not be able to use the {...} syntax.</div><div class="gmail_default" style="font-size:12.8px"><br></div><div class="gmail_default" style="font-size:12.8px">You didn't include the printVarSummary for your "tmp" variable. If it is also rectilinear, then you can subscript it the same way. If it is not, then there are other methods you can use.</div><div class="gmail_default" style="font-size:12.8px"><br></div><div class="gmail_default" style="font-size:12.8px">Please see this page which describes various methods for subsetting your data across lat/lon, given different kinds of grids:</div><div class="gmail_default" style="font-size:12.8px"><br></div><div class="gmail_default"><span style="font-size:12.8px"><a href="http://www.ncl.ucar.edu/Applications/latlon_subset.shtml">http://www.ncl.ucar.edu/Applications/latlon_subset.shtml</a></span><br></div><div class="gmail_default"><span style="font-size:12.8px"><br></span></div><div class="gmail_default"><span style="font-size:12.8px">--Mary</span></div><div class="gmail_default"><span style="font-size:12.8px"><br></span></div><div class="gmail_default"><span style="font-size:12.8px"><br></span></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 21, 2017 at 3:17 AM, dale zuri <span dir="ltr"><<a href="mailto:dalezuri@gmail.com" target="_blank">dalezuri@gmail.com</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><div><div>Hi, <br></div>Please help me to select the same region from the two different files. <br></div>I want to have the region1 selection for both the files. I am not able to set the same region1 for Merra data. <br></div>Thanks<br><div><div><div><div><br>begin<br>;*****************************<wbr>*******************<br>; create pointer to file and read in temperature<br>;*****************************<wbr>*******************<br> diri = "./"<br> fili1 = "MERRA2_300.tavg1_2d_int_Nx.<wbr>20040926.nc4"<br> fili = "<a href="http://tmp2m.2012050218.time.nc" target="_blank">tmp2m.2012050218.time.nc</a>"<br> in = addfile(diri+fili,"r")<br> in1 = addfile(diri+fili1,"r")<br>print(in1)<br>lat_min =-180;;wrong region selection<br> lat_max =180<br> lon_min = 0<br> lon_max = 138<br>;;;;Region1;;;;;;;;;;;;<br>lat_min =15<br>lat_max =40<br> lon_min = 180 <br> lon_max = 255<br> tmp1 = in1->SWNETTOA(:,lat_min:lat_<wbr>max,lon_min:lon_max)<br> tmp = in->TMP_P0_L103_GGA0(:,lat_<wbr>min:lat_max,lon_min:lon_max)<br>printVarSummary(tmp1)<br>; ts = tmp(time|0:0,lat|:,lon|:) <wbr> ; reorder variable<br>;printVarSummary(ts)<br> ts = dim_avg_n_Wrap(tmp, 0 )<br> ts1 = dim_avg_n_Wrap(tmp1, 0 )</div><div><br></div><div>Variable: tmp1<br>Type: float<br>Total Size: 19961856 bytes<br> 4990464 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 24] x [lat | 361] x [lon | 576]<br>Coordinates: <br> time: [0..1380]<br> lat: [ -90.. 90]<br> lon: [-180..179.375]<br>Number Of Attributes: 11<br> long_name : toa_net_downward_shortwave_<wbr>flux<br> units : W m-2<br> _FillValue : 1e+15<br> missing_value : 1e+15<br> fmissing_value : 1e+15<br> scale_factor : 1<br> add_offset : 0<br> standard_name : toa_net_downward_shortwave_<wbr>flux<br> vmax : 1e+15<br> vmin : -1e+15<br> valid_range : ( -1e+15, 1e+15 )<br></div></div></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>