#!/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
print('u200...')
ifiles = '/projects/w/wenchang/data/era5/analysis_wy/plevels/u200/daily/era5.u200.daily.*.nc'
u200 = xr.open_mfdataset(ifiles)['u']

print('u850...')
ifiles = '/projects/w/wenchang/data/era5/analysis_wy/plevels/u850/daily/era5.u850.daily.*.nc'
u850 = xr.open_mfdataset(ifiles)['u']
 
print('v200...')
ifiles = '/projects/w/wenchang/data/era5/analysis_wy/plevels/v200/daily/era5.v200.daily.*.nc'
v200 = xr.open_mfdataset(ifiles)['v']

print('v850...')
ifiles = '/projects/w/wenchang/data/era5/analysis_wy/plevels/v850/daily/era5.v850.daily.*.nc'
v850 = xr.open_mfdataset(ifiles)['v']

#time-mean for mjo phase 1
ofile = f'era5.wind_shear.daily.1979-2020.phase1.nc'
time = [
 '1979-06-01T00:00:00',
 '1979-06-02T00:00:00',
 '1979-06-03T00:00:00',
 '1979-06-04T00:00:00',
 '1979-06-05T00:00:00',
 '1979-06-06T00:00:00',
 '1979-07-12T00:00:00',
 '1979-07-13T00:00:00',
 '1979-07-14T00:00:00',
 '1979-07-16T00:00:00',
 '1979-07-19T00:00:00',
 '1979-07-20T00:00:00',
 '1979-07-21T00:00:00',
 '1979-07-22T00:00:00',
 '1979-07-23T00:00:00',
 '1979-07-24T00:00:00',
 '1979-08-22T00:00:00',
 '1979-08-23T00:00:00',
 '1979-08-24T00:00:00',
 '1979-08-25T00:00:00',
 '1979-08-26T00:00:00',
 '1979-08-27T00:00:00',
 '1979-08-28T00:00:00',
 '1979-08-29T00:00:00',
 '1979-08-30T00:00:00',
 '1979-08-31T00:00:00',
 '1979-09-01T00:00:00',
 '1979-09-02T00:00:00',
 '1979-09-03T00:00:00',
 '1979-09-04T00:00:00',
 '1979-09-05T00:00:00',
 '1979-09-06T00:00:00',
 '1979-09-07T00:00:00',
 '1979-09-10T00:00:00',
 '1979-09-11T00:00:00',
 '1979-09-13T00:00:00',
 '1979-09-14T00:00:00',
 '1979-09-15T00:00:00',
 '1979-10-23T00:00:00',
 '1979-10-24T00:00:00',
 '1979-10-25T00:00:00',
 '1979-10-26T00:00:00',
 '1979-10-27T00:00:00',
 '1979-10-28T00:00:00',
 '1979-10-29T00:00:00',
 '1979-10-30T00:00:00',
 '1980-06-01T00:00:00',
 '1980-06-02T00:00:00',
 '1980-06-03T00:00:00',
 '1980-07-27T00:00:00',
 '1980-07-28T00:00:00',
 '1980-07-29T00:00:00',
 '1980-07-30T00:00:00',
 '1980-07-31T00:00:00',
 '1980-08-01T00:00:00',
 '1980-08-02T00:00:00',
 '1980-08-03T00:00:00',
 '1980-08-07T00:00:00',
 '1980-08-08T00:00:00',
 '1980-08-09T00:00:00',
 '1980-08-10T00:00:00',
 '1980-08-11T00:00:00',
 '1980-08-12T00:00:00',
 '1980-09-19T00:00:00',
 '1980-09-20T00:00:00',
 '1980-09-21T00:00:00',
 '1980-09-22T00:00:00',
 '1980-09-23T00:00:00',
 '1980-09-24T00:00:00',
 '1980-11-10T00:00:00',
 '1980-11-11T00:00:00',
 '1980-11-12T00:00:00',
 '1981-08-24T00:00:00',
 '1981-08-25T00:00:00',
 '1981-08-26T00:00:00',
 '1981-08-27T00:00:00',
 '1981-08-28T00:00:00',
 '1981-08-29T00:00:00',
 '1981-08-30T00:00:00',
 '1981-10-07T00:00:00',
 '1981-10-08T00:00:00',
 '1981-10-09T00:00:00',
 '1982-09-18T00:00:00',
 '1982-09-19T00:00:00',
 '1982-09-20T00:00:00',
 '1982-09-21T00:00:00',
 '1982-09-28T00:00:00',
 '1982-09-29T00:00:00',
 '1982-10-22T00:00:00',
 '1982-10-23T00:00:00',
 '1982-10-24T00:00:00',
 '1982-11-22T00:00:00',
 '1982-11-23T00:00:00',
 '1982-11-24T00:00:00',
 '1983-09-03T00:00:00',
 '1983-09-04T00:00:00',
 '1983-09-05T00:00:00',
 '1983-10-15T00:00:00',
 '1984-11-01T00:00:00',
 '1984-11-02T00:00:00',
 '1984-11-03T00:00:00',
 '1984-11-04T00:00:00',
 '1984-11-05T00:00:00',
 '1984-11-06T00:00:00',
 '1984-11-07T00:00:00',
 '1984-11-08T00:00:00',
 '1984-11-09T00:00:00',
 '1984-11-10T00:00:00',
 '1984-11-11T00:00:00',
 '1985-07-17T00:00:00',
 '1985-07-18T00:00:00',
 '1985-07-19T00:00:00',
 '1985-07-20T00:00:00',
 '1985-07-21T00:00:00',
 '1985-07-22T00:00:00',
 '1985-09-04T00:00:00',
 '1985-10-30T00:00:00',
 '1985-10-31T00:00:00',
 '1985-11-01T00:00:00',
 '1985-11-02T00:00:00',
 '1986-06-01T00:00:00',
 '1986-06-02T00:00:00',
 '1986-08-28T00:00:00',
 '1986-09-25T00:00:00',
 '1986-09-26T00:00:00',
 '1986-09-27T00:00:00',
 '1986-09-28T00:00:00',
 '1986-11-25T00:00:00',
 '1986-11-26T00:00:00',
 '1986-11-27T00:00:00',
 '1986-11-28T00:00:00',
 '1986-11-29T00:00:00',
 '1986-11-30T00:00:00',
 '1987-07-22T00:00:00',
 '1987-07-23T00:00:00',
 '1987-07-24T00:00:00',
 '1987-07-25T00:00:00',
 '1987-07-26T00:00:00',
 '1987-07-27T00:00:00',
 '1987-07-28T00:00:00',
 '1987-07-29T00:00:00',
 '1987-07-30T00:00:00',
 '1987-07-31T00:00:00',
 '1987-08-01T00:00:00',
 '1987-08-02T00:00:00',
 '1987-08-03T00:00:00',
 '1987-08-04T00:00:00',
 '1988-06-19T00:00:00',
 '1988-06-20T00:00:00',
 '1988-06-21T00:00:00',
 '1988-06-22T00:00:00',
 '1988-06-25T00:00:00',
 '1988-06-26T00:00:00',
 '1988-06-27T00:00:00',
 '1988-06-28T00:00:00',
 '1988-06-29T00:00:00',
 '1988-06-30T00:00:00',
 '1988-07-01T00:00:00',
 '1988-07-02T00:00:00',
 '1988-07-03T00:00:00',
 '1988-07-04T00:00:00',
 '1988-07-05T00:00:00',
 '1988-08-11T00:00:00',
 '1988-08-12T00:00:00',
 '1988-08-13T00:00:00',
 '1988-08-14T00:00:00',
 '1988-08-21T00:00:00',
 '1988-08-24T00:00:00',
 '1988-08-25T00:00:00',
 '1988-08-26T00:00:00',
 '1988-08-27T00:00:00',
 '1988-08-28T00:00:00',
 '1988-08-29T00:00:00',
 '1988-08-30T00:00:00',
 '1988-08-31T00:00:00',
 '1988-09-01T00:00:00',
 '1988-09-02T00:00:00',
 '1989-09-02T00:00:00',
 '1989-09-03T00:00:00',
 '1989-09-04T00:00:00',
 '1990-07-02T00:00:00',
 '1990-07-03T00:00:00',
 '1990-07-04T00:00:00',
 '1990-07-05T00:00:00',
 '1990-07-06T00:00:00',
 '1990-07-07T00:00:00',
 '1990-07-08T00:00:00',
 '1990-07-09T00:00:00',
 '1990-07-10T00:00:00',
 '1990-07-11T00:00:00',
 '1990-10-09T00:00:00',
 '1990-10-10T00:00:00',
 '1990-10-11T00:00:00',
 '1990-10-12T00:00:00',
 '1990-10-13T00:00:00',
 '1990-10-14T00:00:00',
 '1990-10-15T00:00:00',
 '1990-10-16T00:00:00',
 '1990-10-17T00:00:00',
 '1990-10-22T00:00:00',
 '1990-10-23T00:00:00',
 '1990-10-24T00:00:00',
 '1991-06-25T00:00:00',
 '1991-06-26T00:00:00',
 '1991-06-27T00:00:00',
 '1991-06-28T00:00:00',
 '1991-06-29T00:00:00',
 '1991-06-30T00:00:00',
 '1991-10-03T00:00:00',
 '1991-10-04T00:00:00',
 '1991-10-05T00:00:00',
 '1991-10-06T00:00:00',
 '1991-10-07T00:00:00',
 '1992-07-02T00:00:00',
 '1992-07-03T00:00:00',
 '1992-07-04T00:00:00',
 '1992-07-05T00:00:00',
 '1992-07-06T00:00:00',
 '1992-07-07T00:00:00',
 '1992-07-08T00:00:00',
 '1992-07-09T00:00:00',
 '1992-09-22T00:00:00',
 '1992-11-09T00:00:00',
 '1992-11-10T00:00:00',
 '1992-11-11T00:00:00',
 '1992-11-12T00:00:00',
 '1992-11-13T00:00:00',
 '1992-11-27T00:00:00',
 '1992-11-28T00:00:00',
 '1992-11-29T00:00:00',
 '1992-11-30T00:00:00',
 '1993-06-24T00:00:00',
 '1993-06-25T00:00:00',
 '1993-06-26T00:00:00',
 '1993-06-27T00:00:00',
 '1993-06-28T00:00:00',
 '1993-06-29T00:00:00',
 '1993-06-30T00:00:00',
 '1993-07-01T00:00:00',
 '1993-07-02T00:00:00',
 '1993-08-20T00:00:00',
 '1993-10-30T00:00:00',
 '1993-10-31T00:00:00',
 '1993-11-01T00:00:00',
 '1993-11-02T00:00:00',
 '1994-08-07T00:00:00',
 '1994-08-08T00:00:00',
 '1994-08-09T00:00:00',
 '1994-08-10T00:00:00',
 '1994-08-11T00:00:00',
 '1994-08-12T00:00:00',
 '1994-08-13T00:00:00',
 '1994-08-20T00:00:00',
 '1994-10-03T00:00:00',
 '1994-10-04T00:00:00',
 '1994-10-05T00:00:00',
 '1994-10-06T00:00:00',
 '1994-10-07T00:00:00',
 '1994-10-08T00:00:00',
 '1994-10-09T00:00:00',
 '1994-10-10T00:00:00',
 '1994-10-11T00:00:00',
 '1994-10-12T00:00:00',
 '1994-10-13T00:00:00',
 '1994-10-14T00:00:00',
 '1994-10-15T00:00:00',
 '1994-10-16T00:00:00',
 '1994-10-17T00:00:00',
 '1994-10-18T00:00:00',
 '1994-10-19T00:00:00',
 '1994-10-20T00:00:00',
 '1994-11-01T00:00:00',
 '1994-11-02T00:00:00',
 '1994-11-03T00:00:00',
 '1994-11-04T00:00:00',
 '1994-11-05T00:00:00',
 '1994-11-06T00:00:00',
 '1994-11-07T00:00:00',
 '1995-06-24T00:00:00',
 '1995-06-25T00:00:00',
 '1995-07-17T00:00:00',
 '1995-07-18T00:00:00',
 '1995-07-19T00:00:00',
 '1995-07-20T00:00:00',
 '1995-07-21T00:00:00',
 '1995-07-22T00:00:00',
 '1995-07-23T00:00:00',
 '1995-07-24T00:00:00',
 '1995-08-09T00:00:00',
 '1995-08-10T00:00:00',
 '1995-08-11T00:00:00',
 '1995-08-12T00:00:00',
 '1995-08-13T00:00:00',
 '1995-08-16T00:00:00',
 '1995-08-17T00:00:00',
 '1995-08-19T00:00:00',
 '1996-06-29T00:00:00',
 '1996-06-30T00:00:00',
 '1996-07-01T00:00:00',
 '1996-07-02T00:00:00',
 '1996-07-03T00:00:00',
 '1996-07-04T00:00:00',
 '1996-08-11T00:00:00',
 '1996-08-12T00:00:00',
 '1996-08-14T00:00:00',
 '1996-08-15T00:00:00',
 '1996-10-05T00:00:00',
 '1996-10-06T00:00:00',
 '1996-10-07T00:00:00',
 '1996-10-08T00:00:00',
 '1996-11-19T00:00:00',
 '1996-11-20T00:00:00',
 '1996-11-21T00:00:00',
 '1997-06-12T00:00:00',
 '1997-06-13T00:00:00',
 '1997-06-14T00:00:00',
 '1997-06-15T00:00:00',
 '1997-09-05T00:00:00',
 '1997-10-13T00:00:00',
 '1997-10-14T00:00:00',
 '1997-10-15T00:00:00',
 '1997-10-16T00:00:00',
 '1997-10-17T00:00:00',
 '1997-11-08T00:00:00',
 '1997-11-09T00:00:00',
 '1997-11-10T00:00:00',
 '1997-11-11T00:00:00',
 '1997-11-30T00:00:00',
 '1998-07-14T00:00:00',
 '1998-07-15T00:00:00',
 '1998-08-24T00:00:00',
 '1998-08-25T00:00:00',
 '1998-08-26T00:00:00',
 '1998-08-27T00:00:00',
 '1999-06-28T00:00:00',
 '1999-06-29T00:00:00',
 '1999-06-30T00:00:00',
 '1999-07-01T00:00:00',
 '1999-07-02T00:00:00',
 '1999-07-03T00:00:00',
 '1999-07-04T00:00:00',
 '1999-07-05T00:00:00',
 '1999-07-06T00:00:00',
 '1999-08-17T00:00:00',
 '1999-08-18T00:00:00',
 '1999-08-19T00:00:00',
 '1999-08-20T00:00:00',
 '1999-08-21T00:00:00',
 '1999-09-19T00:00:00',
 '1999-09-20T00:00:00',
 '1999-09-21T00:00:00',
 '1999-09-22T00:00:00',
 '1999-09-23T00:00:00',
 '1999-09-24T00:00:00',
 '1999-09-25T00:00:00',
 '1999-09-26T00:00:00',
 '1999-09-27T00:00:00',
 '2000-06-18T00:00:00',
 '2000-06-19T00:00:00',
 '2000-06-20T00:00:00',
 '2000-06-21T00:00:00',
 '2000-07-26T00:00:00',
 '2000-07-27T00:00:00',
 '2000-07-28T00:00:00',
 '2000-07-29T00:00:00',
 '2000-07-30T00:00:00',
 '2000-07-31T00:00:00',
 '2000-08-01T00:00:00',
 '2000-08-02T00:00:00',
 '2000-08-03T00:00:00',
 '2000-08-04T00:00:00',
 '2000-08-05T00:00:00',
 '2000-08-06T00:00:00',
 '2000-09-13T00:00:00',
 '2000-09-14T00:00:00',
 '2000-09-17T00:00:00',
 '2000-09-18T00:00:00',
 '2000-09-22T00:00:00',
 '2000-09-23T00:00:00',
 '2000-11-10T00:00:00',
 '2000-11-11T00:00:00',
 '2001-06-01T00:00:00',
 '2001-07-16T00:00:00',
 '2001-07-17T00:00:00',
 '2001-07-18T00:00:00',
 '2001-07-19T00:00:00',
 '2001-07-20T00:00:00',
 '2001-07-21T00:00:00',
 '2001-07-22T00:00:00',
 '2001-07-23T00:00:00',
 '2001-09-06T00:00:00',
 '2001-09-08T00:00:00',
 '2001-09-09T00:00:00',
 '2001-09-11T00:00:00',
 '2001-09-15T00:00:00',
 '2001-10-27T00:00:00',
 '2002-06-01T00:00:00',
 '2002-06-02T00:00:00',
 '2002-06-03T00:00:00',
 '2002-06-04T00:00:00',
 '2002-10-26T00:00:00',
 '2002-10-27T00:00:00',
 '2002-10-28T00:00:00',
 '2003-06-10T00:00:00',
 '2003-06-11T00:00:00',
 '2003-06-12T00:00:00',
 '2003-06-13T00:00:00',
 '2003-06-14T00:00:00',
 '2003-06-15T00:00:00',
 '2003-06-16T00:00:00',
 '2003-06-17T00:00:00',
 '2003-06-18T00:00:00',
 '2003-06-19T00:00:00',
 '2003-06-20T00:00:00',
 '2003-06-21T00:00:00',
 '2003-06-22T00:00:00',
 '2003-06-26T00:00:00',
 '2004-07-06T00:00:00',
 '2004-07-07T00:00:00',
 '2004-07-08T00:00:00',
 '2004-07-09T00:00:00',
 '2004-07-10T00:00:00',
 '2004-07-11T00:00:00',
 '2004-07-12T00:00:00',
 '2004-07-13T00:00:00',
 '2004-07-14T00:00:00',
 '2004-07-15T00:00:00',
 '2004-07-16T00:00:00',
 '2004-07-17T00:00:00',
 '2004-07-18T00:00:00',
 '2004-07-19T00:00:00',
 '2004-07-20T00:00:00',
 '2004-07-21T00:00:00',
 '2004-07-22T00:00:00',
 '2004-07-23T00:00:00',
 '2004-07-24T00:00:00',
 '2004-09-02T00:00:00',
 '2004-09-03T00:00:00',
 '2004-09-04T00:00:00',
 '2004-10-21T00:00:00',
 '2004-10-22T00:00:00',
 '2004-10-23T00:00:00',
 '2004-10-24T00:00:00',
 '2004-10-25T00:00:00',
 '2005-07-10T00:00:00',
 '2005-07-11T00:00:00',
 '2005-07-12T00:00:00',
 '2005-07-13T00:00:00',
 '2005-07-14T00:00:00',
 '2005-07-15T00:00:00',
 '2005-08-13T00:00:00',
 '2005-08-14T00:00:00',
 '2005-08-15T00:00:00',
 '2005-08-16T00:00:00',
 '2005-08-18T00:00:00',
 '2005-08-19T00:00:00',
 '2005-10-05T00:00:00',
 '2005-10-07T00:00:00',
 '2005-10-08T00:00:00',
 '2006-06-03T00:00:00',
 '2006-06-04T00:00:00',
 '2006-06-05T00:00:00',
 '2006-06-06T00:00:00',
 '2006-06-07T00:00:00',
 '2006-06-08T00:00:00',
 '2006-06-09T00:00:00',
 '2006-06-10T00:00:00',
 '2006-06-11T00:00:00',
 '2006-06-12T00:00:00',
 '2006-06-13T00:00:00',
 '2006-06-14T00:00:00',
 '2006-07-25T00:00:00',
 '2006-07-26T00:00:00',
 '2006-07-27T00:00:00',
 '2006-08-23T00:00:00',
 '2006-08-24T00:00:00',
 '2006-08-25T00:00:00',
 '2006-08-26T00:00:00',
 '2006-08-27T00:00:00',
 '2006-09-03T00:00:00',
 '2006-09-06T00:00:00',
 '2006-10-18T00:00:00',
 '2006-10-19T00:00:00',
 '2006-10-20T00:00:00',
 '2006-10-21T00:00:00',
 '2006-10-22T00:00:00',
 '2006-10-24T00:00:00',
 '2006-11-08T00:00:00',
 '2006-11-12T00:00:00',
 '2006-11-13T00:00:00',
 '2006-11-14T00:00:00',
 '2006-11-15T00:00:00',
 '2006-11-16T00:00:00',
 '2006-11-17T00:00:00',
 '2007-06-01T00:00:00',
 '2007-06-02T00:00:00',
 '2007-06-03T00:00:00',
 '2007-06-04T00:00:00',
 '2007-06-05T00:00:00',
 '2007-06-06T00:00:00',
 '2007-07-15T00:00:00',
 '2007-07-16T00:00:00',
 '2007-07-17T00:00:00',
 '2007-07-18T00:00:00',
 '2007-07-19T00:00:00',
 '2007-07-20T00:00:00',
 '2007-07-21T00:00:00',
 '2007-08-24T00:00:00',
 '2007-08-25T00:00:00',
 '2007-08-26T00:00:00',
 '2007-10-18T00:00:00',
 '2007-10-19T00:00:00',
 '2007-10-20T00:00:00',
 '2007-10-21T00:00:00',
 '2008-06-30T00:00:00',
 '2008-07-01T00:00:00',
 '2008-07-02T00:00:00',
 '2008-07-03T00:00:00',
 '2008-07-04T00:00:00',
 '2008-07-05T00:00:00',
 '2008-07-06T00:00:00',
 '2008-07-07T00:00:00',
 '2008-07-08T00:00:00',
 '2008-07-09T00:00:00',
 '2008-07-10T00:00:00',
 '2008-07-11T00:00:00',
 '2008-07-12T00:00:00',
 '2008-07-13T00:00:00',
 '2008-07-14T00:00:00',
 '2009-06-24T00:00:00',
 '2009-06-25T00:00:00',
 '2009-06-26T00:00:00',
 '2009-06-27T00:00:00',
 '2009-08-11T00:00:00',
 '2009-08-12T00:00:00',
 '2009-08-13T00:00:00',
 '2009-08-14T00:00:00',
 '2009-08-19T00:00:00',
 '2009-08-20T00:00:00',
 '2009-08-21T00:00:00',
 '2009-08-22T00:00:00',
 '2009-08-23T00:00:00',
 '2009-10-29T00:00:00',
 '2009-10-30T00:00:00',
 '2009-10-31T00:00:00',
 '2009-11-01T00:00:00',
 '2009-11-02T00:00:00',
 '2009-11-03T00:00:00',
 '2010-06-01T00:00:00',
 '2010-06-02T00:00:00',
 '2010-06-17T00:00:00',
 '2010-06-18T00:00:00',
 '2010-06-19T00:00:00',
 '2010-06-30T00:00:00',
 '2010-07-01T00:00:00',
 '2010-07-02T00:00:00',
 '2010-07-03T00:00:00',
 '2010-07-04T00:00:00',
 '2010-07-05T00:00:00',
 '2010-07-06T00:00:00',
 '2010-07-12T00:00:00',
 '2010-07-13T00:00:00',
 '2010-07-14T00:00:00',
 '2010-07-15T00:00:00',
 '2010-07-16T00:00:00',
 '2010-09-20T00:00:00',
 '2010-09-21T00:00:00',
 '2010-09-22T00:00:00',
 '2010-11-03T00:00:00',
 '2011-07-10T00:00:00',
 '2011-07-11T00:00:00',
 '2011-07-12T00:00:00',
 '2011-07-13T00:00:00',
 '2011-07-14T00:00:00',
 '2011-07-15T00:00:00',
 '2011-08-16T00:00:00',
 '2011-08-21T00:00:00',
 '2011-08-22T00:00:00',
 '2011-08-23T00:00:00',
 '2011-10-15T00:00:00',
 '2011-10-16T00:00:00',
 '2011-10-17T00:00:00',
 '2011-10-18T00:00:00',
 '2011-10-19T00:00:00',
 '2011-11-17T00:00:00',
 '2011-11-18T00:00:00',
 '2011-11-19T00:00:00',
 '2011-11-20T00:00:00',
 '2012-06-17T00:00:00',
 '2012-06-18T00:00:00',
 '2012-06-19T00:00:00',
 '2012-06-20T00:00:00',
 '2012-06-21T00:00:00',
 '2012-06-23T00:00:00',
 '2012-06-24T00:00:00',
 '2012-06-25T00:00:00',
 '2012-06-26T00:00:00',
 '2012-06-27T00:00:00',
 '2012-07-12T00:00:00',
 '2012-07-13T00:00:00',
 '2012-08-09T00:00:00',
 '2012-08-10T00:00:00',
 '2012-08-11T00:00:00',
 '2012-08-12T00:00:00',
 '2012-08-13T00:00:00',
 '2012-08-14T00:00:00',
 '2012-08-15T00:00:00',
 '2012-08-16T00:00:00',
 '2012-10-17T00:00:00',
 '2012-10-18T00:00:00',
 '2012-10-19T00:00:00',
 '2012-10-20T00:00:00',
 '2012-10-21T00:00:00',
 '2013-06-29T00:00:00',
 '2013-06-30T00:00:00',
 '2013-07-01T00:00:00',
 '2013-07-02T00:00:00',
 '2013-07-03T00:00:00',
 '2013-07-04T00:00:00',
 '2013-07-05T00:00:00',
 '2013-07-06T00:00:00',
 '2013-07-07T00:00:00',
 '2013-07-08T00:00:00',
 '2013-07-09T00:00:00',
 '2013-08-26T00:00:00',
 '2013-08-27T00:00:00',
 '2013-08-28T00:00:00',
 '2013-08-29T00:00:00',
 '2013-08-30T00:00:00',
 '2013-08-31T00:00:00',
 '2013-09-01T00:00:00',
 '2013-09-02T00:00:00',
 '2013-09-03T00:00:00',
 '2013-09-04T00:00:00',
 '2013-09-05T00:00:00',
 '2013-11-01T00:00:00',
 '2013-11-02T00:00:00',
 '2014-08-02T00:00:00',
 '2014-08-04T00:00:00',
 '2014-08-05T00:00:00',
 '2014-08-14T00:00:00',
 '2014-08-15T00:00:00',
 '2014-08-16T00:00:00',
 '2014-08-17T00:00:00',
 '2014-08-18T00:00:00',
 '2014-08-19T00:00:00',
 '2014-08-20T00:00:00',
 '2014-08-21T00:00:00',
 '2014-08-22T00:00:00',
 '2014-10-15T00:00:00',
 '2014-10-16T00:00:00',
 '2014-10-17T00:00:00',
 '2014-10-18T00:00:00',
 '2014-10-19T00:00:00',
 '2014-10-20T00:00:00',
 '2014-10-21T00:00:00',
 '2014-10-22T00:00:00',
 '2014-10-23T00:00:00',
 '2014-10-24T00:00:00',
 '2014-10-25T00:00:00',
 '2014-11-11T00:00:00',
 '2014-11-12T00:00:00',
 '2014-11-13T00:00:00',
 '2014-11-14T00:00:00',
 '2014-11-15T00:00:00',
 '2014-11-16T00:00:00',
 '2014-11-17T00:00:00',
 '2015-06-01T00:00:00',
 '2015-06-02T00:00:00',
 '2015-06-03T00:00:00',
 '2015-06-04T00:00:00',
 '2015-06-05T00:00:00',
 '2015-07-18T00:00:00',
 '2015-07-27T00:00:00',
 '2015-07-28T00:00:00',
 '2015-07-29T00:00:00',
 '2015-08-24T00:00:00',
 '2015-08-25T00:00:00',
 '2015-09-02T00:00:00',
 '2015-09-03T00:00:00',
 '2015-09-04T00:00:00',
 '2015-09-05T00:00:00',
 '2015-09-06T00:00:00',
 '2015-10-21T00:00:00',
 '2015-10-22T00:00:00',
 '2016-07-16T00:00:00',
 '2016-07-17T00:00:00',
 '2016-07-18T00:00:00',
 '2016-11-16T00:00:00',
 '2016-11-17T00:00:00',
 '2016-11-18T00:00:00',
 '2017-06-11T00:00:00',
 '2017-06-12T00:00:00',
 '2017-06-13T00:00:00',
 '2017-06-14T00:00:00',
 '2017-06-15T00:00:00',
 '2017-06-16T00:00:00',
 '2017-06-17T00:00:00',
 '2017-06-18T00:00:00',
 '2017-06-19T00:00:00',
 '2018-06-17T00:00:00',
 '2018-06-18T00:00:00',
 '2018-06-19T00:00:00',
 '2018-06-20T00:00:00',
 '2018-06-21T00:00:00',
 '2018-06-22T00:00:00',
 '2018-06-23T00:00:00',
 '2018-06-27T00:00:00',
 '2018-06-28T00:00:00',
 '2018-06-29T00:00:00',
 '2018-10-02T00:00:00',
 '2018-10-03T00:00:00',
 '2018-10-04T00:00:00',
 '2018-10-05T00:00:00',
 '2018-10-06T00:00:00',
 '2018-10-07T00:00:00',
 '2018-10-08T00:00:00',
 '2018-10-09T00:00:00',
 '2018-10-31T00:00:00',
 '2018-11-01T00:00:00',
 '2018-11-02T00:00:00',
 '2019-07-04T00:00:00',
 '2019-07-05T00:00:00',
 '2019-07-06T00:00:00',
 '2019-07-07T00:00:00',
 '2019-07-08T00:00:00',
 '2019-07-09T00:00:00',
 '2019-07-11T00:00:00',
 '2019-07-20T00:00:00',
 '2019-07-21T00:00:00',
 '2019-09-21T00:00:00',
 '2019-09-22T00:00:00',
 '2019-09-23T00:00:00',
 '2019-09-24T00:00:00',
 '2019-09-25T00:00:00',
 '2019-09-26T00:00:00',
 '2019-09-27T00:00:00',
 '2019-09-28T00:00:00',
 '2019-09-29T00:00:00',
 '2019-09-30T00:00:00',
 '2019-10-01T00:00:00',
 '2019-10-02T00:00:00',
 '2019-10-03T00:00:00',
 '2019-10-04T00:00:00',
 '2019-10-05T00:00:00',
 '2019-10-06T00:00:00',
 '2019-10-07T00:00:00',
 '2019-10-08T00:00:00',
 '2019-10-13T00:00:00',
 '2019-10-14T00:00:00',
 '2019-10-15T00:00:00',
 '2019-10-16T00:00:00',
 '2019-10-17T00:00:00',
 '2019-11-20T00:00:00',
 '2019-11-21T00:00:00',
 '2020-06-01T00:00:00',
 '2020-06-02T00:00:00',
 '2020-06-03T00:00:00',
 '2020-06-04T00:00:00',
 '2020-06-05T00:00:00',
 '2020-06-06T00:00:00',
 '2020-06-07T00:00:00',
 '2020-06-23T00:00:00',
 '2020-06-24T00:00:00',
 '2020-06-25T00:00:00',
 '2020-06-26T00:00:00',
 '2020-06-27T00:00:00',
 '2020-07-04T00:00:00',
 '2020-07-05T00:00:00',
 '2020-07-06T00:00:00',
 '2020-07-07T00:00:00',
 '2020-07-20T00:00:00',
 '2020-08-22T00:00:00',
 '2020-08-23T00:00:00',
 '2020-08-24T00:00:00',
 '2020-08-25T00:00:00',
 '2020-08-26T00:00:00',
 '2020-08-27T00:00:00',
 '2020-11-12T00:00:00',
 '2020-11-14T00:00:00']
if os.path.exists(ofile): 
    print('[exists]:', ofile)
else:
    print('wind_shear...')
    wind_shear = ((u200-u850)**2 + (v200 - v850)**2)**0.5
    print('phase1...')
    wind_shear = wind_shear.sel(time=time).load()
    print('wind_shear_mean...')
    wind_shear_mean = wind_shear.mean('time')
    print('wind_shear_mean save...')
    wind_shear_mean.to_dataset(name='wind_shear_mean').to_netcdf(ofile)
    print('[saved]:', ofile)

 
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()