def draw_12month_cycle(cb_min, cb_max, product_n = 3, Ave_or_Std = 'Ave', interval = 15):
xgrid, ygrid, zgrid = subroutine.product_grid_info('t', 'data', product_n)
xlim = [40, 110]
ylim = [-20, 30]
fy = 1990
ly = 2011
Data = np.zeros((12, ygrid.size, xgrid.size))
for i in range(0, 12):
month = i + 1
# データの取得
each_data = np.zeros((ygrid.size, xgrid.size, ly - fy + 1))
for year in range(fy, ly + 1):
k = year - fy
each_data[:, :, k] = load_data_of_npz(year, month, product_n = product_n, data_type = 0)
if Ave_or_Std == 'Ave':
Data[i, :, :] = np.average(each_data, axis = 2)
elif Ave_or_Std == 'Std':
Data[i, :, :] = np.std(each_data, axis = 2)
else:
raise Exception('your Ave_or_Std argument is not valid!')
plta = quick.draw_12month_cycle_with_map(Data, ygrid, xgrid, cb_min, cb_max, xlim = xlim, ylim = ylim, \
clabel = Ave_or_Std + " of T20 Depth(m)", interval = interval)
def draw_12month_cycle(Data, cb_min, cb_max, depthn = 1, product_n = 3, interval = 15, Ekman_or_Geost = 'Ekman', Vector_or_Scalar = 'Scalar'): xlim = [40, 110] ylim = [-20, 30] plta = quick.draw_12month_cycle_with_map(Data, ygrid, xgrid, cb_min, cb_max, xlim = xlim, ylim = ylim, \ clabel = fname + " Current Velocity(m/s) at" + str(zgrid[depthn - 1]) + 'm',\ interval = interval, my_color = 0) return plta
def draw_12month_cycle(cb_min, cb_max, depthn1 = 1, depthn2 = 4, product_n = 3, interval = 15):
_,title_name,_=subroutine.product_n_to_name(product_n)
xgrid, ygrid, zgrid = subroutine.product_grid_info('u', 'data', product_n)
xlim = [40, 110]
ylim = [-20, 30]
fy = 1990
ly = 2011
Data = np.zeros((12, ygrid.size, xgrid.size))
for i in range(0, 12):
month = i + 1
# データの取得
_, strmonth = subroutine.strym(2000, month)
u = CI.load_climatology_or_interannual_variation_as_npz(fy, ly, month, 'u', 'Ave', product_n = product_n)
v = CI.load_climatology_or_interannual_variation_as_npz(fy, ly, month, 'v', 'Ave', product_n = product_n)
data1 = cal_Length_of_Current_Vector(u[:, :, depthn1 - 1], v[:, :, depthn1 - 1])
data2 = cal_Length_of_Current_Vector(u[:, :, depthn2 - 1], v[:, :, depthn2 - 1])
Data[i, :, :] = data1 - data2
plta = quick.draw_12month_cycle_with_map(Data, ygrid, xgrid, cb_min, cb_max, xlim = xlim, ylim = ylim, \
clabel = "Vertical Shear of Current(m/s) " + str(zgrid[depthn1 - 1]) + 'm - ' + str(zgrid[depthn2 - 1]) + 'm',\
interval = interval, my_color = 0)
