parser.error('--%s: wrong argument' % selname) select[selname] = val # Imports from vcmq import DS, map2, MV2, data_sample # The file ncfile = options.ncfile if not os.path.exists(ncfile): ncfile = data_sample(ncfile) if not os.path.exists(ncfile): parser.error('File not found: ' + ncfile) # Read temperature and depth ds = DS(ncfile, options.model, **select) mld = ds.get_mld(mode='deltatemp', squeeze=1) # Time average (if needed) if mld.ndim == 3: mld = MV2.average(mld, axis=0) # Plot map2(mld, proj='merc', figsize=(6, 6), autoresize=0, title=options.title, colorbar_shrink=0.7, right=1, show=options.disp, savefig=options.outputfig,
try: lon = eval(options.lon) except: parser.error('--%s: wrong argument'%selname) select[selname] = val # Imports from vcmq import DS, map2, MV2, data_sample # The file ncfile = options.ncfile if not os.path.exists(ncfile): ncfile = data_sample(ncfile) if not os.path.exists(ncfile): parser.error('File not found: '+ncfile) # Read temperature and depth ds = DS(ncfile, options.model, **select) mld = ds.get_mld(mode='deltatemp', squeeze=1) # Time average (if needed) if mld.ndim==3: mld = MV2.average(mld, axis=0) # Plot map2(mld, proj='merc', figsize=(6, 6), autoresize=0, title=options.title, colorbar_shrink=0.7, right=1, show=options.disp, savefig=options.outputfig, savefig_verbose=True)
print ds.dataset # Logging ds.info('Pratique') # -> ESSAYER WARNING ds.set_loglevel('debug') # Coordonnees lon = ds.get_lon() # -> ESSAYER AU POINT U grid = ds.get_grid() # Variables temp = ds.get_temp(level=slice(-1, None), squeeze=True) sst = ds.get_sst(squeeze=True) # -> VERIFIER ATTRIBUTS ds.info('Plot SST') map2(sst) sal = ds.get_sal(lon=(4., 4.5)) print temp.shape, sal.shape # Generique ds2 = DS(data_sample('mfs.nc'), 'nemo', lon=(4, 5), lat=(42.5, 43.5)) sst2 = ds2.get_sst(squeeze=True) # -> VERIFIER ATTRIBUTS # Plus evolue depth = ds.get_depth() # sigma print depth.shape mld = ds.get_mld() ds.plot_transect('temp', lons=(4.1, 4.9), lats=(43.3, 42.6), figsize=(12, 6)) # -> OUTAXIS=
"""Test :meth:`~vacumm.data.misc.dataset.OceanDataset.get_mld` with deltadens in MENOR""" # Inits ncfile = "menor.nc" # Imports from vcmq import DS, os, map2, data_sample, code_file_name # Read data ds = DS(data_sample(ncfile), 'mars', logger_level='critical') mld = ds.get_mld(mode='deltadens', deltadens=0.06, squeeze=True) # Plot surface density figfile = code_file_name(ext='png') if os.path.exists(figfile): os.remove(figfile) map2(mld, savefig=figfile, close=True, show=False, fill='pcolormesh', cmap='vacumm_previmer', contour=False) # For unittest result = {'files':figfile}
"""Test :meth:`~vacumm.data.misc.dataset.OceanDataset.get_mld` with deltadens in MENOR""" # Inits ncfile = "menor.nc" # Imports from vcmq import DS, os, map2, data_sample, code_file_name # Read data ds = DS(data_sample(ncfile), 'mars', logger_level='critical') mld = ds.get_mld(mode='deltadens', deltadens=0.06, squeeze=True) # Plot surface density figfile = code_file_name(ext='png') if os.path.exists(figfile): os.remove(figfile) map2(mld, savefig=figfile, close=True, show=False, fill='pcolormesh', cmap='vacumm_previmer', contour=False) # For unittest result = {'files': figfile}