2023-01-23 14:57:40,629 [INFO] Getting data from Cyclobs API. Req : https://cyclobs.ifremer.fr/app/api/getData?include_cols=all&filename=20210420t094117-20210420t094300 2023-01-23 14:57:41,122 [INFO] Getting IBTRACS data from Cyclobs API. Req : https://cyclobs.ifremer.fr/app/api/getData?include_cols=all&filename=20210420t094117-20210420t094300&track_source=ibtracs 2023-01-23 14:57:41,599 [INFO] file : /home/datawork-cersat-public/cache/public/ftp/project/sarwing/processings/c39e79a/default/sentinel-1b/L2C/IW/S1B_IW_OWIH_1S/2021/110/S1B_IW_OWIH_CC_20210420T094117_20210420T094300_026549_032B99/post_processing/nclight_L2M/s1b-iw-owi-cm-20210420t094117-20210420t094300-000003-032B99_sw.nc 2023-01-23 14:57:41,601 [INFO] extracting SAR image info for file nclight_L2M... /home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_37/lib/python3.7/site-packages/tc_center/functions.py:790: RuntimeWarning: invalid value encountered in sqrt grad = np.sqrt((np.gradient(wspdcop)[0] ** 2 + np.gradient(wspdcop)[1] ** 2)) 2023-01-23 14:57:44,348 [INFO] Atcf lon/lat 125.9375/16.175 2023-01-23 14:57:49,587 [INFO] Distance to coast of track point is 296058.4630319893 meters. 2023-01-23 14:57:49,590 [INFO] Distance between aquisition border and track point is 27749.064596860604 meters. 2023-01-23 14:57:49,619 [INFO] 32.571608012674176% of area of 100km around atcf track is outside acquisition bounding box 2023-01-23 14:57:49,645 [INFO] 0.0% of area of 100km around atcf track point that is inside acquisition bounding box is land. 2023-01-23 14:57:49,647 [INFO] 32.571608012674176% of area of 100km around track point is not usable (land or outside of acquisition). 2023-01-23 14:57:49,653 [INFO] SAR image has 170291 non-masked points. 2023-01-23 14:58:17,557 [INFO] Starting creation of NetCDF products... 2023-01-23 14:58:17,569 [INFO] Getting data from Cyclobs track API. Req : https://cyclobs.ifremer.fr/app/api/track?include_cols=all&sid=wp022021&freq=15&min_date=2021-04-18&max_date=2021-04-22 2023-01-23 14:58:17,572 [INFO] Acquisition datetime is 2021-04-20 09:42:08.500000 2023-01-23 14:58:17,865 [INFO] Found propagation angle to north is 14.139702782205147°, dir std 0.07510394777928019, speed : 2.9937760415601495 m/s, speed_std 0.2549909103450939 2023-01-23 14:58:17,936 [INFO] Creating aeqd gridded dataset... 2023-01-23 14:58:23,220 [INFO] Done creating aeqd gridded dataset. 2023-01-23 14:58:23,223 [INFO] Converting dataset to polar... 2023-01-23 14:58:24,932 [INFO] Done converting dataset to polar. 2023-01-23 14:58:25,083 [INFO] Rotating polar dataset... /home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_37/lib/python3.7/site-packages/tc_center/polar.py:511: RuntimeWarning: invalid value encountered in remainder polar_ds[var_diff].values + round(direction_angle)), 360) 2023-01-23 14:58:25,394 [INFO] Done rotating polar dataset.. 2023-01-23 14:58:25,546 [INFO] Rotating cartesian dataset... /home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_37/lib/python3.7/site-packages/tc_center/polar.py:573: RuntimeWarning: invalid value encountered in remainder ds_cyclone_gd[var_diff][:] = np.mod(np.where(nan_mask, np.nan, ds_cyclone_gd[var_diff].values + angle), 360) 2023-01-23 14:58:27,083 [INFO] Done rotating cartesian dataset. 2023-01-23 14:58:27,460 [INFO] Adding wind direction related variables... 2023-01-23 14:58:31,127 [INFO] Done adding wind direction related variables. 2023-01-23 14:58:31,359 [INFO] Computing radii using polar north dataset... 2023-01-23 14:58:41,297 [ERROR] Exception: local variable 'closest_rmax' referenced before assignment Traceback (most recent call last): File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_37/lib/python3.7/site-packages/tc_center/__init__.py", line 350, in find_center_save **eye_dict) File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_37/lib/python3.7/site-packages/tc_center/analysis_product.py", line 1551, in center_products_analysis atcf_vmax=atcf_vmax, storm_name=storm_name, sid=sid, eye_ok=eye_ok, sat_file_light=True) File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_37/lib/python3.7/site-packages/tc_center/analysis_product.py", line 1431, in center_products radii_quad, radii_nmi_quad, percent_valid_quad = compute_radii_new(ds_polar_north, rmax_ck22, eye_ok=eye_ok) File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_37/lib/python3.7/site-packages/tc_center/analysis_product.py", line 381, in compute_radii_new range_list=[30, 10]) File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_37/lib/python3.7/site-packages/tc_center/analysis_product.py", line 348, in compute_vmax t_start=t_start, t_stop=t_stop) File "/home/datawork-cersat-public/cache/project/cyclobs/tc_center/tc_center_37/lib/python3.7/site-packages/tc_center/analysis_product.py", line 263, in recursive_vmax_window rmax = closest_rmax.coords["rad"].item() UnboundLocalError: local variable 'closest_rmax' referenced before assignment exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc exc