icepyx icon indicating copy to clipboard operation
icepyx copied to clipboard

Can't download ATL19 / 23 `south_polar` ocean

Open bbuzz31 opened this issue 10 months ago • 2 comments

icepyx doesn't seem to let me download for south_polar region of ATL19 and ATL23 grids, although it works fine for mid_latitude and north_polar.

I'm constructing the query essentially like this:

region_a = ipx.Query(ATL23, [-180,-90,180,90], date_range=['2018-10-01', '2025-04-01'])
region_a.subsetparams(Coverage=region_a.order_vars.wanted)

![Image](https://github.com/user-attachments/assets/e35b93c1-507d-45ff-ab56-e73987baf39b)

region_a.download_granules(self.path_reg)


atl19 = 'longitude latitude delta_time_beg delta_time_end gridcntr_lat gridcntr_lon landmask '\
    'dot_avg_albm dot_avg_uncrtn_albm dot_avgcntr dot_avgcntr_uncrtn ice_conc_albm n_segs_albm depth_avg_albm ' \
    'dot_dfw_albm dot_dfw_uncrtn_albm dot_dfwcntr ssb_avg_albm ssb_avgcntr ssb_dfw_albm ssb_dfwcntr'.split()

dct_prod_var_list = {
                'ATL19':
                {'mid_latitude': atl19,
                'north_polar': atl19[4:] + 'ds_grid_y ds_grid_x'.split(),
                'south_polar': atl19[4:] + 'ds_grid_y ds_grid_x'.split(),
                },
}
dct_prod_var_list['ATL23'] = dct_prod_var_list['ATL19'] # i think they're the same

dct_var_list = dct_prod_var_list[self.prod]
for k, v in dct_var_list.items():
    region_a.order_vars.append(var_list=v, keyword_list=[k])

I've attached a screen shot of a resulting product that's downloaded.

I've also played around with the parameters to subset, especially since there is some weird behavior... if i pass delta_time_*** to north_polar i get an empty group in the resulting hdf5. the same thing happens for south_polar (but at least the group exists then)

Thanks!

bbuzz31 avatar Mar 15 '25 19:03 bbuzz31

Hello @bbuzz31! Thanks for sharing your problem and some code. Could you share exactly what error you're getting? Then I can try and recreate the issue later this week.

JessicaS11 avatar Mar 17 '25 22:03 JessicaS11

Hi @JessicaS11, sorry i forgot to attach the screen shot. But there is no error; the screen shot just shows the south_polar group is entirely missing in the downloaded product

Image

bbuzz31 avatar Mar 18 '25 01:03 bbuzz31

@bbuzz31 I wanted to follow up and see if this was still an issue for you in v2.0.0. If not, feel free to close this issue. Thanks!

JessicaS11 avatar Aug 04 '25 16:08 JessicaS11

@JessicaS11 I can't tell; #690 completely prevents icepyx from working.

bbuzz31 avatar Aug 06 '25 17:08 bbuzz31

@bbuzz31 thanks for your patience - we'll be trying to put through a release with the fix for #690 (making it a valid link) early next week. In the meantime, the update is available on the development branch. Since unfortunately variable subsetting is no longer available, I'm curious if this will continue to be an issue.

JessicaS11 avatar Nov 14 '25 16:11 JessicaS11

HI @JessicaS11 thanks for keeping me updated. I tried on the development branch, but i can't .download_granules:

--> 469 
region_a.download_granules(self.path_reg)
    471 ## create nc
    472 combine_ATL(self.path_reg, self.prod, year)

File ~/Software_SL/icepyx_git/icepyx/core/query.py:1024, in Query.download_granules(self, path, overwrite)
   1022 if hasattr(self, "last_order") is None:
   1023     raise ValueError("No order has been placed yet.")
-> 1024 status = self.last_order.status()
   1026 if status["status"] == "running" or status["status"] == "accepted":

AttributeError: 'Query' object has no attribute 'last_order'
> /home/buzzanga/Software_SL/icepyx_git/icepyx/core/query.py(1024)download_granules()
   1022         if hasattr(self, "last_order") is None:
   1023             raise ValueError("No order has been placed yet.")
-> 1024         status = self.last_order.status()

bbuzz31 avatar Dec 01 '25 01:12 bbuzz31