<div dir="ltr"><div><div><div>[1] area_hi2lores is deprecated. It was originally developed for 30S-30N hi-res precipitation. It should not be used for high latitudes.<br><br></div><div> It would be best to use the ESMF. The first 4 examples illustrate how to interpolate from a hi-res to a low-res grid.<br></div><div> For a smoothly varying variable like temperature, any of the interpolation methods is adequate. However, I would suggest<br></div><div> the conservative method.<br></div><div><br> <a href="https://www.ncl.ucar.edu/Applications/ESMF.shtml">https://www.ncl.ucar.edu/Applications/ESMF.shtml</a><br></div><div><br></div>[2] Why don't you just generate your own ?<br></div><div> You can change the 'LON' to whatever you want.<br></div> <br>LON = fspan(-180,179.75,1440)<br>LON!0 = "LON"<br>LON@units = "degrees_east"<br>LON&LON = LON<br><br>printVarSummary(LON)<br>---<br>Variable: LON<br>Type: float<br>Total Size: 5760 bytes<br> 1440 values<br>Number of Dimensions: 1<br>Dimensions and sizes: [LON | 1440]<br>Coordinates: <br> LON: [-180..179.75]<br>Number Of Attributes: 1<br> units : degrees_east<br>---<br><br></div>Good luck<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 18, 2015 at 11:50 AM, Agnes Lim Huei Ni <span dir="ltr"><<a href="mailto:alim@ssec.wisc.edu" target="_blank">alim@ssec.wisc.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I wish to use area_hi2lores to regrid my high resolution data to 0.25<br>
degree.<br>
<br>
Variable: reorder_int_nr_T<br>
Type: float<br>
Total Size: 1725834240 bytes<br>
431458560 values<br>
Number of Dimensions: 3<br>
Dimensions and sizes: [lv_ISBL3 | 26] x [lat | 2881] x [lon | 5760]<br>
Coordinates:<br>
lat: [ -90.. 90]<br>
lon: [-180..179.9375]<br>
Number Of Attributes: 2<br>
units : K<br>
_FillValue : 1e+15<br>
<br>
Variable: lon_pgbq<br>
Type: float<br>
Total Size: 5760 bytes<br>
1440 values<br>
Number of Dimensions: 1<br>
Dimensions and sizes: [g0_lon_1 | 1440]<br>
Coordinates:<br>
g0_lon_1: [ 0..359.75]<br>
Number Of Attributes: 9<br>
long_name : longitude<br>
GridType : Cylindrical Equidistant Projection Grid<br>
units : degrees_east<br>
Dj : 0.25<br>
Di : 0.25<br>
Lo2 : 359.75<br>
La2 : -90<br>
Lo1 : 0<br>
La1 : 90<br>
<br>
I can't use lonFlip to reorder lon_pgbq. So in this case i use the<br>
following.<br>
<br>
dim=dimsizes(lon_pgbq)<br>
do i=1, dim(0)-1<br>
if(lon_pgbq(i) .gt. 180.0) then<br>
lon+pgbq(i)=lon_pgbq(i)-360<br>
end if<br>
end do<br>
lon_pgbq_1=new((/dim(0)/), "float")<br>
lon_pgbq_1(0:718)=lon_pgbq(721:1439)<br>
lon_pgbq_1(719:1439)=lon_pgbq(0:720)<br>
delete(dim)<br>
<br>
regrid_int_nr_T=area_hi2lores(lon,lat,reorder_int_nr_T,True,1,lon_pgbq_1,lat_pgbq,False<br>
<br>
Variable: regrid_int_nr_T<br>
Type: float<br>
Total Size: 107976960 bytes<br>
26994240 values<br>
Number of Dimensions: 3<br>
Dimensions and sizes: [lv_ISBL3 | 26] x [lat | 721] x [lon | 1440]<br>
Coordinates:<br>
lat: [-90..90]<br>
lon: [-179.75..180]<br>
Number Of Attributes: 1<br>
_FillValue : 1e+15<br>
<br>
<br>
Is there a cleaner way to do it?<br>
<br>
Thanks<br>
Agnes<br>
_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div><br></div>