from scipy.stats import norm
from pandas import DataFrame
from matplotlib import pyplot as plt
from mpl_toolkits.basemap import Basemap
import numpy as np
# Cookie処理
cj = MozillaCookieJar()
account = urlencode({
'email': argv[1],
'password': argv[2],
'action': 'login'
}).encode('utf-8')
opener = build_opener(HTTPCookieProcessor(cj))
login = opener.open("https://rda.ucar.edu/cgi-bin/login", account)
cj.clear_session_cookies()
cj.save("auth.rda.ucar.edu", True, True)
# 海面更正気圧データのダウンロードを行う関数
def prmsl_download(year):
# ダウンロードするファイル
if year < 2014:
_f = ['anl_surf125.002_prmsl.{0}010100_{0}123118'.format(year)]
elif ((year % 4 == 0) & (year % 100 != 0)) | (year % 400 == 0):
_f = [
'anl_surf125.002_prmsl.{0}020100_{0}022918'.format(year),
'anl_surf125.002_prmsl.{0}030100_{0}033118'.format(year),
'anl_surf125.002_prmsl.{0}040100_{0}043018'.format(year),
'anl_surf125.002_prmsl.{0}050100_{0}053118'.format(year),
