<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-&gt;U<br>
              u   = lonFlip(u)<br>
          <br>
          printVarSummary(u) <br>
          <br>
          print("=====================Process=========================")<br>
          <br>
            U =
          linint2(u&amp;longitude,u&amp;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&amp;lat = xlat<br>
            U&amp;lon = xlon<br>
            U&amp;time  = u&amp;time<br>
            U&amp;level = u&amp;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>