#!/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 phase4
time = ['1979-06-20T00:00:00',
 '1979-06-21T00:00:00',
 '1979-06-22T00:00:00',
 '1979-06-23T00:00:00',
 '1979-06-25T00:00:00',
 '1979-09-22T00:00:00',
 '1979-09-23T00:00:00',
 '1979-09-24T00:00:00',
 '1979-09-25T00:00:00',
 '1979-11-26T00:00:00',
 '1979-11-27T00:00:00',
 '1979-11-28T00:00:00',
 '1979-11-29T00:00:00',
 '1980-07-02T00:00:00',
 '1980-08-21T00:00:00',
 '1980-08-22T00:00:00',
 '1980-08-23T00:00:00',
 '1980-09-01T00:00:00',
 '1980-09-02T00:00:00',
 '1980-10-17T00:00:00',
 '1980-10-18T00:00:00',
 '1980-11-24T00:00:00',
 '1981-06-19T00:00:00',
 '1981-06-20T00:00:00',
 '1981-06-21T00:00:00',
 '1981-07-27T00:00:00',
 '1981-07-28T00:00:00',
 '1981-09-16T00:00:00',
 '1981-09-17T00:00:00',
 '1981-09-18T00:00:00',
 '1981-11-03T00:00:00',
 '1981-11-04T00:00:00',
 '1981-11-10T00:00:00',
 '1981-11-11T00:00:00',
 '1982-06-01T00:00:00',
 '1982-06-02T00:00:00',
 '1982-06-03T00:00:00',
 '1982-06-04T00:00:00',
 '1982-11-06T00:00:00',
 '1982-11-07T00:00:00',
 '1983-06-17T00:00:00',
 '1983-06-18T00:00:00',
 '1983-06-19T00:00:00',
 '1983-06-20T00:00:00',
 '1983-06-21T00:00:00',
 '1983-06-22T00:00:00',
 '1983-06-23T00:00:00',
 '1983-06-24T00:00:00',
 '1983-08-11T00:00:00',
 '1983-08-12T00:00:00',
 '1983-08-13T00:00:00',
 '1983-08-14T00:00:00',
 '1983-08-15T00:00:00',
 '1983-08-16T00:00:00',
 '1983-08-17T00:00:00',
 '1983-08-18T00:00:00',
 '1983-08-19T00:00:00',
 '1983-09-18T00:00:00',
 '1983-09-19T00:00:00',
 '1983-09-20T00:00:00',
 '1983-09-21T00:00:00',
 '1983-09-22T00:00:00',
 '1983-09-23T00:00:00',
 '1983-09-24T00:00:00',
 '1985-06-12T00:00:00',
 '1985-06-13T00:00:00',
 '1985-06-14T00:00:00',
 '1985-06-15T00:00:00',
 '1985-11-10T00:00:00',
 '1985-11-11T00:00:00',
 '1985-11-12T00:00:00',
 '1985-11-13T00:00:00',
 '1986-06-13T00:00:00',
 '1986-06-14T00:00:00',
 '1986-06-15T00:00:00',
 '1986-06-16T00:00:00',
 '1986-07-11T00:00:00',
 '1986-07-12T00:00:00',
 '1986-07-13T00:00:00',
 '1986-07-14T00:00:00',
 '1986-07-15T00:00:00',
 '1986-09-09T00:00:00',
 '1986-09-10T00:00:00',
 '1986-09-11T00:00:00',
 '1986-11-04T00:00:00',
 '1986-11-05T00:00:00',
 '1986-11-06T00:00:00',
 '1986-11-07T00:00:00',
 '1986-11-08T00:00:00',
 '1986-11-09T00:00:00',
 '1986-11-10T00:00:00',
 '1987-06-05T00:00:00',
 '1987-06-06T00:00:00',
 '1987-06-07T00:00:00',
 '1987-06-08T00:00:00',
 '1987-08-19T00:00:00',
 '1987-08-20T00:00:00',
 '1987-08-21T00:00:00',
 '1987-08-22T00:00:00',
 '1987-08-23T00:00:00',
 '1987-08-24T00:00:00',
 '1987-08-25T00:00:00',
 '1987-08-26T00:00:00',
 '1987-08-27T00:00:00',
 '1987-09-10T00:00:00',
 '1987-09-11T00:00:00',
 '1987-09-12T00:00:00',
 '1987-09-13T00:00:00',
 '1987-09-14T00:00:00',
 '1987-09-15T00:00:00',
 '1987-11-05T00:00:00',
 '1987-11-06T00:00:00',
 '1987-11-07T00:00:00',
 '1987-11-08T00:00:00',
 '1987-11-09T00:00:00',
 '1987-11-10T00:00:00',
 '1987-11-11T00:00:00',
 '1987-11-12T00:00:00',
 '1988-07-28T00:00:00',
 '1988-09-19T00:00:00',
 '1988-09-20T00:00:00',
 '1988-09-21T00:00:00',
 '1988-09-22T00:00:00',
 '1988-09-23T00:00:00',
 '1988-09-24T00:00:00',
 '1988-09-25T00:00:00',
 '1988-09-26T00:00:00',
 '1988-11-23T00:00:00',
 '1989-07-12T00:00:00',
 '1989-07-13T00:00:00',
 '1989-07-14T00:00:00',
 '1989-10-05T00:00:00',
 '1989-10-06T00:00:00',
 '1989-10-09T00:00:00',
 '1989-10-10T00:00:00',
 '1989-10-11T00:00:00',
 '1989-10-12T00:00:00',
 '1989-10-13T00:00:00',
 '1989-10-14T00:00:00',
 '1989-10-15T00:00:00',
 '1989-10-16T00:00:00',
 '1990-06-02T00:00:00',
 '1990-06-03T00:00:00',
 '1990-07-25T00:00:00',
 '1990-08-29T00:00:00',
 '1991-06-03T00:00:00',
 '1991-07-07T00:00:00',
 '1991-07-08T00:00:00',
 '1991-08-22T00:00:00',
 '1991-10-21T00:00:00',
 '1992-06-12T00:00:00',
 '1992-06-13T00:00:00',
 '1992-06-14T00:00:00',
 '1992-06-15T00:00:00',
 '1992-06-16T00:00:00',
 '1992-06-17T00:00:00',
 '1992-06-18T00:00:00',
 '1992-10-09T00:00:00',
 '1992-10-10T00:00:00',
 '1992-10-11T00:00:00',
 '1992-10-12T00:00:00',
 '1992-10-13T00:00:00',
 '1992-10-14T00:00:00',
 '1992-10-15T00:00:00',
 '1992-10-16T00:00:00',
 '1993-06-04T00:00:00',
 '1993-06-05T00:00:00',
 '1993-06-06T00:00:00',
 '1993-07-23T00:00:00',
 '1993-07-24T00:00:00',
 '1993-11-10T00:00:00',
 '1993-11-11T00:00:00',
 '1993-11-12T00:00:00',
 '1994-06-08T00:00:00',
 '1994-06-09T00:00:00',
 '1994-06-10T00:00:00',
 '1994-06-11T00:00:00',
 '1994-11-30T00:00:00',
 '1995-07-08T00:00:00',
 '1995-07-09T00:00:00',
 '1995-08-29T00:00:00',
 '1995-08-30T00:00:00',
 '1995-08-31T00:00:00',
 '1995-09-01T00:00:00',
 '1995-10-30T00:00:00',
 '1995-10-31T00:00:00',
 '1995-11-01T00:00:00',
 '1995-11-02T00:00:00',
 '1995-11-03T00:00:00',
 '1995-11-04T00:00:00',
 '1995-11-05T00:00:00',
 '1995-11-06T00:00:00',
 '1995-11-07T00:00:00',
 '1995-11-08T00:00:00',
 '1995-11-16T00:00:00',
 '1996-06-13T00:00:00',
 '1996-06-14T00:00:00',
 '1996-06-15T00:00:00',
 '1996-06-16T00:00:00',
 '1996-10-25T00:00:00',
 '1996-10-26T00:00:00',
 '1997-07-09T00:00:00',
 '1997-07-10T00:00:00',
 '1997-07-11T00:00:00',
 '1997-07-12T00:00:00',
 '1997-07-15T00:00:00',
 '1997-07-16T00:00:00',
 '1997-07-17T00:00:00',
 '1997-07-20T00:00:00',
 '1997-07-21T00:00:00',
 '1997-07-22T00:00:00',
 '1997-10-03T00:00:00',
 '1998-06-05T00:00:00',
 '1998-06-06T00:00:00',
 '1998-06-07T00:00:00',
 '1998-06-08T00:00:00',
 '1998-06-09T00:00:00',
 '1998-06-10T00:00:00',
 '1998-09-20T00:00:00',
 '1998-09-21T00:00:00',
 '1998-09-22T00:00:00',
 '1998-09-23T00:00:00',
 '1998-09-24T00:00:00',
 '1998-09-25T00:00:00',
 '1998-09-26T00:00:00',
 '1998-09-27T00:00:00',
 '1998-09-28T00:00:00',
 '1999-07-20T00:00:00',
 '1999-07-21T00:00:00',
 '1999-07-22T00:00:00',
 '1999-07-23T00:00:00',
 '1999-10-20T00:00:00',
 '1999-10-21T00:00:00',
 '1999-10-27T00:00:00',
 '1999-10-28T00:00:00',
 '1999-10-29T00:00:00',
 '1999-10-30T00:00:00',
 '1999-11-30T00:00:00',
 '2000-08-17T00:00:00',
 '2000-10-06T00:00:00',
 '2000-10-07T00:00:00',
 '2000-11-23T00:00:00',
 '2000-11-24T00:00:00',
 '2000-11-25T00:00:00',
 '2000-11-26T00:00:00',
 '2000-11-28T00:00:00',
 '2000-11-29T00:00:00',
 '2001-08-06T00:00:00',
 '2001-08-07T00:00:00',
 '2001-08-08T00:00:00',
 '2001-08-09T00:00:00',
 '2001-10-03T00:00:00',
 '2001-10-04T00:00:00',
 '2001-10-05T00:00:00',
 '2001-10-06T00:00:00',
 '2001-11-08T00:00:00',
 '2001-11-09T00:00:00',
 '2001-11-10T00:00:00',
 '2001-11-11T00:00:00',
 '2001-11-17T00:00:00',
 '2001-11-19T00:00:00',
 '2002-06-16T00:00:00',
 '2002-06-17T00:00:00',
 '2002-06-18T00:00:00',
 '2002-06-19T00:00:00',
 '2002-08-07T00:00:00',
 '2002-08-08T00:00:00',
 '2002-08-09T00:00:00',
 '2002-08-10T00:00:00',
 '2002-10-05T00:00:00',
 '2002-10-06T00:00:00',
 '2002-11-19T00:00:00',
 '2002-11-20T00:00:00',
 '2002-11-21T00:00:00',
 '2002-11-22T00:00:00',
 '2002-11-23T00:00:00',
 '2003-07-19T00:00:00',
 '2003-07-20T00:00:00',
 '2003-07-21T00:00:00',
 '2003-08-24T00:00:00',
 '2003-09-08T00:00:00',
 '2003-09-09T00:00:00',
 '2003-09-10T00:00:00',
 '2003-09-11T00:00:00',
 '2003-09-13T00:00:00',
 '2003-09-15T00:00:00',
 '2003-09-29T00:00:00',
 '2003-09-30T00:00:00',
 '2003-10-09T00:00:00',
 '2003-10-10T00:00:00',
 '2003-10-11T00:00:00',
 '2004-09-30T00:00:00',
 '2004-10-01T00:00:00',
 '2004-10-02T00:00:00',
 '2004-10-03T00:00:00',
 '2005-07-29T00:00:00',
 '2005-07-30T00:00:00',
 '2005-07-31T00:00:00',
 '2005-09-03T00:00:00',
 '2005-09-04T00:00:00',
 '2005-09-05T00:00:00',
 '2005-09-06T00:00:00',
 '2005-09-07T00:00:00',
 '2005-09-08T00:00:00',
 '2005-10-15T00:00:00',
 '2005-10-16T00:00:00',
 '2005-10-27T00:00:00',
 '2005-10-28T00:00:00',
 '2005-10-29T00:00:00',
 '2005-10-30T00:00:00',
 '2005-10-31T00:00:00',
 '2005-11-01T00:00:00',
 '2005-11-02T00:00:00',
 '2005-11-03T00:00:00',
 '2005-11-04T00:00:00',
 '2006-09-21T00:00:00',
 '2006-09-22T00:00:00',
 '2006-09-23T00:00:00',
 '2006-09-24T00:00:00',
 '2006-09-25T00:00:00',
 '2006-09-26T00:00:00',
 '2006-09-27T00:00:00',
 '2006-09-28T00:00:00',
 '2006-09-29T00:00:00',
 '2006-11-27T00:00:00',
 '2006-11-28T00:00:00',
 '2007-06-17T00:00:00',
 '2007-06-18T00:00:00',
 '2007-06-19T00:00:00',
 '2007-06-20T00:00:00',
 '2007-06-21T00:00:00',
 '2007-06-22T00:00:00',
 '2007-09-08T00:00:00',
 '2007-10-29T00:00:00',
 '2008-06-10T00:00:00',
 '2008-06-11T00:00:00',
 '2008-06-12T00:00:00',
 '2008-08-01T00:00:00',
 '2008-08-02T00:00:00',
 '2008-09-06T00:00:00',
 '2008-09-07T00:00:00',
 '2008-09-08T00:00:00',
 '2008-09-09T00:00:00',
 '2008-09-10T00:00:00',
 '2008-10-20T00:00:00',
 '2008-10-21T00:00:00',
 '2008-10-22T00:00:00',
 '2008-10-23T00:00:00',
 '2008-10-24T00:00:00',
 '2008-10-25T00:00:00',
 '2008-10-26T00:00:00',
 '2008-10-27T00:00:00',
 '2008-11-09T00:00:00',
 '2008-11-10T00:00:00',
 '2008-11-11T00:00:00',
 '2008-11-12T00:00:00',
 '2008-11-13T00:00:00',
 '2009-07-13T00:00:00',
 '2009-07-14T00:00:00',
 '2009-09-01T00:00:00',
 '2009-09-02T00:00:00',
 '2009-09-03T00:00:00',
 '2009-09-04T00:00:00',
 '2009-09-05T00:00:00',
 '2009-09-06T00:00:00',
 '2009-09-07T00:00:00',
 '2009-09-08T00:00:00',
 '2009-09-10T00:00:00',
 '2009-09-11T00:00:00',
 '2009-09-17T00:00:00',
 '2009-09-18T00:00:00',
 '2009-09-19T00:00:00',
 '2009-09-20T00:00:00',
 '2009-09-21T00:00:00',
 '2009-09-22T00:00:00',
 '2009-09-23T00:00:00',
 '2009-09-24T00:00:00',
 '2009-09-30T00:00:00',
 '2009-10-01T00:00:00',
 '2009-10-02T00:00:00',
 '2009-10-03T00:00:00',
 '2009-10-04T00:00:00',
 '2009-10-05T00:00:00',
 '2009-10-06T00:00:00',
 '2009-11-16T00:00:00',
 '2009-11-17T00:00:00',
 '2009-11-18T00:00:00',
 '2009-11-19T00:00:00',
 '2009-11-20T00:00:00',
 '2010-10-05T00:00:00',
 '2010-10-06T00:00:00',
 '2010-10-07T00:00:00',
 '2011-06-09T00:00:00',
 '2011-06-10T00:00:00',
 '2011-07-20T00:00:00',
 '2011-08-03T00:00:00',
 '2011-11-02T00:00:00',
 '2011-11-03T00:00:00',
 '2011-11-04T00:00:00',
 '2011-11-05T00:00:00',
 '2012-09-05T00:00:00',
 '2012-09-06T00:00:00',
 '2012-09-07T00:00:00',
 '2012-09-08T00:00:00',
 '2012-09-09T00:00:00',
 '2012-09-10T00:00:00',
 '2012-09-11T00:00:00',
 '2012-09-20T00:00:00',
 '2012-11-06T00:00:00',
 '2013-06-11T00:00:00',
 '2013-06-12T00:00:00',
 '2013-06-13T00:00:00',
 '2013-06-14T00:00:00',
 '2013-09-14T00:00:00',
 '2013-09-16T00:00:00',
 '2013-11-09T00:00:00',
 '2014-06-10T00:00:00',
 '2014-06-11T00:00:00',
 '2014-06-12T00:00:00',
 '2014-06-13T00:00:00',
 '2014-06-14T00:00:00',
 '2014-06-15T00:00:00',
 '2014-06-16T00:00:00',
 '2014-06-17T00:00:00',
 '2014-07-07T00:00:00',
 '2014-07-08T00:00:00',
 '2014-09-01T00:00:00',
 '2014-09-02T00:00:00',
 '2014-09-03T00:00:00',
 '2014-09-04T00:00:00',
 '2014-09-05T00:00:00',
 '2014-09-06T00:00:00',
 '2014-11-28T00:00:00',
 '2014-11-29T00:00:00',
 '2014-11-30T00:00:00',
 '2015-06-16T00:00:00',
 '2015-06-17T00:00:00',
 '2015-06-18T00:00:00',
 '2015-06-19T00:00:00',
 '2015-06-20T00:00:00',
 '2015-06-21T00:00:00',
 '2015-11-30T00:00:00',
 '2016-06-20T00:00:00',
 '2016-06-21T00:00:00',
 '2016-06-22T00:00:00',
 '2016-06-23T00:00:00',
 '2016-06-24T00:00:00',
 '2016-06-25T00:00:00',
 '2016-06-26T00:00:00',
 '2016-09-12T00:00:00',
 '2016-09-13T00:00:00',
 '2016-09-14T00:00:00',
 '2016-09-15T00:00:00',
 '2016-09-16T00:00:00',
 '2016-09-17T00:00:00',
 '2016-09-18T00:00:00',
 '2017-06-01T00:00:00',
 '2017-06-02T00:00:00',
 '2017-07-20T00:00:00',
 '2017-07-21T00:00:00',
 '2017-07-22T00:00:00',
 '2017-08-21T00:00:00',
 '2017-09-14T00:00:00',
 '2017-10-08T00:00:00',
 '2017-10-09T00:00:00',
 '2017-10-10T00:00:00',
 '2017-10-11T00:00:00',
 '2017-10-12T00:00:00',
 '2017-11-25T00:00:00',
 '2017-11-26T00:00:00',
 '2017-11-27T00:00:00',
 '2017-11-29T00:00:00',
 '2017-11-30T00:00:00',
 '2018-07-08T00:00:00',
 '2018-07-09T00:00:00',
 '2018-07-10T00:00:00',
 '2018-11-10T00:00:00',
 '2018-11-11T00:00:00',
 '2018-11-12T00:00:00',
 '2019-06-13T00:00:00',
 '2019-06-14T00:00:00',
 '2019-06-15T00:00:00',
 '2019-06-17T00:00:00',
 '2019-08-06T00:00:00',
 '2019-08-07T00:00:00',
 '2019-08-08T00:00:00',
 '2019-08-09T00:00:00',
 '2019-08-10T00:00:00',
 '2019-08-11T00:00:00',
 '2019-08-25T00:00:00',
 '2019-08-26T00:00:00',
 '2019-08-27T00:00:00',
 '2019-08-30T00:00:00',
 '2019-10-31T00:00:00',
 '2020-07-30T00:00:00',
 '2020-07-31T00:00:00',
 '2020-08-01T00:00:00',
 '2020-08-03T00:00:00',
 '2020-08-04T00:00:00',
 '2020-08-05T00:00:00',
 '2020-08-06T00:00:00',
 '2020-09-02T00:00:00',
 '2020-09-03T00:00:00',
 '2020-09-15T00:00:00',
 '2020-09-16T00:00:00',
 '2020-09-17T00:00:00',
 '2020-09-18T00:00:00',
 '2020-10-08T00:00:00',
 '2020-10-09T00:00:00',
 '2020-10-10T00:00:00',
 '2020-10-11T00:00:00',
 '2020-10-12T00:00:00',
 '2020-10-13T00:00:00',
 '2020-11-25T00:00:00',
 '2020-11-26T00:00:00',
 '2020-11-27T00: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()