# but WITHOUT ANY WARRANTY; without even the implied warranty of # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the # GNU General Public License for more details. # You should have received a copy of the GNU General Public License # along with this program. If not, see <http://www.gnu.org/licenses/>. """Plot positions of Øysteins radiation validation stations on a map """ from pps_matchup import get_radvaldata # ------------------------------------------------------------------------- if __name__ == "__main__": datapoints = get_radvaldata("radval-stlist.txt") lonpos = datapoints['lon'] latpos = datapoints['lat'] from mpl_toolkits.basemap import Basemap import matplotlib.pyplot as plt fig = plt.figure(figsize=(11, 7)) ax = fig.add_axes([0.1, 0.1, 0.8, 0.8]) # m = Basemap(llcrnrlon=-180, llcrnrlat=-80, # urcrnrlon=180, urcrnrlat=80, projection='mill') m = Basemap(projection='stere', lon_0=0, lat_0=90., lat_ts=60, llcrnrlat=40, urcrnrlat=60, llcrnrlon=-30, urcrnrlon=120,
# This program is distributed in the hope that it will be useful, # but WITHOUT ANY WARRANTY; without even the implied warranty of # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the # GNU General Public License for more details. # You should have received a copy of the GNU General Public License # along with this program. If not, see <http://www.gnu.org/licenses/>. """Plot positions of Øysteins radiation validation stations on a map """ from pps_matchup import get_radvaldata # ------------------------------------------------------------------------- if __name__ == "__main__": datapoints = get_radvaldata("radval-stlist.txt") lonpos = datapoints['lon'] latpos = datapoints['lat'] from mpl_toolkits.basemap import Basemap import matplotlib.pyplot as plt fig = plt.figure(figsize=(11, 7)) ax = fig.add_axes([0.1, 0.1, 0.8, 0.8]) # m = Basemap(llcrnrlon=-180, llcrnrlat=-80, # urcrnrlon=180, urcrnrlat=80, projection='mill') m = Basemap(projection='stere', lon_0=0, lat_0=90.,
# along with this program. If not, see <http://www.gnu.org/licenses/>. """Filter/reduce the ascii files with synop-satellite matchups to contain only selected stations """ import pandas as pd #RADVAL_FILE = "./radval-stlist.txt" RADVAL_FILE = "./radval-stlist-wmo.txt" SYNOPFILE = "./matchup_synop_all_npp.txt" OUTPUTFILE = "./matchup_synop_radval_npp.txt" from pps_matchup import get_radvaldata datapoints = get_radvaldata(RADVAL_FILE) pdata = pd.DataFrame(datapoints) ids = [x[1]['id'] for x in pdata.iterrows()] with open(SYNOPFILE, "r") as fpt: lines = fpt.readlines() outlines = [] for line in lines[3:]: sl_ = line.split() if sl_[2] in ids: outlines.append(line) with open(OUTPUTFILE, "w") as fpt: fpt.writelines(lines[0:3] + outlines)
"""Filter/reduce the ascii files with synop-satellite matchups to contain only selected stations """ import pandas as pd #RADVAL_FILE = "./radval-stlist.txt" RADVAL_FILE = "./radval-stlist-wmo.txt" SYNOPFILE = "./matchup_synop_all_npp.txt" OUTPUTFILE = "./matchup_synop_radval_npp.txt" from pps_matchup import get_radvaldata datapoints = get_radvaldata(RADVAL_FILE) pdata = pd.DataFrame(datapoints) ids = [x[1]['id'] for x in pdata.iterrows()] with open(SYNOPFILE, "r") as fpt: lines = fpt.readlines() outlines = [] for line in lines[3:]: sl_ = line.split() if sl_[2] in ids: outlines.append(line) with open(OUTPUTFILE, "w") as fpt: fpt.writelines(lines[0:3] + outlines)