Python script to get a series of surface and sounding stations from a netCDF file. It computes also the different diagnostic variables