def minus_MG_and_argopoint(year,month,var,cb_min,cb_max,area_num,depth): # area_num=0でインド洋全域 import subroutine import MOAA_GPV import WOA01 import AQC import numpy as np import matplotlib.pyplot as plt stryear,strmonth=subroutine.strym(year,month) strym=stryear+strmonth save_dir=subroutine.save_dir() my_color=0 # 青白赤 dataWOA=WOA01.nc_read(month,var,depth+1) # WOAはz=0があるので、1足してMGと一致 dataWOA=subroutine.data_trimming_IO(dataWOA,4) dataMG=MOAA_GPV.nc_read(year,month,var,depth) dataMG=subroutine.data_trimming_IO(dataMG,1) data=dataMG-dataWOA xgrid,ygrid=subroutine.IO_gridvalue() m=subroutine.IO_map(area_num,1,1) x, y = np.meshgrid(xgrid, ygrid) X, Y = m(x, y) title='' plt.title(title,fontsize=25) temp,salt,pres,lon,lat=AQC.get_data(year,month) interval_of_cf=subroutine.cl_res(cb_min,cb_max) data,cmap,cb_interval,plt=subroutine.color(data,my_color,interval_of_cf,plt) CF=m.contourf(X, Y, data,interval_of_cf, cmap=cmap, latlon=True,extend='both') cb=m.colorbar(CF, ticks=cb_interval) SC=plt.scatter(lon,lat,s=45,c='maroon') fnameF = save_dir+strym+var+'_MG-WOA.jpg' plt.savefig(fnameF) plt.clf()
def get_data(year, month, var, depth, title_name = 'ESTOCver03'): import ESTOCver01 import ESTOCver02 import ESTOCver03 import ORAS4 import MOAA_GPV import GECCO2 import WOA01 import HadISST import ERSST import COBESST import RG import ERA if title_name == 'ESTOCver01': data = ESTOCver01.masuda_data_read(year, month, var, depth) elif title_name == 'ESTOCver02': data = ESTOCver02.nc_read(year, month, var, depth) elif title_name == 'ESTOCver03': data = ESTOCver03.masuda_data_read(year, month, var, depth) elif title_name == 'MOAA_GPV': data = MOAA_GPV.nc_read(year, month, var, depth) elif title_name == 'ORAS4': data = ORAS4.nc_read(year, month, var, depth) elif title_name == 'GECCO2': data = GECCO2.grib_read(year, month, var, depth) elif title_name == 'WOA01': data = WOA01.nc_read(month, var, depth) elif title_name == 'COBESST': data = COBESST.nc_read(year, month) elif title_name == 'HadISST': data = HadISST.nc_read(year, month) elif title_name == 'ERSST': data = ERSST.nc_read(year, month) elif title_name == 'Roemmich_Gilson': data = RG.nc_read(year, month, var, depth) elif title_name == 'ERA': data = ERA.get_data(year, month, var) # これで得られるデータはインド洋領域のみなので注意。 else: raise Exception('error! your product_n is not valid!') return data
def get_map_and_argopoint(year, month, var, depth, cb_min, cb_max): import MOAA_GPV import numpy as np import matplotlib.pyplot as plt import netCDF4 my_color = 1 # レインボー product_n = 1 # MOAA_GPV data = MOAA_GPV.nc_read(year, month, var, depth) xgrid, ygrid, d = product_grid_info(var, 'data', product_n) stryear, strmonth = strym(year, month) m = IO_map(6, 1, 1) # 北インド洋全体 x, y = np.meshgrid(xgrid, ygrid) X, Y = m(x, y) title = '' plt.title(title, fontsize = 25) temp, salt, pres, lon, lat = AQC.get_data(year, month) interval_of_cf = cl_res(cb_min, cb_max) data, cmap, cb_interval, plt = color(data, my_color, interval_of_cf, plt) CF = m.contourf(X, Y, data, interval_of_cf, cmap = cmap, latlon = True, extend = 'both') cb = m.colorbar(CF, ticks = cb_interval) SC = plt.scatter(lon, lat, c = salt[:, 0], s = 75, vmin = min(interval_of_cf), vmax = max(interval_of_cf), cmap = cmap) plt.show()