[ncl-talk] lonFlip function

[Sunmin Park]

Thank you Adam

I am trying to do regredding the data using the function linint2_Wrap. 
I got warning 
"warning:linint2: xi, yi, xo, and yo must be monotonically increasing”
therefore I think it is due to LON values. 
How I can remove the warning? 


  diri    = "/Users/spark/vertical/00_DATA/"
  fili    = "precipmonmean.nc"
  f       = addfile (diri+fili+".nc", "r")
  varname = "precip"
  data = f->$varname$

  LON     = f->lon ({0:360})  ;
  LAT     = f->lat ;({-90:90})  ;
;  LON     = lonFlip(LON1)

  printVarSummary(data)
  data1 = lonFlip(data)
  printVarSummary(data1)

  lat = fspan(-85,85,85)
  lon = fspan(-180,180,72)
  lat at units = "degrees_east"
  lon at units = "degrees_north"
  precipreg = linint2_Wrap(LON,LAT,data1,True,lon,lat,0)
printVarSummary(precipreg)

On Feb 20, 2015, at 7:17 AM, Adam Phillips <asphilli at ucar.edu> wrote:

> Hi Sunmin,
> lonFlip flips the attached longitude coordinate variable, and not the array itself. Try:
> printVarSummary(data)
> data = lonFlip(data) 
> printVarSummary(data)
> If you have any further questions please respond to ncl-talk.
> Adam
> 
> On Fri, Feb 20, 2015 at 8:04 AM, Sunmin Park <mireiyue at gmail.com> wrote:
> Dear NCL users
> 
> I am trying to change longitude 0,360 to -180,180. I used the function lonFlip to change longitude but it has wrong values. When I use printVarSummary the coordinate looks fine but when I print it out (print) it doesn’t change LON (the file that I am using is GPCP monthly mean precipitation from NCAR) How can I get correct LON (-180,180)?
> 
> Thank you in advance,
> Sun-
> 
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> 
> begin
> 
> ;==============================================================
> ; Open the file:
> ;================================================================
> 
>   diri    = "/Users/spark/vertical/00_DATA/"
>   fili    = "precipmonmean.nc"
>   f       = addfile (diri+fili+".nc", "r")
>   varname = "precip"
>   data = f->$varname$
> 
>   LON1     = f->lon ;({-180:180})  ;
>   LAT     = f->lat ;({-90:90})  ;
>   LON     = lonFlip(LON1)
> 
>   printVarSummary(data)
>   printVarSummary(LON)
>   print(LON)
> return
> 
> 
> Variable: LON
> Type: float
> Total Size: 576 bytes
>             144 values
> Number of Dimensions: 1
> Dimensions and sizes:   [lon | 144]
> Coordinates:
>             lon: [-178.75..178.75]
> Number Of Attributes: 6
>   units :       degrees_east
>   long_name :   Longitude
>   actual_range :        ( 1.25, 358.75 )
>   standard_name :       longitude
>   axis :        X
>   lonFlip :     longitude coordinate variable has been reordered via lonFlip
> 
> (0)     181.25
> (1)     183.75
> (2)     186.25
> (3)     188.75
> (4)     191.25
> (5)     193.75
> (6)     196.25
> (7)     198.75
> (8)     201.25
> (9)     203.75
> (10)    206.25
> .
> .
> .
> (135)   158.75
> (136)   161.25
> (137)   163.75
> (138)   166.25
> (139)   168.75
> (140)   171.25
> (141)   173.75
> (142)   176.25
> (143)   178.75
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> 
> 
> 
> -- 
> Adam Phillips 
> Associate Scientist,  Climate and Global Dynamics Division, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150220/98672d68/attachment.html