Hello NCL folks,
I am comparing two satellite datasets with different resolutions. Data1 is
a 5x5 global grid. Data2 is .25x.25 degree between 50S and 50N. I am trying
to use the lat/lons of Data1 (the lower resolution dataset) to select the
nearest gridpoint in Data2 (the hi-res dataset).* I keep getting the
following error:
fatal:NclOneDValGetClosestIndex: finish coordinate index out of range,
can't continue
I tried to reduce the latitudinal extent of Data1 to be similar to or even
smaller than that of Data2, but I continue to get the error.
What am I doing wrong? Code snippet and variable info below.
-- Variable: Data1 Type: float Total Size: 653184 bytes 163296 values Number of Dimensions: 3 Dimensions and sizes: [time | 63] x [lat | 36] x [lon | 72] Coordinates: time: [4291..6268] lat: [-87.5..87.5] lon: [-180..175] Number Of Attributes: 2 average_op_ncl : dim_avg_n over dimension(s): lev _FillValue : -999 Variable: Data2 Type: float Total Size: 145152000 bytes 36288000 values Number of Dimensions: 3 Dimensions and sizes: [time | 63] x [lat | 400] x [lon | 1440] Coordinates: time: [2004.666666666667..2009.916666666667] lat: [-49.875..49.875] lon: [-179.875..179.875] Number Of Attributes: 5 _FillValue : -9999.9 missing_value : -9999.9 time_statistic : instantaneous units : mm/hr long_name : precipitation Code snippet: ; ---- Use Data1 (5x5degree) lat/lons to select nearest gridpoints in Data2 (.25x.25degree) ---- latsub1 = -55. ; <-- Still get the error if I choose a smaller number like -49 latsub2 = 55. ; <-- Still get the error if I choose a smaller number like 49 Data1latsub = Data1(:,{latsub1:latsub2},:) data1lat = Data1latsub&lat data1lon = Data1latsub&lon printVarSummary(data1lat) printVarSummary(data1lon) printMinMax(data1lat,True) printMinMax(data1lon,True) printVarSummary(Data2&lat) printVarSummary(Data2&lon) printMinMax(Data2&lat,True) printMinMax(Data2&lon,True) Data2sub = Data2(:,{data1lat},{data1lon}) printVarSummary(Data2sub) exit Output: Variable: data1lat Type: float Total Size: 88 bytes 22 values Number of Dimensions: 1 Dimensions and sizes: [lat | 22] Coordinates: Number Of Attributes: 2 long_name : Latitude units : degrees_north Variable: data1lon Type: float Total Size: 288 bytes 72 values Number of Dimensions: 1 Dimensions and sizes: [lon | 72] Coordinates: Number Of Attributes: 2 units : degrees_east long_name : Longitude (0) (0) Latitude: min=-52.5 max=52.5 (0) (0) Longitude: min=-180 max=175 Variable: lat (coordinate) Type: double Total Size: 3200 bytes 400 values Number of Dimensions: 1 Dimensions and sizes: [lat | 400] Coordinates: Number Of Attributes: 2 long_name : Latitude units : degrees_north Variable: lon (coordinate) Type: double Total Size: 11520 bytes 1440 values Number of Dimensions: 1 Dimensions and sizes: [lon | 1440] Coordinates: Number Of Attributes: 2 long_name : Longitude units : degrees_east (0) (0) Latitude: min=-49.875 max=49.875 (0) (0) Longitude: min=-179.875 max=179.875 fatal:NclOneDValGetClosestIndex: finish coordinate index out of range, can't continue fatal:["Execute.c":8128]:Execute: Error occurred at or near line 160 in file data_corr.ncl -- Many thanks, Bronwen *P.S. A better approach would be to interpolate the hi-res Data2 dataset to match the lat/lons of the low-res Data1 dataset. I looked through the various regridding examples and none seemed appropriate. If anyone has suggestions on the best way to do this, I'd appreciate it! -- Bronwen Konecky, Ph.D. Postdoctoral Fellow School of Earth & Atmospheric Sciences Georgia Institute of Technology 311 Ferst Dr., Atlanta, GA, 30332-0340 (303) 492-5176 http://blkonecky.wordpress.com/
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Mar 17 15:38:08 2014
This archive was generated by hypermail 2.1.8 : Fri Mar 21 2014 - 15:49:21 MDT