<html><head><meta http-equiv="Content-Type" content="text/html charset=gb2312"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Thanks for the reply!<div class=""><br class=""></div><div class=""><div class=""><span style="font-family: Menlo;" class="">xbox=x(:,inds_lat,inds_lon)</span></div></div><div class=""><span style="font-family: Menlo;" class=""><br class=""></span></div><div class="">The above line can not get the appropriate results since inds_lat /inds_lon is 2-d array, the first dimension of inds_lat/inds_lon is</div><div class=""><br class=""></div><div class="">the numbers of grids meeting the selection condition in the former.</div><div class=""><br class=""></div><div class="">Now we got to know the grids numbers in a chosen area ( dimsizes(inds_lat(:,0)) ) and the corresponding location in the original </div><div class=""><br class=""></div><div class="">2-d array of x ( the rightmost dimension of inds_lat/inds_lon denote the location ). However, I <span style="font-family: Menlo;" class="">still find no resolution to</span></div><div class=""><span style="font-family: Menlo;" class=""><br class=""></span></div><div class=""><span style="font-family: Menlo;" class="">rearrange these data </span><font face="Menlo" class="">into another 2d array of y(ygrid,xgrid) or plot 1-d array of x1d with lon1d/lat1d directly.</font></div><div class=""><font face="Menlo" class=""><br class=""></font></div><div class=""><font face="Menlo" class=""><br class=""></font></div><div class=""><font face="Menlo" class="">Besides, the original array is x(ygrid=332,xgrid=316), corresponding to 104912 grids in total. </font></div><div class=""><font face="Menlo" class=""> </font></div><div class=""><font face="Menlo" class=""> box_lon = ind(lon1d.ge.min_lon .AND. lon1d.le.max_lon) ; min_lon = -180 and max_lon =180, and box_lon is a array of 104912*2</font></div><div class=""><font face="Menlo" class=""><br class=""></font></div><div class=""><blockquote type="cite" class=""><div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div class=""><div class=""><font face="Menlo" class=""> box_lat = ind(lat1d.ge.min_lat .AND. lat1d.le.max_lat) ; min_lat = -60. max_lat=-50 and box_lat is in 38444*2</font></div></div></div></blockquote></div><div class=""><div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div class=""><div class=""><font face="Menlo" class=""><br class=""></font></div></div></div></div><div class=""><font face="Menlo" class="">I am puzzled with why 38444 is not divided by 316.</font></div><div class=""><font face="Menlo" class=""><br class=""></font></div><div class=""><br class=""></div><div class="">And this is the webpage of projection</div><div class=""><br class=""></div><div class=""><font face="Menlo" class=""><a href="http://nsidc.org/data/polar-stereo/ps_grids.html" class="">http://nsidc.org/data/polar-stereo/ps_grids.html</a></font></div><div class=""><font face="Menlo" class=""><br class=""></font></div><div class=""><font face="Menlo" class=""><br class=""></font></div><div class=""><font face="Menlo" class=""><br class=""></font></div><div class=""><font face="Menlo" class="">Lin </font></div><div class=""><span style="font-family: Menlo;" class=""><br class=""></span></div><div class=""><div class=""><div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div class=""><div class=""><font face="Menlo" class=""><br class=""></font></div></div></div></div><div class=""><div><blockquote type="cite" class=""><div class="">在 2016年8月23日，14:48，Guido Cioni <<a href="mailto:guidocioni@gmail.com" class="">guidocioni@gmail.com</a>> 写道：</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html charset=utf-8" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">You can use ind_resolve to reshape the indices of the box and then use them to subset. I guess <u class="">something like</u> this would work (untested)<div class=""><br class=""></div><div class=""><div class=""><font face="Menlo" style="font-size: 12px;" class=""> dims_lat = dimsizes(latitude)</font></div><div class=""><font face="Menlo" style="font-size: 12px;" class=""> dims_lon = dimsizes(longitude)</font></div><div class=""><font face="Menlo" style="font-size: 12px;" class=""><br class=""></font></div><div class=""><font face="Menlo" style="font-size: 12px;" class=""> lon1d = ndtooned(longitude)<br class=""> lat1d = ndtooned(latitude) </font></div><div class=""><font face="Menlo" style="font-size: 12px;" class=""><br class=""></font></div><div class=""><font face="Menlo" style="font-size: 12px;" class=""> box_lon = ind(lon1d.ge.min_lon .AND. lon1d.le.max_lon )</font></div><div class=""><font face="Menlo" style="font-size: 12px;" class=""> box_lat = ind(lat1d.ge.min_lat .AND. lat1d.le.max_lat)</font></div><div class=""><font face="Menlo" style="font-size: 12px;" class=""><br class=""></font></div><div class=""><font face="Menlo" style="font-size: 12px;" class=""> inds_lon = ind_resolve(box_lon, dims_lon)</font></div><div class=""><font face="Menlo" style="font-size: 12px;" class=""> inds_lat = ind_resolve(box_lat, dims_lat)</font></div><div class=""><font face="Menlo" style="font-size: 12px;" class=""><br class=""></font></div><div class=""><font face="Menlo" style="font-size: 12px;" class=""> xbox=x(:,inds_lat,inds_lon)</font></div><div class=""><br class=""></div><div class="">
<div style="letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-ligatures: normal; font-variant-position: normal; font-variant-caps: normal; font-variant-numeric: normal; font-variant-alternates: normal; font-variant-east-asian: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class=""><br class="Apple-interchange-newline">Guido Cioni</div><div class=""><a href="http://guidocioni.altervista.org/" class="">http://guidocioni.altervista.org</a> </div><div class=""><br class=""></div></div><br class="Apple-interchange-newline" style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-ligatures: normal; font-variant-position: normal; font-variant-caps: normal; font-variant-numeric: normal; font-variant-alternates: normal; font-variant-east-asian: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><br class="Apple-interchange-newline"></div>
</div>
<br class=""><div class=""><blockquote type="cite" class=""><div class="">Il giorno 23 ago 2016, alle ore 04:53, 林祥 <<a href="mailto:xianglin72@icloud.com" class="">xianglin72@icloud.com</a>> ha scritto:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html charset=gb2312" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hey Cioni,<div class=""><br class=""></div><div class=""> Thanks for your reply. The longitude and latitude array is 2-dimensional, so I rewrite the following scripts</div><div class=""><br class=""></div><div class=""> lon1d = ndtooned(longitude)<div class=""> lat1d = ndtooned(latitude)</div><div class=""> x1d = ndtooned(x)</div><div class=""><br class=""></div><div class=""> box_id = ind(lon1d.ge.min_lon .AND. lon1d.le.max_lon .AND. lat1d.ge.min_lat .AND. lat1d.le.max_lat)</div><div class=""> ;box_id = ind(lat1d.ge.min_lat .AND. lat1d.le.max_lat) ; the result is same as the above line</div><div class=""><br class=""></div><div class=""> xlat = lat1d(box_id)</div><div class=""> xlon = lon1d(box_id)</div><div class=""> xbox = x1d(box_id)</div><div class=""><br class=""></div><div class=""> No I have picked out the data in the chosen area. But I am puzzled by plotting the contouring map for xbox.</div><div class="">Should I have to deal with the array of xbox like the separated station datas?</div><div class=""><br class=""></div><div class="">Thanks</div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">Lin</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><blockquote type="cite" class=""><div class="">在 2016年8月23日，00:12，Guido Cioni <<a href="mailto:guidocioni@gmail.com" class="">guidocioni@gmail.com</a>> 写道：</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html charset=utf-8" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hey Lin,<div class="">when dealing with non-monotonic arrays I usually employ the <font face="Menlo" class="">ind</font> function, so I don’t understand why is not working for you…</div><div class=""><br class=""></div><div class="">If you have, say, <font face="Menlo" class="">lon</font> and <font face="Menlo" class="">lat</font> arrays (<b class="">this assume to have them 1-Dimensio</b>nal, otherwise it gets fishy…).</div><div class=""><br class=""></div><div class=""><font face="Menlo" class="">min_lon= -50.</font></div><div class=""><font face="Menlo" class="">max_lon= -70.</font></div><div class=""><font face="Menlo" class="">min_lat= -180.</font></div><div class=""><font face="Menlo" class="">max_lat= -180.</font></div><div class=""><font face="Menlo" class="">box_lon= ind(lon.ge.min_lon .AND. lon.le.max_lon)</font></div><div class=""><div class=""><font face="Menlo" class="">box_lat= ind(lat.ge.min_lat .AND. lat.le.max_lat)</font></div></div><div class=""><br class=""></div><div class="">should give you exactly all the indices inside the user-prescribed area. Then you can subset using</div><div class=""><br class=""></div><div class=""><font face="Menlo" class="">var_box=var(:,box_lat,box_lon)</font></div><div class=""><br class=""></div><div class="">Please be reminded that box_values will have different dimensions, since it is a subset of the original variable.</div><div class="">Cheers,<br class=""><div class="">
<div style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-position: normal; font-variant-caps: normal; font-variant-numeric: normal; font-variant-alternates: normal; font-variant-east-asian: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br class="Apple-interchange-newline">Guido Cioni</div><div style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-position: normal; font-variant-caps: normal; font-variant-numeric: normal; font-variant-alternates: normal; font-variant-east-asian: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><a href="http://guidocioni.altervista.org" class="">http://guidocioni.altervista.org</a> </div>
</div>
<br class=""><div class=""><blockquote type="cite" class=""><div class="">On 22 Aug 2016, at 01:13, 林祥 <<a href="mailto:xianglin72@icloud.com" class="">xianglin72@icloud.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><span id="cid:E18B13C8-D8C7-4DA5-B89B-AE8032041F39@mpi.zmaw.de" class=""><latitude.csv></span></div></blockquote></div><br class=""></div></div></div></blockquote></div><br class=""></div></div></div></blockquote></div><br class=""></div></div></div></blockquote></div><br class=""></div></div></body></html>