lhenneman / hyspdisp

4 stars 4 forks source link

Unable to apply hyspdisp_zip_link to new units #59

Open ccwright opened 4 years ago

ccwright commented 4 years ago

Hi,

I am trying to link hyspdisp output to ZIPs for the following unit in my data. Lat: 31.4219 Lon: -96.2525 Stack height: 563 Month: 01/2006

I successfully obtain hysplit output for this unit. My problem comes when trying to link to ZIPs using hyspdisp_zip_link(). The function runs, but the returned dataset is empty. Here are the messages in my console: [1] "2019-10-10 01:47:07 Files read and combined" [1] "2019-10-10 01:47:09 PBLs trimmed" [1] "2019-10-10 01:47:20 ZIPs linked" [1] "2019-10-10 01:47:20 No ZIPs to link!" Warning message: In rgdal::rawTransform(projto_int, projfrom, nrow(xy), xy[, 1], : 16860 projected point(s) not finite

I tried to locate in the underlying function where rawTransform() is used, but I cannot find it.

Any help would be much appreciated. Thank you very much.