<div dir="ltr"><div><div><div>Hi Karin,<br></div>What I wanted to do is to interpolate from 144x288 to 73x144.<br></div>It works, but the edges of the data are missing esp. in the north and south poles.<br></div>Thanks<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Apr 2, 2015 at 12:56 PM, Karin Meier-Fleischer <span dir="ltr"><<a href="mailto:meier-fleischer@dkrz.de" target="_blank">meier-fleischer@dkrz.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
Hi Antonio,<br>
<br>
I think you should have a look at you lines<span class=""><br>
<br>
<tt> xlat = fspan(-90.,90,73)</tt><tt><br>
</tt><tt> xlon = fspan(0.,357.5,144)</tt><br>
<br></span>
you're mixing latitude and longitude and/or use wrong dimension
sizes.<br>
<br>
<font color="#3333ff"><tt> xlat = fspan(-90.,90,144)</tt><tt><br>
</tt><tt> xlon = fspan(0.,357.5,288)</tt></font><br>
<br>
Bye,<br>
Karin<br>
<br>
<div>Am 02.04.15 um 12:44 schrieb antonio
avio:<br>
</div>
<blockquote type="cite"><div><div class="h5">
<div dir="ltr">
<div>
<div>
<div>Hallo NCL,<br>
</div>
Im just wondering why the output of the interpolation filed
are missing in the last latitude (south and north poles:
87.5-90). The input data is bellow:<br>
<br>
netcdf MERRA.U.1979.plev {<br>
dimensions:<br>
time = 365 ;<br>
levels = 42 ;<br>
latitude = 144 ;<br>
longitude = 288 ;<br>
variables:<br>
float U(time, levels, latitude, longitude) ;<br>
U:long_name = "Eastward wind component" ;<br>
U:_FillValue = 1.e+15f ;<br>
U:comments = "Unknown1 variable comment" ;<br>
U:grid_name = "grid-1" ;<br>
U:time_statistic = "average" ;<br>
double time(time) ;<br>
time:standard_name = "time" ;<br>
time:units = "hours since 1979-01-01
12:00:00" ;<br>
time:calendar = "standard" ;<br>
double levels(levels) ;<br>
levels:standard_name = "air_pressure" ;<br>
levels:long_name = "Level" ;<br>
levels:units = "hPa" ;<br>
levels:positive = "down" ;<br>
levels:axis = "Z" ;<br>
double latitude(latitude) ;<br>
latitude:standard_name = "latitude" ;<br>
latitude:long_name = "Latitude" ;<br>
latitude:units = "degrees_north" ;<br>
latitude:axis = "Y" ;<br>
double longitude(longitude) ;<br>
longitude:standard_name = "longitude" ;<br>
longitude:long_name = "Longitude" ;<br>
longitude:units = "degrees_east" ;<br>
longitude:axis = "X" ;<br>
<br>
<br>
</div>
My script is:<br>
<br>
oad "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" <br>
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" <br>
<br>
begin<br>
<br>
;===========2.5x2.5=============<br>
<br>
xlat = fspan(-90.,90,73)<br>
xlon = fspan(0.,357.5,144)<br>
xlat@units = "degrees_north"<br>
xlon@units = "degrees_east"<br>
<br>
;===============================<br>
<br>
do y =1979,1979<br>
print(y)<br>
<br>
f = addfile ("MERRA.U."+y+".<a href="http://plev.nc" target="_blank">plev.nc</a>" , "r")<br>
<br>
u = f->U<br>
u = lonFlip(u)<br>
<br>
printVarSummary(u) <br>
<br>
print("=====================Process=========================")<br>
<br>
U =
linint2(u&longitude,u&latitude,u,True,xlon,xlat,0)<br>
<br>
<br>
U!0 = "time"<br>
U!1 = "level"<br>
U!2 = "lat"<br>
U!3 = "lon"<br>
U&lat = xlat<br>
U&lon = xlon<br>
U&time = u&time<br>
U&level = u&levels<br>
<br>
</div>
<div>printVarSummary(U)<br>
<br>
</div>
Thank you<br>
</div>
<br>
<fieldset></fieldset>
<br>
</div></div><pre>_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a>
</pre>
</blockquote>
<br>
<pre cols="72">--
Dipl. Geophys. Karin Meier-Fleischer
Visualization
Application Support
Deutsches Klimarechenzentrum GmbH (DKRZ)
Bundesstrasse 45a - D20146 Hamburg - Germany
Phone: <a href="tel:%2B49%20%280%2940%20460094%20126" value="+4940460094126" target="_blank">+49 (0)40 460094 126</a>
Fax: <a href="tel:%2B49%20%280%2940%20460094%20270" value="+4940460094270" target="_blank">+49 (0)40 460094 270</a>
E-Mail: <a href="mailto:meier-fleischer@dkrz.de" target="_blank">meier-fleischer@dkrz.de</a>
URL: <a href="http://www.dkrz.de" target="_blank">www.dkrz.de</a>
Geschäftsführer: Prof. Dr. Thomas Ludwig
Sitz der Gesellschaft: Hamburg
Amtsgericht Hamburg HRB 39784
</pre>
</div>
<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>
<br></blockquote></div><br></div>