My latest encounter with joys of input format change: With the change to GFSv16 in March 2021 additional height levels were added to the grib output. Especially heights below 1hPa. What is generally good news, is bad news for Flexpart 10.4 (as of 20 Nov 2019). The additional levels cause incompatability with the readwinds_gfs.f90
. Public information on this issue is rare. There is an Issue in Flexparts ticketing system (https://www.flexpart.eu/ticket/305) but thats more related to the mean-wind trajectory module Flextra.
One suggestion is to drop the additional variables with wgrib2, which is somehow dissatisfying (and I could not figure it out on the short run). The only alternative posed ist, to switch to the dev branch (The link in the issue points to a version 9.3.2 legacy branch). So, for 10.4 dev you want to pull git clone https://www.flexpart.eu/gitmob/flexpart --branch dev --single-branch
. With that version, Flexpart switches from grib_api to eccodes, so on ubuntu the libraries apt-get install libeccodes-tools libeccodes-dev
are required.
At least for my docker setup, i needed to adapt the makefile:
INCPATH1 = /usr/include INCPATH2 = /usr/include LIBPATH1 = /usr/lib F90 = gfortran
Then, compiling worked sucessfully and the binary could handle the new GFS files.
Slightly related: natural earth cartopy downloads
Seems like the default natural earth mirror is temporarily unavailable (e.g., http://naciscdn.org/naturalearth/110m/physical), but there is this fix: https://github.com/SciTools/cartopy/issues/1298#issuecomment-920843582