#!/usr/bin/env python # Wenchang Yang (wenchang@princeton.edu) # Tue Feb 18 04:26:56 PM EST 2025 if __name__ == '__main__': import sys,os try: from misc.timer import Timer tt = Timer(f'[{os.getcwd()}] start ' + ' '.join(sys.argv)) except: pass #import sys, os.path, os, glob, datetime #import xarray as xr, numpy as np, pandas as pd, matplotlib.pyplot as plt #more imports #wython = '/tigress/wenchang/wython' #if wython not in sys.path: sys.path.append(wython); print('added to python path:', wython) #from misc import get_kws_from_argv # if __name__ == '__main__': try: tt.check('end import') except: pass # #start from here #mjo phase5 time = ['1979-06-24T00:00:00', '1979-06-26T00:00:00', '1979-06-27T00:00:00', '1979-07-30T00:00:00', '1979-07-31T00:00:00', '1979-08-01T00:00:00', '1979-08-02T00:00:00', '1979-08-03T00:00:00', '1979-08-04T00:00:00', '1979-08-05T00:00:00', '1979-08-06T00:00:00', '1979-08-07T00:00:00', '1979-09-26T00:00:00', '1979-09-27T00:00:00', '1979-09-28T00:00:00', '1979-09-29T00:00:00', '1979-09-30T00:00:00', '1979-10-01T00:00:00', '1979-10-02T00:00:00', '1979-10-03T00:00:00', '1979-10-04T00:00:00', '1979-10-05T00:00:00', '1980-07-03T00:00:00', '1980-07-04T00:00:00', '1980-07-05T00:00:00', '1980-08-24T00:00:00', '1980-08-25T00:00:00', '1980-08-26T00:00:00', '1980-08-27T00:00:00', '1980-08-28T00:00:00', '1980-08-29T00:00:00', '1980-09-05T00:00:00', '1980-09-06T00:00:00', '1980-09-07T00:00:00', '1980-09-08T00:00:00', '1980-09-09T00:00:00', '1980-09-10T00:00:00', '1980-10-19T00:00:00', '1980-10-20T00:00:00', '1980-10-21T00:00:00', '1980-10-22T00:00:00', '1980-10-23T00:00:00', '1980-10-24T00:00:00', '1980-10-25T00:00:00', '1980-10-26T00:00:00', '1981-06-03T00:00:00', '1981-06-12T00:00:00', '1981-06-13T00:00:00', '1981-06-14T00:00:00', '1981-06-22T00:00:00', '1981-07-29T00:00:00', '1981-07-30T00:00:00', '1981-07-31T00:00:00', '1981-08-01T00:00:00', '1981-09-19T00:00:00', '1981-09-20T00:00:00', '1981-09-21T00:00:00', '1981-09-22T00:00:00', '1981-09-23T00:00:00', '1981-09-24T00:00:00', '1981-09-25T00:00:00', '1981-09-26T00:00:00', '1981-09-27T00:00:00', '1981-09-28T00:00:00', '1982-06-06T00:00:00', '1982-09-08T00:00:00', '1982-11-08T00:00:00', '1982-11-09T00:00:00', '1983-06-25T00:00:00', '1983-06-26T00:00:00', '1983-06-27T00:00:00', '1983-06-28T00:00:00', '1983-09-25T00:00:00', '1983-09-26T00:00:00', '1983-09-27T00:00:00', '1983-09-28T00:00:00', '1983-09-29T00:00:00', '1983-09-30T00:00:00', '1983-10-01T00:00:00', '1983-10-02T00:00:00', '1983-10-03T00:00:00', '1983-10-04T00:00:00', '1983-10-05T00:00:00', '1983-10-06T00:00:00', '1983-10-07T00:00:00', '1983-10-08T00:00:00', '1983-10-09T00:00:00', '1983-10-24T00:00:00', '1983-10-25T00:00:00', '1983-10-26T00:00:00', '1983-10-31T00:00:00', '1983-11-01T00:00:00', '1983-11-15T00:00:00', '1983-11-16T00:00:00', '1984-07-19T00:00:00', '1984-07-20T00:00:00', '1984-07-21T00:00:00', '1984-07-22T00:00:00', '1984-07-23T00:00:00', '1984-08-19T00:00:00', '1984-08-20T00:00:00', '1984-08-21T00:00:00', '1984-09-30T00:00:00', '1984-10-10T00:00:00', '1984-10-11T00:00:00', '1984-10-12T00:00:00', '1985-06-16T00:00:00', '1985-06-17T00:00:00', '1985-06-18T00:00:00', '1985-06-19T00:00:00', '1985-06-20T00:00:00', '1985-06-21T00:00:00', '1985-06-22T00:00:00', '1985-06-23T00:00:00', '1985-08-20T00:00:00', '1985-11-14T00:00:00', '1985-11-15T00:00:00', '1985-11-16T00:00:00', '1986-06-17T00:00:00', '1986-06-18T00:00:00', '1986-06-19T00:00:00', '1986-06-20T00:00:00', '1986-06-21T00:00:00', '1986-06-22T00:00:00', '1986-06-23T00:00:00', '1986-08-01T00:00:00', '1986-08-02T00:00:00', '1986-08-03T00:00:00', '1986-08-04T00:00:00', '1986-08-05T00:00:00', '1986-09-12T00:00:00', '1986-09-13T00:00:00', '1986-09-14T00:00:00', '1986-09-15T00:00:00', '1986-11-11T00:00:00', '1986-11-12T00:00:00', '1986-11-13T00:00:00', '1987-06-09T00:00:00', '1987-06-10T00:00:00', '1987-06-11T00:00:00', '1987-06-12T00:00:00', '1987-06-13T00:00:00', '1987-06-14T00:00:00', '1987-06-15T00:00:00', '1987-06-16T00:00:00', '1987-06-17T00:00:00', '1987-06-27T00:00:00', '1987-06-29T00:00:00', '1987-06-30T00:00:00', '1987-07-01T00:00:00', '1987-07-02T00:00:00', '1987-08-28T00:00:00', '1987-08-29T00:00:00', '1987-08-30T00:00:00', '1987-08-31T00:00:00', '1987-09-01T00:00:00', '1987-09-02T00:00:00', '1987-09-03T00:00:00', '1987-09-04T00:00:00', '1987-10-16T00:00:00', '1987-10-17T00:00:00', '1987-10-18T00:00:00', '1987-10-19T00:00:00', '1987-10-20T00:00:00', '1987-10-21T00:00:00', '1987-11-13T00:00:00', '1987-11-14T00:00:00', '1987-11-15T00:00:00', '1988-09-12T00:00:00', '1988-09-27T00:00:00', '1988-09-28T00:00:00', '1988-09-29T00:00:00', '1988-09-30T00:00:00', '1988-10-01T00:00:00', '1988-10-02T00:00:00', '1988-10-03T00:00:00', '1988-10-04T00:00:00', '1988-10-05T00:00:00', '1988-11-08T00:00:00', '1988-11-09T00:00:00', '1988-11-24T00:00:00', '1989-06-08T00:00:00', '1989-06-09T00:00:00', '1989-06-10T00:00:00', '1989-07-15T00:00:00', '1989-07-16T00:00:00', '1989-07-17T00:00:00', '1989-10-19T00:00:00', '1989-10-20T00:00:00', '1989-10-21T00:00:00', '1989-10-22T00:00:00', '1989-10-23T00:00:00', '1989-10-24T00:00:00', '1989-10-25T00:00:00', '1989-11-18T00:00:00', '1989-11-20T00:00:00', '1990-08-14T00:00:00', '1990-08-15T00:00:00', '1990-08-16T00:00:00', '1990-08-17T00:00:00', '1990-08-18T00:00:00', '1990-08-19T00:00:00', '1990-11-10T00:00:00', '1991-06-05T00:00:00', '1991-06-06T00:00:00', '1991-06-07T00:00:00', '1991-06-08T00:00:00', '1991-06-09T00:00:00', '1991-06-10T00:00:00', '1991-06-11T00:00:00', '1991-07-19T00:00:00', '1991-07-20T00:00:00', '1991-07-21T00:00:00', '1991-08-13T00:00:00', '1991-08-14T00:00:00', '1991-08-15T00:00:00', '1991-08-16T00:00:00', '1991-08-17T00:00:00', '1991-08-18T00:00:00', '1991-08-19T00:00:00', '1991-08-20T00:00:00', '1991-08-21T00:00:00', '1991-08-23T00:00:00', '1991-08-24T00:00:00', '1991-08-25T00:00:00', '1991-08-26T00:00:00', '1991-08-27T00:00:00', '1991-08-28T00:00:00', '1991-08-29T00:00:00', '1991-08-30T00:00:00', '1991-09-03T00:00:00', '1991-09-04T00:00:00', '1991-09-05T00:00:00', '1991-09-06T00:00:00', '1992-06-19T00:00:00', '1992-06-20T00:00:00', '1992-06-21T00:00:00', '1992-06-22T00:00:00', '1992-08-17T00:00:00', '1992-08-18T00:00:00', '1992-08-19T00:00:00', '1992-08-28T00:00:00', '1992-08-29T00:00:00', '1992-08-30T00:00:00', '1992-08-31T00:00:00', '1992-09-01T00:00:00', '1992-09-02T00:00:00', '1992-10-17T00:00:00', '1992-10-18T00:00:00', '1992-10-19T00:00:00', '1992-10-20T00:00:00', '1992-10-21T00:00:00', '1992-10-22T00:00:00', '1992-10-23T00:00:00', '1993-06-07T00:00:00', '1993-06-08T00:00:00', '1993-07-25T00:00:00', '1993-07-26T00:00:00', '1993-07-27T00:00:00', '1993-07-28T00:00:00', '1993-07-29T00:00:00', '1993-07-30T00:00:00', '1993-07-31T00:00:00', '1993-08-01T00:00:00', '1993-08-02T00:00:00', '1993-08-03T00:00:00', '1993-08-28T00:00:00', '1993-11-13T00:00:00', '1994-07-08T00:00:00', '1995-09-02T00:00:00', '1995-09-03T00:00:00', '1995-09-04T00:00:00', '1995-09-05T00:00:00', '1995-09-06T00:00:00', '1995-09-19T00:00:00', '1995-09-20T00:00:00', '1995-11-13T00:00:00', '1995-11-14T00:00:00', '1995-11-15T00:00:00', '1995-11-17T00:00:00', '1995-11-18T00:00:00', '1995-11-19T00:00:00', '1995-11-20T00:00:00', '1995-11-21T00:00:00', '1995-11-22T00:00:00', '1995-11-23T00:00:00', '1996-06-17T00:00:00', '1996-06-18T00:00:00', '1996-06-19T00:00:00', '1996-07-11T00:00:00', '1996-07-12T00:00:00', '1996-07-13T00:00:00', '1996-07-14T00:00:00', '1996-09-08T00:00:00', '1996-09-09T00:00:00', '1996-09-10T00:00:00', '1996-09-11T00:00:00', '1996-09-12T00:00:00', '1996-09-13T00:00:00', '1996-09-14T00:00:00', '1996-09-15T00:00:00', '1996-09-17T00:00:00', '1996-10-30T00:00:00', '1996-10-31T00:00:00', '1996-11-01T00:00:00', '1996-11-02T00:00:00', '1996-11-03T00:00:00', '1996-11-04T00:00:00', '1996-11-05T00:00:00', '1996-11-06T00:00:00', '1996-11-07T00:00:00', '1997-07-23T00:00:00', '1997-07-24T00:00:00', '1998-06-11T00:00:00', '1998-09-29T00:00:00', '1998-09-30T00:00:00', '1998-10-01T00:00:00', '1998-10-02T00:00:00', '1998-10-03T00:00:00', '1998-10-04T00:00:00', '1998-10-05T00:00:00', '1998-10-06T00:00:00', '1998-10-07T00:00:00', '1998-10-08T00:00:00', '1998-11-09T00:00:00', '1998-11-10T00:00:00', '1998-11-11T00:00:00', '1998-11-12T00:00:00', '1998-11-13T00:00:00', '1998-11-14T00:00:00', '1998-11-15T00:00:00', '1999-07-27T00:00:00', '1999-07-28T00:00:00', '1999-10-31T00:00:00', '1999-11-01T00:00:00', '1999-11-02T00:00:00', '1999-11-03T00:00:00', '1999-11-04T00:00:00', '1999-11-05T00:00:00', '2000-07-05T00:00:00', '2000-07-06T00:00:00', '2000-07-07T00:00:00', '2000-07-09T00:00:00', '2000-08-15T00:00:00', '2000-08-16T00:00:00', '2000-08-18T00:00:00', '2000-08-19T00:00:00', '2000-08-20T00:00:00', '2000-08-21T00:00:00', '2000-08-22T00:00:00', '2000-08-23T00:00:00', '2000-08-24T00:00:00', '2000-08-25T00:00:00', '2000-10-08T00:00:00', '2000-10-09T00:00:00', '2000-10-10T00:00:00', '2000-10-11T00:00:00', '2000-10-12T00:00:00', '2000-10-13T00:00:00', '2000-10-14T00:00:00', '2000-10-15T00:00:00', '2000-10-16T00:00:00', '2000-10-17T00:00:00', '2000-10-18T00:00:00', '2000-10-19T00:00:00', '2000-10-20T00:00:00', '2000-10-21T00:00:00', '2000-10-22T00:00:00', '2000-10-23T00:00:00', '2000-11-27T00:00:00', '2000-11-30T00:00:00', '2001-06-19T00:00:00', '2001-08-10T00:00:00', '2001-08-11T00:00:00', '2001-08-12T00:00:00', '2001-08-13T00:00:00', '2001-08-14T00:00:00', '2001-08-15T00:00:00', '2001-10-07T00:00:00', '2001-10-08T00:00:00', '2001-10-09T00:00:00', '2001-11-21T00:00:00', '2001-11-22T00:00:00', '2002-06-20T00:00:00', '2002-06-21T00:00:00', '2002-06-22T00:00:00', '2002-06-23T00:00:00', '2002-08-11T00:00:00', '2002-08-12T00:00:00', '2002-08-13T00:00:00', '2002-08-14T00:00:00', '2002-10-07T00:00:00', '2002-10-08T00:00:00', '2002-10-09T00:00:00', '2002-11-24T00:00:00', '2002-11-25T00:00:00', '2002-11-26T00:00:00', '2003-07-22T00:00:00', '2003-07-23T00:00:00', '2003-07-24T00:00:00', '2003-08-23T00:00:00', '2003-09-16T00:00:00', '2003-09-17T00:00:00', '2003-10-12T00:00:00', '2003-10-13T00:00:00', '2003-10-14T00:00:00', '2003-10-15T00:00:00', '2003-10-16T00:00:00', '2003-10-17T00:00:00', '2003-10-18T00:00:00', '2003-10-19T00:00:00', '2004-08-07T00:00:00', '2004-10-04T00:00:00', '2004-10-05T00:00:00', '2004-10-06T00:00:00', '2004-10-07T00:00:00', '2004-10-08T00:00:00', '2004-10-09T00:00:00', '2004-10-10T00:00:00', '2005-07-25T00:00:00', '2005-07-26T00:00:00', '2005-07-27T00:00:00', '2005-08-01T00:00:00', '2005-08-02T00:00:00', '2005-09-09T00:00:00', '2005-09-10T00:00:00', '2005-09-11T00:00:00', '2005-09-12T00:00:00', '2005-09-13T00:00:00', '2005-09-14T00:00:00', '2005-09-15T00:00:00', '2005-09-16T00:00:00', '2005-09-17T00:00:00', '2005-10-17T00:00:00', '2005-10-18T00:00:00', '2005-10-19T00:00:00', '2005-10-20T00:00:00', '2005-11-05T00:00:00', '2006-08-07T00:00:00', '2006-08-08T00:00:00', '2006-08-09T00:00:00', '2006-08-10T00:00:00', '2006-08-11T00:00:00', '2006-09-30T00:00:00', '2006-10-01T00:00:00', '2006-10-02T00:00:00', '2006-10-03T00:00:00', '2006-10-04T00:00:00', '2007-06-25T00:00:00', '2007-06-26T00:00:00', '2007-09-09T00:00:00', '2007-09-10T00:00:00', '2007-09-11T00:00:00', '2007-09-12T00:00:00', '2007-09-13T00:00:00', '2007-09-14T00:00:00', '2007-09-15T00:00:00', '2007-09-16T00:00:00', '2007-09-17T00:00:00', '2007-09-18T00:00:00', '2007-09-19T00:00:00', '2007-09-20T00:00:00', '2007-09-21T00:00:00', '2007-09-22T00:00:00', '2007-09-23T00:00:00', '2007-09-24T00:00:00', '2007-10-27T00:00:00', '2007-10-28T00:00:00', '2007-11-13T00:00:00', '2007-11-14T00:00:00', '2007-11-15T00:00:00', '2007-11-16T00:00:00', '2008-09-11T00:00:00', '2008-09-12T00:00:00', '2008-09-13T00:00:00', '2008-09-14T00:00:00', '2008-09-15T00:00:00', '2008-09-16T00:00:00', '2008-09-17T00:00:00', '2008-09-18T00:00:00', '2008-09-19T00:00:00', '2008-09-20T00:00:00', '2008-09-21T00:00:00', '2008-09-22T00:00:00', '2008-09-23T00:00:00', '2008-10-28T00:00:00', '2008-10-29T00:00:00', '2008-10-30T00:00:00', '2008-10-31T00:00:00', '2008-11-15T00:00:00', '2008-11-16T00:00:00', '2008-11-17T00:00:00', '2008-11-18T00:00:00', '2008-11-19T00:00:00', '2008-11-20T00:00:00', '2008-11-21T00:00:00', '2008-11-22T00:00:00', '2008-11-23T00:00:00', '2008-11-24T00:00:00', '2009-07-28T00:00:00', '2009-09-09T00:00:00', '2009-09-25T00:00:00', '2009-09-26T00:00:00', '2009-09-27T00:00:00', '2009-09-28T00:00:00', '2009-09-29T00:00:00', '2009-10-07T00:00:00', '2009-10-08T00:00:00', '2009-11-21T00:00:00', '2009-11-22T00:00:00', '2009-11-23T00:00:00', '2010-08-23T00:00:00', '2010-08-24T00:00:00', '2010-08-25T00:00:00', '2010-08-26T00:00:00', '2010-08-27T00:00:00', '2010-10-08T00:00:00', '2010-10-09T00:00:00', '2010-10-10T00:00:00', '2010-10-11T00:00:00', '2011-06-11T00:00:00', '2011-07-21T00:00:00', '2011-07-22T00:00:00', '2011-07-23T00:00:00', '2011-07-24T00:00:00', '2011-09-25T00:00:00', '2011-09-26T00:00:00', '2011-09-27T00:00:00', '2011-09-28T00:00:00', '2011-10-02T00:00:00', '2011-11-06T00:00:00', '2011-11-07T00:00:00', '2011-11-08T00:00:00', '2012-06-03T00:00:00', '2012-06-04T00:00:00', '2012-06-05T00:00:00', '2012-06-06T00:00:00', '2012-06-07T00:00:00', '2012-06-08T00:00:00', '2012-09-12T00:00:00', '2012-09-13T00:00:00', '2012-09-14T00:00:00', '2012-09-23T00:00:00', '2012-09-24T00:00:00', '2012-09-27T00:00:00', '2012-09-28T00:00:00', '2012-09-29T00:00:00', '2012-09-30T00:00:00', '2012-10-05T00:00:00', '2013-06-15T00:00:00', '2013-06-16T00:00:00', '2013-06-17T00:00:00', '2013-06-18T00:00:00', '2013-09-15T00:00:00', '2013-09-17T00:00:00', '2013-09-18T00:00:00', '2013-09-19T00:00:00', '2013-09-20T00:00:00', '2013-09-21T00:00:00', '2013-09-22T00:00:00', '2013-09-23T00:00:00', '2013-09-24T00:00:00', '2013-09-29T00:00:00', '2013-09-30T00:00:00', '2013-10-01T00:00:00', '2013-10-02T00:00:00', '2013-10-03T00:00:00', '2013-10-04T00:00:00', '2013-10-05T00:00:00', '2013-10-26T00:00:00', '2014-06-18T00:00:00', '2014-06-19T00:00:00', '2014-06-20T00:00:00', '2014-07-09T00:00:00', '2014-07-10T00:00:00', '2014-07-11T00:00:00', '2014-07-12T00:00:00', '2014-07-13T00:00:00', '2014-07-18T00:00:00', '2014-07-19T00:00:00', '2015-06-22T00:00:00', '2015-06-23T00:00:00', '2015-06-24T00:00:00', '2015-06-25T00:00:00', '2015-06-26T00:00:00', '2016-06-27T00:00:00', '2016-06-28T00:00:00', '2016-08-06T00:00:00', '2016-08-07T00:00:00', '2016-08-08T00:00:00', '2016-08-09T00:00:00', '2016-08-21T00:00:00', '2016-09-19T00:00:00', '2016-09-20T00:00:00', '2016-09-25T00:00:00', '2016-09-26T00:00:00', '2016-09-27T00:00:00', '2016-09-28T00:00:00', '2016-09-29T00:00:00', '2016-09-30T00:00:00', '2016-10-01T00:00:00', '2016-10-02T00:00:00', '2016-10-03T00:00:00', '2016-10-04T00:00:00', '2016-10-05T00:00:00', '2016-10-06T00:00:00', '2016-11-04T00:00:00', '2016-11-05T00:00:00', '2016-11-06T00:00:00', '2017-07-24T00:00:00', '2017-07-25T00:00:00', '2017-07-26T00:00:00', '2017-10-13T00:00:00', '2017-10-14T00:00:00', '2017-10-15T00:00:00', '2017-10-16T00:00:00', '2017-10-17T00:00:00', '2017-10-18T00:00:00', '2017-10-19T00:00:00', '2018-06-01T00:00:00', '2018-06-02T00:00:00', '2018-06-03T00:00:00', '2018-06-04T00:00:00', '2018-07-11T00:00:00', '2018-07-12T00:00:00', '2018-07-13T00:00:00', '2018-07-14T00:00:00', '2018-07-15T00:00:00', '2018-07-16T00:00:00', '2018-07-17T00:00:00', '2018-07-18T00:00:00', '2018-11-13T00:00:00', '2018-11-14T00:00:00', '2018-11-18T00:00:00', '2019-06-16T00:00:00', '2019-06-18T00:00:00', '2019-08-12T00:00:00', '2019-09-03T00:00:00', '2019-09-04T00:00:00', '2019-09-07T00:00:00', '2019-09-08T00:00:00', '2019-09-09T00:00:00', '2019-09-10T00:00:00', '2019-09-11T00:00:00', '2019-11-01T00:00:00', '2019-11-02T00:00:00', '2019-11-03T00:00:00', '2019-11-04T00:00:00', '2019-11-05T00:00:00', '2019-11-06T00:00:00', '2019-11-07T00:00:00', '2020-08-07T00:00:00', '2020-08-08T00:00:00', '2020-08-09T00:00:00', '2020-08-10T00:00:00', '2020-09-22T00:00:00', '2020-09-23T00:00:00', '2020-09-24T00:00:00', '2020-09-25T00:00:00', '2020-09-26T00:00:00', '2020-09-27T00:00:00', '2020-09-28T00:00:00', '2020-09-29T00:00:00', '2020-09-30T00:00:00', '2020-10-01T00:00:00', '2020-10-02T00:00:00', '2020-10-03T00:00:00', '2020-10-05T00:00:00', '2020-10-06T00:00:00', '2020-10-07T00:00:00', '2020-10-14T00:00:00', '2020-10-15T00:00:00', '2020-10-16T00:00:00', '2020-10-17T00:00:00', '2020-10-18T00:00:00', '2020-10-19T00:00:00', '2020-10-20T00:00:00', '2020-10-21T00:00:00', '2020-10-22T00:00:00', '2020-10-23T00:00:00', '2020-10-24T00:00:00', '2020-10-25T00:00:00'] if __name__ == '__main__': #from wyconfig import * #my plot settings #savefig if 'savefig' in sys.argv or 's' in sys.argv: figname = __file__.replace('.py', f'.png') if 'overwritefig' in sys.argv or 'o' in sys.argv: wysavefig(figname, overwritefig=True) else: wysavefig(figname) try: tt.check(f'**Done**') except: pass print() if 'notshowfig' in sys.argv or 'n' in sys.argv: pass else: if 'plt' in globals(): plt.show()