<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi Antonio,<br>
<br>
I think you should have a look at you lines<br>
<br>
<tt> xlat = fspan(-90.,90,73)</tt><tt><br>
</tt><tt> xlon = fspan(0.,357.5,144)</tt><br>
<br>
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 class="moz-cite-prefix">Am 02.04.15 um 12:44 schrieb antonio
avio:<br>
</div>
<blockquote
cite="mid:CACKeJX+CKsmV4C5PawbLht6N3pHQCXMwOxoOVg7tyx28uhmbfg@mail.gmail.com"
type="cite">
<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 moz-do-not-send="true"
href="http://plev.nc">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 class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
<a class="moz-txt-link-freetext" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a>
</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
Dipl. Geophys. Karin Meier-Fleischer
Visualization
Application Support
Deutsches Klimarechenzentrum GmbH (DKRZ)
Bundesstrasse 45a - D20146 Hamburg - Germany
Phone: +49 (0)40 460094 126
Fax: +49 (0)40 460094 270
E-Mail: <a class="moz-txt-link-abbreviated" href="mailto:meier-fleischer@dkrz.de">meier-fleischer@dkrz.de</a>
URL: <a class="moz-txt-link-abbreviated" href="http://www.dkrz.de">www.dkrz.de</a>
Geschäftsführer: Prof. Dr. Thomas Ludwig
Sitz der Gesellschaft: Hamburg
Amtsgericht Hamburg HRB 39784
</pre>
</body>
</html>