Half of Data is FillValue after regridding (linint2_Wrap)

Hi NCL,

I try to regrid Temperature data with a simple linint2_Wrap.
Everything is fine except that in the regrid-plot the half of the data is
filled with the FillValue (uncomment print(newT) in code).
Have anyone an idea what's wrong?

If you compare the plots of original and regridded data the regridded are
obviously moved to the east around half the globe.

Thanks a lot
Chris

See Code:
---------------------------------------------------------------------------------------------------------------------
  in = addfile("T2m.nc", "r") ;
http://www.ncl.ucar.edu/Applications/Data/cdf/T2m.nc
  Ti = in->T(0,:,:)
  print("all data missing Ti? " + all(ismissing(Ti)))
  printVarSummary(Ti)
  ;print(Ti)
  Tv = lonFlip(Ti)
  printVarSummary(Tv)
  print("all data missing Tv? " + all(ismissing(Tv)))
  print("min Tv: " + min(Tv))
  print("max Tv: " + max(Tv))
  ;print(Tv)
  lonv = in->lon
  latv = in->lat
;************************************************
; interpolate to new grid
;***********************************************
  newlat = fspan(-90,90,37)
  ;printVarSummary(newlat)
  ;print(newlat)
  newlon = fspan(-180,180,73)
  ;printVarSummary(newlon)
  ;print(newlon)
  newlat@units = "degrees_north"
  newlon@units = "degrees_east"

  newT = linint2_Wrap(lonv,latv,Tv,False,newlon,newlat,0)

  printVarSummary(newT)
  print("all data missing newT? " + all(ismissing(newT)))
  ;print(newT)
  print("min newT: " + min(newT))
  print("max newT: " + max(newT))

  ;newT!0 = "time" ; must set
  newT!0 = "lat"
  newT!1 = "lon"
  newT&lat = newlat
  newT&lon = newlon

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Jun 16 14:24:21 2012