search

mtalpe / nctoolbox

Automatically exported from code.google.com/p/nctoolbox
0 stars 0 forks source link

For grid_interop, we should only return [time,z,lat,lon] (not extra time coordinates) #25

Closed GoogleCodeExporter closed 8 years ago

GoogleCodeExporter commented 8 years ago 
For grid_interop,  we should only return [time,z,lat,lon] (not extra time 
coordinates).  In the example below, it would be nice if "time_run" was not 
returned.  (nj_tslice just calls grid_interop)

url=['http://michigan.glin.net:8080/thredds/dodsC/',...
'glos/glcfs/michigan/ncasf_his3d_fc/',...
'Lake_Michigan_Forecast_Collection_3D_best.ncd']
nc = ncgeodataset(url);
[t,g]=nj_tslice(nc,'temp',1,1);
>> g

g = 

    time_run: 3
        time: 7.3450e+005
           z: [251x131 double]
         lat: [251x131 single]
         lon: [251x131 single]

Original issue reported on code.google.com by rsignell on 8 Nov 2011 at 3:16

GoogleCodeExporter commented 8 years ago

Original comment by rsignell on 8 Nov 2011 at 3:20

GoogleCodeExporter commented 8 years ago 
I think I would prefer that the method returns it all, instead of ignoring the 
non cdm defined coordinates. We can change the high level functions to just 
filter out the coordinates wants to return.

Original comment by crosb...@gmail.com on 8 Nov 2011 at 4:12

GoogleCodeExporter commented 8 years ago

Original comment by crosb...@gmail.com on 17 Nov 2011 at 5:21

GoogleCodeExporter commented 8 years ago 
I think this is done?

Original comment by crosb...@gmail.com on 17 Apr 2012 at 6:44

GoogleCodeExporter commented 8 years ago 
Nope, the problem is still there:
url=['http://michigan.glin.net:8080/thredds/dodsC/',...
'glos/glcfs/michigan/ncasf_his3d_fc/',...
'Lake_Michigan_Forecast_Collection_3D_best.ncd']
nc = ncgeodataset(url);
[t,g]=nj_tslice(nc,'temp',1,1);

url =

http://michigan.glin.net:8080/thredds/dodsC/glos/glcfs/michigan/ncasf_his3d_fc/L
ake_Michigan_Forecast_Collection_3D_best.ncd

>> g

g = 

    time_run: 3
        time: 7.3487e+005
           z: [251x131 double]
         lat: [251x131 double]
         lon: [251x131 double]

Original comment by rsignell on 17 Apr 2012 at 6:49

GoogleCodeExporter commented 8 years ago 
No fix intended for grid_interop. We can fix it in the nj functions.

Original comment by crosb...@gmail.com on 18 Apr 2012 at 3:05

GoogleCodeExporter commented 8 years ago 
Okay, please make nj_subsetGrid and nj_tslice return just one time coordinate.

Original comment by rsignell on 18 Apr 2012 at 3:39