issues

Table actions

1 row where repo = 13221727 and user = 9618134 sorted by updated_at descending

✎ View and edit SQL

This data as json, CSV (advanced)

Suggested facets: created_at (date), updated_at (date), closed_at (date)

type 1

state 1

repo 1

  • xarray · 1 ✖
id node_id number title user state locked assignee milestone comments created_at updated_at ▲ closed_at author_association active_lock_reason draft pull_request body reactions performed_via_github_app state_reason repo type
159222357 MDU6SXNzdWUxNTkyMjIzNTc= 877 Bug in multiplying datasets? Removes coordinates for lats jormelton 9618134 closed 0     4 2016-06-08T17:10:52Z 2016-06-08T18:31:14Z 2016-06-08T18:31:14Z NONE      

I wasn't sure if this should go on the googlegroup or here, but the xarray docs says bugs go here so here it is.

I am doing a very simple operation. I have two netcdf files that I load in and then multiply. First:

``` Python import xarray as xr ncscd='NCSCDv2_Circumpolar_WGS84_SOCC30_05deg.nc' ncscddata30 = xr.open_dataset(ncscd) ncscd_d30cm = ncscddata30['NCSCDv2'].load() print ncscd_d30cm

xarray.DataArray 'NCSCDv2' (lat: 111, lon: 720) array([[ nan, nan, nan, ..., nan, nan, nan], [ nan, nan, nan, ..., nan, nan, nan], [ nan, nan, nan, ..., nan, nan, nan], ..., [ nan, nan, nan, ..., nan, nan, nan], [ nan, nan, nan, ..., nan, nan, nan], [ nan, nan, nan, ..., nan, nan, nan]]) Coordinates: * lon (lon) float64 -179.7 -179.2 -178.7 -178.2 -177.7 -177.2 -176.7 ... * lat (lat) float64 89.75 89.25 88.75 88.25 87.75 87.25 86.75 86.25 ... Attributes: long_name: NCSCDv2 esri_pe_string: GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]] units: Degree ```

And then the second:

``` Python gridfile ='NCSCDv2_gridarea.nc' gridarea = xr.open_dataset(gridfile).load() print gridarea

xarray.Dataset Dimensions: (lat: 111, lon: 720) Coordinates: * lon (lon) float64 -179.7 -179.2 -178.7 -178.2 -177.7 -177.2 ... * lat (lat) float64 89.75 89.25 88.75 88.25 87.75 87.25 86.75 86.25 ... Data variables: cell_area (lat, lon) float64 1.337e+07 1.337e+07 1.337e+07 1.337e+07 ... Attributes: CDI: Climate Data Interface version 1.6.8 (http://mpimet.mpg.de/cdi) Conventions: CF-1.4 history: Mon May 16 18:13:43 2016: cdo gridarea NCSCDv2_Circumpolar_WGS84_SOCC100_05deg.nc NCSCDv2_gridarea.nc CDO: Climate Data Operators version 1.6.8rc2 (http://mpimet.mpg.de/cdo) ```

Both arrays are the same size. Now I simply multiply them together but note that lats in the output are now 0.

``` Python ncscd_total30cm = gridarea.cell_area * ncscd_d30cm print ncscd_total30cm

xarray.Dataset Dimensions: (lat: 0, lon: 720) Coordinates: * lat (lat) float64 * lon (lon) float64 -179.7 -179.2 -178.7 -178.2 -177.7 -177.2 ... Data variables: cell_area (lat, lon) float64 Attributes: CDI: Climate Data Interface version 1.6.8 (http://mpimet.mpg.de/cdi) Conventions: CF-1.4 history: Mon May 16 18:13:43 2016: cdo gridarea NCSCDv2_Circumpolar_WGS84_SOCC100_05deg.nc NCSCDv2_gridarea.nc CDO: Climate Data Operators version 1.6.8rc2 (http://mpimet.mpg.de/cdo) ```

This doesn't seem like it should be happening. Removing the .load() doesn't help. Bug?

 
{
    "url": "https://api.github.com/repos/pydata/xarray/issues/877/reactions",
    "total_count": 0,
    "+1": 0,
    "-1": 0,
    "laugh": 0,
    "hooray": 0,
    "confused": 0,
    "heart": 0,
    "rocket": 0,
    "eyes": 0
}
   completed xarray 13221727 issue

Advanced export

JSON shape: default, array, newline-delimited, object

CSV options: 

CREATE TABLE [issues] (
   [id] INTEGER PRIMARY KEY,
   [node_id] TEXT,
   [number] INTEGER,
   [title] TEXT,
   [user] INTEGER REFERENCES [users]([id]),
   [state] TEXT,
   [locked] INTEGER,
   [assignee] INTEGER REFERENCES [users]([id]),
   [milestone] INTEGER REFERENCES [milestones]([id]),
   [comments] INTEGER,
   [created_at] TEXT,
   [updated_at] TEXT,
   [closed_at] TEXT,
   [author_association] TEXT,
   [active_lock_reason] TEXT,
   [draft] INTEGER,
   [pull_request] TEXT,
   [body] TEXT,
   [reactions] TEXT,
   [performed_via_github_app] TEXT,
   [state_reason] TEXT,
   [repo] INTEGER REFERENCES [repos]([id]),
   [type] TEXT
);
CREATE INDEX [idx_issues_repo]
    ON [issues] ([repo]);
CREATE INDEX [idx_issues_milestone]
    ON [issues] ([milestone]);
CREATE INDEX [idx_issues_assignee]
    ON [issues] ([assignee]);
CREATE INDEX [idx_issues_user]
    ON [issues] ([user]);