2021-12-08 12:09:49,121 [INFO] Getting data from Cyclobs API. Req : https://cyclobs.ifremer.fr/app/api/getData?include_cols=all&filename=s1b-iw-owi-cm-20191224t213011-20191224t213131-000003-024DE7_ll_gd.nc 2021-12-08 12:09:49,520 [INFO] file : /home/datawork-cersat-public/cache/public/ftp/project/sarwing/processings/c39e79a/default/sentinel-1b/L2C/IW/S1B_IW_OWIH_1S/2019/358/S1B_IW_OWIH_CC_20191224T213011_20191224T213131_019512_024DE7/post_processing/nclight_L2M/s1b-iw-owi-cm-20191224t213011-20191224t213131-000003-024DE7_sw.nc 2021-12-08 12:09:49,786 [INFO] extracting SAR image info for file nclight_L2M... 2021-12-08 12:09:50,200 [INFO] Atcf lon/lat 122.8166666666667/11.75833333333333 2021-12-08 12:10:05,506 [INFO] Distance to coast of track point is 16343.858704168453 meters. 2021-12-08 12:10:05,693 [INFO] Distance between aquisition border and track point is 59969.87502352441 meters. 2021-12-08 12:10:09,274 [INFO] 85.73929234675525% of area of 100km around atcf track is outside acquisition bounding box 2021-12-08 12:10:09,414 [INFO] 0.8827892834445674% of area of 100km around atcf track point that is inside acquisition bounding box is land. 2021-12-08 12:10:09,872 [INFO] 86.62208163019982% of area of 100km around track point is not usable (land or outside of acquisition). 2021-12-08 12:10:09,880 [INFO] SAR image has 78904 non-masked points. /home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_conda/lib/python3.7/site-packages/tc_center/functions.py:811: RuntimeWarning: invalid value encountered in sqrt grad = np.sqrt((np.gradient(wspdcop)[0] ** 2 + np.gradient(wspdcop)[1] ** 2)) 2021-12-08 12:10:13,737 [ERROR] Exception: output array is read-only Traceback (most recent call last): File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_conda/lib/python3.7/site-packages/tc_center/__init__.py", line 231, in find_center_plot_save eye_dict, status = compute_centers(sat_file, resolution, api_data) File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_conda/lib/python3.7/site-packages/tc_center/__init__.py", line 139, in compute_centers atcf_rmw=atcf_rmw) File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_conda/lib/python3.7/site-packages/tc_center/functions.py", line 62, in wrapper result = f(*args, **kwargs) File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_conda/lib/python3.7/site-packages/tc_center/functions.py", line 899, in track_eye_sar var_no_outl[np.where(var >= np.nanquantile(var_no_outl[i_close, j_close], 1. - 200. / np.size(var_no_outl)))] = np.nan File "<__array_function__ internals>", line 6, in nanquantile File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_conda/lib/python3.7/site-packages/numpy/lib/nanfunctions.py", line 1357, in nanquantile a, q, axis, out, overwrite_input, interpolation, keepdims) File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_conda/lib/python3.7/site-packages/numpy/lib/nanfunctions.py", line 1366, in _nanquantile_unchecked return np.nanmean(a, axis, out=out, keepdims=keepdims) File "<__array_function__ internals>", line 6, in nanmean File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_conda/lib/python3.7/site-packages/numpy/lib/nanfunctions.py", line 950, in nanmean avg = _divide_by_count(tot, cnt, out=out) File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_conda/lib/python3.7/site-packages/numpy/lib/nanfunctions.py", line 212, in _divide_by_count return np.divide(a, b, out=a, casting='unsafe') ValueError: output array is read-only