Python利用MJO指数画MJO位相图

 代码由本人拙劣的模仿写成，难免有错漏，欢迎指正

MJO指数数据和样图可在网站下载和查看http://www.bom.gov.au/climate/mjo/

import numpy as np
import datetime as dt

import matplotlib.pyplot as plt
from matplotlib.font_manager import fontManager
from matplotlib.ticker import MultipleLocator
from netCDF4 import Dataset
from netCDF4 import num2date

N=20

fontManager.addfont(r'C:\Users\11620\AppData\Local\Microsoft\Windows\Fonts\SourceHanSansCN-Normal.ttf') #r取消字符串转义

def readtxt(pth): #读取text函数
    rmm1,rmm2=[],[]
    f=open(pth)
    for i,x in enumerate(f.readlines()):
        if(i>2):
            x=x.split()
            rmm1.append(float(x[3]))
            rmm2.append(float(x[4]))
    f.close()
    return rmm1,rmm2

def textin(x,y,tin,size,rot,bg,fontf,ax): #文字输入函数
    if(bg!='none'):
        ax.text(x,y,tin,fontfamily