try: HU_mgrd_xmax = utils_misc.loadvar(path_grd + 'HU_mgrd_xmax') # load from file except: HU_mgrd_xmax = utils_mask.calc_H_mgrd_xmax(ncdat, 'U', path_grd) # ============================================================================= # Streamfunctions # ============================================================================= lat_mgrd = ncdat.TLAT.isel( nlon=0) # mean of LAT for each j #! very inappropriate # ----------------------------------------------------------------------------- # BSF on geographical grid calculated by model fig, map = utils_plt.plot_BSF(BSF.isel(time=0), 'T', nlevels=10) plt.title('BSF model on T grid') #utils_plt.print2pdf(fig, 'testfigures/BSF_model_T') # ----------------------------------------------------------------------------- # MOC on geographical grid calculated by model MOCsel = MOC.isel(time=0) fig, ax = utils_plt.plot_MOC(MOCsel.lat_aux_grid, MOCsel.moc_z, MOCsel, nlevels=40, plttype='pcolor+contour') plt.plot(lat_mgrd, HT_auxgrd_xmax) # plot seafloor plt.title('MOC model') plt.xlim([-36, 90]) #utils_plt.print2pdf(fig, path_figs+'MOC_model')
''' import matplotlib.pyplot as plt import matplotlib as ml import CESM_utils_plt as utils_plt plt.ion() # enable interactive mode path_fig = '../figures/160711/' # ======================================================================================= # CCSM4 representations # ======================================================================================= # ----------------------------------------------------------------------------------------- # BSF on geographical grid calculated by model BSF_model = utils_mask.mask_ATLANTIC(ncdat.BSF.isel(time=0), ncdat.REGION_MASK) fig, map = utils_plt.plot_BSF(BSF_model, 'T', nlevels = 10) plt.title('BSF model on T grid') utils_plt.print2pdf(fig, path_fig+'BSF_model_T') # ----------------------------------------------------------------------------------------- # MOC on geographical grid calculated by model MOC_model = ncdat.MOC.isel(time=0, transport_reg=1, moc_comp=0) #MOC_model = MOC_model - MOC_model[:,-1] # normalisation fig, ax = utils_plt.plot_MOC(MOC_model.lat_aux_grid, MOC_model.moc_z, MOC_model, nlevels=10, plttype='pcolor+contour') plt.plot(lat_auxgrd,HT_auxgrd_xmax) # plot seafloor plt.title('MOC model') plt.xlim([-36,90]) utils_plt.print2pdf(fig, path_fig+'MOC_model') # ======================================================================================= # Calculated on model grid # =======================================================================================
''' import matplotlib.pyplot as plt import matplotlib as ml import CESM_utils_plt as utils_plt plt.ion() # enable interactive mode path_fig = '../figures/160711/' # ======================================================================================= # CCSM4 representations # ======================================================================================= # ----------------------------------------------------------------------------------------- # BSF on geographical grid calculated by model BSF_model = utils_mask.mask_ATLANTIC(ncdat.BSF.isel(time=0), ncdat.REGION_MASK) fig, map = utils_plt.plot_BSF(BSF_model, 'T', nlevels=10) plt.title('BSF model on T grid') utils_plt.print2pdf(fig, path_fig + 'BSF_model_T') # ----------------------------------------------------------------------------------------- # MOC on geographical grid calculated by model MOC_model = ncdat.MOC.isel(time=0, transport_reg=1, moc_comp=0) #MOC_model = MOC_model - MOC_model[:,-1] # normalisation fig, ax = utils_plt.plot_MOC(MOC_model.lat_aux_grid, MOC_model.moc_z, MOC_model, nlevels=10, plttype='pcolor+contour') plt.plot(lat_auxgrd, HT_auxgrd_xmax) # plot seafloor plt.title('MOC model') plt.xlim([-36, 90]) utils_plt.print2pdf(fig, path_fig + 'MOC_model')
# ============================================================================= try: HT_mgrd_xmax = utils_misc.loadvar(path_grd+'HT_mgrd_xmax') # load from file except: HT_mgrd_xmax = utils_mask.calc_H_mgrd_xmax(ncdat, 'T', path_grd) try: HU_mgrd_xmax = utils_misc.loadvar(path_grd+'HU_mgrd_xmax') # load from file except: HU_mgrd_xmax = utils_mask.calc_H_mgrd_xmax(ncdat, 'U', path_grd) # ============================================================================= # Streamfunctions # ============================================================================= lat_mgrd = ncdat.TLAT.isel(nlon=0) # mean of LAT for each j #! very inappropriate # ----------------------------------------------------------------------------- # BSF on geographical grid calculated by model fig, map = utils_plt.plot_BSF(BSF.isel(time=0), 'T', nlevels = 10) plt.title('BSF model on T grid') #utils_plt.print2pdf(fig, 'testfigures/BSF_model_T') # ----------------------------------------------------------------------------- # MOC on geographical grid calculated by model MOCsel = MOC.isel(time=0) fig, ax = utils_plt.plot_MOC(MOCsel.lat_aux_grid, MOCsel.moc_z, MOCsel, nlevels=40, plttype='pcolor+contour') plt.plot(lat_mgrd, HT_auxgrd_xmax) # plot seafloor plt.title('MOC model') plt.xlim([-36,90]) #utils_plt.print2pdf(fig, path_figs+'MOC_model') # ============================================================================= # Streamfunction Indices # =============================================================================